Investigations of the pre-Deep Impact morphology of the dust coma of Comet Tempel-1

The pre-Deep Impact images of Comet Tempel-1 obtained at the Indian Astronomical Observatory are used to investigate the morphology of the dust coma of the comet. We show that the trajectory of a cometary grain under the influence of solar radiation pressure is a reliable diagnostic to estimate its initial velocity. Four main active regions at mean latitudes +45° ± 5°(D), 0° ± 5° (E),−30° ± 5°(A) and−60° ± 5°(F) are found to explain the morphology of the dust coma in the ground-based and published images obtained by the High Resolution Instrument(HRI) cameras aboard the Deep Impact flyby spacecraft. From a χ² fit of the intensity distribution in the observed and the simulated images, we derive the fraction of the productivity of the active vents to the total dust emission of the comet to be 27%. Of this the southern source alone accounts for 19.8%. The grains are found to be ejected with a velocity distribution with an upper limit of 70 ± 7 m s⁻¹. However, the broad region ‘A’ appears to eject slower grains with an upper limit of 24 ± 2.5 m s⁻¹. This source, that is active throughout the cycle is likely to be driven by CO₂ sublimation. We compute the dependence of the percentage contribution of the southern source on the heliocentric distance and show that this ratio varies over the apparition and reaches a maximum at around 260 days before perihelion. The published images of the nucleus of Comet Tempel-1 show significant departure from sphericity. Therefore, the torque exerted by the enhanced activity of the southern region may be significant enough to produce changes in the rotational state of the nucleus before each perihelion passage.     

Ion composition and chemistry in the coma of Comet 1P/Halley—A comparison between Giotto's Ion Mass Spectrometer and our ion-chemical network

In order to understand the cometary plasma environment it is important to track the closely linked chemical reactions that dominate ion evolution. We used a coupled MHD ion-chemistry model to analyze previously unpublished Giotto High Intensity Ion Mass Spectrometer (HIS-IMS) data. In this way we study the major species, but we also try to match some minor species like the CH_x and the NH_x groups. Crucial for this match is the model used for the electrons since they are important for ion-electron recombination. To further improve our results we included an enhanced density of supersonic electrons in the ion pile-up region which increases the local electron impact ionization. In this paper we discuss the results for the following important ions: C⁺, CH⁺, CH⁺₂, CH⁺₃, N⁺, NH⁺, NH⁺₂, NH⁺₃, NH⁺₄, O⁺, OH⁺, H₂O⁺, H₃O⁺, CO⁺, HCO⁺, H₃CO⁺, and CH₃OH⁺₂. We also address the inner shock which is very distinctive in our MHD model as well as in the IMS data. It is located just inside the contact surface at approximately 4550 km. Comparisons of the ion bulk flow directions and velocities from our MHD model with the data measured by the HIS-IMS give indication for a solar wind magnetic field direction different from the standard Parker angle at Halley's position. Our ion-chemical network model results are in a good agreement with the experimental data. In order to achieve the presented results we included an additional short lived inner source for the C⁺, CH⁺, and CH⁺₂ ions. Furthermore we performed our simulations with two different production rates to better match the measurements which is an indication for a change and/or an asymmetric pattern (e.g. jets) in the production rate during Giotto's fly-by at Halley's comet.     

Secular light curve of Comet 28P/Neujmin 1 and of spacecraft target Comets 1P/Halley, 9P/Tempel 1, 19P/Borrelly, 21P/Giacobinni-Zinner, 26P/Grigg-Skjellerup, 67P/Churyumov-Gerasimenko, and 81P/Wild 2

We present the secular light curves of eight comets listed in the title. Two plots per comet are needed to study these objects: a reduced magnitude (to Δ=1 AU = geocentric distance) vs time, and a reduced magnitude vs LogR (R=heliocentric distance). A total of over 16 new parameters, are measured from both plots, and give an unprecedented amount of information to characterize these objects: the onset of sublimation (RON), the offset of sublimation (ROFF), the time lag at perihelion (LAG), the absolute magnitude (m(1,1)), the maximum magnitude at perihelion (mMAX(1,LAG)), the nuclear magnitudes (VN), the amplitude of the secular light curve (ASEC), plus several others, and the photometric functions needed to describe the envelope. The most significant findings of this investigation are: (a) The envelope of the observations is the best representation of the secular light curve. (b) The H10 photometric system is unable to explain the curves and a new set of photometric rules and functions is used. (c) Only four comets exhibit power laws in their secular light curves, and only partially: 1P, 19P, 21P, and 81P. All others have to be described by more complex functions. Of the four, three exhibit a break of the power law, requiring two laws pre-perihelion and one post-perihelion. The reason for this behavior is not understood. (d) We predict the existence of a photometric anomaly in the secular light curve of 67P/Churyumov-Gerasimenko, evidenced by a region of diminished activity from −119 to −6 days before perihelion, that might be interpreted as a topographic effect or the turn off of an active region. (e) We define a photometric parameter (P-AGE) that attempts to measure the relative age of a comet through the activity exhibited in the secular light curve. 81P/Wild 2 (a comet that has recently entered the inner Solar System) is confirmed as a young object, while 28P/Neujmin 1 is confirmed as a very old comet. (f) Arranging the comets by P-AGE also classifies them by shape. A preliminary classification is achieved. (g) The old controversy of what is a nuclear magnitude is clearly resolved.     

Cometary dust properties retrieved from polarization observations: Application to C/1995 O1 Hale-Bopp and 1P/Halley

The analysis of the polarized light scattered by cometary dust particles provides information on the physical properties of the solid component of cometary comae for C/1995 O1 Hale-Bopp and 1P/Halley. A model of light scattering by a size distribution of aggregates of up to 256 submicron-sized grains (spherical or spheroidal) mixed with single spheroidal particles has been developed, with its parameters adjusted to fit the phase angle and wavelength dependence of the polarization observations. The particles are built of two materials: a non-absorbing silicates-type material and a more absorbing organic-type material. The model reproduces accurately the inversion angle and the positive branch of the polarization phase curves from the visible to the near-infrared spectral domains. A negative branch of the polarization phase curves appears in our model, although the negative branch is not deep enough to reproduce accurately the observations. Significant differences are shown between the two comets, with dominance of small grains in the coma of Comet C/1995 O1 Hale-Bopp, well fitted by a distribution of the volume-equivalent diameter, a, following a^−3.0 with a lower cutoff around 0.20 μm and an upper cutoff of at least 40 μm. For 1P/Halley, the size distribution follows a^−2.8 with a lower cutoff around 0.26 μm and an upper cutoff of about 38 μm. The relative amount of organic-type particles is larger for 1P/Halley while the amount of aggregates, significant for both comets, is larger for C/1995 O1 Hale-Bopp.     

Radial distribution of production rates, loss rates and densities corresponding to ion masses ?40 amu in the inner coma of Comet Halley: Composition and chemistry

In this paper we have studied the chemistry of C, H, N, O, and S compounds corresponding to ions of masses ?40 amu in the inner coma of the Comet 1P/Halley. The production rates, loss rates, and ion mass densities are calculated using the Analytical Yield Spectrum approach and solving coupled continuity equation controlled by the steady state photochemical equilibrium condition. The primary ionization sources in the model are solar EUV photons, photoelectrons, and auroral electrons of the solar wind origin. The chemical model couples ion-neutral, electron-neutral, photon-neutral and electron-ion reactions among ions, neutrals, electrons, and photons through over 600 chemical reactions. Of the 46 ions considered in the model the chemistry of 24 important ions (viz., CH₃OH⁺₂, H₃CO⁺, NH⁺₄, H₃S⁺, H₂CN⁺, H₂O⁺, NH⁺₃, CO⁺, C₃H⁺₃, OH⁺, H₃O⁺, CH₃OH⁺, C₃H⁺₄, C₂H⁺₂, C₂H⁺, HCO⁺, S⁺, CH⁺₃, H₂S⁺, O⁺, C⁺, CH⁺₄, C⁺₂, and O⁺₂) are discussed in this paper. At radial distances <1000 km, the electron density is mainly controlled by 6 ions, viz., NH⁺₄, H₃O⁺, CH₃OH⁺₂, H₃S⁺, H₂CN⁺, and H₂O⁺, in the decreasing order of their relative contribution. However, at distances >1000 km, the 6 major ions are H₃O⁺, CH₃OH⁺₂, H₂O⁺, H₃CO⁺, C₂H⁺₂, and NH⁺₄; along with ions CO⁺, OH⁺, and HCO⁺, whose importance increases with further increase in the radial distance. It is found that at radial distances greater than ∼1000 km (±500 km) the major chemical processes that govern the production and loss of several of the important ions in the inner coma are different from those that dominate at distances below this value. The importance of photoelectron impact ionization, and the relative contributions of solar EUV, and auroral and photoelectron ionization sources in the inner coma are clearly revealed by the present study. The calculated ion mass densities are compared with the Giotto Ion Mass Spectrometer (IMS) and Neutral Mass Spectrometer (NMS) data at radial distances 1500, 3500, and 6000 km. There is a reasonable agreement between the model calculation and the Giotto measurements. The nine major peaks in the IMS spectra between masses 10 and 40 amu are reproduced fairly well by the model within a factor of two inside the ionopause. We have presented simple formulae for calculating densities of the nine major ions, which contribute to the nine major peaks in the IMS spectra, throughout the inner coma that will be useful in estimating their densities without running the complex chemical models.     

Organics-glued dust aggregates in Halley's coma and some implications of the snow line inward motion for the formation of comets

In view of the solar nebula models, organics-glued dust aggregates (whose disintegration resulted in the two phenomena found in Halley's coma, the dust boundary and small-scale dust structures) originated due to coagulation of iceless dust particles somewhere within the snow line, and then were incorporated into Halley's nucleus as a consequence of the snow line inward motion. This implies that two types of comets exist: outer comets, formed entirely beyond the snow line, and inner comets, similar to Halley, which are bodies intermediate between outer comets and primitive asteroids. The presence of large iceless dust aggregates in nuclei of inner comets constrains the inward drift velocity of meter-sized dust bodies, which in turn implies that the radial transport of water in the solar nebula was predominantly outward. It is shown that in nuclei of inner comets: both the upper mass limit of iceless dust aggregates and the ice mantle thickness increase with decreasing formation heliocentric distance, while the cumulative mass distribution index decreases; the lower limit of the mass index is ∼0.8, and the upper limit of the ice mantle thickness is ∼10⁻³ cm (∼200 times the interstellar value); the lower limit of the latent heat of organics in organic mantles of submicron particles increases toward small heliocentric distances; the recondensation of organics combined with the growth of dust bodies leads to a fractionation of organics within iceless dust aggregates; last accreted sub-units of an aggregate are always glued by organics with the lowest value of the latent heat, which somewhat exceeds 60 kJ/mol. Based on in situ observations at Halley, the parameters characterizing iceless dust aggregates in that comet are calculated. Finally, feasible observational tests of the conclusions drawn are discussed.     

Non-LTE radiative transfer for sub-millimeter water lines in Comet 67P/Churyumov-Gerasimenko

The European Space Agency (ESA) Rosetta spacecraft (Schulz, R., Alexander, C., Boehnhardt, H., Glassmeier, K.H. (Eds.) [2009]. “ROSETTA - ESA”) will encounter Comet 67P/Churyumov-Gerasimenko in 2014 and spend the next 18 months in the vicinity of the comet, permitting very high spatial and spectral resolution observations of the coma and nucleus. During this time, the heliocentric distance of the comet will change from ∼3.5 AU to ∼1.3 AU, accompanied by an increasing temperature of the nucleus and the development of the coma. The Microwave Instrument for the Rosetta Orbiter (MIRO) will observe the ground-state rotational transition (1₁₀-1₀₁) of H₂¹⁶O at 556.936 GHz, the two isotopologues H₂¹⁷O and H₂¹⁸O and other molecular transitions in the coma during this time (Gulkis, S. et al., [2007]. MIRO: Microwave Instrument for Rosetta Orbiter. Space Sci. Rev. 128, 561-597).The aim of this study is to simulate the water line spectra that could be obtained with the MIRO instrument and to understand how the observed line spectra with various viewing geometries can be used to study the physical conditions of the coma and the water excitation processes throughout the coma. We applied an accelerated Monte Carlo method to compute the excitations of the seven lowest rotational levels (1₀₁, 1₁₀, 2₁₂, 2₂₁, 3₀₃, 3₁₂, and 3₂₁) of ortho-water using a comet model with spherically symmetric water outgassing, density, temperature and expansion velocity at three different heliocentric distances 1.3 AU, 2.5 AU, and 3.5 AU. Mechanisms for the water excitation include water-water collisions, water-electron collisions, and infrared pumping by solar radiation.Synthetic line spectra are calculated at various observational locations and directions using the MIRO instrument parameters. We show that observations at varying viewing distances from the nucleus and directions have the potential to give diagnostic information on the continuum temperature and water outgassing rates at the surface of the nucleus, and the gas density, expansion velocity, and temperature of the coma as a function of distance from the nucleus. The gas expansion velocity and temperature affect the spectral line width and frequency shift of the line from the rest frequency, while the gas density (which is directly related to the outgassing rate) and the line excitation temperature determine the antenna temperature of the absorption and emission signal in the line profile.     

Thermal model of water and CO activity of Comet C/1995 O1 (Hale-Bopp)

An investigation of the activity of Comet C/1995 O1 (Hale-Bopp) with a thermophysical nucleus model that does not rely on the existence of amorphous ice is presented. Our approach incorporates recent observations allowing to constrain important parameters that control cometary activity. The model accounts for heat conduction, heat advection, gas diffusion, sublimation, and condensation in a porous ice-dust matrix with moving boundaries. Erosion due to surface sublimation of water ice leads to a moving boundary. The movement of the boundary is modeled by applying a temperature remapping technique which allows us to account for the loss in the internal energy of the eroded surface material. These kind of problems are commonly referred to as Stefan problems. The model takes into account the diurnal rotation of the nucleus and seasonal effects due to the strong obliquity of Hale-Bopp as reported by Jorda et al. (Jorda, L., Rembor, K., Lecacheux, J., Colom, P., Colas, F., Frappa, E., Lara, L.M. [1997]. Earth Moon Planets 77, 167-180). Only bulk sublimation of water and CO ice are considered without further assumptions such as amorphous ices with certain amount of occluded CO gas. Confined and localized activity patterns are investigated following the reports of Lederer and Campins (Lederer, S.M., Campins, H. [2002]. Earth Moon Planets 90, 381-389) about the chemical heterogeneity of Hale-Bopp and of Bockelée-Morvan et al. (Bockelée-Morvan, D., Henry, F., Biver, N., Boissier, J., Colom, P., Crovisier, J., Despois, D., Moreno, R., Wink, J. [2009]. Astron. Astrophys. 505, 825-843) about a strong CO source at a latitude of 20°. The best fit to the observations of Biver et al. (Biver, N. et al. [2002]. Earth Moon Planets 90, 5-14) is obtained with a low thermal conductivity of 0.01 W m⁻¹ K⁻¹. This is in agreement with recent results of the Deep Impact mission to 9P/Tempel 1 (Groussin, O., A’Hearn, M.F., Li, J.-Y., Thomas, P.C., Sunshine, J.M., Lisse, C.M., Meech, K.J., Farnham, T.L., Feaga, L.M., Delamere, W.A. [2007]. Icarus 187, 16-25) and with previous thermal simulations (Kührt, E. [1999]. Space Sci. Rev. 90, 75-82). The water production curve matches the production rates well from −4 AU pre-perihelion to the outgoing leg while the model does not reproduce so well the water production beyond 4 AU pre-perihelion. The CO production curve is a good fit to the measurements of Biver et al. (2002) over the whole measured heliocentric range from −7 AU pre- to 15 AU post-perihelion.     

Comet classification with new methods for gas and dust spectroscopy

We present the results of a program of comet long-slit spectroscopy with the Kast Dual Spectrograph on the 3-m Shane Telescope at Lick Observatory. A total of 26 comets, from a variety of dynamical families, were observed on 39 different nights from 1996 to 2007. A new statistical method extracted the twilight sky from comet frames, because traditional sky subtraction techniques were inadequate. Because previously published Haser model parent and daughter scale lengths did not fit the data well, unbiased ranges of scale lengths were searched for the best-fitting pairs. Coma gas production rates for OH, CN, C₂, C₃, NH, NH₂, and OH confirmed the widely reported carbon-chain depletion for a sub-class of comets, most notably high-perihelion Jupiter-family comets observed at r_h > 1.5 AU, with different behaviors for C₂ and C₃. Our long-slit spectroscopy data was also adapted for the A(θ)fρ dust production parameter. The assumption that A(θ)fρ is constant throughout the nucleus was not upheld. High dust-to-gas ratios for comets with large perihelia were not a selection effect, and suggest that the dust was released earlier in the formation of the coma than the gas. The dust-to-gas ratio did not exhibit any evolutionary traces between different comet dynamical families. The comet survey illuminates the diversity among comets, including the unusually carbon poor Comet 96P/Machholz.     

Deep Impact's target Comet 9P/Tempel 1 at multiple apparitions: Seasonal and secular variations in gas and dust production

We present results from multi-apparition narrowband photometry of Deep Impact target Comet 9P/Tempel 1. In support of the mission, we obtained data during monthly observing runs between March and September 2005, and these are combined with and compared to observations obtained during the 1983 and 1994 apparitions. A strong seasonal effect is seen, with peak production rates occurring 4-8 weeks before perihelion, with some variation evident among the different species. There is also evidence of a slight systematic shift towards a later time of peak production in 2005 as compared to 1983. Early in the apparition, the radial profile of the dust was much steeper than the canonical 1/ρ, but the slope became progressively smaller until very little departure from 1/ρ remained by late June, a change possibly associated with the general seasonal effects. Unexpectedly, an unprecedented large overall decrease in production rates has taken place since 1983, with water at only about 42% of the 1983 values, CN at about 53%, and dust, based on the proxy A(θ)fρ, at about 77%. Other gas species exhibited declines intermediate between that of CN and of the dust. The large differences in the amount of secular decline among all of the species implies compositional inhomogeneities among source regions on the surface of the nucleus, with one region progressively becoming less active over only a few orbits. While the simplest explanation would invoke either devolatilization or covering up of the ice, no other comet has shown such a rapid change in outgassing unless accompanied by a significant change in its orbit. We, therefore, hypothesize that a change in available solar radiation due to precession of the pole might instead be causing the progressive drop in cometary activity. Given the small obliquity of the rotation axis derived from the Deep Impact observations, and a presumed small rate of precession, the source region would need to be located near the pole to explain both the large secular and seasonal trends.     

Deep Impact, Stardust-NExT and the behavior of Comet 9P/Tempel 1 from 1997 to 2010

We present observational data for Comet 9P/Tempel 1 taken from 1997 through 2010 in an international collaboration in support of the Deep Impact and Stardust-NExT missions. The data were obtained to characterize the nucleus prior to the Deep Impact 2005 encounter, and to enable us to understand the rotation state in order to make a time of arrival adjustment in February 2010 that would allow us to image at least 25% of the nucleus seen by the Deep Impact spacecraft to better than 80 m/pixel, and to image the crater made during the encounter, if possible. In total, ∼500 whole or partial nights were allocated to this project at 14 observatories worldwide, utilizing 25 telescopes. Seventy percent of these nights yielded useful data. The data were used to determine the linear phase coefficient for the comet in the R-band to be 0.045 ± 0.001 mag deg⁻¹ from 1° to 16°. Cometary activity was observed to begin inbound near r ∼ 4.0 AU and the activity ended near r ∼ 4.6 AU as seen from the heliocentric secular light curves, water-sublimation models and from dust dynamical modeling. The light curve exhibits a significant pre- and post-perihelion brightness and activity asymmetry. There was a secular decrease in activity between the 2000 and 2005 perihelion passages of ∼20%. The post-perihelion light curve cannot be easily explained by a simple decrease in solar insolation or observing geometry. CN emission was detected in the comet at 2.43 AU pre-perihelion, and by r = 2.24 AU emission from C₂ and C₃ were evident. In December 2004 the production rate of CN increased from 1.8 × 10²³ mol s⁻¹ to Q_CN = 2.75 × 10²³ mol s⁻¹ in early January 2005 and 9.3 × 10²⁴ mol s⁻¹ on June 6, 2005 at r = 1.53 AU.     

The dust trail of Comet 67P/Churyumov-Gerasimenko between 2004 and 2006

We report on observations of the dust trail of Comet 67P/Churyumov-Gerasimenko (CG) in visible light with the Wide Field Imager at the ESO/MPG 2.2 m telescope at 4.7 AU before aphelion, and at with the MIPS instrument on board the Spitzer Space Telescope at 5.7 AU both before and after aphelion. The comet did not appear to be active during our observations. Our images probe large dust grains emitted from the comet that have a radiation pressure parameter β<0.01. We compare our observations with simulated images generated with a dynamical model of the cometary dust environment and constrain the emission speeds, size distribution, production rate and geometric albedo of the dust. We achieve the best fit to our data with a differential size distribution exponent of −4.1, and emission speeds for a β=0.01 particle of 25 m/s at perihelion and 2 m/s at 3 AU. The dust production rate in our model is on the order of 1000 kg/s at perihelion and 1 kg/s at 3 AU, and we require a dust geometric albedo between 0.022 and 0.044. The production rates of large (>) particles required to reproduce the brightness of the trail are sufficient to also account for the coma brightness observed while the comet was inside 3 AU, and we infer that the cross-section in the coma of CG may be dominated by grains of the order of .     

Nucleus properties of Comet 9P/Tempel 1 estimated from non-gravitational force modeling

The nucleus mass and bulk density of Comet 9P/Tempel 1 have been estimated by utilizing the non-gravitational force modeling technique. Here, the water production rates and non-gravitational perturbations of the orbit are calculated for a large number of model nuclei with different surface ice distribution patterns. By requiring that the empirical water production rate curve is reproduced, a subset of model nuclei are selected, for which masses are calculated by demanding that empirical non-gravitational changes of the orbital period and in the longitude of perihelion (per revolution) are reproduced. We obtain a mass M=5.8(±1.6)×¹⁰¹³ kg, and a bulk density , which compares very well with measurements made by the Deep Impact Science Team. The main goal of the current work is therefore to demonstrate functionality of an indirect method, i.e., mass estimation through non-gravitational force modeling, by comparing such results to ground truth data. Furthermore, the thermal inertia of active areas is estimated as 30-100 MKS, using a comparatively realistic thermophysical model (although a value in the range 100-350 MKS is obtained with a more simple model). An active area fraction of ∼3% is predicted, and these areas are probably confined to the northern hemisphere, being located close to the cometary equator.     

Dust coma morphology in the Deep Impact images of Comet 9P/Tempel 1

We present an overview of the dust coma observations of Comet Tempel 1 that were obtained during the approach and encounter phases of the Deep Impact mission. We use these observations to set constraints on the pre-impact activity of the comet and discuss some preliminary results. The temporal and spatial changes that were observed during approach reveal three distinct jets rotating with a 1.7-day periodicity. The brightest jet produces an arcuate feature that expands outward with a projected velocity of about 12 m s⁻¹, suggesting that the ambient dust coma is dominated by millimeter-sized dust grains. As the spatial resolution improves, more jets and fans are revealed. We use stereo pairs of high-resolution images to put some crude constraints on the source locations of some of the brightest features. We also present a number of interesting coma features that were observed, including surface jets detected at the limb of the nucleus when the exposed ice patches are passing over the horizon, and features that appear to be jets emanating from unilluminated sources near the negative pole. We also provide a list of 10 outbursts of various sizes that were observed in the near-continuous monitoring during the approach phase.     

Imaging polarimetry of the dust coma of Comet Tempel 1 before and after Deep Impact at Haute-Provence Observatory

We have observed the coma of Comet 9P/Tempel 1, the target of the Deep Impact mission, by the polarization imaging technique, before and after the impact event (−32, −7, +43 and +65 h). Our observations were conducted in the red wavelength domain from Haute-Provence Observatory (France), with the 80-cm telescope. The overall polarization of 9P/Tempel 1, as obtained near 41° phase angle, is monitored and compared to data from other (active and less active) comets studied by the same technique. The linear polarization of the dust ejected by the impact is compared to previous observations of dust present in jets, ejected during outbursts or released when comets happen to split. At phase angles of about 41°, the difference in polarization between the comets with a low maximum in polarization and the comets with a high maximum in polarization is about 1%; it may thus be difficult to conclude about the classification. Nevertheless, the overall polarization after the impact rapidly reached a value corresponding to the high polarization class of comets, and later progressively decreased to its initial value. The polarization was measured to be slightly lower (about 1%) before the impact than after it in a 26,000-km aperture. The plume formed from dust ejected by the impact was still present 65 h after it. The variations of the intensity and the polarization in the coma provide some clues to variations of the physical properties of the particles; comparison with other techniques corroborates the presence of large particles and of submicron-sized grains in aggregates.     

Surface temperature of the nucleus of Comet 9P/Tempel 1

The Deep Impact (DI) spacecraft encountered Comet 9P/Tempel 1 on July 4th, 2005 and observed it with several instruments. In particular, we obtained infrared spectra of the nucleus with the HRI-IR spectrometer in the wavelength range of 1.0-4.9 μm. The data were taken before impact, with a maximum resolution of ∼120 m per pixel at the time of observation. From these spectra, we derived the first directly observed temperature map of a comet nucleus. The surface temperature varied from 272±7 to 336±7 K on the sunlit hemisphere, matching the surface topography and incidence angle. The derived thermal inertia is low, most probably <50 W K⁻¹ m⁻² s^1/2. Combined with other arguments, it is consistent with the idea that most of rapidly varying thermal physical processes, in particular the sublimation of volatiles around perihelion, should occur close to the surface. Thermal inertia is sufficient to explain the temperature map of the nucleus of Comet Tempel 1 to first order, but other physical processes like roughness and self-radiation are required to explain the details of the temperature map. Finally, we evaluated that the Standard Thermal Model is a good approximation to derive the effective radius of a cometary nucleus with an uncertainty lower than ∼10% if combined with a thermal infrared light curve.     

HCN and CN in Comet 2P/Encke: Models of the non-isotropic, rotation-modulated coma and CN parent life time

Axisymmetric models of the outgassing of a cometary nucleus have been constructed. Such models can be used to describe a nucleus with a single active region. The models may include a solar zenith angle dependence of the outgassing. They retrieve the outgassing flux at distances from the nucleus where collisions between molecules are unimportant, as function of the angle with respect to the outgassing axis. The observed emissions must be optically thin. Furthermore the models assume that the outflow speed at large distance from the nucleus does not depend on direction. The value of the outflow speed is retrieved. The models are applied to CN images and HCN spectra of Comet 2P/Encke, obtained nearly simultaneously in November 2003 with the 2 m optical telescope on Mount Rozhen, Bulgaria, and with the 10 m Heinrich Hertz Submillimeter Telescope on Mount Graham, Arizona, USA. According to Sekanina (1988), Astron. J. 95, 911-924, at that time a single outgassing source was active. Input parameters to the models like the rotation period of the nucleus and a small correction to Sekanina’s rotation axis are determined from a simpler jet position angle model. The rotation is prograde with a sideric period of 11.056 ± 0.024 h, in agreement with literature values. The best fit model has an outflow speed of 0.95 ± 0.04 km s⁻¹. The same value has been derived from the corkscrew appearing in the CN images. The location of the outgassing axis is at colatitude δ_a = 7.4° ± 2.9° and longitude λ_a = 235° ± 17° (a definition of zero longitude is provided). Comet Encke’s outgassing corresponds approximately to the longitudinally averaged solar input on a spherical nucleus (i.e. very likely comes from deeper layers) but with some deficiency of outgassing at mid-latitudes and non-zero outgassing from the dark polar cap. The presence of gas flow from the dark polar cap is explained as evidence of gas flow across the terminator. The models rely mostly on the CN images. The HCN spectra are more noisy. They provide information how to determine the best fit outflow velocity and the sense of rotation. The model HCN spectra are distinctly non-Gaussian. Within error limits they are consistent with the observations. Models based solely on the HCN spectra are also presented but, because of the lower quality of the data and the unfavorable observing geometry, yield inferior results. As a by-product we determine the CN parent life time from our CN observations. The solar EUV and Lyα radiation field at the time of our observations is taken into account.     

The gas production of Comet 9P/Tempel 1 around the Deep Impact date

The target of the Deep Impact space mission (NASA), Comet 9P/Tempel 1, was observed from two nights before impact to eight nights after impact using the FORS spectrographs at the ESO VLT UT1 and UT2 telescopes. Low resolution optical long-slit spectra were obtained to study the evolution of the gas coma around the Deep Impact event. Following first results of this observing campaign on the CN and dust activity [Rauer, H., Weiler, M., Sterken, C., Jehin, E., Knollenberg, J., Hainaut, O., 2006. Astron. Astrophys. 459, 257-263], this work presents a study of the complete dataset on CN, C₂, C₃, and NH₂ activity of Comet 9P/Tempel 1. An extended impact gas cloud was observed moving radially outwards. No compositional differences between this impact cloud and the undisturbed coma were found as far as the observed radicals are concerned. The gas production rates before and well after impact indicate no change in the cometary activity on an intermediate time scale. Over the observing period, the activity of Comet 9P/Tempel 1 was found to be related to the rotation of the cometary nucleus. The rotational lightcurve for different gaseous species provides indications for compositional differences among different parts of the nucleus surface.     

SOHO/SWAN observations of comets with small perihelia: C/2002 V1 (NEAT), C/2002 X5 (Kudo–Fujikawa), 2006 P1 (McNaught) and 96P/Machholz 1

SWAN, the all-sky hydrogen Lyman-alpha camera on the SOHO spacecraft, designed primarily to image the interplanetary neutral hydrogen around the Sun, also observes comets continuously over large portions of their apparitions to the north and south of the ecliptic and at small solar elongation angles. Because of SOHO’s location at the L1 Lagrange point, analysis of SWAN images provides excellent temporal coverage of water production. We report here our results of observations of some interesting target comets selected from the extensive SWAN archive. These include three Oort Cloud Comets C/2002 V1 (NEAT), C/2002 X5 (Kudo–Fujikawa), C/2006 P1 (McNaught) and three apparitions of atypical short-period Comet 96P/Machholz 1. The common aspect of these four comets is their small perihelion distances, which are 0.19, 0.09, 0.17, and 0.12 AU, respectively. Their water production rates over their whole apparitions can be approximated by power laws in heliocentric distance (r in AU) as follows: 1.3 × 10²⁹ r^−2.1 s⁻¹ for C/2002 V1 (NEAT), 7.5 × 10²⁸ r^−2.0 s⁻¹ for C/2002 X5 (Kudo–Fujikawa), 5.4 × 10²⁹ r^−2.4 s⁻¹ for C/2006 (P1 McNaught) and 4.6 × 10²⁷ r^−2.1 s⁻¹ for 96P/Machholz 1. We also present daily-average water production rates for the long-period comets over long continuous time periods. We examine these results in light of our growing survey of comets that is yielding some interesting comparisons of water production rate variations with heliocentric distance and taxonomic classes.