短周期彗星表面水冰升华分布研究

彗星是太阳系遗留的原始星子,研究彗星彗核的演化对理解太阳系其他天体的形成和演化历史具有重要意义.在太阳的辐射作用下,彗星携带的挥发性成分会发生升华,并带动尘埃运动,造成彗核物质的损失.因此,彗核的升华活动对其表面形貌甚至整体形状演化都会产生影响.从IAU (International Astronomical Union) MPC (Minor Planet Center)获取轨道数据,并考虑了彗核的自转以及进动,利用MONET (Mass lossdriven shape evolution model)形状演化模型对短周期彗星做数值模拟,计算得到了短周期彗星1P/Halley、9P/Tempel 1、 19P/Borrelly、 67P/C-G (Churyumov-Gerasimenko)、 81P/Wild 2和103P/Hartley 2在一个轨道周期内的太阳辐射能量以及表面侵蚀深度的分布,结合其动力学参数讨论了自转、进动和公转等特性对其表面水冰升华分布的影响以及造成南北侵蚀差异的可能性.     

Collisions of comets and meteoroids: The post Stardust‐NExT discussion

The paper considers results of collisions between comets and meteoroids. We re‐discuss the five different approaches to estimate the sizes of holes created during such collisions. The results of the Deep Impact and the Stardust‐NExT missions to comet 9P/Temple 1 are applied to the estimation of these methods. We use the observed amount of ejected mass, the jump of brightness of the comet 9P/Tempel 1 as well as the diameter of the excavated crater. In the paper the simple way of estimation of impact consequences by use of the conception of the fragmentation energy of comet is also discussed. The numerical calculations were carried out for reasonable assumed values of a large range of cometary characteristics. The main conclusion of this paper confirms a general presumption that the main factor which determines the size of the impact crater on the comet 9P/Tempel 1 is the kinetic energy of impactor and strength or fragmentation energy of cometary material. In the considered case the gravitation of a comet has a minor meaning (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The source of energy of the comet 29P/Schwassmann-Wachmann 1 outburst activity: the test of the summary

The comet 29P/Schwassmann-Wachmann 1 is an exceptional comet as far as cometary outbursts are concerned. Despite its large distance from the Sun (about 6 au), it shows quasi-regular outburst activity, usually once or twice a year. Up to now there has not been a generally accepted model that explains this phenomenon. In the first part of this paper, the most well-known hypotheses that attempt to explain the outburst activity of this comet are presented and critically analysed. The main aim of this paper is to present a model for the outburst activity of this comet. The model is based on the global analysis of the internal structure and physical and chemical processes that take place in the cometary nucleus. Numerical calculations were carried out for reasonable assumed values of a large range of cometary characteristics. The obtained results are consistent with observational data.     

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)     

Local subgiants and time-scales of disc formation

This paper concentrates on the relationship between the rate of gas emission from the nucleus of Comet 9P/Tempel 1, the fraction f of the nucleus that is active, and the crater damage inflicted by the recent 2005 July 4 Deep Impact space mission. The cometary nucleus has a surface area of about  1.7 × 10⁸ m²  and a mean radius of about 3700 m. Before the impact it is estimated that only a fraction f = 0.0056 of the nucleus surface was actively producing gas and dust. The active area was about  9.4 × 10⁵ m²  . Absolute magnitudes obtained at recent perihelion passages of this comet indicate that variations in the 0.0074 > f > 0.0039 range can occur from apparition to apparition. Because of the low size of the original active area, the production of a new impact crater in the diameter range 40 to 300 m would lead to a long-term change in the cometary visual magnitude in the range 0.0018 to 0.098 respectively. This is below the limit of detectability. It has been suggested that the cometary dust is in the form of 'talcum powder' not 'beach sand'. We suggest that the dust ejected from the impact site has been broken up by the energetic impact process and thus has a different size distribution from dust locked in the snowy matrix of the nucleus and normally lifted off the nucleus by gentle sublimation processes.     

The search for a cometary outbursts mechanism: a comparison of various theories

In the paper the potential sources of energy of cometary outbursts have been reviewed. Considerations focus on four probable sources of the outbursts' energy. These are the polymerization of hydrogen cyanide HCN, impacts with the meteoroids, destruction of cometary grains in the field of strong solar wind and the transformation of amorphous water Ice into the crystalline one. The values of released energy and jumps of cometary brightness caused by these mechanisms have been discussed. A modern approach to the problem of the thermodynamical evolution of the comet nucleus which includes amorphous water ice is considered as the starting point in the discussion presented in the paper. The main characteristics of an outburst of a hypothetical comet belonging to the Jupiter family comet are calculated. The obtained results are in a good agreement with the characteristics observed during the real outbursts of comets. The main conclusion of this paper confirms a general presumption that the cometary outbursts can have different causes. However, the hypothesis concerning the amorphous water ice transformation appears to be the most probable one. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Large cometary grains – their destruction and changes in the luminosity of comets

The catastrophic thermodynamic destruction of large cometary heterogeneous grains lying on the surface of a comet nucleus is examined. The core–mantle grain-structure model is assumed. Grain fragmentation as an explanation of sudden changes in cometary brightness is proposed. The approach presented to the problem of cometary outbursts is a development of a previous author's paper. The proposed mechanism is based on the idea of thermodynamical destruction of heterogeneous cometary grains. Numerical simulations have been carried out for a wide range of values of physical characteristics of cometary material. The results obtained are consistent with observational data. The main conclusion of this paper is that thermodynamical fragmentation of large grains can explain variations in brightness and also outbursts of comets.     

Cometary outbursts – search of probable mechanisms – case of 29P/Schwassmann‐Wachmann 1

Cometary outbursts, sudden increases in luminosity have not been clearly explained so far and their source is still a mystery. Various possible mechanisms as a source of cometary outbursts at large distances from the Sun have been considered. It has been stated that plausible mechanisms are the polymerization of HCN and the amorphous water ice transformation combined with electrostatic destruction of cometary grains in the head of the comet. The calculations have been carried out for a large range of cometary parameters and it has been shown that the proposed scenario of the outburst gives a jump in the comet brightness which is consistent with the real jump observed during the 29P/Schwassmann‐Wachmann 1 outbursts.     

The Electrostatic Blow-Off of Cometary Grains and Variations and Outbursts of Brightness of Comets at Large Heliocentric Distances

Interaction between the strong solar wind and the dark side of a cometary nucleus was considered. It was calculated that the potential of the dark side of the cometary nucleus could be numerically large and negative. Assuming that the nucleus of the comet has surface mantle which consists of loose, fine dust-ice particles, it was shown that cometary particles could electrostatically levitate over the nucleus. It was examined how this phenomenon affected the changes in the cometary brightness. Calculations were carried out for realistically assumed values of a large range of cometary parameters. It was shown that the considered mechanism could lead to the variations of cometary brightness, sometimes even to the outbursts of brightness.     

The outbursts of the comet 29P/Schwassmann‐Wachmann 1: A new approach to the old problem

As far as outbursts activity is concerned, the 29P/Schwassmann‐Wachmann 1 is the exceptional comet. This Centaur object shows quasi‐regular flares with periodicities of 50 days (Trigo‐Rodriguez et al. 2010). In the introductory part of the presented paper the most well‐known hypotheses which try to explain this cometary behaviour are reviewed. The second, actual part of this paper presents the new model for the outburst activity of this comet. The model is based on the idea of Ipatov (2012), according to which there are large cavities below a considerable fraction of the comet's surface containing material under high gas pressure. In favourite conditions the surface layers over the cavities are thrown away and the interior of these cavities is exposed. Consequently, an outburst of the comet's brightness may be observed. The main characteristics of an outburst of this comet, the brightness jump, is calculated. Numerical simulations were carried out for wide range of possible cometary parameters. The obtained results are in good agreement with the observations. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Short-period Jupiter family comets after Stardust

The NASA Stardust mission has provided for laboratory study an extensive data set of cometary dust of known provenance (from comet 81P/Wild 2) yielding detailed insights into the composition of the comet. Combined with the results of data from other missions to short-period Jupiter family comets (JFC), this has greatly deepened the understanding of such objects. If depressions on the surface of comet 81P/Wild 2 are all taken as evidence of impact cratering, their number suggests a long occupancy in the outer region of the Solar System. The dust from comet 81P/Wild 2 has been shown to be heavily deficient in pre-Solar grains and rich in materials formed at high temperatures in the inner Solar System. Although it is too early to know if this is typical of JFC, it does argue for rapid and thorough mixing of materials in the disk on timescales related to comet formation, and may also suggest outward migration of small icy bodies after their formation. Thus, instead of providing mainly new knowledge of the pre-Solar materials expected to be rich in comets, Stardust and comet 81P/Wild 2 have instead focussed attention on large-scale transport processes during the critical period when cometary parent bodies were forming in the early Solar System.     

The long-term evolution and initial size of comets 46P/Wirtanen and 67P/Churyumov–Gerasimenko

We present a new method to study the long-term evolution of cometary nuclei in order to estimate their original size, and we consider the case of comets 46P/Wirtanen (hereafter 46P) and 67P/Churyumov–Gerasimenko (hereafter 67P). We calculate the past evolution of the orbital elements of both comets over 100 000 yr using a Bulirsch–Stoer integrator and over 450 000 yr using a Radau integrator, and we incorporate a realistic model of the erosion of their nucleus. Their long-term orbital evolution is prominently chaotic, resulting from several close encounters with planets, and this result is independent of the choice of the integrator and of the presence or not of non-gravitational forces. The dynamical lifetime of comet 46P is estimated at ∼133 000 yr and that of comet 67P at ∼105 000 yr. Our erosion model assumes a spherical nucleus composed of a macroscopic mixture of two thermally decoupled components, dust and pure water ice. Erosion strongly depends upon the active fraction and the density of the nucleus. It mainly takes place at heliocentric distances <4 au and lasts for only ∼7 per cent of the lifetime. Assuming a density of 300 kg m⁻³ and an average active fraction over time of 10 per cent, we find an initial radius of ∼1.3 km for 46P and ∼2.8 km for 67P. Upper limit are obtained assuming a density of 100 kg m⁻³ and an active fraction of 100 per cent, and amounts to 21 km for 46P and 25 km for 67P. Erosion acts as a rejuvenating process of the surface so that exposed materials on the surface may only contain very little quantities of primordial materials. However, materials located just under it (a few centimetres to metres) may still be much less evolved. We will apply this method to several other comets in the future.     

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.     

Cometary outbursts – the post-Deep Impact outlook on collisions as possible causes

The possibility of impacts between comets belonging to the Jupiter Family and other small bodies orbiting in the main asteroid belt, and the consequences in relation to cometary activity are discussed. The probability of such events and the jumps in cometary brightness caused by impacts are examined. The results are compared with the results of the Deep Impact mission to Comet 9P/Tempel 1. The main conclusion of this paper is in agreement with previous findings, namely that an impact mechanism cannot be the main cause of the outburst activity of comets.     

The influence of the nucleus shape on the maximum size of grains ejected from a comet by gentle sublimation and jet‐like features

This paper deals with obtaining the maximum size of cometary grains ejected from nuclei of different shapes. Two mechanisms in terms of grain ejection from comets are taken into consideration. The first one is dragging of particles by outflowing gas molecules released by gentle sublimation from the comets. The second one is related with gas jets from the cavities in a nucleus by cometary jet‐like phenomena. We focused on ellipsoidal shapes of cometary nuclei but with different flattening. Calculations have been carried out for a large range of cometary parameters. It has been shown that for fixed mass of the nucleus the maximum size of grains is an increasing function of the nucleus flattening. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Dust in Comet C/1999 S4 (LINEAR) during Its Disintegration: Narrow-Band Images, Color Maps, and Dynamical Models

Comet C/1999 S4 was observed with the 2m-telescopes of the Bulgarian National Observatory and Pik Terskol Observatory, Northern Caucasus, Russia, at the time of its disintegration. Maps of the dust brightness and color were constructed from images obtained in red and blue continuum windows, free from cometary molecular emissions. We analyze the dust environment of Comet C/1999 S4 (LINEAR) taking into account the observed changes apparent in the brightness images and in plots of Afρ profiles as function of the projected distance ρ from the nucleus. We also make use of the syndyne-synchrone formalism and of a Monte Carlo model based on the Finson-Probstein theory of dusty comets. The brightness and color of individual dust particles, which is needed to derive theoretical brightness and color maps of the cometary dust coma from the Monte Carlo model, is determined from calculations of the light scattering properties of randomly oriented oblate spheroids. In general, the dust of Comet C/1999 S4 (LINEAR) is strongly reddened, with reddening values up to 30%/1000 Å in some locations. Often the reddening is higher in envelopes further away from the nucleus. We observed two outbursts of the comet with brightness peaks on July 14 and just before July 24, 2000, when the final disintegration of the comet started. During both outbursts an excess of small particles was released. Shortly after both outbursts the dust coma “turns blue.” After the first outburst, the whole coma was affected; after the second one only a narrow band of reduced color close to the tail axis was formed. This difference is explained by different terminal ejection speeds, which were much lower than normal in case of the second outburst. In particular in the second, final outburst the excess small particles could originate from fragmentation of “fresh” larger particles.     

Periodic variation of Oort Cloud flux and cometary impacts on the Earth and Jupiter

Time variation in impact probability is studied by assuming that the periodic flux of the Oort Cloud comets within 15 au arises from the motion of the Sun with respect to the Galactic mid-plane. The periodic flux clearly shows up in the impact rate of the captured Oort Cloud cometary population, with a phase shift caused by the orbital evolution. Depending on the assumed flux of comets and the size distribution of comets, the impact rate of the Oort Cloud comets of 1 km in diameter or greater is from 5 to 700 impacts Myr⁻¹ on the Earth and from 0.5 to 70 impacts per 1000 yr on Jupiter. The relative fractions of impacts are 0.09, 0.11, 0.26 and 0.54 for long-period comets, Halley type comets, Jupiter family comets and near-Earth objects, respectively. For Jupiter, the corresponding fractions in the first three categories are 0.18, 0.31 and 0.51. If we consider physical fading of comet activity that is compatible with the observations, then the impact rates of active comets are two orders of magnitude smaller than the total impact rates by all kinds of comets and cometary asteroids of size 1 km or greater.     

Cometary masses derived from non-gravitational forces

We compute masses and densities for 10 periodic comets with known sizes: 1P/Halley, 2P/Encke, 6P/d'Arrest, 9P/Tempel 1, 10P/Tempel 2, 19P/Borrelly, 22P/Kopff, 46P/Wirtanen, 67P/Churyumov–Gerasimenko and 81P/Wild 2. The method follows the one developed by Rickman and colleagues, which is based on the gas production curve and on the change in the orbital period due to the non-gravitational force. The gas production curve is inferred from the visual light curve. We found that the computed masses cover more than three orders of magnitude:  ≃(0.3–400) × 10¹²  kg. The computed densities are in all cases very low (≲0.8 g cm⁻³), with an average value of 0.4 g cm⁻³, in agreement with previous results and models of the cometary nucleus depicting it as a very porous object. The computed comet densities turn out to be the lowest among the different populations of Solar system minor bodies, in particular as compared to those of near-Earth asteroids (NEAs). We conclude that the model applied in this paper, in spite of its simplicity (as compared to more sophisticated thermophysical models applied to very few comets), is useful for a statistical approach to the mean density of the cometary nuclei. However, we cannot assess from this simple model if there is a real dispersion among the bulk densities of comets that could tell us about differences in physical structure (porosity) and/or chemical composition.     

A Monte Carlo Code to Compute Energy Fluxes in Cometary Nuclei

A Monte Carlo model designed to compute both the input and output radiation fields from spherical-shell cometary atmospheres has been developed. The code is an improved version of that by H. Salo (1988, Icarus76, 253-269); it includes the computation of the full Stokes vector and can compute both the input fluxes impinging on the nucleus surface and the output radiation. This will have specific applications for the near-nucleus photometry, polarimetry, and imaging data collection planned in the near future from space probes. After carrying out some validation tests of the code, we consider here the effects of including the full 4×4 scattering matrix in the calculations of the radiative flux impinging on cometary nuclei. As input to the code we used realistic phase matrices derived by fitting the observed behavior of the linear polarization as a function of phase angle. The observed single scattering linear polarization phase curves of comets are fairly well represented by a mixture of magnesium-rich olivine particles and small carbonaceous particles. The input matrix of the code is thus given by the phase matrix for olivine as obtained in the laboratory plus a variable scattering fraction phase matrix for absorbing carbonaceous particles. These fractions are 3.5% for Comet Halley and 6% for Comet Hale-Bopp, the comet with the highest percentage of all those observed.The errors in the total input flux impinging on the nucleus surface caused by neglecting polarization are found to be within 10% for the full range of solar zenith angles. Additional tests on the resulting linear polarization of the light emerging from cometary nuclei in near-nucleus observation conditions at a variety of coma optical thicknesses show that the polarization phase curves do not experience any significant changes for optical thicknesses τ?0.25 and Halley-like surface albedo, except near 90° phase angle.     

Are the cometary outbursts caused by impacts with interplanetary vagabonds probable?

The possibility of impacts and their results in relation to the cometary outbursts between comets and other small bodies in the solar system has been investigated. Taking into consideration certain physical features of cometary nuclei and impacting bodies, the probability of impacts of small bodies moving in the main asteroid belt with hypothetical comets which represent three types: Jupiter family comets, Halley family comets and long period comets has been computed. The probability of impacts between comets and meteoroids at large heliocentric distances has also been estimated. Potential consequences of these events in relation to outbursts of the cometary brightness have been discussed. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)