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

彗星是太阳系遗留的原始星子,研究彗星彗核的演化对理解太阳系其他天体的形成和演化历史具有重要意义.在太阳的辐射作用下,彗星携带的挥发性成分会发生升华,并带动尘埃运动,造成彗核物质的损失.因此,彗核的升华活动对其表面形貌甚至整体形状演化都会产生影响.从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在一个轨道周期内的太阳辐射能量以及表面侵蚀深度的分布,结合其动力学参数讨论了自转、进动和公转等特性对其表面水冰升华分布的影响以及造成南北侵蚀差异的可能性.     

Thermal modeling of cometary nuclei

A new model of the sublimation of volatile ices from a cometary nucleus has been developed which includes the effects of diurnal heating and cooling, rotation period and pole orientation, and thermal properties of the ice and subsurface layers. The model also includes the contribution from coma opacity, scattering, and thermal emission, where the properties of the coma are derived from the integrated rate of volatile production by the nucleus. The model is applied to the specific case of the 1986 apparition of Halley's comet. It is found that the generation of a cometary dust coma actually increases the total energy reaching the Halley nucleus. This results because of the significantly greater geometrical cross section of the coma as compared with the bare nucleus, and because the coma provides an essentially isotropic source of multiply scattered sunlight and thermal emission over the entire nucleus surface. For Halley, the calculated coma opacity is approximately 0.2 at 1 AU from the Sun, and 1.2 at perihelion (0.587 AU). At 1 AU this has little effect on dayside temperatures (maximum ≈200°K) but raises nightside temperatures (minimum ≈150°K) by about 40°K. At perihelion the higher opacity results in a nearly isothermal nucleus with only small diurnal and latitudinal temperature variations. The general surface temperature is 205°K with a maximum of 209°K at local noon on the equator. Some possible consequences of the results with respect to the generation of nongravitational forces, observed volatile production rates for comets, and cometary lifetimes against sublimation are discussed.     

Scattering Properties and Composition of Cometary Dust

Composition of the Comet dust obtained by the dust impact analyzer on the Halley probes indicated that the comet dust is a mixture of silicate and carbonaceous material. The collected interplanetary dust particles (IDP's) are fluffy and composite, having grains of several different types stuck together. Using discrete dipole approximation (DDA) we study the scattering properties of composite grains. In particular, we study the angular distribution of the scattered intensity and linear polarization of composite grains. We assume that the composite grains are made up of a host silicate sphere/spheroid with the inclusions of graphite. Results of our calculations on the composite grains show that the angle of maximum polarization shifts, and the degree of polarization varies with the volume fraction of the inclusions. We use these results on the composite grains to interpret the observed scattering in cometary dust.     

Infrared emission from comets

A brief discussion of the infrared observations from 4 to 20 micrometers of seven comets is presented. The observed infrared emission from comets depends primarily on their heliocentric distance. A model based on grain populations composed of a mixture of silicate and amorphous carbon particles in the mass ratio of about 40 to 1, with a power-law size distribution similar to that inferred for comet Halley, is applied to the observations. The model provides a good match to the observed heliocentric variation of both the 10 micrometers feature and the overall thermal emission from comets West and Halley. Matches to the observations of comet IRAS-Araki-Alcock and the antitail of comet Kohoutek require slightly larger grains. While the model does not match the exact profile and position of the 3.4 micrometers feature discovered in comet Halley, it does produce a qualitative fit to the observed variation of the feature's strength as a function of heliocentric distance. The calculations predict that the continuum under the 3.4 micrometers feature is due primarily to thermal emission from the comet dust when the comet is close to the Sun and to scattered solar radiation at large heliocentric distances, as is observed. A brief discussion of the determination of cometary grain temperatures from the observed infrared emission is presented. It is found that the observed shape of the emission curve from about 4 to 8 micrometers provides the best spectral region for estimating the cometary grain temperature distribution.     

Challenging a Paradigm: Do We Need Active and Inactive Areas to Account for Near-Nuclear Jet Activity?

We briefly describe an advanced 3D gas dynamical model developed for the simulation of theenvironment of active cometary nuclei. The model canhandle realistic nucleus shapes and alternative physical models for the gas and dust production mechanism.The inner gas coma structure is computed by solving self-consistently(a) near to the surface the Boltzman Equation(b) outside of it, Euler or Navier-Stokes equations.The dust distribution is computed from multifluid ``zero-temperature' Euler equations,extrapolated with the help of a Keplerian fountain model.The evolution of the coma during the nucleus orbital and spin motion,is computed as a succession of quasi-steady solutions. Earlier versions of the model using simple,``paedagogic' nuclei have demonstrated that the surface orographyand the surface inhomogeneity contribute similarly to structuring the near-nucleusgas and dust coma,casting a shadow on the automatic attribution of such structures to ``active areas'.The model was recently applied to comet P/Halley, for whichthe nucleus shape is available. In the companion paper of this volume,we show that most near-nucleus dust structuresobserved during the 1986 Halley flybys are reproduced, assuming that the nucleus is strictly homogeneous. Here, we investigate the effect of shape perturbations and homogeneityperturbations. We show that the near nucleus gas coma structure is robust vis-a-vissuch effects. In particular, a random distribution of active and inactive areaswould not affect considerably this structure, suggesting that such areas,even if present, could not be easily identified on images of the coma.     

An easy-to-use Model for the Optical Thickness and Ambient Illumination within Cometary Dust Comae

We define a procedure which allows estimation of the optical thickness of a cometarydust coma and the ambient illumination of the nucleus for any given comet, if estimatesof the nucleus radius and the dust activity (Afρ) are available. The calculation isperformed for a singly scattering coma with a cos(ϑ) distribution of dust overits day side. We find that the ambient illumination is of the same order as the incidentsunlight if the optical thickness is of order one. The optical thickness increases, all elseequal, linearly with the nucleus radius. Therefore the effect of the presence of the comamay be neglected for small (≈ 1 km diameter) comets, but is important forcomets such as 1P/Halley and Hale–Bopp.     

The Near-Nuclear Coma of Comet Halley in March 1986

The cameras carried onboard the flyby missions to comet P/Halleyin 1986 imaged the near nuclear jet activity fromseveral spatial directions. The observed, very structured near nucleardust jets were considered at that timeas the result of dust emission from well localized active surface regions(without supporting 3-D model computations, however).Based on the first, recently developed 3-D gas dynamical model ofP/Halley's activity,we have been shown that jet features can be reproduced assuming ahomogeneous dusty icenucleus surface. The dust in the collisional near nuclear comais concentrated along the gas flow discontinuities resulting from thecomplicated surface orography, creating the visual impression ofdust jets. We present here the results of these calculations forthe near nucleus dust distributions,and we compare them with the direct observations made during thethree Halley flybys (Vega 1, Vega 2, and Giotto).     

The composition of very fine dust in the dust shell of comet halley

An analysis of the spectra from the PUMA dust-impact mass spectrometers onboard the Vega-1 and Vega-2 spacecraft shows that a large number of the observed, unidentified small-amplitude peaks are produced by impacts of very-low-mass (from 10^?17 to 10^?20 g) particles. The mass flux of very fine particles accounts for a few percent of the total dust mass flux from comet Halley. The elemental composition of the finest cometary particles is identical to the composition of large particles (10^?12–10^?16 g), in agreement with present views about the nucleus of comet Halley as an aggregate of interstellar dust.     

Hamiltonian Formulation and Perturbations for Dust Motion Around Cometary Nuclei

In this paper we analyze the dynamical behavior of large dust grains in the vicinity of a cometary nucleus. To this end we consider the gravitational field of the irregularly shaped body, as well as its electric and magnetic fields. Without considering the effect of gas friction and solar radiation, we find that there exist grains which are static relative to the cometary nucleus; the positions of these grains are the stable equilibria. There also exist grains in the stable periodic orbits close to the cometary nucleus. The grains in the stable equilibria or the stable periodic orbits won’t escape or impact on the surface of the cometary nucleus. The results are applicable for large charge dusts with small area-mass ratio which are near the cometary nucleus and far from the Solar. It is found that the resonant periodic orbit can be stable, and there exist stable non-resonant periodic orbits, stable resonant periodic orbits and unstable resonant periodic orbits in the potential field of cometary nuclei. The comet gravity force, solar gravity force, electric force, magnetic force, solar radiation pressure, as well as the gas drag force are all considered to analyze the order of magnitude of these forces acting on the grains with different parameters. Let the distance of the dust grain relative to the mass centre of the cometary nucleus, the charge and the mass of the dust grain vary, respectively, fix other parameters, we calculated the strengths of different forces. The motion of the dust grain depends on the area-mass ratio, the charge, and the distance relative to the comet’s mass center. For a large dust grain (> 1 mm) close to the cometary nucleus which has a small value of area-mass ratio, the comet gravity is the largest force acting on the dust grain. For a small dust grain (< 1 mm) close to the cometary nucleus with large value of area-mass ratio, both the solar radiation pressure and the comet gravity are two major forces. If the a small dust grain which is close to the cometary nucleus have the large value of charge, the magnetic force, the solar radiation pressure, and the electric force are all major forces. When the large dust grain is far away from the cometary nucleus, the solar gravity and solar radiation pressure are both major forces.     

All comets are born equal: infrared emission by dust as a key to comet nucleus composition

The infrared emission of various comets can be matched within the framework that all comets are made of aggregated interstellar dust. This is demonstrated by comparing results on Halley (a periodic comet), Borrelly (a Jupiter family short period comet), Hale-Bopp (a long period comet), and extra-solar comets in the β Pictoris disk. Attempts have been made to generalize the chemical composition of comet nuclei based on the observation of cometary dust and volatiles and the interstellar dust model. Finally, we deduce some of the expected dust and surface properties of comet Wirtanen from the interstellar dust model as applied to other comets.     

In situ imaging of μm and sub-μm-sized grains in a cometary environment by atomic force microscopy

The micro-imaging dust analysis system (MIDAS) is an essential element among the scientific payload on the international Rosetta mission to comet 46P/Wirtanen. The MIDAS instrument based on an atomic force microscope (AFM) collects small particles drifting outwards from the nucleus surface. AFM is able to image small structures in 3D at nanometer-scale resolution. These images provide morphological and statistical information like grain size distribution on the dust population. In order to support the development of the flight hardware, optimisation of the control functions and consolidation of a proper scheme of data interpretation, laboratory studies with similar instruments were carried out. The obtained data demonstrate the capabilities of this technique. For the first time an instrument is able to observe the smallest (nm-sized) grains which are predicted by models and were to a certain extent deduced from previous measurements on the Giotto and Vega missions to comet 1P/Halley. On larger (μm-sized) particles the complex morphology will be visualised with high precision in 3D, and if present, within these aggregates crystalline materials with defined crystal faces can be identified.     

The dust distribution within the inner coma of comet P/Halley 1982i: encounter by Giotto's impact detectors

Analysis of the data from Giotto's Dust Impact Detection System experiment (DIDSY) is presented. These data represent measurement of the size of dust grains incident on the Giotto dust shield along its trajectory through the coma of comet P/Halley on 1986 March 13/14. First detection occurred at some 287000 km distance from the nucleus on the inbound leg; the majority of the DIDSY subsystems remained operational after closest approach (604 km) yielding the last detection at about 202000 km from the nucleus. In order to improve the data coverage (and especially for the smallest grains, to approximately 10(-19) kg particle mass), data from the PIA instrument has been combined with DIDSY data. Flux profiles are presented for the various mass channels showing, to a first approximation, a 1/R2 flux dependence, where R is the distance of the detection point from the cometary nucleus, although significant differences are noted. Deviations from this dependence are observed, particularly close to the nucleus. From the flux profiles, mass and geometrical area distributions for the dust grains are derived for the trajectory through the coma. Groundbased CCD imaging of the dust continuum in the inner coma at the time of encounter is also used to derive the area of grains intercepted by Giotto. The results are consistent with the area functions derived by Giotto data and the low albedo of the grains deduced from infrared emission. For the close encounter period (-5 min to +5 min), the cumulative mass distribution function has been investigated, initially in 20 second periods; there is strong evidence from the data for a steepening of the index of the mass distribution for masses greater than 10(-13) kg during passage through dust jets which is not within the error limits of statistical uncertainty. The fluences for dust grains along the entire trajectory is calculated; it is found that extrapolation of the spectrum determined at intermediate masses (cumulative mass index alpha = 0.85) is not able to account for the spacecraft deceleration as observed by the Giotto Radio Science Experiment and by ESOC tracking operations. Data at large masses (>10(-8) kg) recently analysed from the DIDSY data set show clear evidence of a decrease in the mass distribution index at these masses within the coma, and it is shown that such a value of the mass index can provide sufficient mass for consistency with the observed deceleration. The total particulate mass output from the nucleus of comet P/Halley at the time of encounter would be dependent on the maximum mass emitted if this change in slope observed in the coma were also applicable to the emission from the nucleus; this matter is discussed in the text. The flux time profiles have been converted through a simple approach to modeling of the particle trajectories to yield an indication of nucleus surface activity. There is indication of an enhancement in flux at t approximately -29 s corresponding to crossing of the dawn terminator, but the flux detected prior to crossing of the dawn terminator is shown to be higher than predicted by simple modelling. Further enhancements corresponding to jet activity are detected around +190 s and +270 s.     

Plasma environment of Jupiter family comets

As any comet nears the Sun, gas sublimes from the nucleus taking dust with it. Jupiter family comets are no exception. The neutral gas becomes ionized, and the interaction of a comet with the solar wind starts with ion pickup. This key process is also important in other solar system contexts wherever neutral particles become ionized and injected into a flowing plasma such as at Mars, Venus, Io, Titan and interstellar neutrals in the solar wind. At comets, ion pickup removes momentum and energy from the solar wind and puts it into cometary particles, which are then thermalised via plasma waves. Here we review what comets have shown us about how this process operates, and briefly look at how this can be applied in other contexts. We review the processes of pitch angle and energy scattering of the pickup ions, and the boundaries and regions in the comet-solar wind interaction. We use in-situ measurements from the four comets visited to date by spacecraft carrying plasma instrumentation: 21P/Giacobini-Zinner, 1P/Halley, 26P/Grigg-Skjellerup and 19P/Borrelly, to illustrate the process in action. While, of these, comet Halley is not a Jupiter class comet, it has told us the most about cometary plasma environments. The other comets, which are from the Jupiter family, give an interesting comparison as they have lower gas production rates and less-developed interactions. We examine the prospects for Rosetta at comet Churyumov-Gerasimenko, another Jupiter family comet where a wide range of gas production rates will be studied.     

Visible and near-infrared spectrophotometry of the Deep Impact ejecta of Comet 9P/Tempel 1

We have obtained optical spectrophotometry of the evolution of Comet 9P/Tempel 1 after the impact of the Deep Impact probe, using the Supernova Integral Field Spectrograph (SNIFS) at the UH 2.2-m telescope, as well as simultaneous optical and infrared spectra using the Lick Visible-to-Near-Infrared Imaging Spectrograph (VNIRIS). The spatial distribution and temporal evolution of the “violet band” CN (0-0) emission and of the 630 nm [OI] emission was studied. We found that CN emission centered on the nucleus increased in the 2 h after impact, but that this CN emission was delayed compared to the light curve of dust-scattered sunlight. The CN emission also expanded faster than the cloud of scattering dust. The emission of [OI] at 630 nm rose similarly to the scattered light, but then remained nearly constant for several hours after impact. On the day following the impact, both CN and [OI] emission concentrated on the comet nucleus had returned nearly to pre-impact levels. We have also searched for differences in the scattering properties of the dust ejected by the impact compared to the dust released under normal conditions. Compared to the pre-impact state of the comet, we find evidence that the color of the comet was slightly bluer during the post-impact rise in brightness. Long after the impact, in the following nights, the comet colors returned to their pre-impact values. This can be explained by postulating a change to a smaller particle size distribution in the ejecta cloud, in agreement with the findings from mid-infrared observations, or by postulating a large fraction of clean ice particles, or by a combination of these two.     

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.     

Infrared observations of Comet 81P/Wild 2 in 1997

Comet 81P/Wild 2 was observed in the thermal infrared over 6 months during its 1997 perihelion passage. The comet was most active in late February, about 3 months preperihelion; dust production declined by a factor of 3 between February and August. For the GIOTTO Halley dust size distribution, maximum dust production rate was ∼2 × 10⁶ g/s. The comet displayed a 10-μm silicate feature about 25% above the continuum, similar to several other Jupiter-family comets, but much lower than that seen in a number of Oort cloud comets.NASA’s STARDUST sample return mission will encounter P/Wild 2 98 days postperihelion in January 2004. Based on our observations at a similar point in the orbit and the Halley size distribution, we predict that the mass fluence on the spacecraft for a 150 km miss distance will be about 8 × 10⁻⁶ g/cm² for particles up to 1 cm in radius. The corresponding areal coverage will be about 10⁻⁴.     

The destruction of the cometary nucleus – The case of 73P/Schwassmann-Wachmann

The paper presents an analysis of the actual brightness change of comet 73P/Schwassmann-Wachmann, which took place in 1995. The consequence of a cometary outburst is the destruction of a fragment of its surface. This causes the emission of comet material from both the surface and from exposed subsurface layers. Therefore, the calculations take into account the scattering cross-sections that come from ice and dust particles. It was assumed that the dust particles are silicates which are characterized by high irregularity of their structure. This assumption is a consequence of the analysis of the results provided by the Rosetta mission to the comet 67P/Churyumov-Gerasimenko. The main factor determining the amplitude of a cometary outburst is the mass ejected as well as the loss of ice that holds the individual nucleus structures together. Consequently, this phenomenon can significantly contribute to the destruction and even decay of the cometary nucleus.     

The motion of dust and gas in the heads of comets with type II tails

Photographs of Comet Bennett 1969i taken in the dust-scattered continuum reveal that the dust particles, leading to the formation of the type II tail, leave the vicinity of the nucleus only within a certain cone with the aperture in the direction to the Sun. Three parabolic envelopes embracing the nucleus are formed by the dust (vertex always about on the radius vector) reaching distances from the nucleus of 30 000, 60 000 and 100 000 km.There exists no relation between the production and motion of this dust and the production and motion of the neutral coma gases. The cone of expulsion of the dust is identical with the cone of expulsion for the ions leading to the formation of the type I tail. Dust- and ion envelopes have, however, different kinematical properties. The cone of expulsion is identical with Bessel's Ausströmungskegel of visible matter observed by him in Comet Halley 1835.Comet Bennett is compared with Comet Halley 1910; they are related in many respects although Comet Halley had a lower dust production than the Comet Bennett.We ascribe to the dust particles of the tail II from the beginning of the expulsion an electrical charge.     

Large Dust Grains Around Cometary Nuclei

Large amounts of particles ejected from the nucleus surface are present in the vicinity of the cometary nuclei when comets are near the Sun (at heliocentric distances ≤2 AU). The largest dust grains ejected may constitute a hazard for spatial vehicles. We tried to obtain the bounded orbits of those particles and to investigate their stability along several orbital periods. The model includes the solar and the cometary gravitational forces and the solar radiation pressure force. The nucleus is assumed to be spherical. The dust grains are also assumed to be spherical, and radially ejected. We include the effects of centrifugal forces owing to the comet rotation. An expression for the most heavy particles that can be lifted is proposed. Using the usual values adopted for the case of Halley’s comet, the largest grains that can be lifted have a diameter about 5 cm, and the term due to the rotation is negligible. However, that term increases the obtained value for the maximum diameter of the lifted grain in a significant amount when the rotation period is of the order of a few hours.     

Insolation of a cometary crater at the stage of dust-jet formation

An important cause of the activation and development of active processes on the surface of a cometary nucleus is direct solar radiation illuminating a part of the surface that is not shielded by dust. The intensity of solar radiation near the surface of a cometary nucleus depends on the thickness of the dust cloud above the active area. If the size of the dust cloud noticeably changes, the intensity considerably depends on time. In the present paper, we consider the nonlinear equation of radiative transfer in a dust cloud growing towards the incident wave front with a constant velocity. The change in the intensity of direct solar radiation along the dust jet originating from the active surface area of a cometary nucleus has been found. For the sake of comparison, the linear equation of radiative transfer was solved in the framework of this task. It turns out that the linear approach to the solution of the considered problem suggests a noticeable loss in the amount of direct radiation participating in the dust-jet formation. This loss is comparable with the intensity of solar radiation incident to the active area of a cometary nucleus after scattering in the cometary atmosphere.