The Hypothesis of Additional Acceleration in the Motion of Comet Harrington–Abell from Jupiter

By processing 494 observations of Comet Harrington–Abell, we obtained a unified system of elements that includes its turn around the Sun during which it closely approached Jupiter to a minimum distance of 0.037 AU in 1974. A study of the cometary orbit before and after the approach showed that, probably, at the approach of the comet to Jupiter, apart from the well-known gravitational perturbations, its motion was affected by an additional force. An improvement of the cometary orbit by assuming that an additional acceleration inversely proportional to the square of the distance to Jupiter exists in its motion yielded the following values: (4.57 ± 0.42) × 10^–10 and (–7.20 ± 0.42) × 10^–10 AU day^–2 for the radial and transversal acceleration components, respectively. As a plausible explanation of the changes in the cometary orbit, we additionally considered a model based on the hypothesis of partial disintegration of the cometary nucleus. The parameter that characterizes the instant displacement of the center of inertia along the jovicentric radius vector was estimated to be –1.83 ± 0.75 km. Based on a unified numerical theory of cometary motion, we determined the nongravitational parameters using Marsden's model for two periods: A ₁ = (11.68 ± 1.74) × 10^–10 AU day^–2, A ₂ = (0.53 ± 0.0357) × 10^–10 AU day^–2 for 1975–1999 and A ₁ = (5.92 ± 5.86) × 10^–10 AU day^–2, A ₂ = (0.08 ± 0.028) × 10^–10 AU day^–2 for 1955–1969, under the assumption that the nongravitational acceleration changed at the approach of the comet to Jupiter.     

Plasma clouds associated with Comet P/Borrelly dust impacts

The NASA DS1 spacecraft encountered Comet P/Borrelly on September 22, 2001 at a distance of ∼2171 km on the sunward side of the comet. The flyby speed was ∼16.5 km s⁻¹. Using high temporal resolution (50 μs) absolute electric field amplitude measurements from a ∼1 m dipole antenna, new features of plasma clouds created by cometary dust impacts have been detected. The pulses have 1/e exponential decays of ∼650 μs duration, exponentially shaped overshoots with rise times of ∼2 ms, and exponential-shaped overshoot decay times of ∼10 ms. Assuming a plasma temperature of 10⁴ K, these pulse features have been explained as plasma cloud space charge effects from the electron, proton and heavy ion portions of the clouds passing the antenna. Complex pulse shapes were also detected. These are believed to be due to either plasma cloud scattering off of the spacecraft, or to secondary impacts. Small electric pulses of duration 10-15 ms of cometary origin were detected but are presently unexplained. The electric component of the plasma wave spectra at closest approach had an f^−2.4 power law shape from 10 Hz to 1 kHz. The electron cyclotron frequency was approximately 1 kHz. One possible explanation of the wave spectrum is that whistler mode waves associated with phase steepened cometary plasma waves are dispersed, leading to the broad spectrum. Finally, based on the present results, a new type of low-cost, large-area dust detector is proposed.     

Solar wind interactions with Comet 19P/Borrelly

The Plasma Experiment for Planetary Exploration (PEPE) made detailed observations of the plasma environment of Comet 19P/Borrelly during the Deep Space 1 (DS1) flyby on September 22, 2001. Several distinct regions and boundaries have been identified on both inbound and outbound trajectories, including an upstream region of decelerated solar wind plasma and cometary ion pickup, the cometary bow shock, a sheath of heated and mixed solar wind and cometary ions, and a collisional inner coma dominated by cometary ions. All of these features were significantly offset to the north of the nucleus-Sun line, suggesting that the coma itself produces this offset, possibly because of well-collimated large dayside jets directed 8°-10° northward from the nucleus as observed by the DS1 MICAS camera. The maximum observed ion density was 1640 ion/cm³ at a distance of 2650 km from the nucleus while the flow speed dropped from 360 km/s in the solar wind to 8 km/s at closest approach. Preliminary analysis of PEPE mass spectra suggest that the ratio of CO⁺/H₂O⁺ is lower than that observed with Giotto at 1P/Halley.     

On the production of positive molecular ions in cometary comas

Positively charged molecular ions, such as H₂O⁺, which have been observed in cometary. comas, may be efficiently produced by the evaporation of positively charged clathrate grains of radii in the range 10^–6–10^–5 cm. Such grains may be expelled from nuclei of comets, along with gaseous molecules. Grain charging occurs via interaction with solar ultraviolet photons and/or solar wind protons. Observational data on the total quantities as well as the distributions of H₂O and H₂O⁺ in cometary comas are shown to be in accord with detailed model calculations.On leave from: Tata Institute of Fundamental Research, Bombay, India.     

Motion of Dust Structures in the Circumnuclear Region of Comet Hale–Bopp and Rotation of the Cometary Nucleus

The motion of dust structures in the circumnuclear region of comet Hale–Bopp is studied. About 270 envelope images were obtained with the AZT-8 reflector (D = 0.7 m, F = 28 m) and the Filin-3 image intensifier. We carried out our observations at the observational station of the Astronomical Observatory of Shevchenko Kiev State University in the village of Lesniki (near Kiev). The recording from the image-intensifier screen was made on Foto-100 film during 23 nights from March 24 to May 10, 1997. The circumnuclear region was imaged both in white light (without filters) and with IHW CO⁺ (_ef = 426 nm), C₃ (_ef = 496 nm), C₂ (_ef = 514 nm), and RC (red continuum, _ef = 684 nm) narrow-band interference filters. Based on our measurements of the radial expansion of dust structures, we determined the velocities, 0.61–1.99 km s^–1; accelerations, from –18.3 × 10^–3 to 4.0 × 10^–3 m s^–2; and rotation period of the cometary nucleus, 111.41^h ± 0.05^h.     

CoMA — An advanced space experiment forin situ analysis of cometary matter

Exploring one of the most pristine bodies in our solar system — a comet — with a spacecraft will be a great step towards a deeper understanding of our solar system's beginnings. We here present the advanced space experiment CoMA (cometarymatteranalyzer), which will be flown on NASA's cometary rendezvous and asteroid flyby mission CRAF. CoMA is a high resolution time-of-flight secondary ion mass spectrometer. It will analyze m-sized cometary dust grains and cometary gases with an unprecedented mass resolution and will yield data about the elemental, isotopic, and molecular composition.     

On the possibility of H2-near-infrared detection in comets with large gas production

Comets with large gas production offer a unique chance to observe a H₂-flux of about 10⁵ photon cm^–2 s^–1 sr^–1 (1 Rayleigh) at wavelengths 8497.4 Å, 8560.2 Å and 8747.9 Å-i.e., where photon counting methods are still applicable. In the following it will be shown that population of the vibrational levels, giving rise to these quadrupole overtone transitions, is dominated by photodissociation of methane, and that the emission even of quadrupole lines is not attenuated by collisional quenching. Wavelength scanning by ±1 Å is shown to be enough to discriminate between cometary and atmospheric emissions by phase-sensitive subtraction techniques. Solid angle of <10^–7 sr has to be used, whence follows that a large ground-based telescope combined with a tilting Fabry-Perotfilter is best suited for detection of the near-infrared H₂-emissions at reasonable counting rates and sufficient rejection of the atmospheric background.Since H₂ is supposed to be formed mainly by the photolysis of CH₄, the optimum time for detection will be during approach to perihelion when, because of its high vapor pressure, methane will vaporize from the cometary nucleus. Variation of the source strength of both CH₄ and its photolysis product H₂ with time are particularly valuable indicators for the structure of the nucleus, its thermal history and conditions of formation. A high-resolution tilting filter photometer, which allows phase-sensitive background subtraction was used for the first time for near-infrared observations on the dust coma of Comet Kohoutek (Barbieriet al., 1974). The same technique was successfully used for the determination of an upper limit for CH₄ production at 3.3 by airborne observations on the same comet (Cosmoviciet al., 1974).     

Ionization mechanisms of cometary comas

A new approach is considered to the problem of ionization of the inner comas of comets connected with two phenomena: meteor-like process due to cometary molecules — interplanetary meteoroids impacts and explosion-type process due to high-velocity collisions between cometary dust grains and interplanetary meteoroids. It is found that the efficiency of explosive ionization exceeds the efficiency of meteor ionization approximately 100 times. The explosive ionization may be possible mechanism for anomaly ionization of the inner comae of dusty comets like Halley 1986 III with the dust to gas production rate ratio more than 0.1.     

On the origin of cometary nuclei in the presolar nebula

If we assume that the cometary nuclei originated by the gravitational instability of a dust layer, which formed in the equatorial plane of the outer parts of the presolar nebula (PSN) during a period of approximate equilibrium between gravity, centrifugal force, and the pressure gradient, a simple relation is derived between the PSN's temperature and the upper limit to the mass of the planetesimals. It contains, besides the density of the cometary nuclei _p , only the fraction (by mass) of the condensable elements in the PSN, which became part of the dust particle disc, which, on the basis of available observational evidence on the solid particles in interplanetary and interstellar space and of theoretical considerations on the relationship between them and on the sedimentation process, is found to be of the order of ~10%; this estimate will require still further justification. Assuming a temperature in the range 15–20 K, an equatorial diameter of the PSN of 0.1 pc and _p few 0.1 g/cm³, upper limits for the planetesimal's mass of 10¹⁸g and for their radius of 10 km are obtained (on the basis of conservation of circulation, of mass and of angular momentum in the differentially rotating disc), in fair agreement with observation. With the dispersion of those parts of the PSN — of an assumed original mass of 2–3M —, which did not become part of the Sun or the planets, by the young Sun's activity, the planetesimals must have lost a large part of their gravitational binding energy and their orbits must have become so large (semimajor axis several 10⁴ A.U. or more, if not negative), that stellar perturbations produced the distribution in configurational and in velocity space now observed.Paper dedicated to Professor Hannes Alfvén on the occasion of his 70th birthday, 30 May, 1978.The earlier work done since about 1950 in the U.S.S.R. is described in Safronov (1972).     

Surface thermoelastic erosion of atmosphereless solar system bodies under bombardment by multicharge cosmic ray ions

The mechanism of ion-stimulated erosion of atmosphereless solar system bodies is suggested and investigated. A theoretical model for the brittle surface erosion resulting under the effect of multicharge ion cosmic rays is analyzed. It is shown that the thermoelastic waves originated in the energetic track of a very heavy ion can result in the near-surface stresses exceeding the dynamic tensile strength of the surface material for any atmosphereless solar system body. The thermoelastic wave surface arrival yields brittle erosion of the material and ejection of this latter fragments (the track-breaking process). Thus ejected dust grains have plano-oblong shape, average mass on the order of 10^–17 g and velocity up to 400 m/sec providing the surface erosion rate of 10^–1 ÷ 3 · 10² »/year (near the Earth orbit) which depends upon the surface material (rock or ice). Possible track-breaking consequences, in particular, presence of the dust fraction of ultramicron grains and their aggregates on the lunar surface are discussed. Near the bodies with the radii from 10 to 300 km predicted is the existence of extended dust cocoons consisting of ultramicron and submicron grains. Smaller objects (asteroids, comets, smallest satellites of planets, meteoroids, etc.) can serve sources of permanent dust wind of ultramicron and submicron sized grains escaping from their surfaces. The interplanetary dust yield owing to the ion-stimulated erosion of these bodies is not less than 10¹² g/year. Possible interpreting in the frames of track-breaking process some observational data and effects, including existence of dust grains with the mass of 10^–18 ÷ 10^–17 g near the Halley's comet and the nature of 2060 Chiron dust coma is discussed. To prove the theory, observational identification and investigation of dust phenomena complex related to the ion-stimulated erosion of atmosphereless bodies, suggested is employing extreme ultraviolet and far infrared/submillimeter wavelengths, as well as polarimetric methods.     

Expansion velocities of cometary gas

Heating processes are expected to strongly affect the structure and dynamics of cometary comas. A radial expansion velocity of less than 1 km s^–1 in the inner coma is quite compatible with a few km s^–1 in the outer regions of large comets.     

Radioactive heating and layered structure of cometary nuclei

Following the work of Whipple and Stefanik, radioactive heating by uranium, thorium and pottasium of a cometary nucleus is discussed. The assumed composition is that of interstellar medium. If thermal diffusivity is 10^–4 cm²s^–1, the central temperature of a nucleus with radius 10 km can be above 50 K, while if the thermal diffusivity is 5 × 10^–3, the central temperature can be only 25 K or so. Volatile gases such as N₂ and CO will flow toward the outer part of the nucleus and are lost in their first several approaches to the sun. This mechanism appears capable of explaining the depletion of N₂ and CO relative to the interstellar abundance. It is argued that unfamiliar activity of comet Bowell could be explained by sublimation mainly of N₂ and CO.     

Collective plasma effects and the electron number problem in solar hard X-ray bursts

Due to the relatively high stream densities involved, collective interactions with the ambient plasma are likely to be important for the electrons producing solar hard X-ray bursts. In thick- and thin-target bremsstrahlung models the most relevant process is limitation of the invoked electron beams by ion sound wave generation in the neutralizing reverse current established in the atmosphere. For the thick target model it is shown that typical electron fluxes are near the maximum permitted by stability of the reverse current so that ion-sound wave generation may be the process which limits the electron injection rate. On the other hand the chromospheric reverse current is sufficient to supply the large total number of electrons which have to be accelerated in the corona. For the thin target the low density of the corona severely limits the possible reverse current so that the maximum upward flux of fast electrons is probably much too small to explain X-ray bursts but compatible with observations of interplanetary electrons.A distinct class of model postulates a small number of electrons confined by resonant scattering in a dense coronal slab surrounding a current sheet with continuous stochastic acceleration offsetting collisional losses. The energetic aspects of such a situation described by Hoyng (1975) are developed here by addition of equations describing the slab geometry in terms of electron diffusion by whistler scattering and of the collisional damping of the accelerating Langmuir waves. Solution of these equations results in values for the fieldB(70–350 G), densityn ₀(2–5 × 10¹² cm ^–3), slab dimensions (10¹⁸ km² × 0.3–3 km) and relative Langmuir energy density (10^–3 – 10^–2) required to produce the observed range of bursts. It is pointed out, however, that there may be no real gain in electron number requirements since the fast electrons in the emitting slab would be constantly swept out along with the frozen-in plasma as dissipation proceeds so that a large total number of electrons is still required. It could in fact be that just such a coronal region is the injection mechanism for the thick-target model.On leave from Department of Astronomy, University of Glasgow, Scotland.     

The masers of E-type methanol in orion KL and SGR B2

Using a simplified model the statistical equilibrium and radiative transfer equations of E-type-CH₃OH are solved for Orion KL and SgrB2. According to our calculation results and the observation data taken by Matsakiset al. (1980) and Morimotoet al. (1985a, b), the physical conditions of both sources are estimated. In theJ ₂-J ₁ E methanol maser region of Orion KL, the density, kinetic temperature, dust temperature, and the fractional abundance are 0.8–2×10⁶ cm^–3, 150, 30–90 K, 0.8–8×10^–6. In the 4_–1-3₀ E and 5_–1-4₀ E methanol maser region of Sgr B2 the correspondance physical conditions above are 10⁴ cm³, 45, 23 K, and 7×10^–7, respectively.     

On the CN (0, 0) spectrum of Comet Mrkos 1957d

We have carried out an analysis of the (0, 0) vibrational band of the CN molecule in Comet Mrkos 1957d, including the effect of collisions. We found that the sum of the squares of the residuals can be reduced by a factor of ten, if collisions account for 46±3% of the population of the lower level. A rotational temperature can be assigned to the cometary gas. The best value found was 410±40 K. The best fit for the constantR ₁ was (1.07±0.10)×10^–4. The velocity of the comet was left as a free parameter. We found for it a value of 34.38±0.10 km s^–1. This result is in disagreement with the nuclear orbital velocity of 34.74 km s^–1. The discrepancy can be explained, if the CN molecules are ejected from the cometary nucleus preferentially in the sunward direction, with a mean velocity that corresponds to the above temperature.     

Dust outflow velocity and the gas-to-dust ratio in comets

The observational determination of coma outflow velocity for gaseous species is fairly straightforward using high-resolution spectroscopy. The determination of the outflow speed of the dust is much more difficult. Most sources cite Bobrovnikoff (1954). This brief report is not a strictly refereed publication, however, and mixes data from different comets.We present here a simple analysis of some data from the International Halley Watch (IHW) archive. Differences between continuum images from successive nights show dust jets and shells clearly. Their motion is apparent to first order from the edges of the features. The component of the dust outflow velocity perpendicular to the observer's line of sight may thus be determined. This is of course a lower limit on the dust outflow velocity. Many measurements, at different heliocentric distances (R), allow determination of the heliocentric dependence of the dust outflow velocity.We find that the dust outflow velocity in comet P/Halley varied as R ^–0.41 . If data from an outburst at 14 AU (Sekanina et al. 1992) is included in the fit, this dependence becomes R ^–0.55. This confirms the canonical (e.g. Delsemme 1982) inverse-square-root law, and supports the conclusion of Storrs et al. (1992) on the variability of cometary gas-to-dust ratios.Space Telescope Science Institute     

Opposition Effects in Brightness,Color, and Polarization of Comet 1P/Halley: Comparison with Atmosphereless Solar System Bodies

Systematic and uniform sets of photometric and polarimetric observations of comet 1P/Halley have been analyzed. The phase dependence of brightness for comet Halley was obtained at phase angles α ranging from 1.4° ≤ α ≤ 65°. The following parameters were determined: the amplitude of the opposition effect Δm = 0.75^m ± 0.06^m; the half-width at a half-maximum of intensity HWHM = 6.4° ± 1.6°; the linear phase coefficient β = 0.0045 ± 0.0001 mag/deg for α from 30° ≤ α ≤ 65°; and the phase angle at which a nonlinear increase in brightness starts, α _opp ≈ 31°. For the first time, the phase-angle dependence was obtained for the color of the dust of comet Halley: the color index BC-RC systematically decreases with increasing phase angle. Such a phase behavior of the dust color can be caused by the decrease in the mean size of dust particles that occurs when the comet approaches the Sun. For comet Halley, the negative polarization branch is almost symmetric; the minimal value of polarization is P _min = −1.54% at a phase angle α_min = 10.5°, and the inversion angle is α_inv = 21.7°. A comparative analysis of the phase functions of brightness and polarization has been performed for the cometary dust and atmosphereless bodies. Among the latter are low-albedo asteroids of the P and C types (102 Miriam and 47 Aglaja, respectively), as well as Deimos; high-albedo objects, such as the E-type asteroid 64 Angelina and the icy satellite of Jupiter Europa; and the Moon with its intermediate albedo. The possibility of a weak depression in the negative polarization branch of comets Halley and 47P/Ashbrook-Jackson at phase angles smaller than 2° is discussed.__________Translated from Astronomicheskii Vestnik, Vol. 39, No. 4, 2005, pp. 353–363.Original Russian Text Copyright © 2005 by Rosenbush.     

Electromagnetic velocity gradient instabilities in the magnetosphere

The velocity gradients of the contrastreaming electron beams observed in the Earth's magnetosphere can excite three types of ordinary mode instabilities, namely (i) B-resonance electron instability, (ii) ion cyclotron instability, and (iii) unmagnetized ion instability. The B-resonance electron instability occurs at small values of the shear parameter 10^–4<S<10^–3, whereS = [(1/e){dU _o(x)}/(dx)] (U ₀(x) and _e being the streaming velocity of the electron beams and the electron cyclotron frequency, respectively). Near the equatorial plane of the bouncing electron beams region, this instability can generate electromagnetic waves having frequenciesf(0.045–0.2) Hz and wavelentghs (0.5–10)km, and the wave magnetic field is polarised in a radial direction. This instability can also occur in the plasma sheet region during the earthwards and tailwards plasma flows events and can generate waves, with wave magnetic field polarised along north-south direction, in the frequency rangef(0.007–0.02) Hz with (10–100)km nearR=–35R _E . For 10^–3<S<10^–2, the ion cyclotron instability is excited and it can generate waves up to 5th harmonic or so of ion cyclotron frequency. ForS>10^–2, the unmagnetized ion instability is excited which can generate electromagnetic waves having frequences from 5 to 50 Hz and typical wavelengths (0.5–6)km. The growth rates of all the three velocity shear driven instabilities are reduced in the presence of cold background plasma. The turbulence generated by these instabilities may give rise to enhanced effective electron-electron and electron-ion collisions and broaden the bouncing electron beams.     

Characteristics of the dust-plasma interaction near Enceladus’ South Pole

We present RPWS Langmuir probe data from the third Enceladus flyby (E3) showing the presence of dusty plasma near Enceladus’ South Pole. There is a sharp rise in both the electron and ion number densities when the spacecraft traverses through Enceladus plume. The ion density near Enceladus is found to increase abruptly from about 10² cm⁻³ before the closest approach to 10⁵ cm⁻³ just 30 s after the closest approach, an amount two orders of magnitude higher than the electron density. Assuming that the inconsistency between the electron and ion number densities is due to the presence of dust particles that are collecting the missing electron charges, we present dusty plasma characteristics down to sub-micron particle sizes. By assuming a differential dust number density for a range in dust sizes and by making use of Langmuir probe data, the dust densities for certain lower limits in dust size distribution were estimated. In order to achieve the dust densities of micrometer and larger sized grains comparable to the ones reported in the literature, we show that the power law size distribution must hold down to at least 0.03 μm such that the total differential number density is dominated by the smallest sub-micron sized grains. The total dust number density in Enceladus’ plume is of the order of 10² cm⁻³ reducing to 1 cm⁻³ in the E-ring. The dust density for micrometer and larger sized grains is estimated to be about 10⁻⁴ cm⁻³ in the plume while it is about 10⁻⁶-10⁻⁷ cm⁻³ in the E-ring. Dust charge for micron sized grains is estimated to be about eight thousand electron charges reducing to below one hundred electron charges for 0.03 μm sized grains. The effective dusty plasma Debye length is estimated and compared with inter-grain distance as well as the electron Debye length. The maximum dust charging time of 1.4 h is found for 0.03 μm sized grains just 1 min before the closest approach. The charging time decreases substantially in the plume where it is only a fraction of a second for 1 μm sized grains, 1 s for 0.1 μm sized grains and about 10 s for 0.03 μm sized grains.     

The photochemistry and dynamics of a dusty cometary atmosphere

A self-consistent solution of the dynamical and thermal structure of an H₂O-dominated, two-phase, dusty-gas cometary atmosphere has been obtained by solving the simultaneous set of differential equations representing conservation of number density, momentum and energy together with the transfer of solar radiation in the streams responsible for the major photolytic processes and the heating of the nucleus. The validity of the model is restricted to the collision-dominated region where all the gas species are assumed to attain a common velocity and common temperature. Two models are considered for the transfer of solar radiation through the circum-nuclear dust halo. In the first only the direct extinction by the dust is considered. In the second, the finding of some recent models, that the diffuse radiation field due to multiple scattering by the dust halo more or less compensates for radiation removed by direct absorption when the optical depth is near unity, is approximated by neglecting the attenuation of the radiation by the dust altogether.As has been shown earlier, the presence of dust results in a transonic solution, and it is obtained by a two-step iterative procedure which makes use of the asymptotic behaviour of the radiation fields sufficiently far from the nucleus and a regularity condition at the sonic point.The calculations were performed for a medium sized comet (R _n =2.5 km) having a dust to gas production rate ratio of unity, at a heliocentric distance of 1 AU. The dust grains were assumed to be of the same radius (1), of low density (1g cm^–3) and be strongly absorbing (having the optical properties of magnetite).The main effect of the dust on the cometary atmosphere is dynamic. While the dust-gas coupling persists to about 20R _n , the strong throat effect of the dust friction on the gas causes the latter to go supersonic quite rapidly. Consequently the sub-sonic region around the nucleus is very thin, varying between 45 and 85m in the two models considered. On the other hand, while this highly absorbing dust has a temperature substantially above that of the gas in the inner coma, heat exchange between them does not significantly change the temperature profile of the gas. This is because of the predominance of the expansion cooling, and even more importantly, the IR-cooling by H₂O, in the inner coma. Consequently, the gas temperature goes through a strong inversion, as in the dust-free case, achieving a temperature as low as about 6K within about 50km of the nucleus, before increasing to about 700K atr=10⁴km, due to the high efficiency of photolytic heating over the cooling process in the outer coma. The Mach number achieves a maximum value of about 10 at the distance of the temperature minimum, thereafter steadily decreasing to a value of about 2.5 atr10⁴km.It is shown that while the dust attenuation has a strong effect on the production rate of H₂O, it also has an interesting effect on the electron density profile. It increases the electron density in the inner coma over the unattenuated case, while at the same time, decreasing it in the outer coma. In conclusion, the limitations of the present model and the necessity to extend it using a multi-fluid approach are discussed.