Dust Morphology Of Comet Hale-Bopp (C/1995 O1). Ii. Introduction Of A Working Model

A Monte Carlo image simulation code for dust features in comets is applied to comet Hale-Bopp in order to model the object's persistent porcupine-like appearance on high-resolution images taken between May 11 and Nov. 2, 1996. A self-consistent fan model is proposed, with six isolated sources of dust emission assumed at various locations on the surface of the rotating nucleus and with the spin axis undergoing a complex motion in an inertial coordinate system. In the framework of this model, jet pairs represent boundaries of fan-shaped formations described by dust ejected from isolated sources during periods of time when the Sun is above the local horizon. The spin axis is found to have traveled through a field of 10° by 20° during the examined period of nearly six months. Still more successful is a fan model with large diurnal dust-emission fluctuations, which is consistent with an inertially fixed position of the spin axis and requires only three discrete sources. In this scenario, the dust-emission profile is dominated by several brief flare-ups, or “puffs”, in the production of dust from one of the sources. The results are insensitive to the spin rate, but the observed dust coma appearance is more typical of a rapidly rotating comet. This revised version was published online in July 2006 with corrections to the Cover Date.     

Spin Axis Direction of Comet 19p/Borrelly Based on Observations from 2000 and 2001

We calculate the direction of the rotational angular momentum vector,M, of comet 19P/Borrelly based on rotational lightcurve data from 2000, groundbased imaging of the coma during the Deep Space 1encounter, and the basic near-nucleus coma morphology as revealed by the Deep Space 1 spacecraft. For the most likely direction, we derivea family of solutions (with center at RA = 221°, Dec = -7°) if the direction of M is towards the sunward hemisphere during the Deep Space 1 encounter, whereas if the rotation is of opposite sense, the diametrically opposite family of solutions (with center at RA = 41°, Dec = 7°) would result. We argue that the coma morphology in September 2001 is consistent with the nucleus being a principal axis rotator or one observationally indistinguishable from it. Therefore, for all practical purposes, the direction of the rotational angular momentum vector coincides with the spin axis. We also discuss why the determination of the spin axis direction based on observations from the last apparition is in disagreement with the current result.     

Coma Morphology And Constraints On The Rotation Of Comet Hale–Bopp (C/1995 O1)

We present constraints on the spin state of comet Hale-Bopp based on coma morphology. Three cases of rotational states are compatible with near perihelion observations: (1) principal-axis rotation, (2) complex rotational state with a small precessional angle, or (3) complex rotational state with a large ratio between the component periods. For principal axis rotators, images from 1996 (pre-perihelion) are consistent with a rotational angular momentum vector, M, directed at ecliptic longitude and latitude (250°, -5°) while images from late 1997 (post-perihelion) indicate (310°, -40°). This may suggest a change in M. A complex rotational state with small precessional angle requires only a small or no change in M over the active orbital arc. In this case, M is directed near ecliptic longitude and latitude (270°, -20°). A rotationally excited nucleus with a large ratio between component periods requires the nucleus to be nearly spherical. The transformation of dust coma morphology from near-radial jets to bright arcs and then again to near-radial jets is interpreted as a heliocentric and geocentric distance dependent evolutionary sequence. The spiral structures seen in CN filters near perihelion (in contrast to sunward side arcs seen in continuum) can be explained if the precursor of CN molecules (likely sub-micron grains) are emitted from the nucleus at low levels (≈ 10% of the peak daytime emission) during the nighttime. This may be indicative of a nucleus with a CO-rich active area(s). This revised version was published online in July 2006 with corrections to the Cover Date.     

Rotational modeling of Hyperion

Saturn’s moon, Hyperion, is subject to strongly-varying solid body torques from its primary and lacks a stable spin state resonant with its orbital frequency. In fact, its rotation is chaotic, with a Lyapunov timescale on the order of 100 days. In 2005, Cassini made three close passes of Hyperion at intervals of 40 and 67 days, when the moon was imaged extensively and the spin state could be measured. Curiously, the spin axis was observed at the same location within the body, within errors, during all three fly-bys—~ 30° from the long axis of the moon and rotating between 4.2 and 4.5 times faster than the synchronous rate. Our dynamical modeling predicts that the rotation axis should be precessing within the body, with a period of ~ 16 days. If the spin axis retains its orientation during all three fly-bys, then this puts a strong constraint on the in-body precessional period, and thus the moments of inertia. However, the location of the principal axes in our model are derived from the shape model of Hyperion, assuming a uniform composition. This may not be a valid assumption, as Hyperion has significant void space, as shown by its density of 544± 50  kg m⁻³ (Thomas et al. in Nature 448:50, 2007). This paper will examine both a rotation model with principal axes fixed by the shape model, and one with offsets from the shape model. We favor the latter interpretation, which produces a best-fit with principal axes offset of ~ 30° from the shape model, placing the A axis at the spin axis in 2005, but returns a lower reduced χ ² than the best-fit fixed-axes model.     

Planet X and the origins of the shower and steady state flux of short-period comets

In the planet X model periodic comet showers are associated with passages of the planet's perihelion and aphelion points through a primordial disk of comets believed to lie beyond the orbit of Neptune. A strong feature of this model is that the required orbital elements and mass of planet X are consistent with independently predicted values based on the residuals in the motions of Uranus and Neptune. Here we present a more extensive analysis of the model taking into account the fact that only those comets scattered directly into the zones of influence of Saturn and Jupiter can contribute to a shower whose duration is consistent with observation (≲ 15 myr). These requirements impose a minimum planetary inclination of ≈25°, which in turn restricts the semimajor axis to be ≲100 AU. A fraction of the comets scattered directly into the zones of influence of Uranus and Neptune will evolve on time scales of ∼10⁸ years into the steady state flux of short-period comets. We find that the absolute numbers of shower and steady state are comparable and compatible with the known terrestrial cratering rate, assuming the existence of long-lived extinct comet cores. Canonical planet X model parameters, deduced in part from the scattering dynamics analysis, are: semimajor axis ≈80 AU, eccentricity ≈0.3, inclination ≈45°, and mass ≈5m_⊕. An analysis is given which suggests that planet X, in its present orbit, can create the requisite density gradient of comets near perihelion and aphelion during the lifetime of the Solar System. The required inclination of planet X's orbit (≳25°) may explain the failure of previous surveys to discover the planet as its present latitude is not likely to be near the ecliptic. It it exists, the best immediate hope of finding planet X is the ongoing IRAS search in the 100-μm band and the full sky optical survey by Shoemaker and Shoemaker. Independent of the question of periodic comet showers, the existence of planet X and the comet disk can readily explain the origin of the steady state flux of short-period comets over a wide range of parameters.     

New measures of the orientation of the linear structure in the sunward hemisphere of the coma of Comet 19P/Borrelly: apparitions of 1994, 2001, and 2008

The sky‐projected orientation (position angle) of the axis (line of maximum density or maximum brightness) of the long time‐known, linear structure (LS) in the sunward hemisphere of the coma of Comet 19P/Borrelly is measured on 45 photographs taken by different observers under different projection conditions and covering three consecutive apparitions (1994, 2001, and 2008) for a total time interval of 5174 days. The analysis of the results by a tomographic approach yields an LS axis constantly oriented towards a fixed point in the space, at Right Ascension 214°.4 ± 0°.5 and Declination –7°.0 ± 0°.5 (J2000), corresponding to an obliquity of 103°.5 ± 1° and an orbital longitude of 147°.2 ± 1°, throughout the relevant interval. Such coordinates are close to the ones found by other authors for the spatial orientation of the nucleus spin axis during the apparitions of 1994 and 2001. In the hypothesis of an LS orientation aligned with the nucleus spin axis, the new results confirm the previous ones and show that this orientation remained unchanged during the subsequent 2008 apparition (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Rotational dynamics of subsolar sublimating triaxial comets

A model with subsolar water sublimation on a triaxial, ellipsoidal comet nucleus is presented for the calculation of reactive torques. The resulting differential equations describing the comet's rotation are then Hamiltonian, and gravity-gradients are trivial to include. While effects derived from a weak perturbing function are neither able to change the rotational excitation nor the spin magnitude of the nucleus, it is shown how the spin orientation of comets can change significantly over an orbital run. However, of the four comets studied, 1P, 19P, 46P and the Rosetta target 67P, 19P and 46P were the only objects clearly exhibiting this feature, thereby confirming a technique used to derive the consequences of a more elaborate model of sublimation induced torques. In particular, the rotational parameters of 67P were seen to be very stable, indicating that a highly kinematical model of its rotation for the mapping of the comet's gravitational field during the Rosetta mission can be used. The model's hierarchy with 1P/Halley as the object with highest excitation probability, is consistent with observations.     

Tidal gravitational forces: The infall of “new” comets and comet showers

The tidal gravitational field of the Galaxy directed into the galactic plane changes the angular momentum of comets in the Oort cloud. For comet orbits with semimajor axis greater than 2 × 10⁴ AU, the change of angular momentum in one orbit is sufficient to bring comets from the Oort cloud into the visible region, causing the infall of “new” comets. The limiting size orbit is weakly dependent on the angle between the major axis of the comet orbit and the galactic plane. The flux of comets into the inner Solar System caused by the galactic tidal field will be continuous and nearly isotropic. This effect appears to exclude any determination of the trajectories of passing stars by analysis of the angular distribution of new comets. The production of intense comet showers by the tidal field of a solar companion or of an interstellar cloud is considered. We show that the direction of a solar companion cannot be found from the present distribution of observable comets. The frequency of comet showers induced by encounters with interstellar clouds is found to be much lower than that from passing stars, and the tidal fields of interstellar clouds are not strong enough to cause comet showers of sufficient intensity to result in Earth impacts.     

Comet McNaught (260P/2012 K2): spin axis orientation and rotation period

Extensive observations of comet 260P/McNaught were carried out between August 2012 and January 2013. The images obtained were used to analyze the comet’s inner coma morphology at resolutions ranging from 250 to about 1000 km/pixel. A deep investigation of the dust features in the inner coma allowed us to identify only a single main active source on the comet’s nucleus, at an estimated latitude of ?50^°±15^°. A thorough analysis of the appearance and of the motion of the morphological structures, supported by graphic simulations of the geometrical conditions of the observations, allowed us to determine a pole orientation located within a circular spot of a 15^°-radius centered at RA=60^°, Dec=0^°. The rotation of the nucleus seems to occur on a single axis and is not chaotic, furthermore no precession effects could be estimated from our measurements. The comet’s spin axis never reached the plane of the sky from October 2012 to January 2013; during this period it did not change its direction significantly (less than 30^°), thus giving us the opportunity to observe mainly structures such as bow-shaped jets departing from the single active source located on the comet’s nucleus. Only during the months of August 2012 and January 2013 the polar axis was directed towards the Earth at an angle of about 45^° from the plane of the sky; this made it possible to observe the development of faint structures like fragments of shells or spirals. A possible rotation period of 0.340±0.01 days was estimated by means of differential photometric analysis.     

Relationship of comets and planets

We analyze our earlier data on the numerical integration of the equations of motion for 274 short-period comets (with the period P<200 yr) on a time interval of 6000 yr. As many as 54 comets had no close approaches to planets, 13 comets passed through the Saturnian sphere of action, and one comet passed through the Uranian sphere of action. The orbital elements of these 68 comets changed by no more than ±3 percent in a space of 6000 yr. As many as 206 comets passed close to Jupiter. We confirm Everhart’s conclusion that Jupiter can capture long-period comets with q = 4–6 AU and i < 9° into short-period orbits. We show that nearly parabolic comets cross the solar system mainly in the zone of terrestrial planets. No relationship of nearly parabolic comets and terrestrial planets was found for the epoch of the latest apparition of comets. Guliev’s conjecture about two trans-Plutonian planets is based on the illusory excess of cometary nodes at large heliocentric distances. The existence of cometary nodes at the solar system periphery turns out to be a solely geometrical effect.     

Non-gravitational force modeling of Comet 81P/Wild 2: II. Rotational evolution

In this paper, we have studied both the dynamical and the rotational evolution of an 81P/Wild 2-like comet under the effects of the outgassing-induced force and torque. The main aim is to study if it is possible to reproduce the non-gravitational orbital changes observed in this comet, and to establish the likely evolution of both orbital and rotational parameters. To perform this study, a simple thermophysical model has been used to estimate the torque acting on the nucleus. Once the torque is calculated, Euler equations are solved numerically considering a nucleus mass directly estimated from the changes in the orbital elements (as determined from astrometry). According to these simulations, when the water production rate and changes in orbital parameters for 1997, as well as observational rotational parameters for 2004 are imposed as constraints, the change in the orbital period of 81P/Wild 2, , will decrease so that to , which is similar to the actual tendency observed from 1988 up to 1997. This nearly constant decreasing can be explained as due to a slight drift of the spin axis orientation towards larger ecliptic longitudes. After studying the possible spin axis orientations proposed for 1997, simulations suggest that the spin obliquity and argument (I,Φ)=(56°,167°) is the most likely. As for rotational evolution, changes per orbit smaller than 10% of the actual spin velocity are probable, while the most likely value corresponds to a change between 2 and 7% of the spin velocity. Equally, net changes in the spin axis orientation of 4°-8° per orbit are highly expected.     

Aperture Polarimetry And Photometry Of Comet Hale-Bopp

We present results of polarimetric and photometric observations of bright comet C/1995 O1 (Hale-Bopp) obtained at the 0.7 m telescope of Kharkov University Observatory from June 18, 1996 to April 24, 1997. The IHW and HB comet filters were used. The C₂ and C₃ production rates for Hale-Bopp are more than one order of magnitude larger and the dust production rates are more than two orders of magnitude larger than the Halley ones at comparable distances. Hence, Hale-Bopp was one of the most dusty comets. The average UC-BC and BC-RC colours of the dust were −0.02 and 0.13 mag, respectively. The polarization of comet Hale-Bopp at small phase angles of 4.8–13.0° was in good agreement with the date for comet P1/Halley at the same phase angles in spite of the fact that the heliocentric distances of comments differed nearly twice. However, at intermediate phase angles of 34–49° the polarization of comet Hale-Bopp was significantly larger than the polarization of the other dusty comets. It is the first case of such a large difference found in the continuum polarization of comets. The wavelength dependence of polarization for Hale-Bopp was steeper than for other dusty comets. The observed degree of polarization for the anti-sunward side of the coma was permanently higher than that for the sunward shell side. The polarization phase dependence of Hale-Bopp is discussed and compared with the polarization curves for other dusty comets. The peculiar polarimetric properties of comet Hale-Bopp are most likely caused by an over-abundance of small or/and absorbing dust particles in the coma. This revised version was published online in July 2006 with corrections to the Cover Date.     

Evolution of the Rotational State of Irregular Cometary Nuclei

A simplified thermal model has been used to calculate thenon-gravitational forces acting on small irregular nuclei in the orbit ofComet 46P/Wirtanen. The torque of thenon-gravitational force has beencalculated and the Euler equations have been solved in order toinvestigate the rotational evolution of several irregular nuclei duringa single orbital step. Several initial spin axis orientations andactivity patterns on their surfaces have been considered. The nucleiconsidered have a mean radius of 1 km and their inertia moments havebeen calculated assuming a homogeneous bulk density of 500 kg/m³. Inall the simulations, the initial spin period is 6h and the nucleiinitially rotate around their shortest axis. Under these assumptions,significant changes in the angular momentum and in spin period have beenobtained in all the simulations, but the nucleus is found to practicallyremain in its spin state of lowest energy during the entire orbitalperiod.     

The chaotic rotation of Hyperion

A plot of spin rate versus orientation when Hyperion is at the pericenter of its orbit (surface of section) reveals a large chaotic zone surrounding the synchronous spin-orbit state of Hyperion, if the satellite is assumed to be rotating about a principal axis which is normal to its orbit plane. This means that Hyperion's rotation in this zone exhibits large, essentially random variations on a short time scale. The chaotic zone is so large that it surrounds the ½ and 2 states, and libration in the 3/2 state is not possible. Stability analysis shows that for libration in the synchronous and ½ states, the orientation of the spin axis normal to the orbit plane is unstable, whereas rotation in the 2 state is attitude stable. Rotation in the chaotic zone is also attitude unstable. A small deviation of the principal axis from the orbit normal leads to motion through all angles in both the chaotic zone and the attitude unstable libration regions. Measures of the exponential rate of separation of nearby trajectories in phase space (Lyapunov characteristic exponents) for these three-dimensional motions indicate the the tumbling is chaotic and not just a regular motion through large angles. As tidal dissipation drives Hyperion's spin toward a nearly synchronous value, Hyperion necessarily enters the large chaotic zone. At this point Hyperion becomes attitude unstable and begins to tumble. Capture from the chaotic state into the synchronous or ½ state is impossible since they are also attitude unstable. The 3/2 state does not exist. Capture into the stable 2 state is possible, but improbable. It is expected that Hyperion will be found tumbling chaotically.     

Fan-shaped coma,orientation of rotation axis,and surface structure of a cometary nucleus I. Test of a model on four comets

Conspicuous anisotropy in the outgassing from comets, especially from short-period ones, appears to be the factor responsible for a frequent occurrence of a fan-shaped coma, extending in the general direction of the Sun. It is proposed that the pattern of deviations from the sunward direction contains information on the orientation of the spin axis and on the time lag in the sublimation process, which in turn provides insight into the nature of the nuclear surface. A simple model of a spherical rotating nucleus is formulated and a trial-and-error technique devised to determine the axis-orientation constants and a lage angle, a meaasure of the time lag in units of the rotation period. The results of application of this method to periodic comets Encke, Tempel 2, Borrelly, and Schwassmann-Wachmann 3 are presented. It is shown that the sense of rotation determined in this fashion is consistent with the results established for three of the four comets from the transverse component of the nongravitational force affecting their orbital motions. It is found that in general the time lag is strongly time dependent and that lag angles approaching 90° are rather common near perihelion, suggesting a complex surface structure that involves an insulting crust of dust of variable thickness and strength. These results are compared with the observed lightcurves of the four comets and with the calculated distributions of integrated insolation at the nuclear surface as functions of the cometocentric latitude and time. Noticed is a tendency of the comets to turn their spin axes to the Sun near perihelion and to replace, on the outbound leg of orbit, the established fan-orientation pattern by a “late”-tail pattern indicative of old, slowly accelerated particles. It is suggested that the motion of P/Schwassmann-Wachmann 3, which is due for a favorable return in 1979, was affected by a secular deceleration in 1930.     

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.     

Spectropolarimetry of comets Austin and Churyumov-Gerasimenko

Spectropolarimetric observations from 5000 to 8000 Å have been obtained for comets P/Austin (1982g) and P/Churyumov-Gerasimenko (1982f). The observations were spaced over phase angles of 50–125° for comet Austin and 10–40° for comet Churyumov-Gerasimenko. The use of spectropolarimetry allowed an evaluation of continuum polarization without molecular line contamination. Especially for comet Churyumov-Gerasimenko, the curve of polarization versus phase angle resembles curves for asteroids, where the polarization is negative (electric vector maximum parallel to the scattering plane) for phase angles less than 20° and the most negative polarization is from ?1 to ?2%. The negative polarization at backscattering angles may be due to multiple scattering in agglomerated grains, as assumed for asteroids, or to Mie scattering by small dielectric particles. If multiple scattering is important in comet dust, polarization measurements may imply a low albedo, less than 0.08. The polarization of comet Austin remained steady during a large change in the dust production rate. Both comets increased continuum flux by a factor of 2 near perihelion. The continuum of comet Churyumov-Gerasimenko had the shape of the solar spectrum with derivations less than 5%. The equivalent width of spectral features of C₂, NH₂, and O varied as r^?2.     

Stardust-NExT,Deep Impact,and the accelerating spin of 9P/Tempel 1

The evolution of the spin rate of Comet 9P/Tempel 1 through two perihelion passages (in 2000 and 2005) is determined from 1922 Earth-based observations taken over a period of 13 year as part of a World-Wide observing campaign and from 2888 observations taken over a period of 50 days from the Deep Impact spacecraft. We determine the following sidereal spin rates (periods): 209.023 ± 0.025°/dy (41.335 ± 0.005 h) prior to the 2000 perihelion passage, 210.448 ± 0.016°/dy (41.055 ± 0.003 h) for the interval between the 2000 and 2005 perihelion passages, 211.856 ± 0.030°/dy (40.783 ± 0.006 h) from Deep Impact photometry just prior to the 2005 perihelion passage, and 211.625 ± 0.012°/dy (40.827 ± 0.002 h) in the interval 2006–2010 following the 2005 perihelion passage. The period decreased by 16.8 ± 0.3 min during the 2000 passage and by 13.7 ± 0.2 min during the 2005 passage suggesting a secular decrease in the net torque. The change in spin rate is asymmetric with respect to perihelion with the maximum net torque being applied on approach to perihelion. The Deep Impact data alone show that the spin rate was increasing at a rate of 0.024 ± 0.003°/dy/dy at JD2453530.60510 (i.e., 25.134 dy before impact), which provides independent confirmation of the change seen in the Earth-based observations.The rotational phase of the nucleus at times before and after each perihelion and at the Deep Impact encounter is estimated based on the Thomas et al. (Thomas et al. [2007]. Icarus 187, 4–15) pole and longitude system. The possibility of a 180° error in the rotational phase is assessed and found to be significant. Analytical and physical modeling of the behavior of the spin rate through of each perihelion is presented and used as a basis to predict the rotational state of the nucleus at the time of the nominal (i.e., prior to February 2010) Stardust-NExT encounter on 2011 February 14 at 20:42.We find that a net torque in the range of 0.3–2.5 × 10⁷ kg m² s^?2 acts on the nucleus during perihelion passage. The spin rate initially slows down on approach to perihelion and then passes through a minimum. It then accelerates rapidly as it passes through perihelion eventually reaching a maximum post-perihelion. It then decreases to a stable value as the nucleus moves away from the Sun. We find that the pole direction is unlikely to precess by more than ～1° per perihelion passage. The trend of the period with time and the fact that the modeled peak torque occurs before perihelion are in agreement with published accounts of trends in water production rate and suggests that widespread H₂O out-gassing from the surface is largely responsible for the observed spin-up.     

The Influence of Reactive Torques on Comet Nucleus Rotation

Reactive torques, due to anisotropic sublimation on a comet nucleus surface, produce slow variations of its rotation. In this paper the secular effects of this sublimation are studied. The general rotational equations of motion are averaged over unperturbed fast rotation around the mass center (Euler-Poinsot motion) and over the orbital comet motion. We discuss the parameters that define typical properties of the rotational evolution and discover different classifications of the rotational evolution. As an example we discuss some possible scenarios of rotational evolution for the nuclei of the comets Halley and Borrelly.     

Galactic perturbations on nearly-parabolic cometary orbits

The effect of galactic perturbations on long-period comet orbits is examined via numerical and analytical means. Relations are found between a comet's initial perihelion position and the positions of succeeding perihelia. It was found that the galactic effects were strongest on the comets initially at galactic latitudes close to 40°. In such cases the galactic perturbations caused the orbit to become almost circular before becoming nearly parabolic again. This effect allows comets with semimajor axes of about 25 000 AU to make only a few passages through the inner solar system in a time interval of 10⁹yr. Thus the galactic field is an important factor in the evolution of long-period comet orbits. The observed distribution of perihelia of long-period comets indicates that galactic effects have been active.