The Composition of Dust in Jupiter-Family Comets as Inferred from Infrared Spectroscopy

We review the composition of Jupiter-family comet dust as inferred from infrared spectroscopy. We find that Jupiter-family comets have 10 micron silicate emission features with fluxes roughly 20-25% over the dust continuum (emission strength 1.20-1.25), similar to the weakest silicate features in Oort Cloud comets. We discuss the grain properties that change the silicate emission feature strength (composition, size, and structure/shape), and emphasize that thermal emission from the comet nucleus can have significant influence on the derived silicate emission strength. Recent evidence suggests that porosity is the dominant parameter, although more observations and models of silicates in Jupiter-family comets are needed to determine if a consistent set of grain parameters can explain their weak silicate emission features. Models of 8 m telescope and Spitzer Space Telescope observations have shown that Jupiter-family comets have crystalline silicates with abundances similar to or less than those found in Oort Cloud comets, although the crystalline silicate mineralogy of comets 9P/Tempel and C/1995 O1 (Hale-Bopp) differ from each other in Mg and Fe content. The heterogeneity of comet nuclei can also be assessed with mid-infrared spectroscopy, and we review the evidence for heterogeneous dust properties in the nucleus of comet 9P/Tempel. Models of dust formation, mixing in the solar nebula, and comet formation must be able to explain the observed range of Mg and Fe content and the heterogeneity of comet 9P/Tempel, although more work is needed in order to understand to what extent do comets 9P/Tempel and Hale-Bopp represent comets as a whole.Comment: 21 pages, 4 figures, 2 tables. Accepted for publication in Planetary and Space Scienc

Comet Mineralogy as Inferred from Infrared Spectra of Comets

For most comets, infrared (IR) spectroscopy (remote sensing) is the method through which we diagnose the mineralogy and size distribution of dust in their comae. The shape and contrast of the IR spectral features depend on the particle size: optically active minerals (absorbing of visible and near-IR solar photons) and submicron solid grains or highly porous (> 90% vacuum) grains primarily contribute to the shapes of the observed resonances. Comet mineralogies typically are determined by fitting thermal emission models of ensembles of discrete mineral grains to observed IR spectral energy distributions. The absorptivities (Q-abs) and scattering efficiencies (Q-scat) of the discrete mineral grains are computed using Mie scattering, Maxwell-Garnet mixing, Discrete Dipole Approximation, and Multi-Layered Sphere codes. These techniques when applied to crystalline minerals, specifically olivine (Mg_x, Fe_1-x)2 Si04, x>0.9, require the use of ellipsoidal shaped particles with elongated axial ratios or hollow spheres to produce the shapes of the resonances observed both from comet comae and laboratory samples. The wavelength positions of the distinct resonances from submicron-radii crystalline silicates, as well as their thermal equilibrium temperatures, constrain the crystalline olivine to have a relatively high Mg-content (x>0.9, or Fo>90). Only resonances computed for submicron Mg-rich crystalline olivine and crystalline orthopyroxene match the observed IR spectral features. However, this has led to the interpretation that micron-radii and larger crystals are absent from comet comae. Furthermore, the mass fraction of silicate crystals is dependent upon whether just the submicron portion of the size distribution is being compared or the submicron crystals compare to the aggregates of porous amorphous silicates that are computationally tractable as porous spheres. We will discuss the Deep Impact results as examples of these challenges to interpreting mid-IR spectra of short-period comets

Characterization of the High-Albedo NEA 3691 Bede

Characterization of NEAs provides important inputs to models for atmospheric entry, risk assessment and mitigation. Diameter is a key parameter because diameter translates to kinetic energy in atmospheric entry. Diameters can be derived from the absolute magnitude, H(PA=0deg), and from thermal modeling of observed IR fluxes. For both methods, the albedo (pv) is important - high pv surfaces have cooler temperatures, larger diameters for a given Hmag, and shallower phase curves (larger slope parameter G). Thermal model parameters are coupled, however, so that a higher thermal inertia also results in a cooler surface temperature. Multiple parameters contribute to constraining the diameter. Observations made at multiple observing geometries can contribute to understanding the relationships between and potentially breaking some of the degeneracies between parameters. We present data and analyses on NEA 3691 Bede with the aim of best constraining the diameter and pv from a combination of thermal modeling and light curve analyses. We employ our UKIRT+Michelle mid-IR photometric observations of 3691 Bede's thermal emission at 2 phase angles (27&43 deg 2015-03-19 & 04-13), in addition to WISE data (33deg 2010-05-27, Mainzer+2011). Observing geometries differ by solar phase angles and by moderate changes in heliocentric distance (e.g., further distances produce somewhat cooler surface temperatures). With the NEATM model and for a constant IR beaming parameter (eta=constant), there is a family of solutions for (diameter, pv, G, eta) where G is the slope parameter from the H-G Relation. NEATM models employing Pravec+2012's choice of G=0.43, produce D=1.8 km and pv0.4, given that G=0.43 is assumed from studies of main belt asteroids (Warner+2009). We present an analysis of the light curve of 3691 Bede to constrain G from observations. We also investigate fitting thermophysical models (TPM, Rozitis+11) to constrain the coupled parameters of thermal inertia (Gamma) and surface roughness, which in turn affect diameter and pv. Surface composition can be related to pv. This study focuses on understanding and characterizing the dependency of parameters with the aim of constraining diameter, pv and thermal inertia for 3691 Bede

Crystalline Silicates in Comets: Modeling Irregularly-Shaped Forsterite Crystals and Its Implications on Condensation Conditions

Crystalline silicates in comets are a product of the condensation in the hot inner regions (T > or approx. equals 1400 K [1]) of our proto-planetary disk or annealing at somewhat lower temperatures (T > or approx. equals 1000-1200 K) [2, 3, 4] in shocks coupled with disk evolutionary processes that include radial transport of crystals from their formation locations out to the cold outer regions where comet nuclei formed. The grain shape of forsterite (crystals) could be indicative of their formation pathways at high temperatures through vapor-solid condensation or at lower temperatures through vapor-liquid-solid formation and growth [5, 6, 7]. Experiments demonstrate that crystals that formed from a rapidly cooled highly supersaturated silicate vapor are characterized by bulky, platy, columnar/needle and droplet shapes for values of temperature and supersaturation, T and sigma, of 1000-1450 C and 230, respectively [7]. The experimental columnar/needle shapes, which form by vapor-liquid-solid at lower temperatures (<820 C), are extended stacks of plates, where the extension is not correlated with an axial direction: columnar/needles may be extended in the c-axis or a-axis direction, can change directions, and/or are off-kilter or a bit askew extending in a combination of the a- and c-axis direction

Lunar Crater Observation and Sensing Satellite (LCROSS) Instrument Calibration Summary

This document describes the calibration of the LCROSS instruments. It will be released to the public via the Planetary Data System. We need a quick review, if possible, because the data has been delivered to the PDS, and this document is needed to interpret the LCROSS impact data fully. [My mistake [shirley) in not realizing this needed to be treated as a normal publication.] The LCROSS instruments are commercially available units except for one designed and built at Ames. The commercially available instruments don't seem to me to present ITAR issues (Sony video camera, thermal camera from England, and so on.) Also, the internal design details of the instruments are not included in this report, only the process of calibrating them against standard targets. Only very high-level descriptions of the spacecraft are included, comparable to the level of detail included in the public web pages on nasa.gov

Spitzer Observations of Comet 67P/Churyumov-Gerasimenko at 5.5-4.3 AU From the Sun

We report Spitzer Space Telescope observations of comet 67P/Churyumov-Gerasimenko at 5.5 and 4.3 AU from the Sun, post-aphelion. Comet 67P is the primary target of the European Space Agency's Rosetta mission. The Rosetta spacecraft will rendezvous with the nucleus at heliocentric distances similar to our observations. Rotationally resolved observations at 8 and 24 microns (at a heliocentric distance, rh, of 4.8 AU) that sample the size and color-temperature of the nucleus are combined with aphelion R-band light curves observed at the Very Large Telescope (VLT) and yield a mean effective radius of 2.04 +/- 0.11 km, and an R-band geometric albedo of 0.054 +/- 0.006. The amplitudes of the R-band and mid-infrared light curves agree, which suggests that the variability is dominated by the shape of the nucleus. We also detect the dust trail of the comet at 4.8 and 5.5 AU, constrain the grain sizes to be less than or similar to 6 mm, and estimate the impact hazard to Rosetta. We find no evidence for recently ejected dust in our images. If the activity of 67P is consistent from orbit to orbit, then we may expect the Rosetta spacecraft will return images of an inactive or weakly active nucleus as it rendezvous with the comet at rh = 4 AU in 2014.Comment: 19 pages, 2 tables, 10 figures. Accepted for publication in the Astronomical Journa

Comet Dust After Deep Impact

When the Deep Impact Mission hit Jupiter Family comet 9P/Tempel 1, an ejecta crater was formed and an pocket of volatile gases and ices from 10-30 m below the surface was exposed (A Hearn et aI. 2005). This resulted in a gas geyser that persisted for a few hours (Sugita et al, 2005). The gas geyser pushed dust grains into the coma (Sugita et a1. 2005), as well as ice grains (Schulz et al. 2006). The smaller of the dust grains were submicron in radii (0-25.3 micron), and were primarily composed of highly refractory minerals including amorphous (non-graphitic) carbon, and silicate minerals including amorphous (disordered) olivine (Fe,Mg)2SiO4 and pyroxene (Fe,Mg)SiO3 and crystalline Mg-rich olivine. The smaller grains moved faster, as expected from the size-dependent velocity law produced by gas-drag on grains. The mineralogy evolved with time: progressively larger grains persisted in the near nuclear region, having been imparted with slower velocities, and the mineralogies of these larger grains appeared simpler and without crystals. The smaller 0.2-0.3 micron grains reached the coma in about 1.5 hours (1 arc sec = 740 km), were more diverse in mineralogy than the larger grains and contained crystals, and appeared to travel through the coma together. No smaller grains appeared at larger coma distances later (with slower velocities), implying that if grain fragmentation occurred, it happened within the gas acceleration zone. These results of the high spatial resolution spectroscopy (GEMINI+Michelle: Harker et 4. 2005, 2006; Subaru+COMICS: Sugita et al. 2005) revealed that the grains released from the interior were different from the nominally active areas of this comet by their: (a) crystalline content, (b) smaller size, (c) more diverse mineralogy. The temporal changes in the spectra, recorded by GEMIM+Michelle every 7 minutes, indicated that the dust mineralogy is inhomogeneous and, unexpectedly, the portion of the size distribution dominated by smaller grains has a more diverse mineralogy. The lower spatial resolution, high sensitivity Spitzer IRS data reveal resonances of refractory minerals (those seen by GEMINI+Michelle plus ortho-pyroxene)) as well resonances that can be attributed to phillosilicates (layer lattice silicates such as Montmorillonite) (Lisse et al. 2006). Pre- and post-impact, micron to submicron grains were deciphered to be present in the coma by the modeling the high spatial resolution images to account for nucleus plus inner coma fluxes (Wooden et al. 2005, 2006; Harker et al. 2005, 2006a). Note also that crystalline silicates were released from the interior of 73P-B/SW-3 as it disintegrated (Harker et al. 2006b). From the Deep Impact and the disintegration of 73P-B, we are led to ask the questians: Why is the mineralogy of the dust released from a volatile-rich pocket beneath the surface different from the dust that is released from the nominally active areas? Could the most volatile pockets be exhausted quickly? Why would crystalline silicates be associated with more volatile materials? Perhaps the structure of the comet is so inhomogeneous, e.g., the layered pile mode2 of the nucleus (Belton et al. 2006), that a reservoir of crystalline silicate and submicron grains just happens to not be released by the nominally active areas of comet 9P? Perhaps comets lose matter through their mantles from below their surfaces, thus preserving ancient topographic structures and radiation damaged silicates and carbon? We will discuss and ponder different scenarios. We will discuss future directions for coordinated observations of JF comets

Absorption Efficiencies of Forsterite. I: DDA Explorations in Grain Shape and Size

We compute the absorption efficiency (Qabs) of forsterite using the discrete dipole approximation (DDA) in order to identify and describe what characteristics of crystal grain shape and size are important to the shape, peak location, and relative strength of spectral features in the 8-40 {\mu}m wavelength range. Using the DDSCAT code, we compute Qabs for non-spherical polyhedral grain shapes with a_eff = 0.1 {\mu}m. The shape characteristics identified are: 1) elongation/reduction along one of three crystallographic axes; 2) asymmetry, such that all three crystallographic axes are of different lengths; and 3) the presence of crystalline faces that are not parallel to a specific crystallographic axis, e.g., non-rectangular prisms and (di)pyramids. Elongation/reduction dominates the locations and shapes of spectral features near 10, 11, 16, 23.5, 27, and 33.5 {\mu}m, while asymmetry and tips are secondary shape effects. Increasing grain sizes (0.1-1.0 {\mu}m) shifts the 10, 11 {\mu}m features systematically towards longer wavelengths and relative to the 11 {\mu}m feature increases the strengths and slightly broadens the longer wavelength features. Seven spectral shape classes are established for crystallographic a-, b-, and c-axes and include columnar and platelet shapes plus non-elongated or equant grain shapes. The spectral shape classes and the effects of grain size have practical application in identifying or excluding columnar, platelet or equant forsterite grain shapes in astrophysical environs. Identification of the shape characteristics of forsterite from 8-40 {\mu}m spectra provides a potential means to probe the temperatures at which forsterite formed.Comment: 55 pages, 15 figure

Mid-Infrared Spectrophotometric Observations of Fragments B and C of Comet 73P/Schwassmann-Wachmann 3

We present mid-infrared spectra and images from the GEMINI-N (+Michelle) observations of fragments SW3-[B] and SW3-[C] of the ecliptic (Jupiter Family) comet 73P/Schwassmann-Wachmann 3 pre-perihelion. We observed fragment B soon after an outburst event (between 2006 April 16 - 26 UT) and detected crystalline silicates. The mineralogy of both fragments was dominated by amorphous carbon and amorphous pyroxene. The grain size distribution (assuming a Hanner modified power-law) for fragment SW3-[B] has a peak grain radius of a_p ~ 0.5 micron, and for fragment SW3-[C], a_p ~ 0.3 micron; both values larger than the peak grain radius of the size distribution for the dust ejected from ecliptic comet 9P/Tempel 1 during the Deep Impact event (a_p = 0.2 micron. The silicate-to-carbon ratio and the silicate crystalline mass fraction for the submicron to micron-size portion of the grain size distribution on the nucleus of fragment SW3-[B] was 1.341 +0.250 -0.253 and 0.335 +0.089 -0.112, respectively, while on the nucleus of fragment SW3-[C] was 0.671 +0.076 -0.076 and 0.257 +0.039 -0.043, respectively. The similarity in mineralogy and grain properties between the two fragments implies that 73P/Schwassmann-Wachmann 3 is homogeneous in composition. The slight differences in grain size distribution and silicate-to-carbon ratio between the two fragments likely arises because SW3-[B] was actively fragmenting throughout its passage while the activity in SW3-[C] was primarily driven by jets. The lack of diverse mineralogy in the fragments SW3-[B] and SW3-[C] of 73P/Schwassmann-Wachmann 3 along with the relatively larger peak in the coma grain size distribution suggests the parent body of this comet may have formed in a region of the solar nebula with different environmental properties than the natal sites where comet C/1995 O1 (Hale-Bopp) and 9P/Tempel 1 nuclei aggregated.Comment: 31 pages, 5 figure, accepted for publication in A

Comet Dust: The Diversity of "Primitive" Particles and Implications

Comet dust is primitive and shows significant diversity. Our knowledge of the properties of primitive particles has expanded significantly through microscale investigations of cosmic dust samples ( IDP's(Interplanetary Dust Particles) and AMM's (Antarctic Micrometeorites)) and of comet dust samples (Stardust and Rosetta's COSIMA), as well as through remote sensing (spectroscopy and imaging) via Spitzer and via spacecraft encounters with 103P/Hartley 2 and 67P/Churyumov-Gerasimenko. Microscale investigations show that comet dust and cosmic dust are particles of unequilibrated materials, including aggregates of materials unequilibrated at submicron scales. We call unequilibrated materials "primitive" and we deduce they were incorporated into ice-rich (H2O-, CO2-, and CO-ice) parent bodies that remained cold, i.e., into comets, because of the lack of aqueous or thermal alteration since particle aggregation; yet some Stardust olivines suggest mild thermal metamorphism. Primitive particles exhibit a diverse range of: structure and typology; size and size distribution of constituents; concentration and form of carbonaceous and organic matter; D-, N-, and O- isotopic enhancements over solar; Mg-, Fe-contents of the silicate minerals; the compositions and concentrations of sulfides, and of less abundant mineral species such as chondrules, CAIs and carbonates. The uniformity within a group of samples points to: aerodynamic sorting of particles and/or particle constituents; the inclusion of a limited range of oxygen fugacities; the inclusion or exclusion of chondrules; a selection of organics. The properties of primitive particles imply there were disk processes that resulted in different comets having particular selections of primitive materials. The diversity of primitive particles has implications for the diversity of materials in the protoplanetary disk present at the time and in the region where the comets formed