Access the full text.
Sign up today, get DeepDyve free for 14 days.
Loading next page...
/lp/arxiv-cornell-university/the-morphological-elastic-and-electric-properties-of-dust-aggregates-r4IDwaIO5r
References
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
- ISSN
- 0032-0633
- eISSN
- ARCH-3330
- DOI
- 10.1016/j.pss.2019.104825
- Publisher site
- See Article on Publisher Site
Abstract
To appear in Planetary and Space Science (doi:10.1016/j.pss.2019.104825) Preprint typeset using L T X style AASTeX6 v. 1.0 THE MORPHOLOGICAL, ELASTIC, AND ELECTRIC PROPERTIES OF DUST AGGREGATES IN COMETS: A CLOSE LOOK AT COSIMA/ROSETTA'S DATA ON DUST IN COMET 67P/CHURYUMOV-GERASIMENKO Hiroshi Kimura Planetary Exploration Research Center (PERC), Chiba Institute of Technology, Tsudanuma 2-17-1, Narashino, Chiba 275-0016, Japan and Max Planck Institute for Solar System Research, Justus-von-Liebig-Weg 3, 37077 G ottingen, Germany Martin Hilchenbach, Sihane Merouane, John Paquette, and Oliver Stenzel Max Planck Institute for Solar System Research, Justus-von-Liebig-Weg 3, 37077 G ottingen, Germany hiroshi kimura@perc.it-chiba.ac.jp ABSTRACT The Cometary Secondary Ion Mass Analyzer (COSIMA) onboard ESA's Rosetta orbiter has revealed that dust particles in the coma of comet 67P/Churyumov-Gerasimenko are aggregates of small grains. We study the morphological, elastic, and electric properties of dust aggregates in the coma of comet 67P/Churyumov-Gerasimenko using optical microscopic images taken by the COSIMA instrument. Dust aggregates in COSIMA images are well represented as fractals in harmony with morphologi- cal data from MIDAS (Micro-Imaging Dust Analysis System) and GIADA (Grain Impact Analyzer and Dust Accumulator) onboard Rosetta. COSIMA's images, together with the data from the other Rosetta's instruments such as MIDAS and GIADA do not contradict the so-called rainout growth of 10 m-sized particles in the solar nebula. The elastic and electric properties of dust aggregates measured by COSIMA suggest that the surface chemistry of cometary dust is well represented as carbonaceous matter rather than silicates or ices, consistent with the mass spectra, and that organic matter is to some extent carbonized by solar radiation, as inferred from optical and infrared ob- servations of various comets. Electrostatic lofting of cometary dust by intense electric elds at the terminator of its parent comet is unlikely, unless the surface chemistry of the dust changes from a dielectric to a conductor. Our ndings are not in con ict with our current understanding of comet formation and evolution, which begin with the accumulation of condensates in the solar nebula and follow with the formation of a dust mantle in the inner solar system. Keywords: comets: general | comets: individual (67P/Churyumov-Gerasimenko) | meteorites, me- teors, meteoroids | protoplanetary disks 1. INTRODUCTION The formation and evolution of comets are a long-standing issue for planetary scientists, while a study on the physical and chemical properties of cometary dust provides an important clue to correct understanding of comets. In the 20th century, the most popular scenario for the formation of comets was the accumulation of presolar interstellar grains that were preserved in the solar nebula (e.g., Cameron 1975; Greenberg 1998). As a consequence, an expected picture of cometary dust is an aggregate particle of submicrometer-sized elongated interstellar grains consisting of a silicate core, an organic inner mantle, and an ice outer mantle (Greenberg & Hage 1990). Thanks to laboratory analyses of interplanetary dust particles (IDPs) collected in the stratosphere of the Earth, our knowledge of cometary dust has advanced over the decades (Brownlee 1985). There is consensus among experts on IDPs that a chondritic porous (CP) subset of IDPs is of cometary origin and bears compositionally strong resemblance to dust in comet 1P/Halley measured in situ (e.g., Jessberger 1999). The major constituent of CP IDPs is amorphous silicate called GEMS (glass embedded with metal and sul des) that is isotopically homogeneous, compositionally non-solar, and submicrometer- sized (Bradley 1994). GEMS grains were claimed to be of interstellar origin, because they share common properties with interstellar grains such as their submicrometer-size, amorphous structure, and infrared spectra (Bradley et al. arXiv:1912.13198v2 [astro-ph.EP] 6 Jan 2020 2 1999). However, it turned out from thorough laboratory examination of GEMS early this century that GEMS grains have complementary compositions to crystalline grains and thus both formed from the same reservoir of the solar nebula (Keller & Messenger 2011, 2013) . Accordingly, we are now aware that comets are made out of condensates in the solar nebula rather than presolar interstellar grains, contrary to the most popular scenario of comet formation in the 20th century. Therefore, thorough investigation of cometary dust will help us to understand the formation of pristine materials in the solar system and the early stages of planet formation. CP IDPs are aggregate particles consisting of submicrometer-sized Mg-rich crystalline silicates and GEMS glued together by carbonaceous material (Keller et al. 2000; Flynn et al. 2013). Such an aggregate structure is a natural consequence of dust growth by coagulation in the solar nebula, according to a model for the formation of comets (Weidenschilling 1997). It should be noted that the relatively low degree of uness of CP IDPs is not necessarily a typical example of cometary dust, because only compact endmembers of cometary dust are transported to the Earth by the Poynting-Robertson e ect (Kimura et al. 2016). Owing to a stronger radiation pressure force on a u y particle than a compact particle of the same mass, the former tends to have a diculty of staying in a bound orbit around the sun (cf. Kres ak 1976). The Deep Impact mission to comet 9P/Tempel 1 revealed that the radiation pressure on dust particles excavated by a projectile into the comet is strong enough to place the particles in a hyperbolic orbit (see Kobayashi et al. 2013). It is interesting to note that a dust mantle of a short-period comet is most likely composed of relatively compact aggregates with a radius of tens of microns (Li & Greenberg 1998; Kolokolova et al. 2007; Yamamoto et al. 2008; Kobayashi et al. 2013). Therefore, we may expect that in-situ measurements of cometary dust give new insights into not only the composition, but also the morphology of aggregate particles in comets. A large number of dust particles from comet 67P/Churyumov-Gerasimenko have been collected by the Cometary Secondary Ion Mass Analyzer (COSIMA) onboard the ESA's Rosetta orbiter (Langevin et al. 2016; Merouane et al. 2017). COSIMA is equipped with an optical microscope called COSISCOPE that acts as a powerful tool to investigate the morphology of dust particles and their collisional and electrostatic characteristics. We shall explore the data obtained by Rosetta/COSIMA and discuss them together with the other Rosetta's data from a theoretical point of view in terms of physical and chemical properties of cometary dust. 2. THEORETICAL BACKGROUNDS 2.1. Coagulation of Dust Aggregates in the Solar Nebula In the early stage of comet formation, dust aggregates grow under ballistic cluster-cluster aggregation (BCCA) process, because the motion of dust particles is controlled by Brownian motion where dust particles of similar size hit and stick each other (Weidenschilling et al. 1989). Once aggregate particles grow to the size of 10{10 m, they begin to settle toward the central plane of the solar nebula (\rainout") and inevitably sweep up smaller aggregate particles (Weidenschilling 1997). Therefore, we may expect that the rainout particles are well characterized by aggregate particles grown under ballistic particle-cluster aggregation (BPCA) process. A plausible scenario for the formation of planetesimals in the solar nebula requires that dust aggregates grew against mutual collisions even at a relative velocity of 50 m s at 30 au from the Sun (Weidenschilling & Cuzzi 1993; Weidenschilling 1997). It is not impossible for aggregates of silicate grains or aggregates of ice grains to stick each other at a collision velocity of 50 m s , but it is much easier for aggregates of organic grains to grow in the solar nebula (Kimura et al. 2015). As a result, it seems plausible that the surface of submicrometer-sized grains originally consists of organic matter when the grains grew to comets by coagulation in the solar nebula. 2.2. The Concept of Fractal Geometry It has been well known that the structure of aggregate particles grown under the BCCA and BPCA processes is well characterized by the concept of fractals (Meakin & Donn 1988). The number N of constituent grains with radius r and the characteristic radius r of a fractal aggregate de ne a fractal dimension D of the aggregate: N = k ; (1) where k is a proportionality constant of order unity. Hereafter we assume log k = 0:576D + 0:915 (Kimura et al. 0 0 1997). BCCA and BPCA particles are known to have their fractal dimensions D 2 and D 3, respectively, if their There are presolar GEMS grains with isotopic anomalies in CP IDPs, but they occupy only a few percent of total number of GEMS grains present in CP IDPs (Keller & Messenger 2011). This is also called the Hausdor dimension. 3 constituent grains are assumed to hit and stick on contact. In reality, the fractal dimension D of BCCA particles depends on the the relative velocity of mutual collisions where low velocities with no restructuring indicate D 2 and high velocities with the maximum compression result in D 2:5 (Wada et al. 2008). We adopt the following de nition for the porosity po of an aggregate particle proposed by Mukai et al. (1992): D 3 po = 1 k ; (2) with r = r ; (3) c g where r denotes the radius of gyration. The porosity of aggregate particles decreases with their size if D < 3, while the porosity of BPCA particles with D 3 is approximately constant at po 0:85, if their constituent grains hit and stick on contact (Mukai et al. 1992; Kimura et al. 2016). This is in accordance with po = 0:87 derived for the outermost dust mantle of comet 67P/Churyumov-Gerasimenko from the so-called Hapke's modeling of a spectral variation in the re ectance of the nucleus (Fornasier et al. 2015). 2.3. Fractal Dimensions The fractal dimension d in a two-dimensional (2-D) Euclidean plane is related to the fractal dimension D in three- dimensional (3-D) Euclidean space as follows (Meakin 1991): 2 if D 2; d = (4) D if D < 2: Therefore, for BCCA particles with D 2 and BPCA particles with D 3 in 3-D Euclidean space, their projected images onto a two-dimensional Euclidean plane are both characterized by a fractal dimension d 2. It is not straight- forward to derive the fractal dimension D in three-dimensional Euclidean space from two-dimensional projections of fractal aggregates (Maggi & Winterwerp 2004). We are aware that the fractal dimension of an aggregate particle slightly depends on the way how the fractal dimension is measured, although di erent dimensions could be generalized. It should be noted that the correlation dimension d might be smaller than the box-counting dimension d , which is de ned by 2 0 log n(") d = ; (5) log (1=") where " is the box size and n(") is the minimum number of boxes that is required to cover the projected area of the aggregate particle. The box-counting dimension is also called the capacity dimension and usually equivalent to the Hausdor dimension, namely, d = d in 2-D Euclidean space. Accordingly, aggregate particles with d < 2 in a 2-D 0 2 Euclidean plane are not necessarily associated with a low fractal dimension of D < 2 in 3-D Euclidean space. 2.4. Processing of Dust Aggregates in a Dust Mantle When comets evolve under solar radiation, the surface of the comets is processed and form the so-called dust mantle, a layer of dust particles devoid of volatile ices (Prialnik et al. 2004; Yamamoto et al. 2008). Optical and infrared photo-polarimetric observations of dust in comets revealed that dust paricles in short-period comets are physically and chemically processed during the formation of a dust mantle (Kolokolova et al. 2007). Infrared spectra of cometary dust suggest that the organic refractory component of the dust is not intact, but is to some extent carbonized, resulting in the formation of amorphous carbon (Li & Greenberg 1998; Kimura 2014) . Therefore, it is reasonable to assume that the surface of submicrometer-sized grains consists of amorphous carbon to a large extent rather than pristine organic matter. During the formation of a dust mantle on the surface of comets, dust particles with a high porosity, namely, a low fractal dimension are selectively ejected from the surface by gas drag, owing to their high cross-section to mass ratios. Moreover, dust particles with a low porosity, namely, a high fractal dimension tend to fall back on the surface of a comet and thus most likely elevate the fractal dimension of dust particles in the surface of the comet as a result of Here the term \carbonization" is used to indicate the loss of H, N, and O from the organic refractory component of comet dust by photolysis, radiolysis, and chemical reactions (see Roessler & Nebeling 1987; Jenniskens et al. 1993) 4 2.0 2.0 1.9 1.9 1.8 1.8 1.7 1.7 Estelle (COSIMA) Nick (COSIMA) d = 1.984 ± 0.005 d = 1.998 ± 0.003 0 0 1.6 1.6 2 4 6 8 2 4 6 8 2 4 6 8 2 4 6 8 2 4 6 8 2 4 6 8 0.01 0.1 1 10 0.01 0.1 1 10 1/ε 1/ε Figure 1. The box-counting dimension d (") (red circles) of the 100 m-sized aggregates Estelle (left) and Nick (right) versus the inverse of the maximum box size ". The tting curve (blue line) gives the box-counting dimension of Estelle and that of Nick, namely, the asymptotic values of lim d (") = 1:984 0:005 for Estelle and lim d (") = 1:998 0:003 for Nick. "!0 0 "!0 0 2.0 1.9 1.8 1.7 Particle E (MIDAS) d = 1.997 ± 0.012 1.6 2 4 6 8 2 4 6 8 2 4 6 8 0.01 0.1 1 10 1/ε Figure 2. The box-counting dimension d (") (red circles) of the 100 m-sized aggregate particle E (MIDAS) versus the inverse of the maximum box size ". The tting curve (blue line) gives the box-counting dimension of particle E, namely, the asymptotic value of lim d (") = 1:997 0:012. "!0 0 inelastic collisions. Therefore, it is reasonable to assume that dust particles in a dust mantle of a comet are composed of aggregate particles with relatively compact structures, in other words, high fractal dimensions, compared with those in the inner nucleus of the comet (see Yamamoto et al. 2008). 3. INTERPRETATION OF COSIMA'S DATA We shall explore what COSIMA's data as well as the other Rosetta's results tell us the morphological, elastic, and electric properties of cometary dust on the assumption that the morphology of the dust is well represented as a fractal. 3.1. Morphologies COSISCOPE images of dust particles collected on the target of COSIMA revealed that the particles are an assemblage of subunits with an apparent radius of 10 m (Hornung et al. 2016; Hilchenbach et al. 2017). Although COSISCOPE was not able to image submicrometer-sized grains expected as the constituent of dust aggregates in comets because of its resolution of 14 m, the subunits do not seem to have broken up into minute fragments. Therefore, COSIMA's data may indicate that the so-called rainout began at aggregate size of 10 m expected from a model for the formation of comets (see Weidenschilling 1997). By applying a box-counting algorithm to the COSISCOPE image of a 100 m-sized aggregates named Estelle and Nick shown in Langevin et al. (2016), we obtain a box-counting dimension d = 1:984 0:005 for the former and d = 1:998 0:003 for the latter (see Fig. 1) . The fractal dimension d ' 2 0 0 0 in 2-D Euclidean space indicates that at least Estelle and Nick are fractal aggregates whose fractal dimensions in 3-D Euclidean space lie in the range of D = 2{3 (see Eq. (4)). Therefore, COSISCOPE images are still in line with the We use the Java-based image processing program ImageJ to obtain a box-counting dimension with the smallest box size being one pixel of images, namely, 14 m. When we perform curve ttings to box-counting dimensions, we restrict the data points to small box sizes of 1=" 0:1. 0 5 rainout process in which aggregate particles are characterized by D 3, although d 2 does not necessarily mean D 3. Mannel et al. (2016) derived the correlation dimension of d 1:70 0:10 for a 10 m-sized particle E from its 2-D projected image using a topographic image of the particle taken by MIDAS (Micro-Imaging Dust Analysis System) onboard Rosetta. Since a correlation dimension d could be smaller than a Hausdor dimension d, it is not clear whether the particle E has a fractal dimension D < 2 in three-dimensional Euclidean space. They also determined a Hausdor dimension of the particle E to be d 1:76 0:29 in a two-dimensional Euclidean plane, but the large uncertainty prevents justi cation of the claim that the particle E has a smaller fractal dimension than d 2:0. Moreover, the number of grains in the particle E is only 112 in two-dimensional Euclidean space and the fractal dimension might be underestimated due to the small number of grains used in their analysis. We derive a box-counting dimension of d = 2:00 0:01 from the 2-D projected image of the particle E in Bentley et al. (2016), as shown in Fig. 2. Therefore, we may expect that the fractal dimension of the particle E in three-dimensional Euclidean space lies in the range of 2:0 D 3:0. Fulle et al. (2016) claimed that dust particles detected by GIADA (Grain Impact Analyzer and Dust Accumulator) on board Rosetta are extremely u y and characterized by D = 1:87 consistent with d 1:70 0:10 for the particle E of MIDAS. We should emphasize that dust particles detected only by GIADA/GDS are not necessarily characterized by D < 2, because there is no strong argument for the fractal dimension of D = 1:87 in Fulle et al. (2016). According to Fulle et al. (2015), Fulle et al. (2016) assumed N = 10 for GDS-only dust particles whose typical geometric cross-section A amounts to A = 8 10 m . However, we are aware that Fulle et al. (2015) did not provide solid 6 8 evidence for N = 10 and their arguments seem to be also consistent with D 3 on the assumption of N 6 10 and po 0:85. They attributed the fragile nature of the so-called GDS-only showers to the uness of aggregate particles, but instead one could attribute it to the carbonization of organic matter, as will be discussed in the next section. Consequently, the detection of the GDS-only showers by GIADA does not con ict with COSIMA's data as well as our current understandings of coagulation in the solar nebula and processing in a dust mantle. Rotundi et al. (2015) and Fulle et al. (2017) provide GIADA's data on the geometric cross section and the mass of each particle measured by both the GDS and the IS (Impact Sensor) sub-systems . Figure 3 compares the GIADA/GDS+IS data with the relationship between A and m for fractal aggregates grown under hit-and-stick coagulation process without restructuring of constituent grains (dotted line: BPCA particles; dash-dotted line: BCCA particles). The A-m relations for BCCA and BPCA particles formed under hit-and-stick coagulation processes without restructuring have been numerically determined as an analytical function of N (Ossenkopf 1993; Minato et al. 2006). It turns out that the GIADA/GDS+IS data seem not to be inconsistent with the geometrical cross section A of aggregate particles derived from their COSISCOPE images if the particles grow by the ballistic particle-cluster coagulation process (open squares). The slope of the GIADA's data are scattered around the A-m relation for BPCA particles with D 3 and po 0:85, although the GIADA's data tend to exceed the mass of fractal aggregates with the same geometric cross section. Because the slope of BCCA particles with D 2 is much gentler compared to the one for higher fractal dimensions, the GIADA/GDS+IS data do not seem to be in harmony with fractal aggregates with D < 3. The di erence in the slope between GIADA's data and BPCA particles without restructuring might indicate that GIADA's aggregate particles su ered from restructuring either when they grew by coagulation under the ballistic particle-cluster process in the solar nebula or when they formed a dust mantle on the surface of comet 67P/Churyumov-Gerasimenko. 3.2. Chemical Composition The elemental abundances of dust particles in comet 67P/Churyumov-Gerasimenko have been derived from the mass spectra of 30 particles measured by COSIMA onboard Rosetta (Bardyn et al. 2017). Figure 4 compares the elemental abundances of dust particles in comet 67P/Churyumov-Gerasimenko with those in comets 1P/Halley derived from 79 high quality mass spectra measured by PUMA-1 onboard VeGa-1 and those in 81P/Wild 2 from laboratory analyses of cometary matter in aerogel samples and Al foil residues returned by Stardust (Jessberger et al. 1988; Stephan et al. 2008b,a; Leroux et al. 2008; Cody et al. 2008). We also plot the elemental abundances of CP IDPs collected in the stratosphere and the solar photosphere, by considering that CP IDPs are likely of cometary origin and dust particles in comets are originally condensates from the nebular gas of solar composition (Schramm et al. 1989; Arndt et al. 1996; Keller et al. 1995; Asplund et al. 2009). Note that PUMA-1 measured the bulk composition of dust particles in This value was evaluated for showers of u y dust particles detected only by the GDS (Grain Detection System) sub-system of GIADA. While Fulle et al. (2017) analyzed the whole set of GIADA's data inclusive of the limited data analyzed by Rotundi et al. (2015), it turned out that the A and m values derived from the same data set by Fulle et al. (2017) and Rotundi et al. (2015) do not coincide each other, implying that the values are strongly dependent on their calibrations. 6 Characteristic radius (µm) 10 100 6 8 2 4 6 8 2 4 -4 GIADA: Rotundi et al. (2015) -5 GIADA: Fulle et al. (2017) COSIMA -6 -7 -8 BPCA -9 po ≈ 0.846 -10 -11 -12 BCCA po ≈ 0.989–0.999 -13 -10 -9 -8 -7 -6 10 10 10 10 10 Geometric cross section (m ) Figure 3. The relation between the geometric cross section A and the mass m of dust particles detected by both the GDS and the IS sub-systems of GIADA (Rotundi et al. 2015; Fulle et al. 2017). The dash-dotted and dotted lines are the A-m relations for ballistic cluster-cluster aggregates (BCCAs) and ballistic particle-cluster aggregates (BPCAs), respectively, formed under hit-and-stick coagulation process. Also plotted as squares are the A-m relation derived from COSISCOPE images of dust particles collected by COSIMA. The characteristic radius of a particle is intended to give an estimate of the particle size and 1=2 here it is given by (A=) with A being the geometric cross-section of the particle. 67P/C-G 1P/Halley 81P/Wild 2 CP IDPs Sun -1 -2 -3 -4 C N O Na Mg Al Si S K Ca Ti Cr Mn Fe Co Ni element Figure 4. The elemental abundances of dust particles in comet 67P/Churyumov-Gerasimenko normalized to the abundance of Mg ( lled squares) measured by COSIMA (Bardyn et al. 2017). Also plotted as lled circles, triangles, diamonds, and a solid line are the Mg-normalized elemental abundances of dust particles in comets 1P/Halley and 81P/Wild 2, CP IDPs, and the solar photosphere (Jessberger et al. 1988; Stephan et al. 2008b,a; Leroux et al. 2008; Cody et al. 2008; Schramm et al. 1989; Arndt et al. 1996; Keller et al. 1995; Asplund et al. 2009). comet 1P/Halley, most of which seem to consist of a chondritic core and an organic mantle (Kissel & Krueger 1987). In contrast, COSIMA is sensitive only to the surface matter of dust particles down to a depth of a few nanometers and thus the chemical composition of dust aggregates might be biased toward surface materials of their constituent grains. As a result, if chondritic materials are enclosed with carbonaceous matter, then COSIMA tends to trace the elemental abundances of the carbonaceous matter. On the one hand, the elemental abundances of dust particles in comet 67P/Churyumov-Gerasimenko measured by COSIMA have revealed that Si is markedly supra-solar, when normalized to Mg. On the other hand, one may realize elemental abundance relative to Mg Mass (kg) 7 that dust particles in comet 67P/Churyumov-Gerasimenko are depleted of Mg and Ca, when the elemental abundances are normalized to Si or Fe (cf. Bardyn et al. 2017). Here we should emphasize that COSIMA does not necessarily probe the bulk composition of the chondritic materials, even if chondritic materials are not covered by carbonaceous matter. A careful analysis of GEMS and pyroxene grains in CP IDPs has shown that the Si=Mg and Si=Ca atomic ratios are elevated near the surface of the grains, compared with the center of the grains (Bradley 1994). By analogy with CP IDPs, the Mg=Si and Ca=Si ratios in comet 67P/Churyumov-Gerasimenko may be lower at the surfaces of constituent grains in dust aggregates than the centers of the grains (see Appendix A). Therefore, we attribute the supra-solar Si=Mg ratio of dust particles in comet 67P/Churyumov-Gerasimenko to a systematic bias of COSIMA measurements in favor of surface materials. The elemental abundances provide clues about the volume fractions of organic-rich carbonaceous matter and rock- forming components such as Mg-rich silicates, metals, and sul des. Mann et al. (2004) have described how to derive the volume fractions of organic-rich carbonaceous matter and rock-forming components from elemental abundances. They found that the volume fraction of carbonaceous matter, f , and that of rock-forming matters, f , to CHON ROCK be f = 0:66 and f = 0:34 for dust particles in comet 1P/Halley (cf. Kimura et al. 2006). In the same CHON ROCK manner, we obtain f = 0:76 and f = 0:24 for dust particles in comet 67P/Churyumov-Gerasimenko using CHON ROCK the elemental abundances measured by COSIMA. These estimates highlight the overall dominance of carbonaceous matter over rock-forming matters on the surfaces of dust particles from comets 67P/Churyumov-Gerasimenko and 1P/Halley. We may estimate the mass fraction of carbonaceous matter, , and that of rock-forming matters, , from CHON ROCK their volume fractions, if their bulk densities are given. With a bulk density = 1800 kg m for carbonaceous CHON matter and = 3500 kg m for rock-forming matters, one gets 0:5 and 0:5 for dust particles ROCK CHON ROCK in the coma of comet 1P/Halley from their elemental abundances (Mann et al. 2004; Kimura et al. 2006). Fulle et al. (2018) obtained = 0:42 0:04 and = 0:58 0:04 for dust particles in the coma of comet 67P/Churyumov- CHON ROCK Gerasimenko from GIADA's data, by assuming a bulk density of = 1200 kg m for carbonaceous matter. If we CHON adopt = 1800 kg m instead, then we get = 0:48 0:07 and = 0:52 0:07 for GIADA's data, but CHON CHON ROCK = 0:6 and = 0:4 for COSIMA's data. The mass fractions derived from GIADA's data and COSIMA's CHON ROCK data are marginally consistent, but the di erences could be attributed to the fact that COSIMA tends to trace the carbonaceous matter, if it encloses rock-forming matter. As a result, there is no clear evidence that the mass fractions of carbonaceous matter and rock-forming matter in comet 67P/Churyumov-Gerasimenko greatly deviate from those in comet 1P/Halley. 3.3. Responses to a Collision Hereafter we shall demonstrate how the outcome of collision for aggregate particles with the target of the COSIMA instrument depends on the composition of the particles. For the sake of simplicity, we model cometary dust by aggregate particles consisting of N identical spherical grains (monomers) of radius r and volumetric mass density . The critical velocities for the onset of restructuring of an aggregate particle, v , the maximum compression of restr the aggregate, v , the onset of losing single spheres from the aggregate, v , and catastrophic disruption of the comp loss aggregate, v , are given by (Chokshi et al. 1993; Dominik & Tielens 1997; Wada et al. 2007) catastr 1=2 crit v = ; (6) restr r N 1=2 crit v = ; (7) comp " # 1=6 3 2 5 2 62:37 1 v = ; (8) loss 3 2 r E " # 1=6 3 2 5 2 207:9 1 v = ; (9) catastr 3 2 r E If an aggregate loses half the monomers or more upon collision with the target, the collisional outcome of the aggregate is referred to as catastrophic disruption. 8 where , E, and are the surface energy, Youngs modulus, and Poissons ratio, respectively. Here denotes the crit critical displacement that the contact area between the grains starts to move and we take = 0:2 nm (Dominik & crit Tielens 1995). Note that the above formulae for critical velocities not only demonstrate the outcomes of numerical simulations based on the contact mechanics, but also agree with laboratory experiments of aggregate collisions (Kimura et al. 2015) . Table 1 gives the physical properties of polyurethane (PUR), hydrogenated amorphous carbon, amorphous silica, and amorphous water ice, which are used to model organic matter, carbonaceous matter, silicate, and ice, respectively, in comets. Table 1. Physical properties of materials analogous to cometary matter Material E Reference 3 3 2 (10 kg m ) (J m ) (GPa) PUR 1:2 0:073 1:2 10 0:48 1:74 1, 2, 3 a-C:H 1:7 0:034 120 0:30 5:86 3, 4, 5 a-SiO 2:0 0:24 70 0:17 3:63 6, 7, 8 a-H O 1:0 0:11 7:0 0:25 107 7, 9, 10 Note|(1) P otschke et al. (2002); (2) McNicholas & Rankilor (1969); (3) Louh et al. (2005); (4) Piazza & Morell (2009); (5) Marques et al. (2003); (6) Kimura et al. (2015); (7) This study; (8) Henning & Mutschke (1997); (9) Wada et al. (2007); (10) Johari & Whalley (1981). Figure 5 shows the critical velocities v , v , v , and v as a function of the number of constituent grains restr comp loss catastr (monomers) in an aggregate particle consisting of polyurethane, hydrogenated amorphous carbon, amorphous silica, and amorphous water ice. We nd that aggregate particles do not experience any type of fragmentation at impact velocity below 475 m s if they are composed of organic matter, but the carbonization of aggregate particles reduces the impact velocity of fragmentation. It turns out that aggregate particles consisting of amorphous carbon grains with radius r = 0:1 m fragment if an impact velocity of the particles lies in the range of 2.2{4:1 m s . They are catastrophically disrupted at impact velocities above 4:1 m s before they experience maximum compression, because of v > v . Aggregate particles of silicate or water ice grains with radius r = 0:1 m may su er comp catastr 0 from restructuring, but not from fragmentation unless the impact velocity goes well beyond 10 m s . COSIMA's microscope images have revealed that aggregate particles of apparent diameter larger than 100 m fragmented into small parts upon collision with the collection plates of the instrument (Hornung et al. 2016). The impact velocities of the aggregate particles onto the plates have been estimated to be 3:5 1:5 m s (Rotundi et al. 2015). Therefore, our estimates on the outcome of collision support the idea that organic refractory component of cometary dust is carbonized during the formation of a dust mantle. Moreover, Fig. 5 indicates that the fragments of aggregate particles found in the COSIMA instrument cannot be consistent with aggregate particles of submicrometer-sized grains whose surfaces are dominated by silicate nor water ice nor a mixture of both. The carbon-rich nature of dust in the coma of comet 67P/Churyumov-Gerasimenko has been derived from the mass spectra of the dust by the COSIMA instrument (Bardyn et al. 2017). Therefore, we conclude that COSIMA's microscope images provide further evidence for the carbon-rich nature of the dust, indicating the carbonization of organic matter to a certain extent. 3.4. Responses to an Electric Field If an electric eld is applied to aggregate particles, then the particles might be separated into two or more subunits, depending on the strength of the applied eld and the composition of the particles (Kimura et al. 2014). By analogy with the Moon, Mendis et al. (1981) suggested that lofting and hovering of dust particles take place on the surfaces of comets when outgassing activities of comet nuclei cease at heliocentric distances greater than 5 au from the Sun. The rebound of aggregate particles sometimes observed in laboratory experiments has not well been modeled by numerical simulations on mutual collisions between aggregate particles (cf. Wada et al. 2011). It is, however, worthwhile noting that the rebound in the experiments by Ellerbroek et al. (2017) seems to take place only for aggregate particles larger than r 40 m, which exceeds the largest size (r 3 m) c c in the simulations. 9 Characteristic radius of aggregates (µm) Characteristic radius of aggregates (µm) 10 100 10 100 3 4 5 6 7 2 3 4 5 6 7 3 4 5 6 7 2 3 4 5 6 7 3 3 10 10 2 2 10 10 1 1 10 10 0 0 10 10 -1 -1 10 10 PUR a-C:H -2 -2 10 Onset of restructuring 10 Onset of restructuring Maximum compression Maximum compression Onset of losing monomers Onset of losing monomers Catastrophic disruption Catastrophic disruption -3 -3 10 3 4 5 6 7 10 3 4 5 6 7 10 10 10 10 10 10 10 10 10 10 Number of monomers Number of monomers Characteristic radius of aggregates (µm) Characteristic radius of aggregates (µm) 10 100 10 100 3 4 5 6 7 2 3 4 5 6 7 3 4 5 6 7 2 3 4 5 6 7 3 3 10 10 2 2 10 10 1 1 10 10 0 0 10 10 -1 -1 10 10 a-SiO a-H O 2 2 -2 -2 10 Onset of restructuring 10 Onset of restructuring Maximum compression Maximum compression Onset of losing monomers Onset of losing monomers Catastrophic disruption Catastrophic disruption -3 -3 10 10 3 4 5 6 7 3 4 5 6 7 10 10 10 10 10 10 10 10 10 10 Number of monomers Number of monomers Figure 5. The critical velocities for the onset of restructuring (blue), the maximum compression (green), the onset of losing monomers (red), and the catastrophic disruption (orange) of aggregate particles consisting of polyurethane (left top), hydro- genated amorphous carbon (right top), amorphous silica (left bottom), or amorphous water ice (right bottom) as a function of the number of monomers with radius r = 0:1 m. The characteristic radius of the particles given along the horizontal top axis corresponds to that for fractal aggregates with a fractal dimension of D = 2:5. Impact velocity (m/s) Impact velocity (m/s) Impact velocity (m/s) Impact velocity (m/s) 10 Dust aggregates in comets consist of submicrometer-sized grains whose surfaces are usually dominated by organic- rich carbonaceous material ( 50 wt%), in contrast to lunar agglutinates and Itokawa's samples (Greenberg 1982; Kissel & Krueger 1987; Kimura et al. 2003). This is also the case for CP IDPs of supposedly cometary origin where submicrometer-sized GEMS grains are typically encased in organic-rich carbonaceous matter (Keller et al. 2000). The ratio of organic to silicate components varies among the grains and the aggregates, but no single grain was found to consist of pure silicate nor organic components in dust from comet 1P/Halley (Jessberger et al. 1988). It is worthwhile noting that a force acting on conductive particles in an electric eld is repulsive, while a force acting on dielectric particles in an electric eld is attractive (Arp & Mason 1977; Nakajima & Matsuyama 2002). This dependence of electrostatic force on the grain material was con rmed by laboratory experiments on granular materials of 100 m sized grains (Holstein-Rathlou et al. 2012). Therefore, all forces acting on dust grains on the surfaces of airless bodies become attraction, if the outer layers of dust grains are covered by organic substrate, which is the most insulating material (cf. McCarty & Whitesides 2008). Nevertheless, electric charges on the surface of dust particles in an electric eld exert additional forces on the particles and thus might be able to lift up or disrupt the particles. The electrostatic repulsive force F on aggregate particles consisting of N dielectric spheres with radius r in an applied eld E is given el 0 by (cf. Nakajima & Sato 1999; Sow et al. 2013) 2 2 1 N Q 1 1 3 1 2 2 F = + 1 + NQE " r E ; (10) el 0 16" r 2 + 2 2 + 2 where Q is the average electric charge on each grain, is the static dielectric constant, and " is the permittivity of free space. Note that only the second term inside the braces of Eq. (10) corresponds to the repulsive force due to the electric eld and the rst and the third terms are attractive forces due to image charges and induced dipoles, respectively. The electrostatic repulsive force must overcome the adhesive force given by F = r (1 ) n ; (11) co 0 r c where (1 ) is de ned as the ratio of cohesive forces for a rough surface to a smooth one and n is the number of r c contacts for the grain at the point of disruption. We nd that numerical results of for a charged dielectric sphere on a 1:2 dielectric at surface by Nakajima & Sato (1999) are reproduced to within approximately 20% by 0:15 . As far 2 2 as aggregates consisting of identical dielectric spheres of N & 10 are concerned, we may assume 10 (Nakajima & Matsuyama 2002). The threshold of electric eld strength E , below which dust particles cannot be detached from th the surface by electric eld, can be determined by the balance between Eqs. (10) and (11): 8 9 " # 2 2=D 2 < = 2N ( + 2) ( + 1) Q 1 k " (1 ) n Q 0 0 r c E = 1 1 + : (12) th 2 2 2 2 4r : + 1 6 N 4r N 4r ; " ( 1) 0 0 0 0 Note that E in Eq. (12) may be underestimated, because the gravitational attractive force acting on the particles th from the surface must elevate the value of E . Figure 6 shows the dependence of E on the number N of dielectric th th spheres in an aggregate particle with a fractal dimension of D = 2:5. The particles are assumed to consist of organic matter (left top), hydrogenated amorphous carbon (right top), amorphous silica (left bottom), or amorphous water ice (right bottom). If the outer layer of cometary dust is covered by water ice and contains one elementary charge per grain on average or more, we expect that the dust cannot be lofted o the surfaces of short-period comets by electrostatic forces even at the terminator nor does hovering of grains take place. COSIMA's experiments on collected particles have observed fracture of an aggregate particle in an electric eld of 1:5 MV m after sucient exposure to indium ion beams (Hilchenbach et al. 2017). The experimental result was interpreted as evidence for accumulation of positive charges on the particle on the order of 75{300 C m by sequential irradiation of indium ion beams. To investigate whether the electrostatic repulsive force on insulators in the COSIMA's electric eld is strong enough to overcome the cohesive force, we shall estimate the minimum strength E of electric eld for the electrostatic lofting of aggregate particles. Because COSIMA collected dust particles by th their impacts on metal targets, we may assume = 1 in the rst term of Eq. (10). Figure 7 shows the dependence of the threshold electric eld E on the number N of dielectric spheres in an aggregate particle with D = 2:5. It turns th out that the results are insensitive to the choices of , , and because the image charge force (the rst term of Eq. (10)) is the most dominant attractive force. If positive charges on the order of 75{300 C m are accumulated on the surface of aggregate particles, then the attractive force is too large to lift up the particles from the metal targets. Photo-processing of carbonaceous material due to exposure to ultraviolet (UV) radiation reduces the band gap and 11 Characteristic radius of aggregates (µm) Characteristic radius of aggregates (µm) 10 100 10 100 3 4 5 6 7 2 3 4 5 6 7 3 4 5 6 7 2 3 4 5 6 7 5 5 10 10 4 4 10 10 Q = 10 e 3 3 Q = 10 e 10 10 terminator terminator 2 2 Q = 1 e 10 10 PUR a-C:H Q = 1 e φ = 0 φ = 0 r r φ = 0.99 φ = 0.99 r r 1 1 10 10 3 4 5 6 7 3 4 5 6 7 10 10 10 10 10 10 10 10 10 10 Number of monomers Number of monomers Characteristic radius of aggregates (µm) Characteristic radius of aggregates (µm) 10 100 10 100 3 4 5 6 7 2 3 4 5 6 7 3 4 5 6 7 2 3 4 5 6 7 5 5 10 10 Q = 10 e 4 4 10 10 Q = 10 e Q = 1 e 3 3 10 10 terminator terminator 2 Q = 1 e 2 10 10 a-SiO a-H O 2 2 φ = 0 φ = 0 r r φ = 0.99 φ = 0.99 r r 1 1 10 10 3 4 5 6 7 3 4 5 6 7 10 10 10 10 10 10 10 10 10 10 Number of monomers Number of monomers Figure 6. The threshold electric eld for lofting of aggregate particles consisting of dielectric grains as a function of the total number of grains (monomers) in an aggregate particle with fractal dimension D = 2:5. Left top: polyurethane; right top: hydrogenated amorphous carbon; left bottom: amorphous silica; right bottom: amorphous water ice. Here (1 ) is a reduction factor of cohesive forces due to the roughness of grain surfaces (solid line: = 0; long dashed line: = 0:99) and r r Q is the average electric charge per grain (red line: Q = 10 e; blue line: Q = 1 e). Also plotted as a dotted green line is the maximum strength of an electric eld that could be achieved at the terminator. 2 3 enhances the sp =sp ratio (Jones 2012). Consequently, dielectric properties may be transformed into semimetallic properties by graphitization, as the carbonization of organic matter proceeds with exposure to UV radiation. Therefore, an external electric eld might be able to lift up aggregate particles of highly carbonized grains, but not aggregate particles of less carbonized grains. COSIMA's experiments on collected particles have shown electrostatic lofting of some particles when an electric eld of 1:5 MV m is applied to the particles (Hilchenbach et al. 2017). COSIMA's results may indicate that the lofted particles were highly carbonized by solar UV radiation on the surface of the comet or even graphitized to some extent. To investigate whether the electrostatic repulsive force on conductors in the COSIMA's electric eld is strong enough to overcomes the cohesive force, we shall estimate the minimum strength E th of electric eld for the electrostatic lofting of conducting aggregate particles. If the aggregate particles consist of N spherical conducting monomers having radius r , then E is given by (Kimura et al. 2014) 0 th " s # r (1 ) n 3 (1 ) n 0 r c r c E = min ; ; (13) th 1 1 4e (3) + 8" (3) + r 0 0 6 6 Threshold electric field (kV/m) Threshold electric field (kV/m) Threshold electric field (kV/m) Threshold electric field (kV/m) 12 Characteristic radius of aggregates (µm) Characteristic radius of aggregates (µm) 10 100 10 100 3 4 5 6 7 2 3 4 5 6 7 3 4 5 6 7 2 3 4 5 6 7 5 5 10 10 –2 4 4 –2 300 µC m 10 10 300 µC m –2 –2 75 µC m 75 µC m 3 3 10 10 COSIMA COSIMA 2 2 10 10 PUR a-C:H φ = 0 φ = 0 r r φ = 0.99 φ = 0.99 r r 1 1 10 10 3 4 5 6 7 3 4 5 6 7 10 10 10 10 10 10 10 10 10 10 Number of monomers Number of monomers 2 2 Figure 7. The same as Fig. 6, but in the case for the surface charge densities of 75 C m (blue line) and 300 C m (red line) on the surfaces of polyurethane (left) and hydrogenated amorphous carbon (right). Also plotted as a dotted line is the strength of an electric eld that was applied to cometary dust by COSIMA's experiments. Characteristic radius of aggregates (µm) 10 100 3 4 5 6 7 2 3 4 5 6 7 COSIMA graphite φ = 0 φ = 0.99 3 4 5 6 7 10 10 10 10 10 Number of monomers Figure 8. The threshold electric eld for lofting of aggregate particles consisting of graphite (representing amorphous carbon with zero band gap) grains with surface roughness of = 0 (solid line) and = 0:99 (long dashed line) as a function of the r r total number of grains (monomers) in an aggregate particle with fractal dimension D = 2:5. Also plotted as a dotted line is the strength of an electric eld that was applied to cometary dust by COSIMA's experiments. where e is the elementary charge and is Riemann's zeta function. Figure 8 depicts the dependence of E on the size of th the particles, on the assumption that the surface energy of amorphous carbon with zero band gap is equivalent to that 2 1 of graphite (i.e., = 0:165 J m ) (cf. Abrahamson 1973). We nd that the strength of electric eld of 1:5 MV m is sucient to loft aggregates with characteristic radius r & 9 m. However, it should be noted that E in Eq. (13) c th might be overestimated by more than an order of magnitude, because the rugosity of each monomer enhances the induction charge and reduces the cohesion. Therefore, we expect that the electrostatic lofting of aggregate particles could be observed in COSISCOPE even for the smallest size of the particles, if their surfaces are heavily carbonized. This may, however, not be the case for most of the particles, because the electrostatic lofting of aggregate particles in COSIMA's experiments was assisted with indium ion beams. Because the molecular structure of a surface is easily modi ed by adsorption of atoms and molecules on the surface, Threshold electric field (kV/m) Threshold electric field (kV/m) Threshold electric field (kV/m) 13 Characteristic radius of aggregates (µm) 10 100 3 4 5 6 7 2 3 4 5 6 7 COSIMA a-C:H+In φ = 0 φ = 0.99 3 4 5 6 7 10 10 10 10 10 Number of monomers Figure 9. The same as Fig. 8, but for hydrogenated amorphous carbon covered with a thin layer of indium. an insulating surface could behave like a conducting one, depending on the adsorbed atoms and molecules on the surface (Voorthuyzen et al. 1987; Kimura et al. 2014). This is also the case for implantation of indium ions on the surface of an insulator that results in the formation of a thin metallic lm on the surface of the insulator (Yoshimura et al. 2010). Therefore, it is most likely that once a dose of indium ion radiation exceeds a threshold, then even aggregate particles of less carbonized grains could behave like conductors and be lofted o by an applied electric eld. In Figure 9, we plot the size dependence of E for aggregate particles of hydrogenated amorphous carbon grains with th indium coating, using = 0:37 J m (see Appendix C). It turns out that E is low enough for aggregate particles th of hydrogenated amorphous carbon grains with indium coating to be lofted o the COSIMA's target by an electric eld applied in COSIMA. Therefore, we expect that the electrostatic lofting of aggregate particles could be observed in COSISCOPE, irrespective of their size, if their surfaces are exposed to suciently high doses of indium ions. Eq. (13) indicates that the larger the size of aggregate particles is, the easier the electrostatic lofting of the particles is, in the case of conducting monomers. Figure 10 depicts the size distribution of dust particles derived from COSISCOPE images, the number of the particles lofted o the target surface in an applied electric eld, and the fraction of the lofted particles in each size bin. Dotted lines show tting lines in logarithmic scales and the error bars for the fraction of lofted particles are computed from square roots of the numbers for all particles and lofted particles. The fraction of lofted particles tend to increase with the size of the particles, in accordance with responses of electrically conductive particles to an electric eld given in Eq. (13). This proves that the surface layer of dust particles imaged by COSISCOPE is at least in part composed of conductive material rather than dielectric material. 4. SUMMARY We have examined the morphological, elastic, and electric properties of cometary dust using optical microscopic images of the dust taken by COSISCOPE of the COSIMA instrument onboard Rosetta. Cometary dust is an aggregate of subunits, the spatial arrangement of which is characterized by the concept of fractals and consistent with the so- called rainout growth of 10 m-sized particles in the solar nebula. Electric responses of aggregates to an electric eld in COSIMA indicate that the surfaces of the aggregates are dominated by dielectric materials, while their elastic responses to a collision onto the COSIMA target favor the surfaces covered by carbonaceous matter rather than silicates and ices. Since the subunits of aggregates were most likely glued together by organic matter during the coagulation growth of aggregates in the solar nebula, COSIMA's experiments imply that organic matter was carbonized during the formation of a dust mantle. Consequently, COSISCOPE images of cometary dust have embodied the most plausible scenario for the formation and evolution of comets, in which solar nebula condensates conglomerated into comets and solar radiation alters the surface structure and composition of comets in the inner solar system. We would like to thank Marco Fulle for useful discussion on the data for the cross section and mass of dust particles taken by GIADA onboard Rosetta. We also thank both the reviewers for their careful readings and bene cial comments Threshold electric field (kV/m) 14 0.1 0.01 all lofted 2 3 4 5 6 7 8 9 10 100 Characteristic radius (µm) Figure 10. The size distribution of dust particles derived from COSISCOPE images (red) and the number of the particles lofted o the target surface in an applied electric eld (blue) in the lower panel, and the fraction of the lofted particles (green) in the upper panel. on the manuscript. We acknowledge support from the Faculty of the European Space Astronomy Centre (ESAC), contract number PO 5001018751. H.K. is grateful to a stipend from the Max Planck Institute for Solar System Research and JSPS's Grants-in-Aid for Scienti c Research (KAKENHI #23244027, #26400230, #15K05273, #19H05085). Number of particles Fraction of lofted particles 15 67P/C-G GEMS (B) GEMS (C) Pyroxene Sun -1 -2 -3 -4 O Mg Al Si Ca Cr Mn Fe element Figure A1. The elemental abundances of dust particles in comet 67P/Churyumov-Gerasimenko normalized to the abundance of Si ( lled squares) measured by COSIMA (Bardyn et al. 2017). Also plotted as lled circles, triangles, and diamonds are the Si-normalized elemental abundances of two GEMS grains and pyroxenes at their surfaces (Bradley 1994). The solar photospheric abundances are depicted as a solid line, after they are normalized to Si (Asplund et al. 2009). APPENDIX A. DEPLETION OF MAGNESIUM AND CALCIUM FROM THE SURFACE OF CP IDPS The compositional pro les of two GEMS grains in CP IDPs revealed a gradual variation in the chemical composition from the core to the surface (Bradley 1994). Figure A1 compares the abundances of major GEMS-forming elements between dust particles in comet 67P/Churyumov-Gerasimenko ( lled squares) and the surfaces of two GEMS grains ( lled circles and triangles). Also plotted as lled diamonds and solid lines are the elemental abundances measured at the surface of a low-Fe pyroxene crystal in CP IDPs and the solar photospheric abundances, respectively (Bradley 1994; Asplund et al. 2009). If we take geometric means of the Mg=Si and Ca=Si atomic ratios over the GEMS and the +0:55 +0:0141 pyroxene, we obtain Mg=Si = 0:12 and Ca=Si = 0:0087 for the mineral surfaces. Since the corresponding 0:09 0:0054 +0:17 +0:0072 values of 67P/Churyumov-Gerasimenko are Mg=Si = 0:11 and Ca=Si = 0:0053 , the depletions of Mg and 0:08 0:0035 Ca relative to Si found in comet 67P/Churyumov-Gerasimenko by COSIMA seem to be consistent with the chemical compositions of IDP's mineral surfaces. B. ESTIMATES OF SURFACE ENERGY FROM SURFACE TENSION In spite of its importance, the surface energy of a solid that is composed of a material of astronomical interest is not always available in the literature. It is not often that we can obtain experimental data on the surface energy of a solid, if the solid matter is a liquid or a gas under Earth's ambient pressure and temperature. When experimental data on the surface energy are unavailable, we may derive the surface energy of a solid from the surface tension of a liquid using the following relationship (Ip & Toguri 1994) = T ; (B1) d T where T is the temperature of the liquid. The surface tension usually decreases with the temperature (i.e., d =d T < 0) and vanishes at the critical temperature T . As long as the surface tension is expressed by a linear equation = a +b T with a and b being t coecients, we have = at T = 0 (E otv os 1886; Palit 1956). To check the validity of Eq. (B1), we plot the temperature dependences of surface tension measured (solid lines) for (a) SiO and (b) H O 2 2 as well as their linear extrapolations (dashed lines) to lower temperature in Fig. B2 (Janz et al. 1969; International Association for the Properties of Water Steam 1994). Interestingly, the temperature coecient of SiO is positive (i.e., d =d T > 0) in the range of 1773 K < T < 2080 K, but the surface tension is expected to decrease with temperature at elemental abundance relative to Si 16 0.32 0.12 SiO 0.11445 ± 0.00007 0.30 0.10 0.28 0.08 0.26 0.06 0.24308 ± 0.00012 0.24 0.04 0.22 0.02 H O 0.20 0.00 0 500 1000 1500 2000 2500 0 200 400 600 Temperature (K) Temperature (K) Figure B2. The dependence of surface tension on the temperature for amorphous SiO (left) and H O (right). Open circles: 2 2 experimental data by Parikh (1958); closed circles: experimental data by Kingery (1959); dotted line: a model by Kraus et al. (2012); solid curves: compilations of data by Janz et al. (1969) for amorphous SiO (left) and International Association for the Properties of Water Steam (1994) for H O (right); dashed lines: linear extrapolations of the compiled data. higher temperature (T & 2000 K) and vanish at the critical temperature T (Kraus et al. 2012) . The surface energies 2 2 = 243:08 0:12 mJ m for SiO and = 114:45 0:07 mJ m for H O ice derived from the extrapolations SiO 2 H O 2 2 2 2 2 agree with the values in the literature: 250 mJ m for amorphous silica (Kimura et al. 2015); 110 mJ m for water ice (Israelachvili 2011). C. EFFECT OF ADSORPTION ON SURFACE ENERGY When the surface of a bulk is covered by a thin layer of material, then the surface energy can be approximately given by the equation that is expressed as (Israelachvili 1972, 2011): p p = ( ) ; (C2) 1 2 where and are the surface energies of the bulk and the thin layer, respectively. To check the validity of Eq. (C2), 1 2 we consider the case that the surface of amorphous silica is covered by water molecules under atmospheric condition. If 2 2 we insert = and = into Eq. (C2), we obtain = 0:024 J m , which is consistent with = 0:025 J m 1 SiO 2 H O 2 2 measured for hydrophilic amorphous silica under ambient conditions (Kendall et al. 1987). Accordingly, we may apply Eq. (C2) to the determination of the surface energy for carbonaceous matter covered by a thin layer of indium. On the assumption that the surface energy of indium is = 0:631 J m , Eq. (C2) with = and = gives In 1 a-C:H 2 In = 0:37 J m for indium-covered hydrogenated amorphous carbon (cf. Alchagirov et al. 2001, 2014). Kraus et al. (2012) adopted T = 5130 K for SiO , but the critical temperature of SiO is still open to debate in the range of c 2 2 T = 4862{13500 K (Iosilevskiy et al. 2013). Surface tension (N/m) Surface tension (N/m) 17 REFERENCES Abrahamson, J. 1973, Carbon, 11, 337. Ip, S. W., & Toguri, J. M. 1994, J. Mater. Sci., 29, 688. doi:10.1016/0008-6223(73)90075-4. doi:10.1007/BF00445980. Alchagirov, A. B., Alchagirov, B. B., Taova, T. M., Khokonov, Israelachvili, J. N. 1972, Proc. R. Soc. Lond. A, 331, 39. K. B.. 2001, Transactions of JWRI, 30, 287. doi:10.1098/rspa.1972.0163. Alchagirov, B. B., Dadashev, R. K., Dyshekova, F. F., & Israelachvili, J. N. 2011, Intermolecular and Surface Forces (3rd Elimkhanov, D. Z. 2014, High Temp., 52, 920. ed.; Boston, MA: Academic). doi:10.1134/S0018151X14060017. Janz, G. J., Dijkhuis, C. G. M., Lakshminarayanan, G. R., Arndt, P., Bohsung, J., Maetz, M., & Jessberger, E. K. 1996, Tomkins, R. P. T., & Wong, J. 1969, Molten Salts Meteorit. Planet. Sci., 31, 817. (NSRDS-NBS 28), (Washington, D.C.: NIST). doi:10.1111/j.1945-5100.1996.tb02116.x. Jenniskens, P., Baratta, G. A., Kouchi, A., de Groot, M. S., Arp, P. A., & Mason, S. G. 1977, Colloid. Polym. Sci., 255, 566. Greenberg, J. M., & Strazzulla, G. 1993, A&A, 273, 583. doi:10.1007/BF01549744. Jessberger, E. K. 1999, Space Science Reviews, 90, 91. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, doi:10.1023/A:1005233727874. ARA&A, 47, 481. doi:10.1146/annurev.astro.46.060407.145222. Jessberger, E. K., Christoforidis, A., & Kissel, J. 1988, Nature, Bardyn, A., Baklouti, D., Cottin, H., et al. 2017, MNRAS, 469, 332, 691. doi:10.1038/332691a0. S712. doi:10.1093/mnras/stx2640. Johari, G. P., & Whalley, E. 1981, JChPh, 75, 1333. Bentley, M. S., Schmied, R., Mannel, T., et al. 2016, Nature, doi:10.1063/1.442139. 537, 73. doi:10.1038/nature19091. Jones, A. P. 2012, A&A, 540, A2. Bradley, J. P. 1994, Science, 265, 925. doi:10.1051/0004-6361/201117624. doi:10.1126/science.265.5174.925. Keller, L. P., & Messenger, S. 2011, GeoCoA, 75, 5336. Bradley, J. P., Keller, L. P., Snow, T. P., et al. 1999, Science, doi:10.1016/j.gca.2011.06.040. 285, 1716. doi:10.1126/science.285.5434.1716. Keller, L. P., & Messenger, S. 2013, GeoCoA, 107, 341. Brownlee, D. E. 1985, Annual Review of Earth and Planetary doi:10.1016/j.gca.2012.11.012. Sciences, 13, 147. doi:10.1146/annurev.ea.13.050185.001051. Keller, L. P., Thomas, K. L., Bradley, J. P., & McKay, D. S. Cameron, A. G. W. 1975, Icarus, 24, 128. 1995, Meteoritics, 30, 526. doi:10.1016/0019-1035(75)90167-0. Keller, L. P., Messenger, S., & Bradley, J. P. 2000, Chokshi, A., Tielens, A. G. G. M., & Hollenbach, D. 1993, ApJ, J. Geophys. Res., 105, 10397. doi:10.1029/1999JA900395. 407, 806. doi:10.1086/172562. Kendall, K., Alford, N. M., & Birchall, J. D. 1987, Nature, 325, Cody, G. D., Ade, H., Alexander, C. M. O' D., et al. 2008, 794. doi:10.1038/325794a0. Meteorit. Planet. Sci., 43, 353. Kimura, H. 2014, Icarus, 232, 133. doi:10.1111/j.1945-5100.2008.tb00627.x. doi:10.1016/j.icarus.2014.01.009. Dominik, C., Tielens, A. G. G. M. 1995, PMagA, 72, 783. Kimura, H., Ishimoto, H., & Mukai, T. 1997, A&A, 326, 263. doi:10.1080/01418619508243800. Kimura, H., Kolokolova, L., & Mann, I. 2003, A&A, 407, L5. Dominik, C., & Tielens, A. G. G. M. 1997, ApJ, 480, 647. doi:10.1051/0004-6361:20030967. doi:10.1086/303996. Kimura, H., Kolokolova, L., & Mann, I. 2006. A&A, 449, 1243. Ellerbroek, L. E., Gundlach, B., Landeck, A., et al. 2017, doi:10.1051/0004-6361:20041783. MNRAS, 469, S204. doi:10.1093/mnras/stx1257. Kimura, H., Senshu, H., & Wada, K. 2014, Planet. Space Sci., E otv os, R. 1886, Ann. Phys., 263, 448. 100, 64. doi:10.1016/j.pss.2014.03.017. doi:10.1002/andp.18862630309. Kimura, H., Wada, K., Senshu, H., & Kobayashi, H. 2015, ApJ, Flynn, G. J., Wirick, S., & Keller, L. P. 2013, Earth Planet 812, 67. doi:10.1088/0004-637X/812/1/67. Space, 65, 1159. doi:10.5047/eps.2013.05.007. Kimura, H., Kolokolova, L., Li, A., & Lebreton, J. 2016, in Light Fornasier, S., Hasselmann, P. H., Barucci, M. A., et al. 2015, Scattering Reviews Vol. 11, ed. A. Kokhanovsky (Berlin, A&A, 583, A30. doi:10.1051/0004-6361/201525901. Springer), 363, doi:10.1007/978-3-662-49538-4 8. Fulle, M., Della Corte, V., Rotundi, A., et al. 2015, ApJ, 802, Kingery, W. D. 1959, J. Am. Ceram. Soc., 42, 6. L12. doi:10.1088/2041-8205/802/1/L12. doi:10.1111/j.1151-2916.1959.tb09134.x. Fulle, M., Altobelli, N., Buratti, B., et al. 2016, MNRAS, 462, Kissel, J., & Krueger, F. R. 1987, Nature, 326, 755. S2. doi:10.1093/mnras/stw1663. doi:10.1038/326755a0. Fulle, M., Della Corte, V., Rotundi, A., et al. 2017, MNRAS, Kobayashi, H., Kimura, H., & Yamamoto, S. 2013, A&A, 550, 469, S45. doi:10.1093/mnras/stx983. A72. doi:10.1051/0004-6361/201220464. Fulle, M., Bertini, I., Corte, Della, V., et al. 2018, MNRAS, 476, Kolokolova, L., Kimura, H., Kiselev, N., & Rosenbush, V. 2007, 2835. doi:10.1093/mnras/sty464. A&A, 463, 1189. doi:10.1051/0004-6361:20065069. Greenberg, J. M. 1982, in Comets, ed. L. L. Wilkening (Tuscon, Kraus, R. G., Stewart, S. T., Swift, D. C., et al. 2012, AZ: Univ. Arizona Press), 131. J. Geophys. Res., 117, E09009. doi:10.1029/2012JE004082. Greenberg, J. M. 1998, A&A, 330, 375. Kres ak, L. 1976, Bulletin of the Astronomical Institutes of Greenberg, J. M., & Hage, J. I. 1990, ApJ, 361, 260. Czechoslovakia, 27, 35. doi:10.1086/169191. Langevin, Y., Hilchenbach, M., Ligier, N., et al. 2016, Icarus 271, Henning, T., & Mutschke, H. 1997, A&A, 327, 743. 76. doi:10.1016/j.icarus.2016.01.027. Hilchenbach, M., Fischer, H., Langevin, Y., et al. 2017, Phil. Leroux, H., Stroud, R. M., Dai, Z. R., et al. 2008, Meteorit. Trans. R. Soc. A, 375, 20160255. doi:10.1098/rsta.2016.0255. Planet. Sci., 43, 143. doi:10.1111/j.1945-5100.2008.tb00615.x. Holstein-Rathlou, von, C., Merrison, J. P., Brdstrup, C. F., & Li, A., & Greenberg, J. M. 1998, A&A, 338, 364. N rnberg, P. 2012, Icarus, 220, 1. Louh, S. P., Leu, I. C., & Hon, M. H. 2005, Diamond and Related doi:10.1016/j.icarus.2012.04.011. Materials, 14, 1000. doi:10.1016/j.diamond.2004.12.004. Hornung, K., Merouane, S., Hilchenbach, M., et al. 2016, Maggi, F., & Winterwerp, J. C. 2004, PhRvE, 69, 011405. Planet. Space Sci., 133, 63. doi:10.1016/j.pss.2016.07.003. doi:10.1103/PhysRevE.69.011405. International Association for the Properties of Water Steam Mann, I., Kimura, H., & Kolokolova, L. 2004, JQSRT, 89, 291. 1994, Proceedings of the 12th International Conference on the doi:10.1016/j.jqsrt.2004.05.029. Properties of Water and Steam, ed. H. J. White, Jr., J. V. Mannel, T., Bentley, M. S., Schmied, R., et al. 2016, MNRAS, Sengers, D. B. Neumann, & J. C. Bellows (Orlando: Begell 462, S304. doi:10.1093/mnras/stw2898. House), a107. Marques, F. C., Lacerda, R. G., Champi, A., et al. 2003, Appl. Iosilevskiy, I., Gryaznov, V., & Solov'ev, A. 2013, High Phys. Lett., 83, 3099. doi:10.1063/1.1619557. Temp.-High Press., 43, 227. 18 McCarty, L. S., & Whitesides, G. M. 2008, Angew. Chem. Int. Rotundi, A., Sierks, H., Della Corte, V., et al. 2015, Science, 347, Ed., 47, 2188. doi:10.1002/anie.200701812. aaa3905. doi:10.1126/science.aaa3905. McNicholas, J. B., & Rankilor, P. R. 1969, Strain, 5, 74. Sow, M., Widenor, R., Akande, A. R., et al. 2013, J. Braz. doi:10.1111/j.1475-1305.1969.tb01582.x. Chem. Soc., 24, 273. Meakin, P. 1991, Rev. Geophys., 29, 317. Schramm, L. S., Brownlee, D. E., Wheelock, M. M. 1989, doi:10.1029/91RG00688. Meteoritics, 24, 99. doi:10.1111/j.1945-5100.1989.tb00950.x. Meakin, P., & Donn, B. 1988, ApJ, 329, L39. doi:10.1086/185172. Stephan, T., Flynn, G. J., Sandford, S. A., Zolensky, M. E. Mendis, D. A., Hill, J. R., Houpis, H. L. F., & Whipple, E. C., 2008a, Meteorit. Planet. Sci., 43, 285. Jr. 1981, ApJ, 249, 787. doi:10.1086/159337. doi:10.1111/j.1945-5100.2008.tb00623.x. Merouane, S., Stenzel, O., Hilchenbach, M., et al. 2017, MNRAS, Stephan, T., Rost, D., Vicenzi, E. P., et al. 2008b, Meteorit. 469, S459. doi:10.1093/mnras/stx2018. Planet. Sci., 43, 233. doi:10.1111/j.1945-5100.2008.tb00619.x. Minato, T., K ohler, M., Kimura, H., Mann, I.,& Yamamoto, T. Voorthuyzen, J. A., Keskin, K., & Bergveld, P. 1987, Surf. Sci., 2006, A&A, 452, 701. doi:10.1051/0004-6361:20054774. 187, 201. doi:10.1016/S0039-6028(87)80132-2. Mukai, T., Ishimoto, H., Kozasa, T., Blum, J., & Greenberg, Wada, K., Tanaka, H., Suyama, T., Kimura, H., & Yamamoto, J. M. 1992, A&A, 262, 315. T. 2007, ApJ, 661, 320. doi:10.1086/514332. Nakajima, Y., & Sato, T. 1999, J. Electrostat., 45, 213. Wada, K., Tanaka, H., Suyama, T., Kimura, H., & Yamamoto, doi:10.1016/S0304-3886(98)00051-5. T. 2008, ApJ, 677, 1296. doi:10.1086/529511. Nakajima, Y., & Matsuyama, T. 2002, J. Electrostat., 55, 203. Wada, K., Tanaka, H., Suyama, T., Kimura, H., & Yamamoto, doi:10.1016/S0304-3886(01)00198-X. T. 2011, ApJ, 737, 36. doi:10.1088/0004-637X/737/1/36. Ossenkopf, V. 1993, A&A, 280, 617. Weidenschilling, S. J. 1997, Icarus, 127, 290. Palit, S. R. 1956, Nature, 177, 1180. doi:10.1038/1771180a0. doi:10.1006/icar.1997.5712. Parikh, N. M. 1958, J. Am. Ceram. Soc., 41, 18. Weidenschilling, S. J., & Cuzzi, J. N. 1993, in Protostars and doi:10.1111/j.1151-2916.1958.tb13497.x. Planets III, ed. E. H. Levy & J. I. Lunine (Tuscon, AZ: Univ. Piazza, F., & Morell, G. 2009, Diamond Relat. Mater., 18, 43. Arizona Press), 1031. doi:10.1016/j.diamond.2008.09.023. Weidenschilling, S. J., Donn, B., & Meakin, P. 1989, in The P otschke, P., Pionteck, J., & Stutz, H. 2002, Polymer, 43, 6965. Formation and Evolution of Planetary Systems, ed. H. A. doi:10.1016/S0032-3861(02)00612-2. Weaver, L. Danly, & S. Fall (Cambridge: Cambridge Univ. Prialnik, D., Benkho , J., & Podolak, M. 2004, in Comets II, ed. Press,), 131. M. C. Festou, H. U. Keller, & H. A. Weaver (Tucson: Univ. Yamamoto, S., Kimura, H., Zubko, E., Kobayashi, H., Wada, K., Arizona Press), 359. Ishiguro, M., & Matsui, T. 2008, ApJ, 673, L199. Roessler, K., & Nebeling, B. 1987, Proceedings of the doi:10.1086/527558. International Symposium on the Diversity and Similarity of Yoshimura, S., Hine, K., Kiuchi, M., et al. 2010, Appl. Surf. Sci., Comets, ESA SP-278, 509. 257, 192. doi:10.1016/j.apsusc.2010.06.063.
Journal
Astrophysics – arXiv (Cornell University)
Published: Dec 31, 2019