We present experimental phase functions of three types of millimeter-sized dust grains consisting of enstatite, quartz and volcanic material from Mount Etna, respectively. The three grains present similar sizes but different absorbing properties. The measurements are performed at 527 nm covering the scattering angle range from 3 to 170 degrees. The measured phase functions show two well defined regions i) soft forward peaks and ii) a continuous increase with the scattering angle at side- and back-scattering regions. This behavior at side- and back-scattering regions are in agreement with the observed phase functions for the Fomalhaut and HR 4796A dust rings. Further computations and measurements (including polarization) for millimeter sized-grains are needed to draw some conclusions about the fluffy or compact structure of the dust grains.

Looking for the primary rainbow in starlight that is reflected by exoplanets appears to be a promising method to search for liquid water clouds in exoplanetary atmospheres. Ice water clouds, that consist of water crystals instead of water droplets, could potentially mask the rainbow feature in the planetary signal by covering liquid water clouds. Here, we investigate the strength of the rainbow feature for exoplanets that have liquid and icy water clouds in their atmosphere, and calculate the rainbow feature for a realistic cloud coverage of Earth. We calculate flux and polarization signals of starlight that is reflected by horizontally and vertically inhomogeneous Earth--like exoplanets, covered by patchy clouds consisting of liquid water droplets or water ice crystals. The planetary surfaces are black. On a planet with a significant coverage of liquid water clouds only, the total flux signal shows a weak rainbow feature. Any coverage of the liquid water clouds by ice clouds, however, dampens the rainbow feature in the total flux, and thus the discovery of liquid water in the atmosphere. On the other hand, detecting the primary rainbow in the polarization signal of exoplanets appears to be a powerful tool for detecting liquid water in exoplanetary atmospheres, even when these clouds are partially covered by ice clouds. In particular, liquid water clouds covering as little as 10%-20% of the planetary surface, with more than half of these covered by ice clouds, still create a polarized rainbow feature in the planetary signal. Indeed, calculations of flux and polarization signals of an exoplanet with a realistic Earth--like cloud coverage, show a strong polarized rainbow feature.

In this paper we present two efficient implementations of the diffusion approximation to be employed in Monte Carlo computations of radiative transfer in dusty media of massive circumstellar disks. The aim is to improve the accuracy of the computed temperature structure and to decrease the computation time. The accuracy, efficiency and applicability of the methods in various corners of parameter space are investigated. The effects of using these methods on the vertical structure of the circumstellar disk as obtained from hydrostatic equilibrium computations are also addressed. Two methods are presented. First, an energy diffusion approximation is used to improve the accuracy of the temperature structure in highly obscured regions of the disk, where photon counts are low. Second, a modified random walk approximation is employed to decrease the computation time. This modified random walk ensures that the photons that end up in the high-density regions can quickly escape to the lower density regions, while the energy deposited by these photons in the disk is still computed accurately. A new radiative transfer code, MCMax, is presented in which both these diffusion approximations are implemented. These can be used simultaneously to increase both computational speed and decrease statistical noise. We conclude that the diffusion approximations allow for fast and accurate computations of the temperature structure, vertical disk structure and observables of very optically thick circumstellar disks.

We present measurements of ratios of elements of the scattering matrix of Martian analogue palagonite particles for scattering angles ranging from 3 to 174 degrees and a wavelength of 632.8 nm. To facilitate the use of these measurements in radiative transfer calculations we have devised a method that enables us to obtain, from these measurements, a normalized synthetic scattering matrix covering the complete scattering angle range from 0 to 180 degrees. Our method is based on employing the coefficients of the expansions of scattering matrix elements into generalized spherical functions. The synthetic scattering matrix elements and/or the expansion coefficients obtained in this way, can be used to include multiple scattering by these irregularly shaped particles in (polarized) radiative transfer calculations, such as calculations of sunlight that is scattered in the dusty Martian atmosphere.

In order to deduce properties of dust in astrophysical environments where dust growth through aggregation is important, knowledge of the way aggregated particles interact with radiation, and what information is encoded in the thermal radiation they emit, is needed. The emission characteristics are determined by the size and structure of the aggregate and the composition and shape of the constituents. We thus aim at performing computations of compositionally inhomogeneous aggregates composed of irregularly shaped constituents. In addition we aim at developing an empirical recipe to compute the optical properties of such aggregates in a fast and accurate manner. We performed CDA computations for aggregates of irregularly shaped particles with various compositions. The constituents of the aggregate are assumed to be in the Rayleigh regime, and in addition we assume that the dominant interaction of the aggregate constituents is through dipole-dipole interactions. We computed the spectral structure of the emission efficiency in the 10 micron region for aggregates with 30% amorphous carbon and 70% silicates by volume with various fractions of crystalline and amorphous components. We find that the spectral appearance of the various components of the aggregate are very different and depend on their abundances. Most notably, materials that have a very low abundance appear spectroscopically as if they were in very small grains, while more abundant materials appear, spectroscopically to reside in larger grains. We construct a fast empirical approximate method, based on the idea of an effective medium approximation, to construct the spectra for these aggregates which almost perfectly reproduces the more exact computations. This new method is fast enough to be easily implemented in fitting procedures.

We investigate the composition and shape distribution of silicate dust grains in the interstellar medium. The effect of the amount of magnesium in the silicate lattice is studied. We fit the spectral shape of the interstellar 10 mu extinction feature as observed towards the galactic center. We use very irregularly shaped coated and non-coated porous Gaussian Random Field particles as well as a statistical approach to model shape effects. For the dust materials we use amorphous and crystalline silicates with various composition and SiC. The results of our analysis of the 10 mu feature are used to compute the shape of the 20 mu silicate feature and to compare this with observations. By using realistic particle shapes we are, for the first time, able to derive the magnesium fraction in interstellar silicates. We find that the interstellar silicates are highly magnesium rich (Mg/(Fe+Mg)>0.9) and that the stoichiometry lies between pyroxene and olivine type silicates. This composition is not consistent with that of the glassy material found in GEMS in interplanetary dust particles indicating that these are, in general, not unprocessed remnants from the interstellar medium. Also, we find a significant fraction of SiC (~3%). We discuss the implications of our results for the formation and evolutionary history of cometary and circumstellar dust. We argue that the fact that crystalline silicates in cometary and circumstellar grains are almost purely magnesium silicates is a natural consequence of our findings that the amorphous silicates from which they were formed were already magnesium rich.

We model the 10 micron absorption spectra of nonspherical particles composed of amorphous silicate. We consider two classes of particles, compact ones and fractal aggregates composed of homogeneous spheres. For the compact particles we consider Gaussian random spheres with various degrees of non-sphericity. For the fractal aggregates we compute the absorption spectra for various fractal dimensions. The 10 micron spectra are computed for ensembles of these particles in random orientation using the well-known Discrete Dipole Approximation. We compare our results to spectra obtained when using volume equivalent homogeneous spheres and to those computed using a porous sphere approximation. We conclude that, in general, nonspherical particles show a spectral signature that is similar to that of homogeneous spheres with a smaller material volume. This effect is overestimated when approximating the particles by porous spheres with the same volume filling fraction. For aggregates with fractal dimensions typically predicted for cosmic dust, we show that the spectral signature characteristic of very small homogeneous spheres (with a volume equivalent radius r_V<0.5 micron) can be detected even in very large particles. We conclude that particle sizes are underestimated when using homogeneous spheres to model the emission spectra of astronomical sources. In contrast, the particle sizes are severely overestimated when using equivalent porous spheres to fit observations of 10 micron silicate emission.

We discuss the composition and size distribution of the dust in the coma of comet Hale-Bopp. We do this by fitting simultaneously the infrared emission spectrum measured by the infrared space observatory (ISO) and the measured degree of linear polarization of scattered light at various phase angles and 12 different wavelengths. The effects of particle shape on the modeled optical properties of the dust grains are taken into account. We constrain our fit by forcing the abundances of the major rock forming chemical elements to be solar. The infrared spectrum at long wavelengths reveals that large grains are needed in order to fit the spectral slope. The size and shape distribution we employ allows us to estimate the sizes of the crystalline silicates. The ratios of the strength of various forsterite features show that the crystalline silicate grains in Hale-Bopp must be submicron sized. We exclude the presence of large crystalline silicate grains in the coma. Because of this lack of large crystalline grains combined with the fact that we do need large amorphous grains to fit the emission spectrum at long wavelengths, we need only approximately 4% of crystalline silicates by mass. After correcting for possible hidden crystalline material included in large amorphous grains, our best estimate of the total mass fraction of crystalline material is approximately 7.5%, significantly lower than deduced in previous studies in which the typical derived crystallinity is 20-30%. The implications of this on the possible origin and evolution of the comet are discussed. The crystallinity we observe in Hale-Bopp is consistent with the production of crystalline silicates in the inner solar system by thermal annealing and subsequent radial mixing to the comet forming region.

We provide a theoretical foundation for the statistical approach for computing the absorption properties of particles in the Rayleigh domain. We present a general method based on the Discrete Dipole Approximation (DDA) to compute the absorption and scattering properties of particles in the Rayleigh domain. The method allows to separate the geometrical aspects of a particle from its material properties. Doing the computation of the optical properties of a particle once, provides them for any set of refractive indices, wavelengths and orientations. This allows for fast computations of e.g. absorption spectra of arbitrarily shaped particles. Other practical applications of the method are in the interpretation of atmospheric and radar measurements as well as computations of the scattering matrix of small particles as a function of the scattering angle. In the statistical approach, the optical properties of irregularly shaped particles are represented by the average properties of an ensemble of particles with simple shapes. We show that the absorption cross section of an ensemble of arbitrarily shaped particles with arbitrary orientations can always be uniquely represented by the average absorption cross section of an ensemble of spheroidal particles with the same composition and fixed orientation. This proves for the first time that the statistical approach is generally viable in the Rayleigh domain.

In this paper we study the combined effects of size and shape of small solid state particles on the absorption, emission and scattering characteristics. We use the statistical approach to calculate these optical properties. In this approach the average optical properties of an ensemble of particles in random orientation are represented by the average optical properties of an ensemble of simple shapes. The validity of this approach is studied in detail for a uniform distribution of hollow spheres where the fractional volume of the central inclusion is varied. We apply the results to two different areas of interest, i) infrared spectroscopy and ii) polarization of scattered light. The effects of particle size and shape on the optical characteristics are discussed. We compare the results using the distribution of hollow spheres with those obtained by using randomly oriented spheroids. Also we compare the results with observations and laboratory measurements. The distribution of hollow spheres is very successful in reproducing laboratory measurements of the scattering angle distribution of the degree of linear polarization for incident unpolarized light of randomly oriented irregular quartz particles. Furthermore, we show that we are able to derive the size distribution of dust grains by fitting the measured degree of linear polarization using computational result for hollow spheres. It is shown that the distribution of hollow spheres is a powerful tool for studying light scattering, absorption and emission by cosmic dust grains and in particular when large numbers of particle parameters need to be considered since the computational demand of the distribution of hollow spheres is small.