A&A 428, 99-108 (2004)
Research School of Astronomy & Astrophysics, Institute of Advanced Studies, The Australian National University, Cotter Road, Weston Creek, ACT 2611, Australia
Received 31 July 2003 / Accepted 28 July 2004
The optical properties of small spherical iron grains are derived using a Kramers-Kronig-consistent model of the dielectric function including its dependence on temperature and size. Especially discussed is the effect of the size dependence, which results from the limitation of the free path of the free electrons in the metal by the size of the grain, on the absorption behaviour of small iron spheres and spheroids. The estimated absorption properties are applied to study the temperature behaviour of spherical and spheroidal grains which are heated by the interstellar radiation field.
Key words: ISM: dust, extinction - ISM: general
According to the cosmic abundances given by Anders & Grevesse (1989), iron is the ninth most abundant element in the interstellar medium (ISM) and should therefore be an essential constituent of interstellar dust grains. Iron particles are expected to form independently of the oxygen-to-carbon ratio in circumstellar shells (see Jones 1990, and references therein) and have been supposed to be, apart from silicates, the main condensates in the metal-rich core region of supernova ejecta (Hoyle & Wickramasinghe 1970). SN 1987A was the first supernova where the condensation of new grains in ejected material could be observed. The composition of these grains could not be clarified but indications, summarised by Wooden (1997), point to predominantly iron-rich grains.
Even though the ISM will probably be enriched in iron particles it is uncertain whether they can survive and play an important role as a dust component of interstellar grains (Duley 1980; Jones 1990). If they are present, iron particles may affect the interstellar extinction (Wickramasinghe & Nandy 1971), the chemistry in the ISM due to catalytic processes on grain surfaces (Tabak 1978), and may contribute to the diffuse emission from interstellar dust grains (Chlewicki & Laureijs 1988). Further, iron grains may be an important grain species to explain the Rosseland mean opacity in molecular clouds cores (above 575 K) and accretion disks (above 425 K) (Pollack et al. 1994). On the other hand it has been claimed, based on observations at 90 GHz, that not more than 5% of interstellar iron is in form of metallic iron grains or inclusions (Draine & Lazarian 1999). This would support the hypothesis that the survival time of a relatively pure iron grain is quite short. However, even though overall iron may be not an important constituent of interstellar grains it does not exclude the possibility of local enhancements where these grains may form.
Mathis et al. (1977) found, from analysis of the extinction curve, that the interstellar medium is mainly composed of grains much smaller than 0.25 m. Very small grains are furthermore thought to be responsible for the diffuse IR emission at wavelengths shorter than 60 m (see e.g. Draine & Anderson 1985) due to the stochastic heating process, where very small grains reach temperatures much higher than the corresponding equilibrium temperature. To model their emission or their optical properties, it is usually assumed that the dielectric function derived for solids can also be applied for very small grains. At least for small metallic grains, as shown by experiments with very small silver spheres, it is known that the dielectric function has to be corrected for grain size (Kreibig & Fragstein 1969; Kreibig 1974). The influence of the grain size on the dielectric function can be even stronger than its dependence on temperature. One expects that the optical properties of small metallic spheres should be much less affected by temperature than those of larger grains as shown, for example, by calculations of the absorption behaviour of spherical graphite grains (Draine & Lee 1984).
Here, we present calculations for the temperature behaviour of pure iron grains in the ISM heated by the interstellar radiation field (ISRF) where we have taken into account both the temperature and the size dependence of the dielectric function. It has been suggested that a part of the interstellar grains may have a non-spherical shape (see Voshchinnikov et al. 1999). Therefore, apart from simple iron spheres, spheroidal iron grains will also be considered. The influence of the shape of small grains of different compositions, including iron, on temperature has also been discussed by Voshchinnokov et al. (1999). Their results are based on the dielectric function of solids using a sophisticated code to derive the optical properties of spheroidal grains. Here a much simpler approach is used demonstrating the effect of the size dependence of the dielectric function on the cooling behaviour of spheroidal iron grains.
In Sect. 2 we describe the model used to derive the optical properties of small iron spheres and analyse in more detail the absorption behaviour of small spherical and spheroidal iron grains. In Sect. 3 the optical properties are used to derive the temperature behaviour of spherical and spheroidal small iron grains in the ISM which are heated by the ISRF. The derived temperature distributions are used to obtain the emission spectrum of spherical iron grains assuming a grain size distribution as proposed for the ISM including the stochastic heating process of very small grains. The results are discussed in Sect. 4 and summarised in Sect. 5.
The interaction of an electromagnetic wave of a given energy with a grain is generally determined by the response of the corresponding solid, expressed by the complex dielectric function , the shape, and the size of the grain.
As for non-metallic grains the magnetic permeability is taken to be . As will also be shown this assumption should be sufficient to describe the cooling of iron grains in the ISM, which is the main purpose of this paper. The contribution of the magnetizability of metallic grains to the optical properties has been studied by Draine & Lazarian (1999) and is only significant at wavelengths larger than 1000 m.
To derive the optical properties of spherical iron grains we used a self consistent simple model for the complex dielectric function based on optical data derived for bulk iron at a temperature of .
We followed basically the method used by Draine & Lee (1984) to determine the optical properties of small graphite grains and separated the dielectric function in
which correspond to the free and the bound electrons of the grain:
The parameters of the Drude term of bulk iron are determined at long wavelengths where the free electron gas dominates the optical behaviour of metals. To derive a self consistent dielectric function we have chosen the imaginary part as given by and used the Kramers-Kronig relation (Landau & Lifschitz 1985) to estimate the real part .
The main dependence of on size and temperature is given by the collision rate which can be taken as the sum of the collision rates and due to scattering events in the solid and at the outer surface of the grain (Kreibig & Fragstein 1969; Kreibig 1974). A further temperature dependence is assumed to be caused by the thermal expansion of the solid which alters the electron density and therefore ( ).
Inside the solid, electrons are scattered at impurities, defects and phonons, but only the phonon scattering events lead to a strong temperature dependence. In the case of bulk iron, the scattering at impurities and defects is negligible with the exception of very low temperatures where they are the dominant scattering process inside the grain.
As a first approximation the scattering rate at phonons may be assumed to be independent of the amount of impurities and defects and the scattering rate at impurities and defects independent on temperature so that the effects can be easily de-coupled. Under these assumptions the collision rate due to impurities and defects can be scaled by a factor to allow consideration of the cases where the iron grains are either less pure or have more defects in comparison to pure iron grains.
The total collision rate is given by
|Figure 1: Critical grain sizes of iron spheres where the collision rate of the free electrons in the material is equal to the collision rate at the surface shown for , 10, and 100. Below the curves the collision rate is dominated by surface scattering.|
The collision rate at the surface becomes important when the free path of the electrons is greater than the size of the grain. Where the collision rate in the solid is equal to the collision rate at the surface is shown in Fig. 1 for different values . In the case of pure very cold iron grains () the scattering should be dominated by surface scattering for all grains with radii smaller than 20 m. If the grains are less pure and have more defects the effect is limited to smaller grain sizes.
As the purity of "astrophysical'' iron grains is not known we considered pure iron grains. The collision rate assumed in this paper can therefore be regarded as a lower limit of the collision rate of actual iron grains. However, this assumption should not affect interstellar iron grains which, if they do exist, are believed to be quite small with maximum sizes not much larger than 0.1 m. As seen in Fig. 1, their mean free path should be determined by the size of the grain unless the grains are highly disordered or impure with . The effect of the purity and defects on the temperature behaviour of larger iron grains will be discussed in Sect. 4.1.1.
The optical data used to derive the optical properties can be taken from Table 1. We interpolated data from Lynch & Hunter (1996) and Moravec et al. (1976) in the energy range 51.57 eV to 26 eV using . A smooth transition was made in the range 3.5 to 5 eV between data from Moravec et al. and the data published in the "CRC Handbook for Chemistry and Physics'' (1998). Above 0.8 eV we used the averaged data from Ordal et al. (1988).
|Energy intervall [eV]||Reference|
|10 000-51.57||Lynch & Hunter (1996)|
|51.57-26.0||interpolation (see text)|
|26.-5||Moravec et al. (1976)|
|1.55-0.0124||Ordal et al. (1988)|
|0.0124-0.00124||extrapolation (see text)|
From the optical data in the range of 1 to we determined the plasma frequency and the internal collision rate of bulk iron at room temperature which are given in Table 2. To prevent amplification of the electromagnetic wave inside the solid we have taken care that in the considered energy range. The derived parameters of the Drude term ( , ) are slightly smaller than the ones used by Ordal et al. for extrapolation above ( and ) but very close to those used by Draine & Lazarian (1999).
The DC conductivity of iron at room temperature is found to be which is only about a factor of 1.2 lower than the measured value of (CRC 1998). For comparison in the case of graphite the DC conductivity given by the parameters of the Drude term is a factor of 3 smaller than the measured value (Draine & Lee 1984).
|Plasma frequency (298 K)a||[s-1]|
|Internal collision rate (298 K)a||[s-1]|
The derived dielectric function is shown in Figs. 2 and 3.
At long wavelengths (Fig. 2) the functional form of
is qualitatively similar to the dielectric function
of graphite perpendicular
to the main axis (Draine & Lee 1984), which shows metal-like behaviour.
The imaginary part
with wavelength and reaches its highest values at 60 m. Beyond
assumed that the values of
decrease with a power law.
The explanation of the behaviour of the dielectric function
in the IR
is out of the scope of this paper and it cannot be excluded that this behaviour is
a result of the simplification in modelling the free electron regime with a single Drude term.
|Figure 2: Real and imaginary part of the dielectric function for individual pure small iron spheres at energies above 1 eV. For comparison also the real and imaginary parts of are shown as thin solid lines. The part of shown as thin dotted line are extrapolated values used for the Kramers-Kronig relation to derive . The real part of the dielectric function is always negative except for the three branches shown. The short and the long dashed thin lines in the right hand figure refer to the dielectric function of an iron sphere with and , respectively.|
|Figure 3: Real and imaginary part, given as conductivity , of the dielectric function above 1eV at a temperature of 298 K. The chosen values of the imaginary part and the result for the real part after Kramers-Kronig analysis are shown as solid lines. The values differ somewhat from the published values. The derived values of the real part of the dielectric function lie between those from Johnson & Cristy (1974) and those from Moravec et al. (1976). The real part (shown as dashed-dotted line) is the direct result derived with the Kramers-Kronig relation from the imaginary part .|
The temperature dependent DC conductivities of bulk iron are taken from the "CRC Handbook of Chemistry and Physics'' (1998). The expansion coefficients of iron for correcting the plasma frequency with temperature was taken from the "American Institute of Physics Handbook'' (1972). For the collision rate at the surface we used for the free electron gas of iron a Fermi velocity (Ashcroft & Mermin 1976).
We estimated the optical properties of spheres above using the MIE-formalism from Bohren & Huffman (1983). At longer wavelengths we used the dipole approximation for spheres presented by Draine & Lee (1984). The results for iron grains at a temperature of 25 K are shown in Fig. 4. At this temperature the mean free path of the electrons would be larger than the grains so that the collision rate is given by the collision rate at the grain surface and therefore through . While scattering ( , ) is almost insensitive to the grain temperature, the absorption of larger iron grains can have large temperature variations.
|Figure 4: Derived optical properties of spherical iron grains. Shown are the absorption coefficients , scattering coefficients , and the scattering parameter at 25 K for various grain radii from to with a difference in radius of . Thick lines are labeled with the corresponding grain radii in microns.|
The dependence of the absorption behaviour of iron grains on temperature is shown in Fig. 5 for two different grain sizes with radii of 0.01 and . For comparison we also show the absorption cross section due to magnetisation effects inside iron spheres well below the Curie temperature as derived by Draine & Lazarian (1999). If correct this effect clearly dominates the absorption of cold iron grains at wavelengths longer than 1000 m. Above the Curie temperature the iron spheres become non magnetic and the absorption at long wavelength should be determined by the electronic response of the material. As the magnetisation does not affect the cooling behaviour of the grains this effect will be ignored in the following.
Due to the high Fermi velocity only the optical properties of larger iron grains with are affected by grain temperature. Even a variation of grain temperature over a huge range (10-1200 K) causes only a slight change in the absorption behaviour of grains with . Grains with a radius of on the other hand have, at a wavelength of , an order of 10 higher absorption at 300 K than at 10 K. In this wavelength region the absorption coefficients increase with temperature. At longer wavelengths it is possible that colder iron grains emit more readily than warm grains do.
|Figure 5: Effect of the temperature on the absorption behaviour of spherical iron grains. For various temperatures, the absorption coefficients of iron grains with radii 1.0 and for (thick lines) are shown. For comparison the absorption coefficients due to magnetisation effects for magnetic single and multi-domain iron grains well below the Curie temperature ( ) are given (thin dashed and solid line).|
However, as seen in Figs. 4 and 5, inside certain regions of the parameter space defined by the wavelength and the grain size a the absorption coefficient can still be described by a simple approximation with fixed values for s and t. The absorption properties at long wavelengths can be explained using the electric and the magnetic dipole approximation for the cross sections and are discussed in some detail in the Appendices A and B.
As a special case of non-spherical iron particles we will consider spheroidal iron grains. In adopting the model of the dielectric function described in Sect. 2.1 to
non-spherical grains one has to choose the appropriate value of the collision rate in Eq. (3). As an approximation of the mean free path due to surface scattering we have taken the averaged distance from the center to the surface:
In general for the calculation of the optical properties of non-spherical grains, even for spheroidal grains, this can be very complicated (see e.g. Voshchinnikov & Farafonov 1993). Here we are interested in the temperature behaviour of grains in the interstellar medium where the radiation is almost isotropic and unpolarised. In such cases it is sufficient to use only the cross section averaged over the orientations of the spheroids.
|Figure 6: Angle-averaged absorption cross section of prolate (a>b) and oblate (a<b) iron grains with volumes equivalent to the volume of a sphere with radius . For comparison also the absorption of a sphere is given. Solid and dashed lines correspond to absorption where the size and temperature dependence of the dielectric function has been taken into account and has been ignored, respectively. The dotted line shows the electric dipole absorption of a sphere for the case of temperature and size independent dielectric function.|
For simplicity, we considered spheroidal grains small in comparison to the wavelength of the absorbed photon so that the optical properties can be described by the dipole approximation. For the electric dipole absorption we used the approximation for ellipsoidal grains and for the magnetic dipole absorption we adopted the expression given for spheres (Draine & Lee 1984).
The angle averaged absorption cross sections of spheroidal iron grains for a number of different ratios a/b are shown in Fig. 6 (solid lines). The volume of the spheroids is chosen to be that of a sphere with a radius of . For comparison also the cross sections for temperature and size independent dielectric function (dashed lines) is given.
The absorption behaviour at long wavelengths of a spherical iron grain with is dominated by the electric dipole absorption if the size dependence is taken into account. This is also true for the spheroids as the contribution of the electric dipole absorption only increases with axis ratio while that of the magnetic dipole absorption is assumed to be constant.
As can be seen in the figure, by considering the size dependence of the dielectric function, the absorption at long wavelengths is clearly enhanced in comparison to the calculation where this dependence is ignored. On the other hand a resonance absorption in the IR due to surface modes (best seen for a prolate grain with a/b=100) is reduced.
The absorption properties at long wavelengths is qualitatively different if the magnetic dipole absorption becomes important, which is the case of larger iron grains. The behaviour should be similar to the case shown in Fig. 6 where the size and temperature dependence of the dielectric function is ignored. As can be seen in the figure, at long wavelengths the absorption behaviour of the sphere is then dominated by the magnetic dipole absorption. Because of the assumption of magnetic dipole absorption, the absorption efficiency at those wavelengths is almost insensitive to the axis ratio a/b unless the electric dipole absorption becomes important.
The absorption behaviour in the optical on the other hand increases with increasing axis ratio as long as the ratio is small. Therefore, slightly prolate or oblate small iron grains with a volume corresponding to a sphere with can absorb optical light more efficiently than spheres but show a similar cooling behaviour. If the grains are heated by the ISRF for example this results in an increase of the grain temperature with increasing axis ratio (Sect. 3.1.2).
As an example we considered spherical and small spheroidal iron grains which are heated by the ISRF. For the radiation field we have taken the energy density in the solar neighbourbood as derived by Mathis et al. (MMP 1982), scaled by a factor :
First we ignore the temperature variation in the case of very small grains and characterise all iron grains with single temperatures
where the heating is equal to the cooling of the grain:
For the photon density we assumed three different strengths and . For comparison we also derived the equilibrium temperatures for spheres where we assumed a dielectric function which is independent of temperature and size. The derived equilibrium temperatures of iron grains with radii ranging from 10 Å to 10 are shown in Fig. 7.
As mentioned by Chlewicki & Laureijs (1988), the equilibrium temperature of iron grains in the ISM is a strong function of grain size. Where the magnetic dipole effect becomes important at the temperature falls from relatively high to low values which are close to temperatures of silicate or graphite grains. The strong temperature variation with grain size is also expected for collisionally heated iron spheres (see Fischera et al. 2002).
|Figure 7: Equilibrium temperatures of spherical iron grains in the ISM for three different radiation strengths (, and ). Dashed and dotted lines correspond to calculations in which we included and ignored the temperature and size dependence of the dielectric function.|
The grain temperatures are so low that the collision rate of the free electrons is dominated by surface scattering (see Fig. 1). As a consequence of the absorption behaviour in the IR (Sect. 2.2.1) the cooling of the smallest grains increases while the cooling of larger grains becomes less efficient. Therefore, if the size dependence is taken into account, the temperature of smallest iron spheres is lower and the temperature of larger iron spheres higher compared with the case where the size dependence of the dielectric function has been ignored. We mention that the cooling of the smallest grains could be even stronger if the dielectric function in the IR were lower than assumed here.
For comparison we considered spheroids with a radius and . The derived temperatures of these iron grains are shown in Fig. 8. We also give the equilibrium temperatures of iron spheroids for a dielectric function ignoring the temperature or size dependence. The temperatures for extremely elongated or flattened iron grains, in particular for the grain with , should be taken as approximate values because of the fact that in this case the electric dipole approximation in the UV gives only a poor representation of the actual absorption behaviour.
Taking the size and temperature dependence of the dielectric function into account, the two sizes considered show qualitatively a very different temperature behaviour. We found that the temperature of a spheroid with a radius becomes smaller with an increasing axis ratio. This behaviour is therefore the same found for spheroids of other grain compositions (Voshchinnikov et al. 1999). On the other hand the temperature of an iron grain with increases first before the temperature decreases for larger axis ratios. This is expected from the discussion at the end of Sect. 2.3.
|Figure 8: Approximate equilibrium temperatures of spheroidal iron grains with a volume equal to a volume of a sphere with radius and as function of the ratio a/b. The solid and the dashed line correspond to calculations where the dependence of the dielectric function on grain size and temperature has been taken into account and has been ignored, respectively.|
At extreme values of the ratio a/b the equilibrium temperatures of spheroids with are higher than for spheroids with . This can be explained by the higher collision rate due to surface scattering in the smaller grain leading to a higher electric dipole absorption.
By including the size and temperature dependence of the dielectric function the temperatures of spheroids with are up to 10 K lower than the temperatures where this dependence is ignored. If the dependences are not taken into account the temperature behaviour would qualitatively be similar to the larger grains showing an increase of the temperature for intermediate values of the ratio a/b first and a decreasing temperature for more oblate and prolate shapes.
|Figure 9: Temperature distribution of spherical iron grains with radius 0.001 and heated by the ISRF (MMP). Solid and dashed lines correspond to calculations in which we included and ignored the dependence of the dielectric function on temperature and grain size.|
The temperature distributions of iron grains with radius a=10 Å and a=40 Å heated by the ISRF () are shown in Fig. 9 as . For comparison we also included in this figure the theoretical temperature variation of small spheres with a dielectric function of bulk iron at room temperature. If the size dependence of the collision rate is taken into account the probabilities of high temperatures are partly strongly reduced while the probabilities of low temperatures is clearly increased. Iron spheres with a=10 Å for example should have temperatures lower than 10 K most of the time. In contrast a model where the size dependence of the dielectric function is ignored predicts for the same grain size temperatures higher than 10 K with a most probable temperature of 30 K.
The maximum temperature very small grains reach is roughly the temperature where the thermal energy is equal to the maximum photon energy of 13.6 eV. Because of the high heat capacity of iron grains the maximum temperature is clearly lower in comparison with other grain species like silicate or graphite grains. By using the heat capacity given by Dwek (1986) graphite grains with a radius of 10 Å attain a maximum temperature of above 400 K (Fischera 2000). In contrast the maximum temperature of iron grains with the same size is less than 300 K.
From the analysis of the extinction curve Mathis et al. (1977) proposed a power law with for the grain size distribution in the ISM. The minimum and maximum grain size as the power of the size distribution depends somewhat on the assumed grain composition. For the iron grains we considered a maximum grain size of . To study the effect of stochastic heated grains we assumed as minimum grain size different radii 10, 20 and 40 Å.
The mean intensity emitted at the surface of a flickering grain with radius a is given by
|Figure 10: Emission spectra of individual small iron grains. The sizes are chosen to be 10, 20, and 40 Å. The spectra shown as solid lines correspond to calculations in which the size and temperature dependence of the dielectric function has been taken into account. The spectra shown as dashed lines correspond to calculations in which this dependence has been ignored.|
The total flux density of a distribution of spherical iron grains at distance D is given by:
The SED of iron grains of at 1 kpc is shown in Fig. 11. For comparison we give the SEDs for the calculation including and ignoring the size and temperature dependence of the dielectric function. Due to the large variation of the equilibrium temperatures the SED is relatively broad. In comparison to graphite or silicate grains the SED is less effected by flickering grains. A clear increase at the blue side of the spectrum requires grains with radii smaller than Å. The inclusion of the size and temperature dependence of the dielectric function reduces the emission at the blue side which is contributed by stochastically heated small grains. The peak of the emission on the other hand is shifted towards shorter wavelengths. At long wavelengths we obtained for the two different cases of the dielectric function almost the same emission spectrum.
|Figure 11: Spectral energy distribution (SEDs) of spherical iron grains with a mass of at a distance of 1kpc heated by the ISRF (MMP). The grains are assumed to have a power law with a maximum grain size and different minimum grain sizes of 10, 20 and . SEDs in which calculation the temperature and size dependence of the dielectric function was included are shown as solid lines. Dashed lines correspond to SEDs in which this dependence has been ignored.|
The study of the temperature behaviour of metallic grains in the ISM requires a realistic model of their optical properties. The grains emit mainly in the IR, in which range the absorption is strongly affected by the free electrons. This leads to an absorption that depends on temperature and, in cases where the free path of the electrons is larger than the grain size, especially on the grain size.
A simplification can easily cause bigger uncertainties in the equilibrium temperatures of spherical and spheroidal grains or in the probability distribution of temperatures from flickering small grains.
To derive the cooling behaviour of iron grains Tabak and Straitiff (1983) and Tabak (1987) only considered the temperature but no size dependence of the optical properties. In addition it was assumed that the grain absorption is sufficiently described by the electric dipole approximation. Both simplifications lead to absorption coefficients which are too low. Therefore, equilibrium temperatures in the range of 100 K for iron grains heated by the ISRF as presented by Tabak & Straitiff (1983) and also Tabak (1987) should be in general too high.
The temperatures of iron spheroids with as given by Voshchinnicov et al. (1999) show qualitatively the peculiar behaviour we derived only for larger spheriods where the absorption behaviour at long wavelengths of the corresponding sphere is dominated by the magnetic dipole absorption. The peculiar curve of spheroidal grains with should disappear if the size dependence of the dielectric function is taken into account.
We found by considering iron spheres heated by the ISRF that the grain temperature as a function of radius shows a minimum at 0.3 m and that the temperature increases towards even larger iron grains. This is caused by the fact, as is shown in Sects. A and B, that the absorption behaviour of big iron grains is dominated by the magnetic dipole absorption (Eq. (A.7)) which decreases towards larger grains if the scattering rate is dominated by scattering events at the grain surface and therefore given by . In the case of less pure and disordered iron grains the collision rate will be dominated by scattering events at impurities and defects so that the absorption coefficient becomes essentially constant with grain size. It is therefore possible that the increase in grain temperature will disappear if iron grains in astrophysical environments are considered.
In this case, as can be seen in Fig. 2, the IR-values of the dielectric function of very small grains could be lower and their absorption efficiency therefore higher than assumed here. If the part of the bounded electrons would be negligible so that the dielectric function is determined by the free electron part, the absorption or cooling efficiency of an iron grain with could increase by almost one order of magnitude. It is therefore possible that the probabilities at high temperatures of flickering small iron grains is even smaller than the probabilities presented here.
The optical properties of small spherical and spheroidal iron grains were derived using a self consistent model of the dielectric function including its dependence on size and temperature. We have shown that the limitation of the free path of the electrons due to the size of the grain has a non-negligible effect on the absorption behaviour of small spheres and spheroids. We applied the estimated absorption coefficients for small spherical and spheroidal iron grains in the ISM which are heated by the ISRF and derived the emission spectrum of stochastically heated small iron spheres assuming as grain size distribution a power law as proposed for grains in the ISM and chosing the radius of the largest grain as . The main results for the temperature behaviour of small iron grains are:
The temperature behaviour of larger iron grains is likely to be effected by impurities and defects. If these grains are relatively pure the temperature should increase with size. Otherwise their temperature should be size independent.
By considering the size dependence of the dielectric function we found in comparison to calculations where this dependence is ignored:
I am thankful for the support provided by the Max-Planck-Institut für Kernphysik in Heidelberg and the Research School of Astronomy & Astrophysics of the Australian National University in Canberra. I acknowledge the financial support through the ARC Discovery project DP0208445. I would like to thank Prof. Krätschmer for fruitful discussions and Prof. Dopita, Roberto de Propris, and Helmut Jerjen for useful comments that improved the manuscript.
Here we will discuss in some detail the absorption behaviour of metallic grains using the electric and magnetic dipole approximation. To derive approximations it is assumed that the dielectric function is dominated by the free electron gas and therefore by
The cross section at wavelengths with
is given by (see e.g. Draine & Lee 1984):
If the contribution of the magnetic dipole absorption is negligible the absorption coefficient is known to be given by:
With increasing grain size the magnetic dipole absorption begins to dominate the absorption behaviour at long wavelengths. In the region
the absorption coefficient due to induced eddy currents in the metal spheres can be simplified to the expression:
In the region
the absorption can be described by:
The plateau with constant absorption coefficients at long wavelengths is caused by the magnetic dipole absorption in the region
In this range the absorption coefficients are approximately given by:
For completeness we also give the absorption behaviour for the region
|Figure B.1: Absorption coefficients of metallic iron spheres derived approximating the dielectric function through the free electron part . In the left hand figure the collision rate of the free electrons is taken to be the same as for bulk iron at 298 K. In the right hand figure it is assumed that the collision rate is dominated by scattering events at the grain surface. The absorption coefficients are shown as contours in logarithmic scale (thin solid lines) where the difference between the lines is taken to be . The labels correspond to . The hatched region shows where the absorption behaviour deviates strongly from the dipole absorption ( ). The thick solid line visualises where the electric dipole absorption is equal to the magnetic dipole absorption. Where the frequency of the photons is equal to the collision rate of the free electrons in the grain is shown as long dashed line. The short dashed and dashed dotted line correspond to absorption coefficients where |y|=1 and |y|=3, respectively.|
To visualise the effect of the limitation of the mean free path lengths of the free electrons on the cooling behaviour of iron spheres it might be useful to approximate the dielectric function by a single Drude term using the parameters as given in Table 2. The derived absorption coefficients are shown in Fig. B.1 as function of wavelength and grain size. For comparison the absorption coefficients have been calculated for two different assumptions for the collision rate . In the left hand figure no correction to the measured dielectric function (at room temperature) has been made. In the right hand figure it is assumed that the collision rate is dominated by surface scattering which may be appropriate to describe the situation of cold pure iron grains in the ISM which are heated by the ISRF.
If no correction is applied, the magnetic dipole absorption would be the dominant absorption at long wavelengths for grains with radii larger than 3 nm. The limitation of the free path of the free electrons by the grain size shifts this value to grains with radii larger than 0.015 m. The absorption increases almost for all wavelengths with grain size when no size dependence of the dielectric function has been considered (Fig. B.1, left hand side). When this dependence is taken into account, the absorption in the IR shows a maximum for grains with a radius in the range of 0.3 m (Fig. B.1, right hand side), which is the reason for the lowest temperatures of these grains in the ISM as derived in Sect. 3.1.1. The decrease of the absorption with grain size for grains larger than extends also into the region where the absorption has to be described using the Mie-formalism (hatched region).
The approximations of the magnetic dipole absorption given in Sect. 2.2.1 are valid in the regions left of and right of (shown in Fig. B.1 as short dashed and dashed dotted lines) considering and .
The absorption behaviour is shown for non-magnetic iron grains to demonstrate the different trends. As stated in the text, above the absorption of iron grains might be dominated by magnetisation effects (Fig. 5). At short wavelengths on the other hand the absorption behaviour of the iron grains additionally may be effected by the response of the bound electrons.