It is well known that the dust properties of the diffuse interstellar medium exhibit variations towards different sight-lines on a large scale. We have investigated the variability of the dust characteristics on a small scale, and from cloud-to-cloud. We use low-resolution spectro-polarimetric data obtained in the context of the Large Interstellar Polarisation Survey (LIPS) towards 59 sight-lines in the Southern Hemisphere, and we fit these data using a dust model composed of silicate and carbon particles with sizes from the molecular to the sub-micrometre domain. Large (> 6 nm) silicates of prolate shape account for the observed polarisation. For 32 sight-lines we complement our data set with UVES archive high-resolution spectra, which enable us to establish the presence of single-cloud or multiple-clouds towards individual sight-lines. We find that the majority of these 35 sight-lines intersect two or more clouds, while eight of them are dominated by a single absorbing cloud. We confirm several correlations between extinction and parameters of the Serkowski law with dust parameters, but we also find previously undetected correlations between these parameters that are valid only in single-cloud sight-lines. We find that interstellar polarisation from multiple-clouds is smaller than from single-cloud sight-lines, showing that the presence of a second or more clouds depolarises the incoming radiation. We find large variations of the dust characteristics from cloud-to-cloud. However, when we average a sufficiently large number of clouds in single-cloud or multiple-cloud sight-lines, we always retrieve similar mean dust parameters. The typical dust abundances of the single-cloud cases are [C]/[H] = 92 ppm and [Si]/[H] = 20 ppm.

Polarimetric studies of light transmitted through interstellar clouds may give constraints on the properties of the interstellar dust grains. Traditionally, broadband linear polarisation (BBLP) measurements have been considered an important diagnostic tool for the study of the interstellar dust, while comparatively less attention has been paid to spectropolarimetric measurements. However, spectropolarimetry may offer stronger constraints than BBLP, for example by revealing narrowband features, and by allowing us to distinguish the contribution of dust from the contribution of interstellar gas. Therefore, we have decided to carry out a Large Interstellar Polarisation Survey (LIPS) using spectropolarimetric facilities in both hemispheres. Here we present the results obtained in the Southern Hemisphere with the FORS2 instrument of the ESO Very Large Telescope. Our spectra cover the wavelength range 380--950\,nm at a spectral resolving power of about 880. We have produced a publicly available catalogue of 127 linear polarisation spectra of 101 targets. We also provide the Serkowski-curve parameters, as well as the wavelength gradient of the polarisation position angle for the interstellar polarisation along 76 different lines of sight. In agreement with previous literature, we found that the best-fit parameters of the Serkowski-curve are not independent of each other. However, the relationships that we obtained are not always consistent with what was found in previous studies.

We present polarisation spectra of seven stars in the lines-of-sight towards the Sco OB1 association. Our spectra were obtained within the framework of the Large Interstellar Polarization Survey carried out with the FORS instrument of the ESO VLT. We have modelled the wavelength-dependence of extinction and linear polarisation with a dust model for the diffuse interstellar medium which consists of a mixture of particles with size ranging from the molecular domain of 0.5 nm up to 350 nm. We have included stochastically heated small dust grains with radii between 0.5 and 6 nm made of graphite and silicate, as well as polycyclic aromatic hydrocarbon molecules (PAHs), and we have assumed that larger particles are prolate spheroids made of amorphous carbon and silicate. Overall, a dust model with eight free parameters best reproduces the observations. Reducing the number of free parameters leads to results that are inconsistent with cosmic abundance constraints. We found that aligned silicates are the dominant contributor to the observed polarisation, and that the polarisation spectra are best-fit by a lower limit of the equivolume sphere radius of aligned grains of 70 - 200nm.

We interpret the interstellar extinction observed towards the Galactic Center (GC) in the wavelength range $\lambda = 1 - 20\,\mkm$. Its main feature is the flat extinction at $3 - 8\,\mkm$ whose explanation is still a problem for the cosmic dust models. We search for structure and chemical composition of dust grains that could explain the observed extinction. In contrast to earlier works we use laboratory measured optical constants and consider particles of different structure. We show that a mixture of compact grains of aromatic carbon and of some silicate is better suited for reproducing the flat extinction in comparison with essentially porous grains or aliphatic carbon particles. Metallic iron should be located inside the particle, i.e. cannot form layers on silicate grains as the extinction curves become then very peculiar. We find a model including aromatic carbonaceous particles and three-layered particles with an olivine-type silicate core, a thin very porous layer and a thin envelope of magnetite that provides a good (but still not perfect) fit to the observational data. We suggest that such silicate dust should be fresh, i.e. recently formed in the atmospheres of late-type stars in the central region of the Galaxy. We assume that this region has a radius of about 1 kpc and produces about a half of the observed extinction. The remaining part of extinction is caused by a "foreground" material being practically transparent at $\lambda = 4 - 8\,\mkm$.

We interpret the relation between the polarizing efficiency $P_{\max}/E(B-V)$ and the wavelength of the maximum polarization $\lambda_{\max}$ observed for 17 objects (including 243 stars) separated into two groups: "dark clouds" and "open clusters". The objects are assigned to one of the groups according to the distribution of the parameter $\lambda_{\max}$. We use the model of homogeneous silicate and carbonaceous spheroidal particles with the imperfect alignment and a time-evolving size distribution. The polarization is assumed to be mainly produced by large silicate particles with the sizes $r_{V} \ga r_{V,\rm cut}$. The models with the initial size distribution reproducing the average curve of the interstellar extinction fail to explain the values of $\lambda_{\max} \ga 0.65\,\mkm$ observed for several dark clouds. We assume that the grain size distribution is modified due to accretion and coagulation, according to the model of Hirashita \& Voshchinnikov (2014). After including the evolutionary effects, $\lambda_{\max}$ shifts to longer wavelengths on time-scales $\sim 20 (n_\mathrm{H}/10^3 \mathrm{cm}^{-3})^{-1}$ Myr where $n_\mathrm{H}$ is the hydrogen density in molecular clouds where dust processing occurs. The ratio $P_{\max}/E(B-V)$ goes down dramatically when the size of polarizing grains grows. The variations of the degree and direction of particle orientation influence this ratio only moderately. We have also found that the aspect ratio of prolate grains does not affect significantly the polarizing efficiency. For oblate particles, the shape effect is stronger but in most cases the polarization curves produced are too narrow in comparison with the observed ones.

Massive stars live short lives, losing large amounts of mass through their stellar wind. Their mass is a key factor determining how and when they explode as supernovae, enriching the interstellar medium with heavy elements and dust. During the red supergiant phase, mass-loss rates increase prodigiously, but the driving mechanism has proven elusive. Here we present high-contrast optical polarimetric-imaging observations of the extreme red supergiant VY Canis Majoris and its clumpy, dusty, mass-loss envelope, using the new extreme-adaptive-optics instrument SPHERE at the VLT. These observations allow us to make the first direct and unambiguous detection of submicron dust grains in the ejecta; we derive an average grain radius $\sim$ 0.5 $\mu$m, 50 times larger than in the diffuse ISM, large enough to receive significant radiation pressure by photon scattering. We find evidence for varying grain sizes throughout the ejecta, highlighting the dynamical nature of the envelope. Grains with 0.5 $\mu$m sizes are likely to reach a safe distance from the eventual explosion of VY Canis Majoris; hence it may inject upwards of 10$^{-2}$ M$_\odot$ of dust into the ISM.

We apply the time evolution of grain size distributions by accretion and coagulation found in our previous work to the modelling of the wavelength dependence of interstellar linear polarization. We especially focus on the parameters of the Serkowski curve $K$ and $\lambda_{\max}$ characterizing the width and the maximum wavelength of this curve, respectively. We use aligned silicate and non-aligned carbonaceous spheroidal particles with different aspect ratios $a/b$. The imperfect alignment of grains with sizes larger than a cut-off size $r_{V,\rm cut}$ is considered. We find that the evolutionary effects on the polarization curve are negligible in the original model with commonly used material parameters (hydrogen number density $n_\mathrm{H}=10^3$ cm$^{-3}$, gas temperature $T_\mathrm{gas}=10$~K, and the sticking probability for accretion $S_\mathrm{acc}=0.3$). Therefore, we apply the tuned model where the coagulation threshold of silicate is removed. In this model, $\lambda_{\max}$ displaces to the longer wavelengths and the polarization curve becomes wider ($K$ reduces) on time-scales $\sim (30 - 50) (n_\mathrm{H}/10^3 \mathrm{cm}^{-3})^{-1}$ Myr. The tuned models at $T < 30 (n_\mathrm{H}/10^3 \mathrm{cm}^{-3})^{-1} $ Myr and different values of the parameters $r_{V,\rm cut}$ can also explain the observed trend between $K$ and $\lambda_{\max}$. It is significant that the evolutionary effect appears in the perpendicular direction to the effect of $r_{V,\rm cut}$ on the $K$ - $\lambda_{\max}$ diagram. Very narrow polarization curves can be reproduced if we change the type of particles (prolate/oblate) and/or vary $a/b$.

Polarimetric observations of about 50 stars located in a close vicinity of the Bok globule CB67 having significantly nonspherical shape and rapid rotation are performed. The data obtained are compared with the available observations of this globule at radio and submillimeter wavelengths as well as some theoretical calculations. It is found that the elongation and the rotation moment of CB67 are oriented rather perpendicular to the magnetic fields, which is unusual for Bok globules and is difficult to be explained from the theoretical point of view.

Dust grains grow their sizes in the interstellar clouds (especially in molecular clouds) by accretion and coagulation. Here we model and test these processes by examining the consistency with the observed variation of the extinction curves in the Milky Way. We find that, if we simply use the parameters used in previous studies, the model fails to explain the flattening of far-UV extinction curve for large $R_V$ (flatness of optical extinction curve) and the existence of carbon bump even in flat extinction curves. This discrepancy is resolved by adopting a `tuned' model, in which coagulation of carbonaceous dust is less efficient (by a factor of 2) and that of silicate is more efficient with the coagulation threshold removed. The tuned model is also consistent with the relation between silicon depletion (indicator of accretion) and $R_V$ if the duration of accretion and coagulation is >100(n_H/10^3 cm^{-3})^{-1} Myr, where n_H is the number density of hydrogen nuclei in the cloud. We also examine the relations between each of the extinction curve features (UV slope, far-UV curvature, and carbon bump strength) and $R_V$. The correlation between UV slope and $R_V$, which is the strongest among the three correlations, is well reproduced by the tuned model. For far-UV curvature and carbon bump strength, the observational data are located between the tuned model and the original model without tuning, implying that the large scatters in the observational data can be explained by the sensitive response to the coagulation efficiency. The overall success of the tuned model indicates that accretion and coagulation are promising mechanisms of producing the variation of extinction curves in the Milky Way, although we do not exclude possibilities of other dust-processing mechanisms changing extinction curves.

We present a model for the diffuse interstellar dust that explains the observed wavelength-dependence of extinction, emission, linear and circular polarisation of light. The model is set-up with a small number of parameters. It consists of a mixture of amorphous carbon and silicate grains with sizes from the molecular domain of 0.5 up to about 500nm. Dust grains with radii larger than 6nm are spheroids. Spheroidal dust particles have a factor 1.5 - 3 larger absorption cross section in the far IR than spherical grains of the same volume. Mass estimates derived from submillimeter observations that ignore this effect are overestimated by the same amount. In the presence of a magnetic field, spheroids may be partly aligned and polarise light. We find that polarisation spectra help to determine the upper particle radius of the otherwise rather unconstrained dust size distribution. Stochastically heated small grains of graphite, silicates and polycyclic aromatic hydrocarbons (PAHs) are included. We tabulate parameters for PAH emission bands in various environments. They show a trend with the hardness of the radiation field that can be explained by the ionisation state or hydrogenation coverage of the molecules. For each dust component its relative weight is specified, so that absolute element abundances are not direct input parameters. The model is confronted with the average properties of the Milky Way, which seems to represent dust in the solar neighbourhood. It is then applied to four specific sight lines including the reflection nebula NGC2023. For these sight lines, we present linear and circular spectro-polarimetric observations obtained with FORS/VLT. Using prolate rather than oblate grains gives a better fit to observed spectra; the axial ratio of the spheroids is typically two and aligned silicates are the dominant contributor to the polarisation.

A quantitative interpretation of the observed relation between the interstellar linear polarization curve parameters $K$ and $\lambda_{\max}$ characterizing the width and the wavelength of a polarization maximum, respectively, is given. The observational data available for 57 stars located in the dark clouds in Taurus, Chamaeleon, around the stars $\rho$ Oph and R CrA are considered. The spheroidal particle model of interstellar dust grains earlier applied to simultaneously interpret the interstellar extinction and polarization curves in a wide spectral range is utilized. The observed trend $K \approx 1.7 \lambda_{\max}$ is shown to be most likely related to a growth of dust grains due to coagulation rather than mantle accretion. The relation of the parameters $K$ and $\lambda_{\max}$ with an average size of silicate dust grains is discussed.

The review contains an analysis of the observed and model curves of the interstellar extinction and polarization. The observations mainly give information on dust in diffuse and translucent interstellar clouds. The features of various dust grain models including spherical/non-spherical, homogeneous/inhomogeneous particles are discussed. A special attention is devoted to the analysis of the grain size distributions, alignment mechanisms and magnetic field structure in interstellar clouds. It is concluded that the interpretation of interstellar extinction and polarization is not yet complete.

We compiled the polarimetric data for a sample of lines of sight with known abundances of Mg, Si, and Fe. We correlated the degree of interstellar polarization $P$ and polarization efficiency (the ratio of $P$ to the colour excess $E(B-V)$ or extinction $A_V$) with dust phase abundances. We detect an anticorrelation between $P$ and the dust phase abundance of iron in non silicate - containing grains $< [\rm Fe(rest)/H > ]_\rm d$, a correlation between $P$ and the abundance of Si, and no correlation between $P/E(B-V)$ or $P/A_V$ and dust phase abundances. These findings can be explained if mainly the silicate grains aligned by the radiative mechanism are responsible for the observed interstellar linear polarization.

The light scattering problem for a confocal multilayered spheroid has been solved by the extended boundary condition method (EBCM) with a corresponding spheroidal basis. The solution preserves the advantages of the approach applied previously to homogeneous and core-mantle spheroids, i.e. the separation of the radiation fields into two parts and a special choice of scalar potentials for each of the parts. The method is known to be useful in a wide range of the particle parameters. It is particularly efficient for strongly prolate and oblate spheroids. Numerical tests are described. Illustrative calculations have shown that the extinction factors to converge to average values with a growing number of layers and how the extinction vary with a growth of particle porosity.

We analyse observational correlations for three elements entering into the composition of interstellar silicate and oxide grains. Using current solar abundances (Asplund et al. 2009), we convert the gas-phase abundances into dust-phase abundances for 196 sightlines. We deduce a sharp difference in abundances for sightlines located at low ($|b|<30\degr$) and high ($|b|>30\degr$) galactic latitudes. For high-latitude stars the ratios Mg/Si and Fe/Si in dust are close to 1.5. For disk stars they are reduced to ${\rm Mg/Si} \sim 1.2$ and ${\rm Fe/Si} \sim 1.05$. The derived numbers indicate that 1) the dust grains cannot be the mixture of silicates with olivine and pyroxene composition only and some amount of magnesium or iron (or both) should be in another population and 2) the destruction of Mg-rich grains in the warm medium is more effective than of Fe-rich grains. We reveal a decrease of dust-phase abundances and correspondingly an increase of gas-phase abundances with distance $D$ for stars with $D\ga 400$\,pc. We attribute this fact to an observational selection effect: a systematic trend toward smaller observed hydrogen column density for distant stars. We find differences in abundances for disk stars with low ($E({\rm B-V}) \la 0.2$) and high ($E({\rm B-V}) \ga 0.2$) reddenings which reflect the distinction between the sightlines passing through diffuse and translucent interstellar clouds. For Scorpius-Ophiuchus we detect an uniform increase of dust-phase abundances of Mg and Si with an increase of the ratio of total to selective extinction $R_{\rm V}$ and a decrease of the strength of the far-UV extinction. This is the first evidence for a growth of Mg-Si grains due to accretion in the interstellar medium.

We extend and investigate the spheroidal model of interstellar dust grains used to simultaneously interpret the observed interstellar extinction and polarization curves. We compare our model with similar models recently suggested by other authors, study its properties and apply it to fit the normalized extinction $A(\lambda)/A_{\rm V}$ and the polarizing efficiency $P(\lambda)/A(\lambda)$ measured in the near IR to far UV region for several stars seen through one large cloud. We conclude that the model parameter $\Omega$ being the angle between the line of sight and the magnetic field direction can be more or less reliably determined from comparison of the theory and observations. This opens a way to study the spatial structure of interstellar magnetic fields by using multi-wavelength photometric and polarimetric observations.

The flattening of the 10mu silicate emission feature observed in the spectra of T Tauri and Herbig Ae/Be stars is usually interpreted as an indicator of grain growth. We show in this paper that a similar behaviour of the feature shape occurs when the porosity of composite grains varies. The fluffy aggregates, having inclusions of different sizes, were modeled by multi-layered spheres consisting of amorphous carbon, amorphous silicate and vacuum. It is also found that the inclusion of crystalline silicates in composite porous particles can lead to a shift of the known resonances and production of new ones.

We calculate the wavelength dependence of the ratio of the linear polarization degree to extinction (polarizing efficiency) $P(\lambda)/A(\lambda)$ from the ultraviolet to near-infrared. The prolate and oblate particles with aspect ratios from $a/b=1.1$ up to 10 are assumed to be rotating and partially aligned with the mechanism of paramagnetic relaxation (Davis--Greenstein). Size/shape/orientation effects are analyzed. It is found that the wavelength dependence of $P(\lambda)/A(\lambda)$ is mainly determined by the particle composition and size whereas the values of $P(\lambda)/A(\lambda)$ depend on the particle shape, degree and direction of alignment.

We study the extinction efficiencies as well as scattering properties of particles of different porosity. Calculations are performed for porous pseudospheres with small size (Rayleigh) inclusions using the discrete dipole approximation. Five refractive indices of materials covering the range from $1.20+0.00i$ to $1.75+0.58i$ were selected. They correspond to biological particles, dirty ice, silicate, amorphous carbon and soot in the visual part of spectrum. We attempt to describe the optical properties of such particles using Lorenz-Mie theory and a refractive index found from some effective medium theory (EMT) assuming the particle is homogeneous. We refer to this as the effective model. It is found that the deviations are minimal when utilizing the EMT based on the Bruggeman mixing rule. Usually the deviations in extinction factor do not exceed $\sim 5%$ for particle porosity ${\cal P}=0 - 0.9$ and size parameters $x_{\rm porous} = 2 \pi r_{\rm s, porous}/\lambda \la 25$. The deviations are larger for scattering and absorption efficiencies and smaller for particle albedo and asymmetry parameter. Our calculations made for spheroids confirm these conclusions. Preliminary consideration shows that the effective model represents the intensity and polarization of radiation scattered by fluffy aggregates quite well. Thus, the effective models of spherical and non-spherical particles can be used to significantly simplify computations of the optical properties of aggregates containing only Rayleigh inclusions.

We present a compositional analysis of 8-13um spectra of 32 young stellar objects (YSOs). Our sample consists of 5 intermediate-mass stars and 27 low-mass stars. While the spectra and first scientific results have already been published by Przygodda et al. (2003) and Kessler-Silacci et al. (2004) we perform a more detailed analysis of the 10um silicate feature. In our analysis we assume that this emission feature can be represented by a linear superposition of the wavelength-dependent opacity $\kappa_{\rm abs}(\lambda)$ describing the optical properties of silicate grains with different chemical composition, structure, and grain size. The determination of an adequate fitting equation is another goal of this study. Using a restricted number of fitting parameters we investigate which silicate species are necessary for the compositional fitting. Particles with radii of 0.1um- and 1.5um consisting of amorphous olivine and pyroxene, forsterite, enstatite, and quartz have been considered. Only compact, homogeneous dust grains have been used in the presented fitting procedures. In this context we show that acceptable fitting results can also be achieved if emission properties of porous silicate grains are considered instead. Although some previous studies give reasons for the similarity between the dust in circumstellar disks of TTauri stars and Herbig Ae/Be stars, a quantitative comparison has been missing, so far. Therefore, we conclude with a discussion of the results of a 10um spectroscopic survey of van Boekel et al. (2005) who focus on Herbig Ae/Be stars, the higher mass counterparts of T Tauri stars and draw comparisons to this and other studies. We find that the results of our study of T Tauri systems partly agree with previous studies of Herbig Ae/Be stars.

In many models of dusty objects in space the grains are assumed to be composite or fluffy. However, the computation of the optical properties of such particles is still a very difficult problem. We analyze how the increase of grain porosity influences basic features of cosmic dust -- interstellar extinction, dust temperature, infrared bands and millimeter opacity. Porous grains can reproduce the flat extinction across the $3 - 8 \mkm$ wavelength range measured for several lines of sight by {\it ISO} and {\it Spitzer}. Porous grains are generally cooler than compact grains. At the same time, the temperature of very porous grains becomes slightly larger in the case of the EMT-Mie calculations in comparison with the results found from the layered-sphere model. The layered-sphere model predicts a broadening of infrared bands and a shift of the peak position to larger wavelengths as porosity grows. In the case of the EMT-Mie model variations of the feature profile are less significant. It is also shown that the millimeter mass absorption coefficients grow as porosity increases with a faster growth occurring for particles with Rayleigh/non-Rayleigh inclusions. As a result, for very porous particles the coefficients given by two models can differ by a factor of about 3. It is found that an increase of porosity leads to an increase of extinction cross sections at some wavelengths and a decrease at others depending on the grain model. However, this behaviour is sufficient to reproduce the extinction curve in the direction of the star $\sigma$ Sco using current solar abundances. In the case of the star $\zeta$ Oph our model requires larger amounts of carbon and iron in the dust-phase than is available.

There are indications that interstellar and interplanetary dust grains have an inhomogeneous and fluffy structure. We investigate different methods to describe light scattering by such composite particles. Both a model of layered particles and discrete dipole calculations for particles with Rayleigh and non-Rayleigh inclusions are used. The calculations demonstrate that porosity is a key parameter for determining light scattering. We find that the optical properties of the layered particles depend on the number and position of layers if the number of layers is small ($\la 15$). For a larger number of layers the scattering characteristics become independent of the layer sequence. The optical properties of particles with inclusions depend on the size of inclusions provided the porosity is large. The scattering characteristics of very porous particles with inclusions of different sizes are found to be close to those of multi-layered spheres. We compare the results of these calculations with the predictions of the effective medium theories (EMT) which are often used in astronomy as a tool to calculate the optical properties of composite particles. The results of our analysis show that the internal structure of grains (layers versus inclusions) only slightly affects the optics of particles provided the porosity does not exceed 50%. It is also demonstrated that in this case the optical properties of composite grains calculated with EMT agree with the results of the exact method for layered particles. For larger porosity, the standard EMT rules (i.e., Garnett and Bruggeman rules) give reliable results for particles with Rayleigh inclusions only.

We present a computer program for the simulation of Mie scattering in case of arbitrarily large size parameters. The elements of the scattering matrix, efficiency factors as well as the corresponding cross sections, the albedo and the scattering asymmetry parameter are calculated. Single particles as well as particle ensembles consisting of several components and particle size distributions can be considered.

We present a database containing information on different aspects of calculation and usage of the optical properties of small non-spherical particles -- cosmic dust analogues. The main parts of the DOP are a review of available methods of the light scattering theory, collection of light scattering codes, special computational tools, a graphic library of the optical properties (efficiencies, albedo, asymmetry factors, etc.), a database of the optical constants for astronomy, a bibliographic database of light scattering works, and links to related Internet resources. The general purpose of the DOP having the address http://www.astro.spbu.ru/DOP is twofold -- to help scientists to apply the light scattering theory in astronomy and to give students and beginners a possibility quickly to get necessary knowledge on the subject.

Dust grains have been detected in various astronomical objects. Interpretation of observations of dusty objects includes three components: 1) determination of the materials which can exist in the solid phase and the measurements or acquisition of their optical constants; 2) selection of the light scattering theory in order to convert the optical constants into the optical properties of particles; and 3) the proper choice of the object model which includes, in particular, the correct treatment of the radiative transfer effects. The current state of the components of dust modeling and the reliability of information obtained on the cosmic dust from transmitted, scattered and emitted radiation are discussed.

The approach to model composite interstellar dust grains, using the exact solution to the light scattering problem for multi-layered spheres as suggested by Voshchinnikov & Mathis (1999), is further developed. Heterogeneous scatteres are represented by particles with very large number of shells, each including a homogeneous layer per material considered (amorphous carbon, astronomical silicate and vacuum). The applicability of the effective medium theory (EMT) mostly utilized earlier to approximate inhomogeneous interstellar grains is examined on the basis of the new model. It is shown that the EMT rules generally have an accuracy of several percent in the whole range of particle sizes provided the porosity does not exceed about 50%. For larger porosity, the EMT rules give wrong results. Using the model, we reanalyze various basic features of cosmic dust -- interstellar extinction, scattered radiation, infrared radiation, radiation pressure, etc. As an example of the potential of the model, it is applied to reproduce the extinction curves in the directions to $\zeta$ Oph and $\sigma$ Sco using subsolar cosmic abundances. We also conclude that metallic iron even in negligible amount ($\la 1$ % by the volume fraction) is unlikely to form a layer on or inside a grain because of peculiar absorption of radiation by such particles.

We present the first numerical radiative transfer simulation of multiple light scattering in dust configurations containing aligned non-spherical (spheroidal) dust grains. Such models are especially important if one wants to explain the circular polarization of light, observed in a variety of astronomical objects. The radiative transfer problem is solved on the basis of the Monte Carlo method. Test simulations, confirming the correct numerical implementation of the scattering mechanism, are presented. As a first application, we investigate the linear and circular polarization of light coming from a spherical circumstellar shell. This shell contains perfectly aligned prolate or oblate spheroidal grains. The most remarkable features of the simulated linear polarization maps are so-called polarization null points where the reversal of polarization occurs. They appear in the case when the grain alignment axis is perpendicular to the line of sight. The maps of circular polarization have a sector-like structure with maxima at the ends of lines inclined to the grain alignment axis by $\pm 45\degr$.

Using the method of separation of variables and a new approach to calculations of the prolate spheroidal wave functions, we study the optical properties of very elongated (cigar-like) spheroidal particles. A comparison of extinction efficiency factors of prolate spheroids and infinitely long circular cylinders is made. For the normal and oblique incidence of radiation, the efficiency factors for spheroids converge to some limiting values with an increasing aspect ratio a/b provided particles of the same thickness are considered. These values are close to, but do not coincide with the factors for infinite cylinders. The relative difference between factors for infinite cylinders and elongated spheroids ($a/b \ga 5$) usually does not exceed 20 % if the following approximate relation between the angle of incidence $\alpha (in degrees)$ and the particle refractive index m=n+ki takes the place: $\alpha \ga 50 |m-1| + 5$ where $1.2 \la n \la 2.0$ and $k \la 0.1$. We show that the quasistatic approximation can be well used for very elongated optically soft spheroids of large sizes.

The albedo and dichroic polarization of spheroidal particles are studied. For the description of the phase function two asymmetry parameters $g_{||}$ and $g_{\bot}$ characterizing the anisotropy in forward/backward and left/right directions are introduced and calculated. The consideration is based on the solution to the light scattering problem by the Separation of Variables Method (Voshchinnikov and Farafonov, 1993).

The temperatures of prolate and oblate spheroidal dust grains in the envelopes of stars of various spectral types are calculated. Homogeneous particles with aspect ratios {\small $a/b \le 10$} composed of amorphous carbon, iron, dirty ice, various silicates, and other materials are considered. The temperatures of spherical and spheroidal particles were found to vary similarly with particle size, distance to the star, and stellar temperature. The temperature ratio {\small $T_{\rm d}({\rm spheroid})/T_{\rm d}({\rm sphere}) $} depends most strongly on the grain chemical composition and shape. Spheroidal grains are generally colder than spherical particles of the same volume; only iron spheroids can be slightly hotter than iron spheres. At {\small $a/b \approx 2$}, the temperature differences do not exceed 10 %. If {\small $a/b \ga 4$}, the temperatures can differ by 30-40 %. For a fixed dust mass in the medium, the fluxes at wavelengths {\small $\lambda \ga 100 \mkm$} are higher if the grains are nonspherical, which gives over estimated dust masses from millimeter observations. The effect of grain shape should also be taken into account when modeling Galactic-dust emission properties, which are calculated when searching for fluctuations of the cosmic microwave background radiation in its Wien wing.

Observational data for dusty comets are summarized and systemized. The dependence of the linear polarization on various parameters (the phase angle, wavelength, etc.) is analyzed. An ensemble of aligned spheroidal particles having various sizes and chemical compositions is used to fit the linear and circular polarization observed for the comet Halley.

We model the positional dependence of the optical polarization (BVRI-bands) in the circumstellar disc of beta Pictoris as observed by Gledhill et al. (1991) and Wolstencroft et al. (1995). The particles are spherical, have a size distribution n(a) \sim a^{-q} and their number density decreases with distance from the star as r^{-s}. We consider both compact and porous grains. Varying the grain size and the exponent, s, of the density distribution as well as the refractive index of the grain material m, we find that the measured polarization, colours and brightness distribution in the disc are best reproduced by a model in which the grains are larger than interstellar grains (the minimum grain size a_min=0.15 mkm. The value of the maximum size is ill determined because it has little influence and was taken to be a_max=100 mkm. The best-fit exponents of the power laws are q=3.2 and s=3. The grains have in the R-band a refractive index m_R=1.152-0.005i. Such a value is roughly appropriate for porous grains where half of the volume is ice and the other half vacuum, or where 24% of the volume consists of silicate and the remaining 76% of vacuum.

Interstellar grains may be composite collections of particles of distinct materials, including voids, agglomerated together. We determine the various optical cross sections of such composite grains, given the optical properties of each constituent, using an approximate model of the composite grain. We assume it consists of many concentric spherical layers of the various materials, each with a specified volume fraction. In such a case the usual Mie theory can be generalized and the extinction, scattering, and other cross sections determined exactly. We find that the ordering of the materials in the layering makes some difference to the derived cross sections, but averaging over the various permutations of the order of the materials provides rapid convergence as the number of shells (each of which is filled by all of the materials proportionately to their volume fractions) is increased. Three shells, each with one layer of a particular constituent material, give a very satisfactory estimate of the average cross section produced by larger numbers of shells. We give the formulae for the Rayleigh limit (small size parameter) for multi-layered spheres and use it to propose an ``Effective Medium Theory'' (EMT), in which an average optical constant is taken to represent the ensemble of materials. Multi-layered models are used to compare the accuracies of several EMTs already in the literature.

We use the separation of variables and T-matrix methods to calculate the optical properties of homogeneous spheroids with refractive indices from m = 1.3+0.0i up to 3+4i. It is found that the extinction cross-sections for highly absorbing spheroids are normally 1.5 -- 2 times larger than those for spheres of the same volume. The albedo of the non-spherical particles rather slightly depends on the particle shape and is mainly determined by the imaginary part of the refractive index. Beginning at some size, the spheroidal particles do not polarize the transmitted radiation independent of their shape. We also suggest a new approach for axisymmetric particles which would combine the strong aspects of both methods mentioned above and give several values of the cross-sections as benchmarks in tabular form.

A model of spheroidal particles is used to calculate the steady-state temperature of dust grains immersed in the interstellar radiation field. It is found that the temperature of non-spherical grains with the aspect ratios a/b <= 2 deviates from that of spheres less than 10%. More elongated or flattened particles are usually cooler than spheres of the same mass and in some cases the temperatures may differ by even about a factor of 2. The shape effects increase with the infrared absorptivity of the grain material and seem to be more important in dark interstellar clouds.