Dust tail of the active distant Comet C/2003 WT42 (LINEAR) studied with photometric and spectroscopic observations

We present the study of dust environment of dynamically new Comet C/2003 WT42 (LINEAR) based on spectroscopic and photometric observations. The comet was observed before and after the perihelion passage at heliocentric distances from 5.2 to 9.5 AU. Although the comet moved beyond the zone where water ice sublimation could be significant, its bright coma and extended dust tail evidenced the high level of physical activity. Afρ values exceeded 3000 cm likely reaching its maximum before the perihelion passage. At the same time, the spectrum of the comet did not reveal molecular emission features above the reflected continuum. Reddening of the continuum derived from the cometary spectrum is nonlinear along the dispersion with the steeper slop in the blue region. The pair of the blue and red continuum images was analyzed to estimate a color of the comet. The mean normalized reflectivity gradient derived from the innermost part of the cometary coma equals to 8% per 1000 Å that is typical for Oort cloud objects. However, the color map shows that the reddening of the cometary dust varies over the coma increasing to 15% per 1000 Å along the tail axis. The photometric images were fitted with a Monte Carlo model to construct the theoretical brightness distribution of the cometary coma and tail and to investigate the development of the cometary activity along the orbit. As the dust particles of distant comets are expected to be icy, we propose here the model, which describes the tail formation taking into account sublimation of grains along their orbits. The chemical composition and structure of these particles are assumed to correspond with Greenberg’s interstellar dust model of comet dust. All images were fitted with the close values of the model parameters. According to the results of the modeling, the physical activity of the comet is mainly determined by two active areas with outflows into the wide cones. The obliquity of the rotation axis of the nucleus equals to 20° relative to the comet’s orbital plane. The grains occupying the coma and tail are rather large amounting to 1 mm in size, with the exponential size distribution of a^−4.5. The outflow velocities of the dust particles vary from a few centimeters to tens of meters per second depending on their sizes. Our observations and the model findings evidence that the activity of the nucleus decreased sharply to a low-level phase at the end of April–beginning of May 2007. About 190 days later, in the first half of November 2007 the nucleus stopped any activity, however, the remnant tail did not disappear for more than 1.5 years at least.     

Charge-coupled device photometry of Comet P/Halley

Comet P/Halley has been observed during its approach to perihelion at heliocentric distancesR = 11.0 AU and R = 8.2 AU. No extended coma is seen and limits can be placed on the fraction of the total light contributed by coma. The brightness of the comet varies on a short time scale. The variations may be due to transient activity or to rotation of the irregular nucleus.     

A thermal modelling of cometary activity with a crystallized water ice nucleus

Arguments are presented to suggest that surface layers of the nuclei of periodic comets consist of crystallized rather than amorphous water ice and thermal modelling of such nuclei is presented. The rate of sublimation of water from a rotating nucleus is found to be greater than that from a uniformly heated nucleus. When the model is applied to P/Halley, the sublimation rate at perihelion is found to be 8.1 × 10²⁹ mol s⁻¹ for a nucleus rotating with a period of 50 hours and 7.6 × 10²⁹ for a uniformly heated nucleus on the premise that the effective radius of the nucleus is 2.5 km. The total sublimation of water per revolution is 5.38 × 10³⁶ molecules for P/Halley and 3.91 × 10³⁶ molecules per P/Crommelin. The result so obtained is discussed in relation to the observational data.     

The dust coma of Comet P/Giacobini-Zinner in the infrared

We present 1-20 micrometers photometry of P/Giacobini-Zinner obtained at the NASA Infrared Telescope Facility, during 1985 June-September (r = 1.57-1.03 AU). A broad, weak 10 micrometers silicate emission feature was detected on August 26.6; a similar weak emission feature could have been hidden in the broadband photometry on other dates. The total scattering and emitting cross section of dust in the inner coma was similar to that in other short-period comets, but a factor of 10 (r = 1.56 AU) to 100 (r = 1.03 AU) lower than the amount of dust in Comet Halley. The thermal emission continuum can be fit with models weighted toward either small or large absorbing grains. The dust production rate near perihelion was approximately 10(5) g/s (small-grain model) to approximately 10(6) g/s (large-grain model). The corresponding dust/gas mass ratio on August 26 was approximately 0.1-1. A silicate-rich heterogeneous grain model with an excess of large particles is compatible with the observed spectrum of Giacobini-Zinner on August 26. Thus, weak or absent silicate emission does not necessarily imply an absence of silicates in the dust, although the abundance of silicate particles < or = 1 micrometer radius must have been lower than in Comet Halley.     

On the visibility of nuclei of dusty comets

The problem of visibility of a cometary nucleus discussed in general terms for single scattering by dust grains. The ratio of radiatio scattered in the dust column above the surface and that reflected from the nucleus determines the visibility of features on the nuclear surface. A contrast parameter characterizing the ration of radiation foming from the nuclear surface and that of the nuclear vicinity describes the visibility of the full nucleus against the dust fore- and background. These quantities and the intensity distribution of scattered solar radiation across the nucleus and its vicinity are calculated for the case of comet Halley at a heliocentric distance of 0.9 AU after perihelion (Giotto encounter). The scattering calculations are based on an isotropic dust distribution derived from hydrodynamics gas-dust interactions resulting in a steep densiity increase right above the surface. For Newburn's nominal model of comet Halley, an optical depth of about 0.5 impairs the visibility of the nucleus somewhat.     

Narrow-Band Photometry of Comet C/1995 O1 (Hale-Bopp)

Photometric observations of comet C/1995 O1 (Hale-Bopp) carried out at the Stará Lesná Observatory since February to April 1997 are analyzed and discussed. Emission band fluxes and continuum fluxes are presented, from which the total numbers of molecules in the columns of the coma encircled by diaphragms are calculated. The production rates are estimated from the conventional Haser model. We found that the photometric exponent of dust contribution two months prior perihelion was n = 5.2. The photometric exponent n of the cometary magnitude solely to the C₂ emission alone equals 3.3 and that of CN equals 2.5. These values can be explained by a fact that the maximums of production rates of the gases were reached between March 2and 12 and not at the perihelion as it is valid for dust. These results are compared with the values of 1P/Halley (1986 III) under the similar conditions, obtained with the same method and instrument. C/Hale-Bopp exhibited 4.1 times more molecules radiating the CN-emission than 1P/Halley in the same column of the coma. The continuum flux of C/Hale-Bopp was also very strong. The ratios (to 1P/Halley) are 94:1 (Cont. 484.5) and 74:1 (Cont. 365.0). The cometary colour was the same as that of the Sun. This revised version was published online in July 2006 with corrections to the Cover Date.     

The loss and depth of CO2 ice in comet nuclei

An analytical model has been developed to simulate the chemical differentiation of a homogeneous, initially unmantled cometary nucleus composed of water ice, putative unclathrated CO₂ ice, and silicate dust in specified proportions. Selective sublimation of any free CO₂ ice present in a new comet should produce a surface layer of water ice and dust overlying the undifferentiated core. This surface layer modifies the temperature of buried CO₂ ice and restricts the outflow of gaseous CO₂. On each orbit, water sublimation closer to perihelion temporarily reduces the thickness of the water ice and dust layer and liberates dust. Most of the dust is blown off the nucleus, but a small amount of residual dust remains on the surface (cf. H. L. F. Houpis, W. H. Ip, and D. A. Mendis, 1986, Astrophys. J., in press). Our model includes the effects of nucleus rotation, arbitrary orientation of the rotation axis, latitude, heat conduction into the interior of the nucleus, restriction of CO₂ gas outflow by the water ice and dust layer, and the use of thermal conductivities for both amorphous and crystalline water ice as appropriate, featuresthat were not included in the Houpis et al. model. The model also accounts for the erosion of the water ice surface, which Houpis et al. appear to have accounted for and which is an important effect. Specifically, we investigate the effects of varying the permeability of the surface water ice layer, the mass fraction of CO₂, the orbit and the latitude, using the orbital parameters of Comets Halley and Tempel 2. It is found that CO₂ gas production should exceed H₂O gas production beyond ∼3 AU, and at 1 AU CO₂ gas production should be between 20 to 25% of H₂O gas production. The depth of CO₂ ice and the variation in the depth of CO₂ ice throughout an orbit are affected significantly by the perihelion of the orbit. The effects due to water ice permeability are significant but much less than expected on the basis of flow area. Latitude and CO₂ concentration produce relatively small effects. Under all conditions considered here, CO₂ ice should always be found within ∼1 m from the surface of comet nuclei if it is present as a free species to begin with. This result is probably generally valid for unmantled portions of most comets and qualitatively simulates the behavior of an abundant, highly volatile component in an H₂O/silicate matrix. Comparison of these and similar results with observations could yield information regarding the permeability and chemical composition of cometary material and suggest sampling strategies to minimize fractionation effects. The method is applicable to other nonwater ices.     

C/1995 O1 (Hale-Bopp): Gas Production Curves And Their Interpretation

The discovery of C/1995 O1 (Hale-Bopp) at 7 AU from the Sun provided the first opportunity to follow the activity of a bright comet over a large range of heliocentric distances r_h. Production rates of a number of parent molecules and daughter species have been monitored both pre- and postperihelion. CO was found to be the major driver of the activity far from the Sun, surpassed by water within 3 AU whose production rate reached 10³¹ s⁻¹ at perihelion. Gas production curves obtained for various species show several behaviours with r_h. Gas production curves contain important information concerning the physical state of cometary ices, the structure of the nucleus and all the processes taking place inside the nucleus leading to outgassing. They are relevant to the study of several other phenomena such as the sublimation from icy grains, dust mantling or seasonal effects. For some species, such as H₂CO or HNC, they permit to constrain their origin in the coma. We discuss models of subsurface gas production in distant comets and predictions of how such a source may vary as the comet moves along its orbit, approaching perihelion and receding again. Features in the observed gas production curves of comet Hale-Bopp are generally interpretable in terms of either subsurface production (typical example: CO at large r_h) or free sublimation (typical example: H₂O). Possible implications for the vertical stratification of the cometary ices are reviewed, and preference is found for a model with crystallization of amorphous ice close to the nuclear surface. This revised version was published online in July 2006 with corrections to the Cover Date.     

Have stellar and supernovae radiations altered the pristine nature of comets?

Effect of stellar and supernova radiations on cometary nuclei in the Oort cloud is investigated. Radiation dose received by a comet is calculated and compared with the one which Halley's comet receives by one perihelion passage. Stellar radiation provides 10 to 50% of Halley unit over 4 billion years. Inclusion of sublimation of volatile molecules such as CO or N₂ does not allow the temperature to rise to 30 K by irradiation of bright OB stars, contrary to the claim of Stern and Shull. A chance encounter with a SN provides radiation dose which is just sufficient to raise to 30 K the surface layer which is 1 m thick on the assumption that the radiation is wholly communicated to the interior. Thus, the comets remain pristine under the effect of stellar and SN radiations.     

Heat content and evolution of cometary nuclei

In continuation of an earlier study of the influence of phase transitions on the thermal behavior of cometary nuclei, the heat flux into nuclei at various distances from the Sun before and after perihelion has been investigated for the isothermal case and for the fixed subsolar point. It turns out that this heat flux may be a large fraction of the incident solar heat input, so that the surface temperature and the associated rate of evaporation are lower than usually calculated. The effect is strongly dependent on the porosity of the nucleus. The surface temperature of the nucleus reaches a maximum after perihelion, as does the size of the coma, in agreement with several observations. The denser surface layers made either of ice or of dust may break away. An ideal, initially homogeneous and spherical nucleus cannot remain isothermal so that it must gradually develop considerable surface nonuniformities through localized phase changes, evaporation, and break-away. An explanation of the splitting of comets as far as 9 AU from the Sun is suggested in terms of heating of a CO₂-rich inclusion in a nucleus.     

An idealized short-period comet model: Surface insolation,H2O flux,dust flux,and mantle evolution

A model of cometary activity is developed which integrates the feedback processes involving heat, gas, and dust transport, and dust mantle development. The model includes the effects of latitude, rotation, and spin axis orientation. Results are obtained for various grain size distributions, dust-to-ice ratios, and spin axis orientations. Attention is focused on the development, change of structure and distribution of dust mantles and their mutual interaction with ice surface temperature and gas and dust production. In this model the dust mantle controls the mechanism of gas transport not onlu by its effect on the temperature but, more importantly, by its own dynamic stability. Results suggest that an initially homogeneous short-period comet with a “cosmic” dust-to-water ice ratio, typical orbit, rotation rate, and grain size distribution would develop at most only a thin (<1 mm) cyclic mantle at all points on the nucleus. Such a fully developed temporary mantle would exist throughout the diurnal cycle only beyond ～4AU. Thus, cyclic behavior would be expected for such an idealized comet, at least for most of its lifetime. Long-term irreversible mantle development on comets with typical rotation rates was not found except regionally on Encke and also on objects with perihelia ?1.5 AU. Even in these cases, free silicate exists, after a few cycles, only as relatively rare large grains and agglomerates with radii ～1 cm scattered over exposed ice. Full mantle development would require hundreds to thousands of cycles. In the case of an initially homogeneous comet Encke, this slow incipient mantle development is shown to be the direct result of its peculiar axial orientation. High obliquity appears required for long-term mantle development for typical rotation rates and perihelia ?1.5 AU. Heat conduction into the nucleus for an incompletely mantled or bald comet has been found to be very important in maintaining relatively higher ice surface temperatures, and hence fluxes, during those portions of the diurnal and orbital cycles which would otherwise be cooler. It is also shown to be at least one cause of post perihelion brightness asymmetries, especially in lower obliquity comets. Maximum heliocentric distances at which 1-μm dust, sand, pebbles, cobbles, and boulders can be permanently ejected from the subsolar point by H₂O (CO₂) are (in AU): 6.9 (16.8), 5.2 (11.5), 1.8 (3.0), 0.21 (0.34) and 0.07 (0.11), respectively. A detailed anatomy of temperature, gas and dust fluxes vs latitude and longitude for a homogeneous rotating comet with fixed axis is given for comparison with future observations. Most H₂O flux histories deduced from brightness data are found to be in reasonable agreement with the model, allowing for uncertainty in radius and albedo. A clear exception is Encke. It is shown that the large discrepancy between Encke's observed and model predicted fluxes, based on radar cross section, can be used to evaluate the extent of exposed ice (<10%). The model is then used to place an active area so as to explain a reported sharp drop in flux on approach to the Sun at 0.78 AU. An active area or areas, <10% of the comet's surface, centered near 65°N latitude appears indicated. Although cyclic mantles are generally indicated for the set of parameters we used, our results show that a global mantle only 1 to 3 cm thick (depending on the orbit) consisting of a full range of grain sizes can cause irresversible evolution to a noncometary body. We investigated the long-term evolution of such a postulated initially thinly mantled cometary object. It was found that after the first few passes and until the end of its dynamic lifetime the object averaged <3 × 10^?12 g cm^?1 sec^?1 H₂O flux. Therefore, if cometary objects evolve into Apollo asteroids, ice should always be accessible within 10 m of the surface despite numerous close perihelion passages. The possible impact of factors not included in the model, such as initial inhomogeneities, coma scattering of radiation, and global redistribution of ejected silicate around the nucleus, are discussed.     

The 1995–2002 Long-Term Monitoring of Comet C/1995 O1 (HALE–BOPP) at Radio Wavelength

The bright comet Hale–Bopp provided the first opportunity to follow the outgassing rates of a number of molecular species over a large range of heliocentric distances. We present the results of our observing campaign at radio wavelengths which began in August 1995 and ended in January 2002. The observations were carried out with the telescopes of Nançay, IRAM, JCMT, CSO and, since September 1997, SEST. The lines of nine molecules (OH, CO, HCN, CH₃OH, H₂CO, H₂S, CS, CH₃CN and HNC) were monitored. CS, H₂S, H₂CO, CH₃CN were detected up to r_h= 3–4 AU from the Sun, while HCN and CH₃OH were detected up to 6 AU. CO, which is the main driver of cometary activity at heliocentric distances larger than 3–4 AU, was last detected in August 2001, at r_h= 14 AU. The gas production rates obtained from this programme contain important information on the nature of cometary ices, their thermal properties and sublimation mechanisms.Line shapes allow to measure gas expansion velocities, which, at large heliocentric distances, might be directly connected to the temperature of the nucleus surface. Inferred expansion velocity of the gas varied as r_h ^-0.4 within 7 AU from the Sun, but remained close to 0.4 km s^-1 further away. The CO spectra obtained at large r_hare strongly blueshifted and indicative of an important day-to-night asymmetry in outgassing and expansion velocity. The kinetic temperature of the coma, estimated from the relative intensities of the CH₃OH and CO lines, increased with decreasing r_h, from about 10 K at 7 AU to 110 K around perihelion.     

Estimating the nucleus density of Comet 19P/Borrelly

The nucleus bulk density of Comet 19P/Borrelly has been estimated by modeling the sublimation-induced non-gravitational force acting upon the orbital motion, thereby reproducing the empirical perihelion advance (i.e., the shortening of the orbital period). The nucleus has been modeled as a prolate ellipsoid, covered by various surface activity maps which reproduce the observed water production rate. The theoretical water production rate of active areas has been obtained by applying a sophisticated thermophysical model. This model takes into account net sublimation of ice and thermal reradiation from the surface, solid state conductivity, sub-surface sublimation and recondensation, mass and heat transport by diffusing gas, layer absorption of solar energy, a full treatment of local time-dependent illumination conditions, and a detailed consideration of nucleus/coma interaction mechanisms. The outgassing properties of the modeled nucleus are physically consistent with the gas kinetic structure of the innermost coma since the molecular backflux and surface gas density required in the thermophysical model (as functions of the nucleus surface temperature and the sub-surface temperature profile) have been obtained from Direct Simulation Monte Carlo modeling of inelastic intermolecular collisions in the cometary Knudsen layer. The calculation of local normal forces acting on the nucleus due to outgassing has been made within the same framework—recoil and/or impact momentum transfer to the nucleus caused by sublimating molecules and by recondensing and/or scattered coma molecules is therefore evaluated in accordance with local nucleus/coma conditions. According to this model, the density is found to be 100-300 kg, depending on the applied spin axis orientation and surface activity map. This range can be narrowed down to 180-300 kg by also requiring that the empirical changes (per orbital revolution) of the argument of perihelion and the longitude of the ascending node are reproduced.     

The rotation period of comet 29P/Schwassmann-Wachmann 1 determined from the dust structures (Jets) in the coma

The results of the photometric observations of comet 29P/Schwassmann-Wachmann 1 are analyzed. The comet demonstrates substantial activity at heliocentric distances larger than 5 AU, i.e., outside the water ice sublimation zone. The CCD images of the comet were obtained in wideband R filters at the 6-meter telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences (SAO RAS) and at the 2-meter Zeiss-2000 telescope of the Peak Terskol Observatory. The processing of the images with special digital filters allowed the active structures (jets) to be distinguished in the dust coma of the comet. With the cross-correlation method, the rotation period of the cometary nucleus was determined as 12.1 ± 1.2 and 11.7 ± 1.5 days for the observations taken in December 2008, and February 2009, respectively. The probable causes of the difference in the estimates of the rotation period of the cometary nucleus obtained by different authors are discussed.     

Isophot Observations Of Comet Hale-Bopp

Comet Hale-Bopp has been observed five times with ISOPHOT, the photometer on board the Infrared Space Observatory (ISO), four times before its perihelion passage at heliocentric distances of 4.92, 4.58, 2.93 and 2.81 AU, and at 3.91 AU postperihelion. Each time, multi-filter photometry covering the range between 3.6–175 μm with eight to ten filters was performed to sample the spectral energy distribution of the comet. These measurements were used to determine dust temperatures for the cometary coma. The evolution of the strength of the silicate feature can be followed in the data as well as the flux deficit at longer wavelengths. This revised version was published online in July 2006 with corrections to the Cover Date.     

Photoelectron fluxes in cometary atmospheres

The photoelectron fluxes in cometary atmospheres are calculated by a Monte Carlo method. This is the first quantitative model calculation of this kind. A pure H₂O atmosphere is assumed with a sublimation rate of 10³⁰ molecules sec^?1 at 1 AU. Discussions of the energetics of electron gas and of the elementary collisional processes in determining the fluxes largely concern this water atmosphere. Influences on the photoelectron fluxes are also investigated for CO, another possible constituent. The excitation rate of the ${}^1\text{D}$ level of atomic oxygen in electron impacts is evaluated. It is highly improbable that the photoelectrons are responsible for the observed 6300 Å emission of the order of 1 kR at a heliocentric distance of 1 AU. The structure of the heat equation for thermal electrons is analyzed and a drastic change of the plasma behavior within the coma region is expected.     

Origin of cometary extended sources from degradation of refractory organics on grains: polyoxymethylene as formaldehyde parent molecule

The radial distribution of some molecules (CO, H₂CO, HNC, …) observed in the coma of some comets cannot be explained only by a direct sublimation from the nucleus, or by the photolysis of a detected parent compound. Such molecules present a so-called extended source in comae. We show in this paper that extended sources can be explained by refractory organic material slowly releasing gas from grains ejected from the cometary nucleus, due to solar UV photons or heat. The degradation products are produced throughout the coma and therefore are presenting an extended distribution. To model this multiphase chemistry we derive new equations, which are applied to Comet 1P/Halley for the case of the production of formaldehyde from polyoxymethylene (POM), the polymer of formaldehyde (-CH₂-O-)_n. We show that the presence of a few percent of POM on cometary grains (a nominal value of ∼4% in mass of grains is derived from our calculations) is in good agreement with the observed distribution, which so far were not interpreted by the presence of any gaseous parent molecule.     

Thermal Emission From The Dust Coma Of Comet Hale-Bopp And The Composition Of The Silicate Grains

The dust coma of comet Hale-Bopp was observed in the thermal infrared over a wide range in solar heating (R = 4.9–0.9 AU) and over the full wavelength range from 3 μm to 160 μm. Unusual early activity produced an extensive coma containing small warm refractory grains; already at 4.9 AU, the 10 μm silicate emission feature was strong and the color temperature was 30% above the equilibrium blackbody temperature. Near perihelion the high color temperature, strong silicate feature, and high albedo indicated a smaller mean grain size than in other comets. The 8–13 μm spectra revealed a silicate emission feature similar in shape to that seen in P/Halley and several new and long period comets. Detailed spectral structure in the feature was consistent over time and with different instruments; the main peaks occur at 9.3, 10.0 and 11.2 μm. These peaks can be identified with olivine and pyroxene minerals, linking the comet dust to the anhydrous chondritic aggregate interplanetary dust particles. Spectra at 16–40 μm taken with the ISO SWS displayed pronounced emission peaks due to Mg-rich crystalline olivine, consistent with the 11.2 μm peak. This revised version was published online in July 2006 with corrections to the Cover Date.     

A comparison of modeled and observed intensity profiles for C2, CN, and the blue continuum for Comet Bradfield (1987s)

Intensity profiles were obtained for the C₂ and CN emission and blue continuum of Comet Bradfield (1987s), from observations obtained over a 10 week period starting shortly before perihelion. Model intensity profiles were produced and then fitted to the observed profiles, and used to put constraints on some of the dust and gas parameters. Most of these parameters, including the gas and dust outflow speeds from the cometary nucleus and the molecular lifetimes, were consistent with expected values. The best fitting models incorporate significant dust particle fragmentation and extended emission of CN from dust, both occurring in the inner coma. In addition, although there may have been enhancement of gas and dust emission on the sunward side of the cometary nucleus, it appears that the tailward side maintained a significant level of activity.     

On the evolution and activity of cometary nuclei

The thermal evolution of a spherical cometary nucleus (initial radius of 2.5 km), composed initially of very cold amorphous ice and moving in comet Halley's orbit, is simulated numerically for 280 revolutions. It is found that the phase transition from amorphous to crystalline ice constitutes a major internal heat source. The transition does not occur continuously, but in five distinct rounds, during the following revolutions: 1, 7, 40-41, 110-112, and 248-252. Due to the (slow) heating of the amorphous ice between crystallization rounds, the phase transition front advances into the nucleus to progressively greater depths: 36 m on the first round, and then 91 m, 193 m, 381 m, and 605 m respectively. Each round of crystallization starts when when the boundary between amorphous and crystalline ice is brought to approximately 15 m below the surface, as the nucleus radius decreases due to sublimation. At the time of crystallization, the temperature of the transformed ice rises to 180 K. According to experimental studies of gas-laden amorphous ice, a large fraction of the gas trapped in the ice at low temperatures is released. Whereas some of the released gas may find its way out through cracks in the crystalline ice layer, the rest is expected to accumulate in gas pockets that may eventually explode, forming "volcanic calderas." The gas-laden amorphous ice thus exposed may be a major source of gas and dust jets into the coma, such as those observed on comet Halley by the Giotto spacecraft. The activity of new comets and, possibly, cometary outbursts and splits may also be explained in terms of explosive gas release following the transition from amorphous to crystalline ice.