Thermal and photophoretic properties of dust mantled chondrules and sorting in the solar nebula
Faculty of physics, University of DuisburgEssen, Lotharstr. 1, 47057 Duisburg, Germany
email: christoph.loesche@unidue.de
Received: 8 February 2012
Accepted: 23 June 2012
Many chondrules are found to be surrounded by a fine grained rim. Supposedly, these rims were acquired from dust accreted to the chondrules in a protoplanetary disk. In numerical simulations we study the heat transfer in illuminated bare and dust mantled chondrules. The calculations consider the chondrule size, the dust mantle size, and the thermal conductivities of both components as parameters. We calculate the photophoretic force and compare the numerical results to analytical approximations. We give an expression to quantify the photophoretic force on a spherical particle in the free molecular regime to better than 2%. We describe the influence of a dust mantle on the photophoretic strength by an effective thermal conductivity of the coremantle particle. The effective thermal conductivity significantly depends on the size ratio between mantle and chondrule but not on the absolute sizes. It also strongly depends on the thermal conductivity of the mantle with minor influence of the thermal conductivity of the chondrule. The size ratio between rim and chondrule in meteorites is found to vary systematically with overall size by other authors. Based on this, our calculations show that a photophoretic size sorting can occur for dust mantled chondrules in optically thin disks or at the evolving inner edge of the solar nebula.
Key words: acceleration of particles / conduction / radiation mechanisms: thermal / radiative transfer / methods: numerical / protoplanetary disks
© ESO, 2012
1. Introduction
The work behind this paper is based on two observations concerning the mmsize particles, called chondrules, which are a major component in primitive meteorites.

Chondrules often have specific but different size in different meteorites (Kuebleret al. 1999; Scottet al. 1996; Hughes 1978).

Chondrules are often surrounded by fine grained rims (Metzler et al. 1991).
The origin of these facts is still debated. The rims might originate in the solar nebula. As bare chondrules move through a dusty disk they would collect the dust. This dust is compressed to the observed rims, eventually (Bland et al. 2011; Metzler et al. 1992; Ormel et al. 2008; Morfill et al. 1998; Carballido 2011). From an experimentalists point of view this is very plausible as it is almost impossible to prevent that a chondrule in a dust rich environment will enshroud itself in a mantle of small grains (Blum & Wurm 2008). Specific experiments to study accretion of dust rims and further evolution have recently been carried out by Beitz et al. (2011). As far as the size sorting is concerned also a number of different scenarios are sketched in the literature. The local reservoirs of matter in the solar nebula, which created chondrules and which fed a parent body, might just have produced different size distributions in different locations. More often though it is discussed that certain mechanisms might have acted on a particle reservoir, which separated chondrules of different size (Liffman 2005). Aerodynamic aspects are clearly tempting here as particles of different size couple differently to the embedding gas (Cuzzi et al. 1996). All the explanations for size sorting and rim formation seem to have their advantages and shortcomings. We cannot dismiss one or the other as within their set of assumptions all might currently be viable scenarios. In this paper we assume that the rims set off as highly porous dust mantles on chondrules of nebula origin and add just one more assumption, namely that the particles are illuminated by the young sun. Adding a light source might then naturally lead to a (size) sorting due to photophoresis.
Photophoresis was first introduced to the physics of protoplanetary disks by Krauss & Wurm (2005). First attempts to explain size sorting and concentration of chondrules by photophoresis were made by Wurm & Krauss (2006). They showed that in an optically thin solar nebula chondrule aggregates are sorted due to the size of the constituents in different locations. Key factor here is the thermal conductivity of the aggregates which determines the temperature gradient upon illumination. First experiments to measure photophoresis on chondrules have been carried out by Wurm et al. (2010). The measurements support the idea that photophoresis acts sufficiently strong on chondrules to provide sorting. The assumption of a directed radiation source is vital to the problem and has two very clear sides. On one side, without a directed radiation source (visible or infrared) photophoresis does not exist. Therefore, photophoresis does not exist in the midst of an optical thick protoplanetary disk. However, on the other side, given a light source it is also inevitable that photophoresis acts as significant force on the particles. In the inner parts of an optical thin protoplanetary disks it is orders of magnitudes larger than a stars gravity and cannot be neglected (Moudens et al. 2011; Krauss et al. 2007; Mousis et al. 2007; Wurm & Haack 2009; Takeuchi & Krauss 2008; Herrmann & Krivov 2007). The assumption that the solar nebula was optically thin at certain times and places is supported by observations of protoplanetary disks. There are a number of disks observed now which have large inner clearings, which are optically thin (Calvet et al. 2002; SiciliaAguilar et al. 2008). In some of these disks the star is still accreting matter so that the inner region still must contain gas. Clearly the conditions for photophoresis are given here. Also, disks in general have an inner edge. At the inner edge of a disk sorting of chondrules is possible as some particles might be pushed into the dark a bit stronger than others. This idea has been suggested by Haack & Wurm (2007) and Wurm et al. (2006). To evaluate the potential of size sorting and photophoretic transport in further detail in the future we carried out numerical heat transfer calculations on bare chondrules and dust mantled chondrules as prerequisite here. These are used to quantify the photophoretic strengths on particles more accurately. The structure of this paper is as follows. Section 2 describes our heat transfer simulations, which result in temperature fields at the surface of a (dust mantled) chondrule that is illuminated. In Sect. 3 the resulting temperature field at the surface of a particle is used to calculate photophoretic forces. These forces are compared to values, which are calculated from common equations found in the literature. To facilitate analytical treatment of photophoretic forces in future work we introduce correction factors on the order of 20% for homogeneous particles of given size and thermal conductivity. These results are then used in Sect. 4 to define average thermal conductivities for dust mantled chondrules depending on the chondrule size, dust mantle size and the thermal conductivities of both components. In Sect. 5 we sketch how chondrules might be size sorted due to photophoresis based on our calculations.
2. Heat transfer calculations
Our basic model setup consists of a spherical particle, which is made up of a core (bare chondrule) and a dust mantle (rim) (Fig. 1).
Fig. 1 Model setup of a chondrule core (high thermal conductivity) and a dust mantle (low thermal conductivity). The light source is directed downwards. Due to rotational symmetry only the temperature distribution along one line from top (180°) to the bottom (0°) has to be considered. 

Open with DEXTER 
Both components have specific thermal conductivities. The core, being compact and consisting of silicates, has a thermal conductivity on the order of k_{core} = 1 W/(mK) (Opeil et al. 2010). The dust mantle being highly porous has a very low thermal conductivity. We vary k_{dust} between 0.01 W/(mK) and 0.5 W/(mK) (Presley & Christensen 1997a,b; Krause et al. 2011). This coremantle sphere is illuminated from the top (Fig. 1), absorbs radiation and heats up until in equilibrium it emits as much thermal radiation as it absorbs visible radiation. Due to the directed illumination a temperature gradient over the particle surface is established. The problem to solve is the heat transfer equation (1)where k is the local thermal conductivity, T is the local temperature and Q is a volumetric heat source. As boundary condition we have cooling due to thermal emission and heating due to illumination. We consider a complete absorption of incident light at the surface (i.e. Q = 0). We further assume an emissivity for thermal radiation of 1. For a plane wave light source of light flux I ∗ e_{I} we then get a von Neumann boundary condition (2)where n is normal to the surface pointing outwards and the first term on the righthand side is only considered for the upper illuminated half sphere of the particle. The light flux in all calculations is fixed at I = 20 kW/m^{2}. Equations (1) and (2) are solved numerically using COMSOL, which is based on a finite element method to solve differential equations (COMSOL v4.1 from COMSOL AB, Sweden). In case of no rotation the temperature distribution is symmetric around the axis given by e_{I}. The maximum temperature T_{max} is at the top surface element facing the light source (at 180°). The minimum temperature T_{min} is at the opposite side (bottom or 0°). The environment is assumed to have T_{gas} = 293 K. This is far below the equilibrium temperature for an illuminated particle but reflects experimental settings, e.g. in Wurm et al. (2010). It will also allow us to specify dependencies on temperatures more detailed and not just refer to a single average gas temperature, i.e. we will not assume one average temperature for the gaseous environment and the particle. A sample of calculations is shown in Fig. 2 for a homogeneous 1.12 mm particle (total diameter). To calculate homogeneous particles we use the same thermal conductivity for the dust rim and the chondrule (core).
Fig. 2 Temperature distributions for homogeneous spheres of 1.12 mm diameter with varying thermal conductivities. 

Open with DEXTER 
All areas of validity and restrictions of the model are discussed in Sect. 6.
3. Photophoretic forces
3.1. general case at low gas pressure
The photophoretic force depends on the temperature distribution over the surface of the chondrule. We consider the case of free molecular flow or low gas pressure here, where the mean free path of the gas molecules is much larger than the particle size. In that case the interaction between particle and gas can be treated as individual collisions of the gas molecules with the surface. A gas molecule impinging the surface of the particle is either reflected specular or diffusely. In the first case, due to the larger mass of the solid particle the absolute momentum of the gas molecule does not change and does not depend on the temperature of the particle. If this first fraction of particles, 1 − α, is the same over the surface of the particle, then due to symmetry, specular reflected gas molecules on average do not transfer momentum to the particle. The second fraction of gas molecules, α, which are reflected diffusely, stick to the surface for a short while, accommodate to the surface temperature and the momentum of the reemitted gas molecules differs from the momentum of the incoming gas molecules. On a hot surface part the gas molecules leave faster than on colder parts. Therefore the gas carries a net momentum after the interactions. The total force on the particle is the integral over the total particle surface of the local pressure induced by the interaction at the local temperature and in the free molecular flow regime is given as (Rohatschek & Zulehner 1985; Rohatschek 1995; Hidy & Brock 1970) (3)where denotes the statistical mean temperature and α is the thermal accommodation coefficient describing the fraction of gas particles being adsorbed by the solid particle’s surface, T is the local temperature of the particle surface, T_{gas} and p are the gas temperature and pressure far away from the particle, respectively. As T_{gas} in our calculations we use 293 K.
Fig. 3 Ratio for a given particle size and thermal conductivity k. For thermal conductivities larger than k = 0.2 the factor is essentially constant z^{(1)} = 0.7 (α = 1, this plot) (resp. z^{(1)} = 0.8 for α = 0.8). 

Open with DEXTER 
3.2. rotationalsymmetric temperature field
In case of a spherical rotationalsymmetric temperature field T = T(θ), e.g. if a sphere is illuminated with an intensity I in direction e_{I}, without loss of generality let e_{I}: = −e_{z} (see Fig. 1) in order to develop T in Legendre Polynomials P_{n}: (4)with A_{i} as the ith expansion coefficient. Rohatschek (1995) further linearizes the square root to its first order. Due to the Legendre Polynomial’s orthogonality relation F (see Eq. (3)) is only a function of A_{1}; all other A_{i} do not show up since the corresponding integrals yield zero. Therefore, Eq. (3) yields (Rohatschek 1995) (5)where r denotes the sphere’s radius.
Since A_{1} is usually unknown there are three well known approximations for this coefficient. The first one is (indication (i) refers to the ith approximation presented) (6)with ΔT being the maximum temperature difference between top and bottom ΔT = T_{max} − T_{min}. Thus Eq. (5) yields an approximated force for a nonrotating illuminated sphere: (7)Our first goal here is to calculate the photophoretic force for a homogeneous sphere (core and rim share the same thermal conductivity k) based on Eq. (3) exactly but describe the results in a way similar to the simple analytic expression of Eq. (7). In the following we therefore determine a correction factor from simulations. We define as correction (8)Varying size r and thermal conductivity k we find that in general, the approximation overestimates the photophoretic force. The correction factor can be considered a constant of z^{(1)} = 0.7 for thermal conductivities larger than about k = 0.2 W/mK independent of the particle size (Fig. 3). For smaller thermal conductivities there is a dependency of z^{(1)} on both parameters. As the thermal conductivity of dust falls in the low range a nonconstant correction factor has to be considered. The parameters studied range from r = 0.1...11 mm, k = 0.01...8 W/mK, and α = 0.1...1, the latter is explicitly part of Eq. (5). To give an analytic equation we tried different simple expressions. One expression that fits the correction factor over the whole range with an accuracy being better than 2% is (9)It should be noted, that in Eqs. (9), (13), (18) and (19) r is in units of mm and k in units of W/(mK).
The next step is expressing the temperature difference in terms of an illuminating flux. For a homogeneous sphere of thermal conductivity kRohatschek (1995) finds (10)where I is the radiative flux. As Rohatschek (1995) notes, this is an approximation for good conductors. With low thermal conductivity dust mantles on high conductivity chondrules, it is not clear that the approximation of Eq. (10) is sufficient to study the details of photophoretic particle transport and sorting. However, Eqs. (5) and (10) are useful as they provide a rather simple analytical treatment of photophoresis including all important quantities, i.e. the gas properties p, T_{gas}, the radiative flux I and accommodation coefficient α, and the particle properties r, k. To keep this simplicity but to be accurate enough to apply the equations to problems of photophoretic transport of chondrules we chose again a correction term. We extend our correction factor including deviations introduced by the approximation of Eq. (10) and define (11)with (12)
Fig. 4 Ratio for a given particle size and thermal conductivity k. For thermal conductivities larger than 2 the factor is essentially constant. The accommodation coefficient α is 1. 

Open with DEXTER 
As before, the approximation (Eq. (12)) overestimates the photophoretic force following the same pattern as before. In the same parameter set (r = 0.1...11 mm, k = 0.01...8 W/mK, α = 0.1...1) z^{(2)} can be well fitted again to better than 2% by the following expression (13)This correction factor (Fig. 4) can be rather far from unity especially for low thermal conductivity particles. Clearly, particles with thermal conductivities much lower than 1 W/mK can no longer be regarded as good conductors in the sense of Eq. (10). Doing so would usually overestimate the photophoretic force significantly. This is due to the fact that the approximation of Eq. (10) neglects the effect of thermal radiation. This aspect can be improved by adding a term to Eq. (10) which accounts for this (e.g. cmp. Beresnev et al. 1993) (14)Usually, particles are considered in thermal equilibrium with their surroundings, photophoretic temperature gradients are small and illumination is not significantly changing the average temperature. In that situation the particle temperature T_{par} is approximated by T_{gas}. For our calculations and often for experimental settings this is not the case, i.e. the particle temperature can easily be a factor 2 higher than the surrounding gas temperature. As the radiative cooling is depending on the particle temperature it is the particle temperature and not the gas temperature, which enters in Eq. (14). As seen in Fig. 2 the particle temperature varies strongly over its surface and is still not well defined.
Fig. 5 Ratio for a given particle radius and thermal conductivity k. For thermal conductivities larger than 0.2 the factor can be regarded as constant. The accommodation coefficient α is 1. 

Open with DEXTER 
We use the average black body temperature here and define (15)Based on Eqs. (14) and (15) we finally consider as analytical approximation that captures all physical parameters of the system but is accurate enough to treat dust mantled chondrules (16)and quantify the photophoretic force adding a correction factor z^{(3)}(17)Similar in structure, z^{(3)}, can well be approximated in the whole range of parameters to better than about 2% by the following terms (18)As summary, if the temperature gradient is known the first approximation (Eqs. (7) and (9)) might be used. More general and used further in this paper is the photophoretic force on a homogeneous, perfectly absorbing spherical particle as given by Eqs. (16) to (18) or combined as (19)
4. Dust mantled chondrules
In view of Eq. (19), thermal conductivity of homogeneous spheres is a major particle property that determines photophoresis. For dust mantled chondrules there are two distinct different thermal conductivities for the core and the mantle and it is not clear, a priori, what determines the photophoretic strength of such particles. A core mantle particle can be described by five parameters, which we vary within the following range:

the coreradius, r_{core}: 0.2...1.0 mm,

the mantle thickness, d_{dust}: 0.1...1 in fractions of r_{core} (the total radius then is ρ: = r_{core}^{(}1 + d_{dust}^{)}),

the thermal conductivity of the core, k_{core}: 1...4 W/mK,

the thermal conductivity of the mantle, k_{dust}: 0.01...0.5 W/mK,

and the accommodation coefficient, α: 0.1...1.
For each combination of these parameters a photophoretic strength can be calculated. In principle we could argue on the photophoretic strength depending on these parameters at this point. We note that different configurations of dust coated chondrules might result into the same photophoretic strength, even if the temperature field across the mantle is absolutely different. So two configurations γ_{1},γ_{2} are related to each other if (20)Nevertheless, it is possible to choose a representative γ(ρ,κ) of all those configurations obeying Eq. (20) with a (total) radius ρ and constant thermal conductivity κ. We call κ effective thermal conductivity.
Fig. 6 κ over d_{dust} and k_{dust} for a 1 mm sized core with a thermal conductivity of 4 W/mK. The accommodation coefficient α is 1; white lines indicate constant values and show multiple configurations can lead to the same κ. 

Open with DEXTER 
In view of quantifying thermal properties of dust mantled spherical chondrules with a total radius we use the results of the photophoretic strength calculations and Eqs. (16) and (18) to calculate . It is of course not injective, so multiple configurations can lead to a same κ (see Fig. 6). In the following the dependency of κ on the different parameters is considered. We therefore numerically solved (21)for κ, where is the sphere’s surface temperature (depending on all four parameters for the heat equation system).
Accommodation coefficient: the accommodation coefficient is not connected to κ (by construction). Numerically, variations only influence the result on the 1% level. This is in agreement with the idea that thermal conductivity is independent of the photophoretic property.
Core and mantle size: in view of estimating the potential for size sorting the most important dependencies to explore are the influence of the total size of the chondrule and the influence on the ratio between mantle thickness and core radius. We find that only for low thermal conductivities of the dust mantle does the effective thermal conductivity slightly depend on the total size. The extreme for a dust mantle of the lowest thermal conductivity studied of k_{dust} = 0.01 W/mK is seen in Fig. 7. The thermal conductivity of the core is 1 W/mK. All three lines vary very little. This shows that the dependency on total size for a given mantle to core size ratio is very weak. For given thermal conductivities of core and mantle only the ratio between the sizes of both components determines the effective thermal conductivity. This trend can be understood as the smaller the mantle the more effective is the high thermal conductivity of the core. In the extreme of a vanishing mantle the thermal conductivity has to reach the core value, which is consistently confirmed by the calculations. For thick mantles the low thermal conductivity of the mantle should dominate, which is observed in Fig. 7. However, it should be noted that even for equal size components the effective thermal conductivity is still a factor 2 larger than the thermal conductivity of the pure dust mantle.
Fig. 7 Effective thermal conductivity κ over dust mantle to core size ratio. The 3 different lines represent 3 different core sizes: from bottom to top r_{core} = 0.2,0.5,1.0 mm. The thermal conductivity of the dust mantle is 0.01 W/mK, the core’s thermal conductivity is 1 W/mK and the accommodation coefficient α is 1. 

Open with DEXTER 
Thermal conductivity of core and mantle: from Fig. 7 it is obvious that the thermal conductivity of the mantle is a dominating factor for the effective thermal conductivity. In Fig. 8 the dependency on the size ratio between rim and core and the thermal conductivity of the mantle is shown for two core thermal conductivities. The two different plots are based on the extreme values used for the thermal conductivity of the core in our calculations. The comparison illustrates that the core has a significant influence only for thin mantles and then only if these mantles have high thermal conductivities.
Fig. 8 κ over d_{dust} and k_{dust} for a 1 mm sized core with a thermal conductivity of 1 W/mK (bottom) and 4 W/mK (top). The accommodation coefficient α is 1. 

Open with DEXTER 
Figure 9 shows the dependency in higher resolution again by plotting κ over d_{dust} but for two fixed thermal conductivities of both the mantle (0.01 W/mK, 0.1 W/mK) and core (1 W/mK, 4 W/mK). Even for thin dust mantles, if the thermal conductivity has a reasonable value attributed to dust, then the effective thermal conductivity does not significantly depend on the thermal conductivity of the core.
Conclusion Sect. 4: in total we find that the effective thermal conductivity significantly depends only on three of the four parameters, the relative thickness of the dust mantle and the thermal conductivities of core and mantle but the total size plays no major role for thermal conductivity. The thermal conductivity of the core is visible in the values of the effective thermal conductivity, i.e. by increasing the value of the mantle by a factor of two or more. However, the effective thermal conductivity is insensitive to variations of the core thermal conductivity in the range of reasonable values of k_{dust}. Assuming that the solid chondrules (core) have a fixed thermal conductivity on the order of 1 W/mK, this leaves two variables in the context of photophoretic forces on chondrules, the thermal conductivity of the mantle and the dust mantle to core size ratio.
5. Size sorting of chondrules
The dependency of the effective thermal conductivity on the different parameters provides various options by which dust mantled chondrules can be sorted due to photophoresis. We only give a qualitative account here. Detailed treatments would require specific disk models. This would be a problem based on the results of this paper but otherwise being beyond its scope.
Fig. 9 κ for different pairs . Two bottom lines: k_{dust} = 0.01 W/mK, two top lines: k_{dust} = 0.1 W/mK 

Open with DEXTER 
The sorting mechanism in Wurm & Krauss (2006) was based on a disk density profile which decreases with distance to the star. This way, the gas is supported against the stellar gravity and rotates slower than with Keplerian speed in order to rotate on a stable orbit. Small solid particles like chondrules embedded in the disk would couple to the gas rapidly, would rotate with the same slower orbital frequency but would not feel the support by the pressure gradient. Such particles drift inward (Weidenschilling 1977). The particles are subject to a radial force (22)then which can be considered as residual gravity. Here a is the radial distance to the star, μ is the molar mass of the gas, n is the power of the radial distance dependency of the gas density (n = 11/4 in the minimum mass solar nebula, see Eq. (26) Hayashi et al. 1985), R_{G} is the gas constant and m_{par} is the particle mass. This includes gas drag in the way that the gas grain coupling time is much shorter than the orbital timescale and the problem reduces to a one dimensional radial problem.
Wurm & Krauss (2006) consider the simple photophoretic force estimate given by Eq. (12). and calculate the ratio (23)If this ratio is 1 inward and outward drift balance and chondrules are concentrated. This concentration distance depends on the strength of the photophoretic force. Therefore, chondrules with different photophoretic properties – as outlined in this paper – would concentrate at different locations.
Particle sorting due to photophoresis might also take place at the inner edge of the disk, which might move over the asteroid belt (Haack & Wurm 2007). The idea is sketched in Fig. 10. At the inner edge of a disk which is optically thick the sorting might work differently than described above. Per definition the optical depth decreases further away from the edge. Particles are therefore pushed outward as long as the decreasing photophoretic outward motion is more effective than any other inward directed effect like turbulent diffusion or any other remaining radial inward drift and the particles pushed more vigorously by photophoresis will move outward into the dark the furthest. Here, the sorting distances are strongly disk dependent but the principle is independent of the details.
Wurm & Krauss (2006) considered that particles with thermal conductivities varying by a factor 2 could easily be concentrated at different positions in a solar nebula almost 1 AU apart of each other. A rough linear scaling would imply that a difference in effective thermal conductivity on the order of 10% would separate particles on the order of 0.1 AU or 15 × 10^{6} km. This is certainly sufficient to separate particles far enough to be incorporated into different parent bodies, eventually.
Fig. 10 Principles of sorting dust mantled chondrules by photophoresis. 

Open with DEXTER 
Direct size sorting: as calculated above the effective thermal conductivity does not change with total size. The photophoretic force has an intrinsic dependence on size though as r^{3} (Eq. (19)) with the correction factor z^{(3)} included. For a homogeneous particle the gravitational pull by a star also depends on r^{3} and the ratio of both forces is essentially independent of the total size. Without correction factor photophoresis would not sort by size directly. For the lowest thermal conductivities there is a size dependence in the correction factor which varies the photophoretic force by 10% for the smallest and largest chondrule (see Fig. 5). In that case direct size sorting of dust mantled chondrules with plausible mantle thickness by photophoresis is possible. It might also be considered here that the dust mantle has a different porosity as the core, which might vary by about a factor of 3. Nevertheless, for a constant size ratio between mantle and core a change in total size would not change the average density of the particle and this is not important.
Mantle to core size ratio sorting: it is interesting to note that for a constant dust mantle to core size ratio no sorting due to the effective thermal conductivity would occur. However, considering the growth of dust mantles by accretion, Carballido (2011) shows that the size ratio between mantle and core changes with overall size. This introduces an absolute size dependence of the effective thermal conductivities then: κ(d_{dust}(ρ)). Values estimated from Fig. 9 in Carballido (2011) are e.g. a 0.3 mm core with a 0.2 mm mantle (d_{dust} = 0.6) but a 1 mm core only has 0.4 mm mantle (d_{dust} = 0.4). Corresponding effective thermal conductivities from our simulations assuming the core to have, k_{core} = 1 W/mK and the rim to have k_{dust} = 0.01 W/mK are κ = 0.02 W/mK and κ = 0.03 W/mK, respectively. This is a difference of a factor 1.5 which is certainly sufficient to separate these particles by photophoresis. A thinner mantle also corresponds to a slightly higher average mass density and gravitational attraction. This works in the same direction as decreasing the strength of the photophoretic force.
Mantle thermal conductivity: the exact thermal conductivity of the mantle is unknown as it depends on many parameters. Assuming that the material properties of the dust are the same on average, porosity of the dust rim is still a free parameter. The work by Bland et al. (2011) suggests that mantle porosities are on the order of 0.7 or more. This is consistently on the same order of porosities to be expected from dust experiments (Teiser et al. 2011). For this porosity and micrometer size dust Krause et al. (2011) find thermal conductivities on the order of 0.01 W/mK and below. We assumed this value for dust mantles above. Typically smaller particles collide with dust at lower velocities than larger ones. The experiments by Teiser et al. (2011) and Kothe et al. (2011) show that such differences in collision velocities produce a difference in porosity. In that case the dust mantle porosity would systematically change (decrease) with total size. Typical differences are several percent. Assuming a porosity variation e.g. of 0.05 would lead to a difference in thermal conductivity on the order of 1.5 for the dust mantle (Krause et al. 2011). Varying the dust mantle thermal conductivity by this factor, e.g. between 0.01 W/mK and 0.015 W/mK for a core with 1 W/mK changes the effective thermal conductivity by a similar factor of 1.54 (r_{core} = 1 mm,d_{dust} = 1) or 1.46 (r_{core} = 1 mm,d_{dust} = 0.5). This is on the order of the changes induced by the different mantle to core size ratios discussed in the section before. This provides the basis for a photophoretic size sorting. The trend would go in the same direction as the size ratio variations, i.e. larger particles have a lower effective thermal conductivity and would be found closer to the sun.
Core thermal conductivity: the third parameter influencing the effective thermal conductivity is the thermal conductivity of the core. As it is plausible that dust mantles have rather low thermal conductivity, the dependency on this parameter is not strong. As found in the simulations effective thermal conductivities for thin (20% of core radius) dust mantles vary by 10%. For more conductive rims of k_{dust} = 0.1 W/mK this might be somewhat larger but compared to the dependency on the size ratio between mantle and core the effect is small. We conclude that for coremantle particles the thermal conductivity of the core will not be responsible for size sorting.
Wurm & Krauss (2006) estimate the drift velocity of chondrules due to photophoresis at 1 AU which would be on the order of 0.1 m/s. Separating the particles by 10 × 10^{6} km would require about 10^{7} s or 1 year or 1 orbit then which is fast compared to evolutionary time scales. If all mass is in chondrules and if chondrules do not evolve e.g. by growing to larger bodies then a maximum mass of separated chondrules can be estimated by the number of chondrules which would make the disk optically thick themselves. Assuming a disk height of 0.1 AU at 1 AU distance to the star the area of the inner edge would be 2π × 1 AU × 0.1 AU. Completely covered with chondrules of 1 mm diameter, the number of chondrules is 6.35 × 10^{27} and assuming these chondrules have a mass density of 3 g/cm^{3} one chondrule has a mass of 1.6 mg. The total number then is 10^{22} kg which is about 10 times the mass of Ceres. Being that large the estimate shows that the mechanism provides more than sufficient separation rates even if conditions are less perfect than assumed for this estimate.
6. Caveats
Radiation pressure: photophoresis is driven by radiation. Other movements of particles induced by radiation and more commonly considered are radiation pressure and the PoyntingRobertson effect (Klacka & Kocifaj 2002). The difference is that radiation pressure is based on the direct momentum transfer of photons while photophoresis is based on momentum transfer via gas molecules. For chondrule size particles and plausible gas densities photophoresis is much larger as seen in the following estimate. Radiation pressure can be approximated by (24)Compared to the classical photophoretic force estimate of Eq. (12) which is sufficient for this estimate as photophoretic force it is (25)Taking a typical value of k = 0.1 W/mK, a worst case of a small chondrule of r = 0.2 mm and a worst case temperature of T_{gas} = 280 K we get F_{phot}/F_{rad} ≃ 357 p. The minimum mass solar nebula by Hayashi et al. (1985) has an ambient pressure of about 1.4 Pa at 1 AU. Here, photophoresis is a factor of 500 stronger than radiation pressure. The pressure is 0.04 Pa at 3 AU or distances related to the outer asteroid belt. Here, photophoresis is 14 times stronger than radiation pressure using the above conservative estimates. Therefore, photophoresis dominates as radiative force.
Rotation: particle rotation was not considered here though it is important in the context of photophoresis. If particles would rotate more rapidly than the typical heat transfer time then the temperature gradient would not reach equilibrium and be reduced which leads to a decrease in the photophoretic force. We did not consider rotations for the following reasons. Any random rotation is damped on a time scale of the gasgrain coupling time. In the free molecular flow regime which is relevant for mmparticles in protoplanetary disks the coupling time is given as (Blum et al. 1996) (26)with ρ_{gas} being the gas density, σ_{par} the particle’s geometrical cross section, and the mean thermal velocity . At a = 1 AU in a minimum mass nebula it is (Hayashi et al. 1985) Then, τ_{c} is on the order of 10^{3} s for a 1 mm particle. Random rotation is therefore decaying rapidly. Random rotations would only be important in decreasing the photophoretic strength effectively if they would continuously be excited on similar time scales. A way to generate particle rotation is interparticle collisions. Krauss et al. (2007) estimated that the time between collisions of mmsize particles of massdensity ρ_{par} in a “typical” disk is given by (30)which then is about a factor of 10^{4} longer than the gas grain coupling time (v_{par,rel} = 1 m/s). Therefore, even assuming strongly varying disk conditions, rotation is damped much more rapidly and the total time of potentially decreased photophoretic strength is not significant.
Another source of random rotation to consider is Brownian motion. The thermal rotation energy of a spherical particle is given by , which leads to (31)with J_{par} being the sphere’s moment of inertia. To rotate by 180° and change front and backside it is (32)For chondrules this is on the order of 10^{4} s. The time scale for thermal conduction through the particle is given by (33)Here, c_{par} is the particle’s heat capacity at constant pressure. For chondrules τ_{cond} is on the order of 1 s. Therefore, a temperature gradient reestablishes itself much faster than thermal rotation can disalign it and hence Brownian Rotation can be neglected, too.
Besides random rotation exciting a continuous rotation is possible by photophoresis itself. As particles are never ideal spheres torques are induced. One can argue though that only torques around the direction of light will induce a systematic rotation around the axis aligned to the direction of light. Any other torque will only lead to the stable alignment of the particle. To prove this van Eymeren & Wurm (2012) carried out a detailed study of particle rotation. In their setup particles were levitated by means of thermophoresis and photophoresis (Kelling et al. 2011). van Eymeren & Wurm (2012) found that at least more than 95% of all particles rotate around the vertical direction which coincides with the direction of illumination. Such a rotation does not change front and back with respect to the illumination and is not important for the photophoretic strength along the direction of light. For this reason rotation is no issue in particle transport by stellar radiation. It might be noted that the situation is different for radiation pressure acting in the current solar system. However, here no gas and therefore no damping is present.
Heat exchange with the ambient gas: in our heat transfer calculations we did not include heat exchange with the ambient gas. Heat transfer to the ambient gas for large Knudsen numbers can be quantified as (Stoffels et al. 1996) (34)where I is the heat flux (power per surface area).
In the case of molecular hydrogen, which consists of two atoms, c_{p}/c_{v} can by approximated by 9/7. ΔT denotes the temperature difference between the chondrule and the surrounding gas. For a chondrule at 2 AU in a minimum mass nebula by Hayashi et al. (1985) we have a radiative flux I = 341 W/m^{2}, p = 0.15 Pa,T_{gas} = 262 K. Maximum temperature differences are ΔT = 10 K. This leads to I_{cond} = 35 W/m^{2}. This is an order of magnitude less than the incoming radiative flux and can be neglected in the context of photophoretic sorting of chondrules in protoplanetary disks.
Typical laboratory conditions are I = 20 kW/m^{2} (maximum ΔT = 400 K, see Fig. 2), p = 1 Pa, m_{gas} = 29 u (air), and T_{gas} = 293 K which yields I_{cond} = 315 W/m^{2}. This is also much lower than the incoming light flux and cooling by gas interaction still can be neglected. In case of higher pressures (p = 10 Pa...100 Pa) heat exchange with the gas becomes important, but at the same time the assumption of photophoresis in the free flow regime is also no longer valid. For transitional disks assumed to be of low pressure our treatment holds. On the other hand higher disk (gas) densities lead to a reduced light flux as they are usually assumed to go with higher absorption. This also changes the photophoretic strength and transport. In the extreme case of an optically thick disk, photophoresis does not work. Therefore, high density regions of a disk require an extended treatment of photophoresis, but this is beyond the scope of this paper.
7. Conclusion
In summary there are two options to size sort chondrules with dust mantles by photophoresis. The mantle thickness variation with size (Carballido 2011) naturally induces different effective thermal conductivities. Plausible porosity variations with size act in the same direction and on the same magnitude. Therefore, if dust covered chondrules are subject to a directed radiation they likely get size sorted by photophoresis. In the original work by Wurm & Krauss (2006) the solar nebula was assumed to be optically thin throughout. This is not likely to be the case for the very young nebula. However, an increasing number of transition disks is found, which have an optical thin inner clearing up to tens of AU (Najita et al. 2007; SiciliaAguilar et al. 2008). Some of these show active signs of accretion and therefore the optical thin part is still gaseous. Within such a clearing the idea by Wurm & Krauss (2006) can naturally be applied. The effective thermal conductivities of dust mantled chondrules are 1 to 2 orders of magnitude lower than the thermal conductivity of bare chondrules assumed by Wurm & Krauss (2006). This increases the photophoretic strength by a similar factor. Therefore, dust mantled chondrules can concentrate in the asteroid belt region even if the inner clearing has a significantly reduced gas density. Creating the inner clearing takes some time. This would also fit well with the fact that chondrules are forming their parent bodies only a few million years after the beginning of the solar nebula. Haack & Wurm (2007) also considered the idea that the inner edge of the disk might be important. This would shift the focus from the spatially extended inner clearing to the actual edge that is moving outwards. Here the disk is supposed to be of its original density. Particles exposed to the starlight at the edge would feel a photophoretic force and be pushed outwards until the optical depth is too large. The final position with respect to the edge will depend on the susceptibility to photophoresis and particles with lower effective thermal conductivity will move further out. This will introduce a size sorting at the edge of the disk as it moves inside out. If parent bodies form in this increased density environment chondrules will be size sorted. A detailed treatment of this is beyond the scope of this paper as it would have to focus on disk models, while this work is dedicated to the photophoretic properties needed as input. Here, we evaluated the effective thermal conductivity of a dust mantled chondrule in the context of photophoretic forces. It strongly depends on the size ratio between mantle and core and the thermal conductivity of the mantle. Based on these results a size sorting of dust mantle chondrules by photophoresis is possible.
Acknowledgments
This work is funded by the Deutsche Forschungsgemeinschaft within the priority program SPP 1385. We also thank the referee for a very constructive review of the manuscript.
References
 Beitz, E., Blum, J., Mathieu, R., Pack, A., & Hezel, D. C. 2011 [arXiv:1109.4018] [Google Scholar]
 Beresnev, S., Chernyak, V., & Fomyagin, G. 1993, Phys. Fluids, 5, 2043 [NASA ADS] [CrossRef] [Google Scholar]
 Bland, P. A., Howard, L. E., Prior, D. J., et al. 2011, Nat. Geosci., 4, 244 [NASA ADS] [CrossRef] [Google Scholar]
 Blum, J., & Wurm, G. 2008, ARA&A, 46, 21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Blum, J., Wurm, G., Kempf, S., & Henning, T. 1996, Icarus, 124, 441 [NASA ADS] [CrossRef] [Google Scholar]
 Calvet, N., D’Alessio, P., Hartmann, L., et al. 2002, ApJ, 568, 1008 [NASA ADS] [CrossRef] [Google Scholar]
 Carballido, A. 2011, Icarus, 211, 876 [NASA ADS] [CrossRef] [Google Scholar]
 Cuzzi, J. N., Dobrovolskis, A. R., & Hogan, R. C. 1996, in Chondrules and the Protoplanetary Disk, ed. R. H. Hewins, R. H. Jones, & E. R. D. Scott, 35 [Google Scholar]
 Haack, H., & Wurm, G. 2007, Meteorit. Planet. Sci. Suppl., 42, 5157 [NASA ADS] [Google Scholar]
 Hayashi, C., Nakazawa, K., & Nakagawa, Y. 1985, in Protostars and planets II (Tucson, AZ, University of Arizona Press), eds. D. C. Black, & M. S. Matthews, 1100 [Google Scholar]
 Herrmann, F., & Krivov, A. V. 2007, A&A, 476, 829 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hidy, G., & Brock, J. 1970, The Dynamics of Aerocolloidal Systems (Oxford: Pergamon Press) [Google Scholar]
 Hughes, D. W. 1978, Earth Planet. Sci. Lett., 38, 391 [NASA ADS] [CrossRef] [Google Scholar]
 Kelling, T., Wurm, G., & Dürmann, C. 2011, Review of Scientific Instruments, 82, 115105 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Klacka, J., & Kocifaj, M. 2002, in Proceedings of Asteroids, Comets, Meteors – ACM, ed. B. Warmbein (Noordwijk, Netherlands: ESA Publications Division), ESA SP, 500, 169 [Google Scholar]
 Kothe, S., Güttler, C., & Blum, J. 2011, in EPSCDPS Joint Meeting 2011, Nantes, France, 412 [Google Scholar]
 Krause, M., Blum, J., Skorov, Y. V., & Trieloff, M. 2011, Icarus, 214, 286 [NASA ADS] [CrossRef] [Google Scholar]
 Krauss, O., & Wurm, G. 2005, ApJ, 630, 1088 [NASA ADS] [CrossRef] [Google Scholar]
 Krauss, O., Wurm, G., Mousis, O., et al. 2007, A&A, 462, 977987 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kuebler, K. E., McSween, H. Y., Carlson, W. D., & Hirsch, D. 1999, Icarus, 141, 96 [NASA ADS] [CrossRef] [Google Scholar]
 Liffman, K. 2005, Meteorit. Planet. Sci., 40, 123 [NASA ADS] [CrossRef] [Google Scholar]
 Metzler, K., Bischoff, A., & Morfill, G. 1991, LPI Contributions, 766, 155 [NASA ADS] [Google Scholar]
 Metzler, K., Bischoff, A., & Stoeffler, D. 1992, Geochim. Cosmochim. Acta, 56, 2873 [NASA ADS] [CrossRef] [Google Scholar]
 Morfill, G. E., Durisen, R. H., & Turner, G. W. 1998, Icarus, 134, 180 [NASA ADS] [CrossRef] [Google Scholar]
 Moudens, A., Mousis, O., Petit, J.M., et al. 2011, in Lunar and Planetary Inst. Technical Report, Lunar and Planetary Institute Science Conference Abstracts, 42, 1409 [NASA ADS] [Google Scholar]
 Mousis, O., Petit, J.M., Wurm, G., et al. 2007, A&A, 466, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Najita, J. R., Carr, J. S., Glassgold, A. E., & Valenti, J. A. 2007, Protostars and Planets V, 951, 507 [NASA ADS] [Google Scholar]
 Opeil, C. P., Consolmagno, G. J., & Britt, D. T. 2010, Icarus, 208, 449 [NASA ADS] [CrossRef] [Google Scholar]
 Ormel, C. W., Cuzzi, J. N., & Tielens, A. G. G. M. 2008, ApJ, 679, 1588 [NASA ADS] [CrossRef] [Google Scholar]
 Presley, M. A., & Christensen, P. R. 1997a, J. Geophys. Res., 102, 6535 [NASA ADS] [CrossRef] [Google Scholar]
 Presley, M. A., & Christensen, P. R. 1997b, J. Geophys. Res., 102, 6551 [NASA ADS] [CrossRef] [Google Scholar]
 Rohatschek, H. 1995, J. O. Aerosol Sci., 26, 717 [CrossRef] [Google Scholar]
 Rohatschek, H., & Zulehner, W. 1985, J. O. Colloid a. Interface Sci., 10, 457 [CrossRef] [Google Scholar]
 Scott, E. R. D., Love, S. G., & Krot, A. N. 1996, in Chondrules and the Protoplanetary Disk, ed. R. H. Hewins, R. H. Jones, & E. R. D. Scott, 87 [Google Scholar]
 SiciliaAguilar, A., Henning, T., Juhász, A., et al. 2008, ApJ, 687, 1145 [NASA ADS] [CrossRef] [Google Scholar]
 Stoffels, W. W., Stoffels, E., Kroesen, M. W., & de Hoog, F. J. 1996, J. Vac. Sci. Technol. A, 14, 588 [NASA ADS] [CrossRef] [Google Scholar]
 Takeuchi, T., & Krauss, O. 2008, ApJ, 677, 1309 [NASA ADS] [CrossRef] [Google Scholar]
 Teiser, J., Engelhardt, I., & Wurm, G. 2011, ApJ, 742, 5 [NASA ADS] [CrossRef] [Google Scholar]
 van Eymeren, J., & Wurm, G. 2012, MNRAS, 420, 183 [NASA ADS] [CrossRef] [Google Scholar]
 Weidenschilling, S. J. 1977, Ap&SS, 51, 153 [NASA ADS] [CrossRef] [Google Scholar]
 Wurm, G., & Haack, H. 2009, Meteorit. Planet. Sci., 44, 689 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Wurm, G., & Krauss, O. 2006, Icarus, 180, 487 [NASA ADS] [CrossRef] [Google Scholar]
 Wurm, G., Krauss, O., & Haack, H. 2006, Meteorit. Planet. Sci. Suppl., 41, 5115 [NASA ADS] [Google Scholar]
 Wurm, G., Teiser, J., Bischoff, A., Haack, H., & Roszjar, J. 2010, Icarus, 208, 482 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1 Model setup of a chondrule core (high thermal conductivity) and a dust mantle (low thermal conductivity). The light source is directed downwards. Due to rotational symmetry only the temperature distribution along one line from top (180°) to the bottom (0°) has to be considered. 

Open with DEXTER  
In the text 
Fig. 2 Temperature distributions for homogeneous spheres of 1.12 mm diameter with varying thermal conductivities. 

Open with DEXTER  
In the text 
Fig. 3 Ratio for a given particle size and thermal conductivity k. For thermal conductivities larger than k = 0.2 the factor is essentially constant z^{(1)} = 0.7 (α = 1, this plot) (resp. z^{(1)} = 0.8 for α = 0.8). 

Open with DEXTER  
In the text 
Fig. 4 Ratio for a given particle size and thermal conductivity k. For thermal conductivities larger than 2 the factor is essentially constant. The accommodation coefficient α is 1. 

Open with DEXTER  
In the text 
Fig. 5 Ratio for a given particle radius and thermal conductivity k. For thermal conductivities larger than 0.2 the factor can be regarded as constant. The accommodation coefficient α is 1. 

Open with DEXTER  
In the text 
Fig. 6 κ over d_{dust} and k_{dust} for a 1 mm sized core with a thermal conductivity of 4 W/mK. The accommodation coefficient α is 1; white lines indicate constant values and show multiple configurations can lead to the same κ. 

Open with DEXTER  
In the text 
Fig. 7 Effective thermal conductivity κ over dust mantle to core size ratio. The 3 different lines represent 3 different core sizes: from bottom to top r_{core} = 0.2,0.5,1.0 mm. The thermal conductivity of the dust mantle is 0.01 W/mK, the core’s thermal conductivity is 1 W/mK and the accommodation coefficient α is 1. 

Open with DEXTER  
In the text 
Fig. 8 κ over d_{dust} and k_{dust} for a 1 mm sized core with a thermal conductivity of 1 W/mK (bottom) and 4 W/mK (top). The accommodation coefficient α is 1. 

Open with DEXTER  
In the text 
Fig. 9 κ for different pairs . Two bottom lines: k_{dust} = 0.01 W/mK, two top lines: k_{dust} = 0.1 W/mK 

Open with DEXTER  
In the text 
Fig. 10 Principles of sorting dust mantled chondrules by photophoresis. 

Open with DEXTER  
In the text 