Interpretation of HMC images by a combined thermal and gasdynamic model

Images of comet Halley's nucleus taken by the HMC camera during the GIOT-TO encounter in 1986 show that a major part of the total dust production is localized in a few active areas which are the sources of gas-dust jets. The global dust distribution in the inner coma is dominated by two main jets roughly directed to the sun. A combination of a 1D thermal nucleus model with an axisymmetric continuum model of the jet outflow was used to investigate the properties of the inner coma. Detailed investigations show that the characteristics of the observed jets can be reproduced by outgassing from free sublimating active areas of a few km in diameter, a dust to gas ratio of 1–2.5 and a size distribution dominated by the larger grains. It is further shown that most of the observational constraints provided by the HMC data can be met simultaneously by a model of three jets superimposed on a weak background.     

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.     

Interpretation of HMC images by a combined thermal and gasdynamic model

Images of comet Halley's nucleus taken by the HMC camera during the GIOT-TO encounter in 1986 show that a major part of the total dust production is localized in a few active areas which are the sources of gas-dust jets. The global dust distribution in the inner coma is dominated by two main jets roughly directed to the sun. A combination of a 1D thermal nucleus model with an axisymmetric continuum model of the jet outflow was used to investigate the properties of the inner coma. Detailed investigations show that the characteristics of the observed jets can be reproduced by outgassing from free sublimating active areas of a few km in diameter, a dust to gas ratio of 1–2.5 and a size distribution dominated by the larger grains. It is further shown that most of the observational constraints provided by the HMC data can be met simultaneously by a model of three jets superimposed on a weak background.     

Models of P/Wirtanen nucleus: active regions versus non-active regions

From the current understanding we know that comet nuclei have heterogeneous compositions and complex structures. It is believed that cometary activity is the result of a combination of physical processes in the nucleus, like sublimation and recondensation of volatile ices, dust grains release, phase transition of water ice, depletion of the most volatile components in the outer layers and interior differentiation.The evolution of the comet depends on the sublimation of ices and the release of different gases and dust grains: the formation of a dust crust, the surface erosion and the development of the coma are related to the gas fluxes escaping from the nucleus. New observations, laboratory experiments and numerical simulations suggest that the gas and dust emissions are locally generated, in the so-called active regions. This localized activity is probably superimposed to the global nucleus activity. The differences between active and inactive regions can be attributed to differences in texture and refractory material content of the different areas.In this paper we present the results of numerical models of cometary nucleus evolution, developed in order to understand which are the processes leading to the formation of active and non-active regions on the cometary surface. The used numerical code solves the equations of heat transport and gas diffusion within a porous nucleus composed of different ices—such as water (the dominant constituent), CO₂, CO- and of dust grains embedded in the ice matrix.By varying the set of physical parameters describing the initial properties of comet P/Wirtanen, the different behaviour of the icy and dusty areas can be followed.Comet P/Wirtanen is the target of the international ROSETTA mission, the cornerstone ESA mission to a cometary nucleus. The successful design of ROSETTA requires some knowledge of comet status and activity: surface temperatures, amount of active and inactive surface areas, gas production rate and dust flux.     

Dust tail streamers of a comet with rotating nucleus

Synthetic images of the dust tail are presented for a comet which has a rotating nucleus with one predominant dust source fixed to it. The images have been generated using a new computer model which, unlike similar models, allows for the study of dust tails caused by a rotating nucleus with an anisotropic distribution of sources.The dust tail is studied in the post-perihelion phase of a parabolic comet with a perihelion distance of 0.5 AU. One finds that in the case of a rotating nucleus with anisotropic emission characteristics streamers caused solely by the dynamics of the dust particles are forming in the dust tail even if there is no dependence between the solar irradiation angle of the source and the amount of dust emitted. If the dust emission depends on the solar irradiation angle of the dust source, then the brightest tail regions do not necessarily coincide with the synchrones for the times of maximum dust emission.As a consequence, a thorough analysis of streamer patterns in a cometary dust tail requires assumptions on the rotational state and the dust source distribution of the nucleus. Otherwise, it seems not possible to discern between streamers which are caused dynamically by nucleus rotation and others which reflect variations in the emission activity.     

Thermal destruction of cometary materials – application to activity of comets

Destruction mechanisms connected with thermodynamical behaviour of cometary material are reviewed with a special consideration of their effects on activity of comets. Consequences of thermal stresses which occur in the interior of a comet are discussed with reference to changes in the cometary brightness. Moreover, thermal destruction of grains placed in the head of the comet as well as on the surface of the nucleus is considered. It has been shown that the destruction of the cometary material can lead to an essential increase in the activity of the comet. Calculations have been carried out for a large assumed range of cometary parameters. The obtained simulated changes in the brightness of comets are consistent with the ones observed during the real variations and outbursts of brightness. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Splitting of Comet Halley

In combination with the authors' previous observation about the splitting of Comet Halley in March 1986,the events involving the sharp,straight feature in the antisolar direction observed in the head of Cornet halley in 1910 (such as those occurring on May 14,25 and 31,and June 2) are redis-cussed.The analysis leads to the following scenario;When comet Halley explodes and splits,a fragment jettisoned or thrown off from the nucleus will,after moving in the direction of its tail,develop into a mini-comet.Although not weel developed or permanent,it has its own plasma tail and,sometimes,a dust tail.If Bobrovnikoff's definition of a secondary nucleus is assumed,then the fragment should be considered as a real secondary nucleus.It seems that the current idea of a tailward jet suggested by Sekanina and Larson is a wrong explanation for the plasma tail of a mini-comet and hence the rotation period of 52-53h for Comet Halley is doubtful.     

Search for circum‐planetary material and orbital period variations of short‐period Kepler exoplanet candidates

A unique short‐period (P = 0.65356(1) d) Mercury‐size Kepler exoplanet candidate KIC012557548b has been discovered recently by Rappaport et al. (2012). This object is a transiting disintegrating exoplanet with a circum‐planetary material–comet‐like tail. Close‐in exoplanets, like KIC012557548b, are subjected to the greatest planet‐star interactions. This interaction may have various forms. In certain cases it may cause formation of the comet‐like tail. Strong interaction with the host star, and/or presence of an additional planet may lead to variations in the orbital period of the planet. Our main aim is to search for comet‐like tails similar to KIC012557548b and for long‐term orbital period variations. We are curious about frequency of comet‐like tail formation among short‐period Kepler exoplanet candidates. We concentrate on a sample of 20 close‐in candidates with a period similar to KIC012557548b from the Kepler mission. We first improved the preliminary orbital periods and obtained the transit light curves. Subsequently we searched for the signatures of a circum‐planetary material in these light curves. For this purpose the final transit light curve of each planet was fitted with a theoretical light curve, and the residuals were examined for abnormalities. We then searched for possible long‐term changes of the orbital periods using the method of phase dispersion minimization. In 8 cases out of 20 we found some interesting peculiarities, but none of the exoplanet candidates showed signs of a comet‐like tail. It seems that the frequency of comet‐like tail formation among short‐period Kepler exoplanet candidates is very low. We searched for comet‐like tails based on the period criterion. Based on our results we can conclude that the short‐period criterion is not enough to cause comet‐like tail formation. This result is in agreement with the theory of the thermal wind and planet evaporation (Perez‐Becker & Chiang 2013). We also found 3 cases of candidates which showed some changes of the orbital period. Based on our results we can see that orbital period changes are not caused by comet‐like tail disintegration processes, but rather by possible massive outer companions. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nucleus and tail studies of comet P/Swift-tuttle

CCD images of comet P/Swift-Tuttle, obtained in April 1994 with the 2.2m telescope at ESO La Silla/Chile, showed a comaless stellar nucleus. From absolute photometry we estimated the equivalent radius of the cometary nucleus to be about 11 km (assuming an albedo of 0.04 as for P/Halley) for two rotation phase angles which differ by about 75 deg. From that we conclude that the nucleus is either of rather spherical shape or that the viewing geometry was almost pole-on during our observations.An analysis of the plasma tail and inner coma of the comet by means of photographic plates and CCD images through IHW and BVR filters, obtained with the 80cm Schmidt camera and the 1.2m telescope at Calar Alto/Spain in November 1992, revealed several tail rays, head streamers and substructures in brightness excess areas in the coma. While some of the tail rays extended to several million km nuclear distance, most of them can be traced to starting points which lie in a region just 20000–35000 km projected distance tailward from the nucleus.     

Indication for precession of comet Halley's nucleus from dust tail analysis

The distinct structures ta comet Halley's dust tail around the beginning of March 1986 are analyzed by means of a computer simulation based on nucleus data obtained by the Giotto mission. It is shown that the assumption of a considerable free precession is required to understand the ground based dust tail observations from that time supposing a rotational period of some 50 hr. But a precession-free rotation with a period of about 7 days does not contradict an analysis of the dust tail structure. In both cases, an asymmetric distribution of the relevant emission sources is required.     

Helical waves in type-1 comet tails

Oscillations of type-1 comet tails with plasma compressibility taken into account are studied. A comet tail is treated as a plasma cylinder separated by a tangential discontinuity surface from the solar wind. The dispersion equation obtained in the linear approximation is solved numerically with typical plasma parameters. A sufficient condition for instability of the cylindrical tangential discontinuity in the compressible fluid is obtained. The phase velocity of helical waves is shown to be approximately coincident with Alfvén speed in the tail in the reference system moving with the bulk velocity of the plasma outflow in the tail. The instability growth rate is calculated.This theory is shown to be in good agreement with observations in the tails of Comets Kohoutek, Morehouse and Arend-Roland. Hence we conclude that helical waves observed in type-1 comet tails are produced due to the Kelvin-Helmholtz instability, and the model under consideration is justified. If so, one may estimate comet tail magnetic field from the pressure balance at the tangential discontinuity; it turns out to be of the order of the interplanetary magnetic field.     

Shape and obliquity effects on the thermal evolution of the Rosetta target 67P/Churyumov-Gerasimenko cometary nucleus

This work is dedicated to the application to 67P/Churyumov-Gerasimenko of a new quasi-3D approach for non-spherically shaped comet nuclei with the aim to interpret the current activity of the comet in terms of initial characteristics and to predict shape and internal stratification evolution of the nucleus. The model is applied to differently shaped nuclei taking into account the characteristics of Comet 67P/Churyumov-Gerasimenko deduced from observations. We focus our attention on the combined effects that shapes and obliquity have on the comet surface and sub-surface evolution. We discuss the results in terms of activity, local dust mantle formation and disruption, erosion of the surface and internal stratigraphy.The results show that differently shaped nuclei can have different internal structures leading to different activity patterns and behaviors. Our calculations have shown that local variations in the dust and gas fluxes can be induced by the nucleus shape. The distribution of “active” areas on Comet 67P/Churyumov-Gerasimenko is different because of different shapes, reflecting the illumination conditions on the surface. These shapes can influence the structure of the inner coma, but the coma far away from the nucleus is only marginally affected by the nucleus shape. However, different comet behaviors can arise from differently shaped comet nuclei, especially in terms of local activity, surface and sub-surface characteristics and properties. The water flux local distribution is the most influenced by the shape as it is directly linked to the illumination. Irregular shapes have large shadowing effects that can result in activity patterns on the comet surface.The effects of different pole directions are discussed to see the relations with the nucleus activity and internal structure. It is shown that the orientation of the rotation axis plays a strong role on the surface evolution of 67P/Churyumov-Gerasimenko, determining seasonal effects on the fluxes. The activity of the comet changes greatly with the nucleus obliquity leading to pre-post-perihelion differences in the activity and seasonal effects. The effects of the dust deposition and crust formation on the cometary activity have also been simulated and are discussed with respect to 67P/Churyumov-Gerasimenko observations. The dust mantling is also strongly obliquity dependent, with different surface distributions of the dust-covered regions according to the different comet pole orientations. Finally, we show that our model can reproduce the fluxes behavior near perihelion in terms of amplitude and asymmetry, and we estimate 20% of the illuminated surface to be active.     

Evolution of Comet Nucleus Rotation

The secular evolution of comet nucleus rotation states subject to outgassing torques is studied. The dynamical model assumes that the nucleus inertia ellipsoid is axially symmetric. The outgassing torques acting on the surface are modeled using standard cometary activity formulae. The general rotational equations of motion are derived and separately averaged over the fast rotational dynamics terms and the comet orbit. Special cases where the averaging assumptions cannot be applied are evaluated separately. The modification of the comet orbit due to comet outgassing is neglected. Resulting from this analysis is a system of secular differential equations that describes the dynamics of the comet nucleus angular momentum and rotation state. We find that the qualitative secular evolution of the rotation state is controlled by a single parameter that combines parameters related to the comet orbit and parameters related to the nucleus surface geometry and activity. From this solution, we find qualitatively different evolutionary paths for comet nuclei whose entire surface is active, as compared to nuclei with only a single active region. For surface activity models between these extremes, we show that certain evolutionary paths are more likely than others. Additionally, our solution indicates that a comet nucleus' rotational angular momentum will tend to increase over time, potentially contributing to the observed phenomenon of comet nucleus splitting.     

Model analysis of the dust tail of comet Hale-Bopp

The paper presents the results of dynamic simulation for the dust tail formation of comet C/1995 O1 (Hale-Bopp). To simulate the dust tail, the trajectories of 2 × 10⁶ dust particles were traced. The sizes, ejection moments, outflow directions and velocities of the dust particles were defined by the Monte Carlo algorithm. The obtained three-dimensional tail was projected on the sky plane to compare it with the observed images. The brightness distribution in the comet tail was fitted to similar model parameters for three different dates. According to our model experiments, the observed tails could be formed by particles with sizes from 0.3 to 8.0 μm, ejection velocities from 0.155 to 0.670 km/s, and power index of the exponential size distribution from −3.6 to −3.7. It is shown that the inclusion of the particles fragmentation processes leads to a noticeable improvement of the simulation results.     

Rosetta—one comet rendezvous and two asteroid fly-bys

One of the two planetary cornerstone missions of the European Space Agency is the Rosetta mission to comet 67P/Churyumov-Gerasimenko. Rosetta is a rendezvous mission with a comet nucleus, which combines an Orbiter with a Lander. It will monitor the evolution of the comet nucleus and the coma as a function of increasing and decreasing solar flux input along the comet’s pre- and post-perihelion orbit. Different instrumentations will be used in parallel, from multi-wavelength spectrometry to in-situ measurements of coma and nucleus composition and physical properties. Rosetta will go in orbit around the nucleus of its target comet 67P/Churyumov-Gerasimenko, when it is still far from the Sun and accompany the comet along its way to perihelion and beyond. In addition the Rosetta Lander Philae will land on the nucleus surface, before the comet is too active to permit such a landing (i.e. at around r = 3 AU) and examine the surface and subsurface composition of the comet nucleus as well as its physical properties.     

Activity and nucleus properties of 46 Pz.urule;Wirtanen

Astronomical observations of size and of related outgassing rates seem not to be compatible for the nucleus of comet 46 Pz.urule;Wirtanen, the target comet of the ROSETTA mission. This possible disagreement has caused speculations about peculiar properties of this comet nucleus. It is shown by model calculations which also takes into account vertical heat fluxes into the nucleus that there is a possibility to combine the results of astronomical observations within a model of a freely sublimating ice surface of this comet with an outgassing area of about half the dayside surface. The resulting half-size parameter (i.e. the radius of an equivalent sphere) can be shown to be of about R ≈ (725±230) m, and the nucleus is shown to have an active area of about half of the dayside surface, i.e. of about 25% of the total surface.     

On the observed polarization of comet Ikeya-Seki (1965 VIII)

It is shown that the grain material should be of the silicate type in order to produce the observed polarization reversal in the comet Ikeya-Seki. The polarization reversal is shown to arise as a result of the segregation of different sizes of grains in the tail of the comet due to the effect of radiation pressure.     

Physical characteristics of the plasma tail of comet C/2006 M4 (SWAN)

Wide-angle images with narrow-band filters isolating the radiation of CO⁺ ions and continuum of comet C/2006 M4 (SWAN) were acquired. The brightness distribution in the plasma tail shortly after the outburst of the comet on October 26, 2006, has been derived. The model brightness distribution was calculated, and the parameters of the diffusion model of the plasma tail have been obtained. For the date of observations, the magnetic-field induction in the plasma tail, as well as the lifetime of luminous particles and the lengthwise and transverse diffusion coefficients of ions have been estimated. Possible association between the change in the tail structure and that in the physical parameters along the plasma tail of the comet is discussed.     

Lineaments on the Surface of the Consolidated Material of the Comet 67P/Churyumov−Gerasimenko Nucleus

Morphologic analysis of the fragments (500 × 500 pixels) of images of comet 67P/Churyumov? Gerasimenko obtained with the OSIRIS camera for 18 regions, where the consolidated material is exposed on the surface, has been carried out. In terms of resolution, the images form a series from 0.04 to 1.29 m/pixel; consequently, the areas covered by these fragments vary from 400 to ~400000 m². In all of the regions, lineaments resembling tension fractures—several dozens of structures per region—are seen; and their number scarcely changes when passing from high-resolution images to lower-resolution ones. It is clear that relatively small lineaments cease to be reliably distinguished when the image resolution worsens, but the largest and well-defined ones remain observable, while the number of larger lineaments grows proportionally to the survey area increasing. Undoubtedly, this is an observational effect, which was demonstrated by examples with an artificially worsened resolution. However, on the other hand, this means that the tension fractures of the consolidated nucleus material represent a hierarchic population of smaller (meters long) to larger (decameters and longer) features. The lineaments seen in the analyzed fragments of images were counted, their lengths were measured, and the spacings (the mean distances between lineaments) and, from them, the depths of penetration of fractures into the nucleus material were estimated. It has been shown that the mean length of lineaments within each studied region depends on its area (which here correlates with the image resolution) and the depths of fracture penetration into the nucleus body depends on the mean length of lineaments. Both dependences are close to a power law. In the images of four regions covering the areas from 100000 to 400000 m² with a resolution of 0.66?1.29 m/pixel, the structures that look like layering or sheet jointing are seen in addition to fracture lineaments. The tension fractures are apparently formed due to seasonal and diurnal variations of the temperature, while the nature of the formation of the layer-like structures is not yet fully understood.     

On the cometary nature of near-Earth asteroid 2003 EH1

Differential evolution of the orbits of near-Earth asteroid (NEA) 2003 EH1 and comet 96P/Machholz 1 under perturbing action of planets was investigated for the time interval of 28000 years. The similarity of the orbits was analyzed with the Southworth–Hawkins criterion D _SH. It has been shown that both the comet and the asteroid can be fragments of a nucleus of the same larger comet being a progenitor of the Quadrantid complex. A break-up of the parent comet apparently occurred approximately 9500 years ago. NEA 2003 EH1 is actually a dormant fragment of a nucleus of the parent comet. It was concluded that comet 96P/Machholz 1, NEA (186256) 2003 EH1 of the Amore group, and the Quadrantid meteorite swarm form a family of related objects.