(47171) 1999 TC36, A transneptunian triple

We present new analysis of HST images of (47171) 1999 TC₃₆ that confirm it as a triple system. Fits to the point-spread function (PSF) consistently show that the apparent primary is itself composed of two similar-sized components. The two central components, A1 and A2, can be consistently identified in each of nine epochs spread over 7 years of time. In each instance, the component separation, ranging from 0.023 ± 0.002 to 0.031 ± 0.003 arcsec, is roughly one half of the Hubble Space Telescope’s diffraction limit at 606 nm. The orbit of the central pair has a semi-major axis of a ∼ 867 km with a period of P ∼ 1.9 days. These orbital parameters yield a system mass that is consistent with M_sys = 12.75 ± 0.06 × 10¹⁸ kg derived from the orbit of the more distant secondary, component B. The diameters of the three components are . The relative sizes of these components are more similar than in any other known multiple in the Solar System. Taken together, the diameters and system mass yield a bulk density of . HST photometry shows that component B is variable with an amplitude of ?0.17 ± 0.05 magnitudes. Components A1 and A2 do not show variability larger than 0.08 ± 0.03 magnitudes approximately consistent with the orientation of the mutual orbit plane and tidally distorted equilibrium shapes. The system has high specific angular momentum of J/J′ = 0.93, comparable to most of the known transneptunian binaries.     

A global solution for the Mars static and seasonal gravity, Mars orientation, Phobos and Deimos masses, and Mars ephemeris

With the collection of six years of MGS tracking data and three years of Mars Odyssey tracking data, there has been a continual improvement in the JPL Mars gravity field determination. This includes the measurement of the seasonal changes in the gravity coefficients (e.g., , , , , , ) caused by the mass exchange between the polar ice caps and atmosphere. This paper describes the latest gravity field MGS95J to degree and order 95. The improvement comes from additional tracking data and the adoption of a more complete Mars orientation model with nutation, instead of the IAU 2000 model. Free wobble of the Mars' spin axis, i.e. polar motion, has been constrained to be less than 10 mas by looking at the temporal history of and . A strong annual signature is observed in , and this is a mixture of polar motion and ice mass redistribution. The Love number solution with a subset of Odyssey tracking data is consistent with the previous liquid outer core determination from MGS tracking data [Yoder et al., 2003. Science 300, 299-303], giving a combined solution of k₂=0.152±0.009 using MGS and Odyssey tracking data. The solutions for the masses of the Mars' moons show consistency between MGS, Odyssey, and Viking data sets; Phobos GM=(7.16±0.005)×¹⁰⁻⁴ km³/s² and Deimos GM=(0.98±0.07)×¹⁰⁻⁴ km³/s². Average MGS orbit errors, determined from differences in the overlaps of orbit solutions, have been reduced to 10-cm in the radial direction and 1.5 m along the spacecraft velocity and normal to the orbit plane. Hence, the ranging to the MGS and Odyssey spacecraft has resulted in position measurements of the Mars system center-of-mass relative to the Earth to an accuracy of one meter, greatly reducing the Mars ephemeris errors by several orders of magnitude, and providing mass estimates for Asteroids 1 Ceres, 2 Pallas, 3 Juno, 4 Vesta, and 324 Bamberga.     

Optical and thermal infrared observations of six near-Earth asteroids in 2002

We present thermal infrared photometry and spectrophotometry of six Near-Earth Asteroids (NEAs) using the 3.8 m United Kingdom Infrared Telescope (UKIRT) together with quasi-simultaneous optical observations of five NEAs taken at the 1.0 m Jacobus Kapteyn Telescope (JKT). For Asteroid (6455) 1992 HE we derive a rotational period P=2.736±0.002 h, and an absolute visual magnitude H=14.32±0.24. For Asteroid 2002 HK₁₂ we derive . The Standard Thermal Model (STM), the Fast Rotating Model (FRM) and the Near-Earth Asteroid Thermal Model (NEATM) have been fitted to the measured fluxes to derive albedos and effective diameters. The derived geometric albedos and effective diameters are (6455) 1992 HE: pv=0.26±0.08, Deff=3.55±0.5 km; 1999 HF₁: pv=0.18±0.07, ; 2000 ED₁₀₄: pv=0.18±0.05, Deff=1.21±0.2 km; 2002 HK₁₂: , Deff=0.62±0.2 km; 2002 NX₁₈: pv=0.031±0.009, Deff=2.24±0.3 km; 2002 QE₁₅: , Deff=1.94±0.4 km. The limitations of using the NEATM to observe NEAs at high phase angles are discussed.     

(42355) Typhon-Echidna: Scheduling observations for binary orbit determination

We describe a strategy for scheduling astrometric observations to minimize the number required to determine the mutual orbits of binary transneptunian systems. The method is illustrated by application to Hubble Space Telescope observations of (42355) Typhon-Echidna, revealing that Typhon and Echidna orbit one another with a period of 18.971±0.006 days and a semimajor axis of 1628±29 km, implying a system mass of (9.49±0.52)×¹⁰¹⁷ kg. The eccentricity of the orbit is 0.526±0.015. Combined with a radiometric size determined from Spitzer Space Telescope data and the assumption that Typhon and Echidna both have the same albedo, we estimate that their radii are and , respectively. These numbers give an average bulk density of only , consistent with very low bulk densities recently reported for two other small transneptunian binaries.     

Titania's radius and an upper limit on its atmosphere from the September 8, 2001 stellar occultation

On September 8, 2001 around 2 h UT, the largest uranian moon, Titania, occulted Hipparcos star 106829 (alias SAO 164538, a V=7.2, K0 III star). This was the first-ever observed occultation by this satellite, a rare event as Titania subtends only 0.11 arcsec on the sky. The star's unusual brightness allowed many observers, both amateurs or professionals, to monitor this unique event, providing fifty-seven occultations chords over three continents, all reported here. Selecting the best 27 occultation chords, and assuming a circular limb, we derive Titania's radius: (1-σ error bar). This implies a density of using the value derived by Taylor [Taylor, D.B., 1998. Astron. Astrophys. 330, 362-374]. We do not detect any significant difference between equatorial and polar radii, in the limit , in agreement with Voyager limb image retrieval during the 1986 flyby. Titania's offset with respect to the DE405 + URA027 (based on GUST86 theory) ephemeris is derived: ΔαTcos(δT)=−108±13 mas and ΔδT=−62±7 mas (ICRF J2000.0 system). Most of this offset is attributable to a Uranus' barycentric offset with respect to DE405, that we estimate to be: and ΔδU=−85±25 mas at the moment of occultation. This offset is confirmed by another Titania stellar occultation observed on August 1st, 2003, which provides an offset of ΔαTcos(δT)=−127±20 mas and ΔδT=−97±13 mas for the satellite. The combined ingress and egress data do not show any significant hint for atmospheric refraction, allowing us to set surface pressure limits at the level of 10-20 nbar. More specifically, we find an upper limit of 13 nbar (1-σ level) at 70 K and 17 nbar at 80 K, for a putative isothermal CO₂ atmosphere. We also provide an upper limit of 8 nbar for a possible CH₄ atmosphere, and 22 nbar for pure N₂, again at the 1-σ level. We finally constrain the stellar size using the time-resolved star disappearance and reappearance at ingress and egress. We find an angular diameter of 0.54±0.03 mas (corresponding to projected at Titania). With a distance of 170±25 parsecs, this corresponds to a radius of 9.8±0.2 solar radii for HIP 106829, typical of a K0 III giant.     

On the population, physical decay and orbital distribution of Jupiter family comets: Numerical simulations

We study the Jupiter family comet (JFC) population assumed to come from the Scattered Disk and transferred to the Jupiter’s zone through gravitational interactions with the Jovian planets. We shall define as JFCs those with orbital periods and Tisserand parameters in the range 2<T?3.1, while those comets coming from the same source, but that do not fulfill the previous criteria (mainly because they have periods ) will be called ‘non-JFCs’. We performed a series of numerical simulations of fictitious comets with a purely dynamical model and also with a more complete dynamical-physical model that includes besides nongravitational forces, sublimation and splitting mechanisms. With the dynamical model, we obtain a poor match between the computed distributions of orbital elements and the observed ones. However with the inclusion of physical effects in the complete model we are able to obtain good fits to observations. The best fits are attained with four splitting models with a relative weak dependence on q, and a mass loss in every splitting event that is less when the frequency is high and vice versa. The mean lifetime of JFCs with radii and is found to be of about 150-200 revolutions (∼. The total population of JFCs with radii within Jupiter’s zone is found to be of 450±50. Yet, the population of non-JFCs with radii in Jupiter-crossing orbits may be ∼4 times greater, thus leading to a whole population of JFCs + non-JFCs of ∼2250±250. Most of these comets have perihelia close to Jupiter’s orbit. On the other hand, very few non-JFCs reach the Earth’s vicinity (perihelion distances ) which gives additional support to the idea that JFCs and Halley-type comets have different dynamical origins. Our model allows us to define the zones of the orbital element space in which we would expect to find a large number of JFCs. This is the first time, to our knowledge, that a physico-dynamical model is presented that includes sublimation and different splitting laws. Our work helps to understand the role played by these erosion effects in the distribution of the orbital elements and lifetimes of JFCs.     

Thermal infrared and optical observations of four near-Earth asteroids

We present thermal infrared photometry and spectrophotometry of four near-Earth asteroids (NEAs), namely (433) Eros, (66063) 1998 RO₁, (137032) 1998 UO₁, and (138258) 2000 GD₂, using the United Kingdom Infrared Telescope (UKIRT) in 2002. For two objects, i.e. (433) Eros and (137032) 1998 UO₁, quasi-simultaneous optical observations were also obtained, using the Jacobus Kapteyn Telescope (JKT). For (127032) 1998 UO₁, we obtain a rotation period P=3.0±0.1 h and an absolute visual magnitude HV=16.7±0.4. The Standard Thermal Model (STM), Fast Rotating Model (FRM) and near-Earth asteroid Thermal Model (NEATM) have been fitted to the IR fluxes to determine effective diameters Deff, geometric albedos pv, and beaming parameters η. The derived values are (433) Eros: Deff=23.3±3.5 km (at lightcurve maximum), pv=0.24±0.07, η=0.95±0.19; (66063) 1998 RO₁: , ; (137032) 1998 UO₁: Deff<1.13 km, pv>0.29; (138258) 2000 GD₂: Deff=0.27±0.04 km, , η=0.74±0.15. (66063) 1998 RO₁ is a binary asteroid from lightcurve characteristics [Pravec, P., and 56 colleagues, 2006. Icarus 181, 63-93] and we estimate the effective diameter of the primary (Dp) and secondary (Ds) components: and . The diameter and albedo of (138258) 2000 GD₂ are consistent with the trend of decreasing diameter for S- and Q-type asteroids found by Delbó et al. [Delbó, M., Harris, A.W., Binzel, R.P., Pravec, P., Davies, J.K., 2003. Icarus 166, 116-130]. A possible trend of increasing beaming parameter with diameter for small (less than about 3 km) S- and Q-type asteroids is found.     

Low-speed impacts between rubble piles modeled as collections of polyhedra, 2

We present the results of additional calculations involving the collisions of km-scale rubble piles. In new work, we used the Open Dynamics Engine (ODE), an open-source library for the simulation of rigid-body dynamics that incorporates a sophisticated collision-detection and resolution routine. We found that using ODE resulted in a speed-up of approximately a factor of 30 compared with previous code. In this paper we report on the results of almost 1200 separate runs, the bulk of which were carried out with 1000-2000 elements. We carried out calculations with three different combinations of the coefficients of friction η and (normal) restitution ?: low (η=0,?=0.8), medium (η=0,?=0.5), and high (η=0.5,?=0.5) dissipation.For target objects of ∼1 km in radius, we found reduced critical disruption energy values in head-on collisions from 2 to 100 J kg⁻¹ depending on dissipation and impactor/target mass ratio. Monodisperse objects disrupted somewhat more easily than power-law objects in general. For oblique collisions of equal-mass objects, mildly off-center collisions (b/b₀=0.5) seemed to be as efficient or possibly more efficient at collisional disruption as head-on collisions. More oblique collisions were less efficient and the most oblique collisions we tried (b/b₀=0.866) required up to ∼200 J kg⁻¹ for high-dissipation power-law objects. For calculations with smaller numbers of elements (total impactor or 200 elements) we found that collisions were more efficient for smaller numbers of more massive elements, with values as low as for low-dissipation cases. We also analyzed our results in terms of the relations proposed by Stewart and Leinhardt [Stewart, S.T., Leinhardt, Z.M., 2009. Astrophys. J. 691, L133-L137] where where Q_R is the impact kinetic energy per unit total mass m_i+m_T. Although there is a significant amount of scatter, our results generally bear out the suggested relation.     

Embedded star clusters and the formation of the Oort Cloud

Observations suggest most stars originate in clusters embedded in giant molecular clouds [Lada, C.J., Lada, E.A., 2003. Annu. Rev. Astron. Astrophys. 41, 57-115]. Our Solar System likely spent 1-5 Myrs in such regions just after it formed. Thus the Oort Cloud (OC) possibly retains evidence of the Sun's early dynamical history and of the stellar and tidal influence of the cluster. Indeed, the newly found objects (90377) Sedna and 2000 CR₁₀₅ may have been put on their present orbits by such processes [Morbidelli, A., Levison, H.F., 2004. Astron. J. 128, 2564-2576]. Results are presented here of numerical simulations of the orbital evolution of comets subject to the influence of the Sun, Jupiter and Saturn (with their current masses on orbits appropriate to the period before the Late Heavy Bombardment (LHB) [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461]), passing stars and tidal force associated with the gas and stars of an embedded star cluster. The cluster was taken to be a Plummer model with 200-400 stars, with a range of initial central densities. The Sun's orbit was integrated in the cluster potential together with Jupiter and Saturn and the test particles. Stellar encounters were incorporated by directly integrating the effects of stars passing within a sphere centred on the Sun of radius equal to the Plummer radius for low-density clusters and half a Plummer radius for high-density clusters. The gravitational influence of the gas was modeled using the tidal force of the cluster potential. For a given solar orbit, the mean density, 〈ρ〉, was computed by orbit-averaging the density of material encountered. This parameter proved to be a good measure for predicting the properties of the OC. On average 2-18% of our initial sample of comets end up in the OC after 1-3 Myr. A comet is defined to be part of the OC if it is bound and has q>35 AU. Our models show that the median distance of an object in the OC scales approximately as 〈ρ〉^−1/2 when . Our models easily produce objects on orbits like that of (90377) Sedna [Brown, M.E., Trujillo, C., Rabinowitz, D., 2004. Astrophys. J. 617, 645-649] within ∼1 Myr in cases where the mean density is or higher; one needs mean densities of order to create objects like 2000 CR₁₀₅ by this mechanism, which are reasonable (see, e.g., Guthermuth, R.A., Megeath, S.T., Pipher, J.L., Williams, J.P., Allen, L.E., Myers, P.C., Raines, S.N., 2005. Astrophys. J. 632, 397-420). Thus the latter object may also be part of the OC. Close stellar passages can stir the primordial Kuiper Belt to sufficiently high eccentricities (e?0.05; Kenyon, S.J., Bromley, B.C., 2002. Astron. J. 123, 1757-1775) that collisions become destructive. From the simulations performed it is determined that there is a 50% or better chance to stir the primordial Kuiper Belt to eccentricities e?0.05 at 50 AU when . The orbit of the new object 2003 UB₃₁₃ [Brown, M.E., Trujillo, C.A., Rabinowitz, D.L., 2005. Astrophys. J. 635, L97-L100] is only reproduced for mean cluster densities of the order of , but in the simulations it could not come to be on its current orbit by this mechanism without disrupting the formation of bodies in the primordial Kuiper Belt down to 20 AU. It is therefore improbable that the latter object is created by this mechanism.     

Radar and infrared observations of binary near-Earth Asteroid 2002 CE26

We observed near-Earth Asteroid (NEA) 2002 CE26 in August and September 2004 using the Arecibo S-band (2380-MHz, 12.6-cm) radar and NASA's Infrared Telescope Facility (IRTF). Shape models obtained based on inversion of our delay-Doppler images show the asteroid to be 3.5±0.4 km in diameter and spheroidal; our corresponding nominal estimates of its visual and radar albedos are 0.07 and 0.24, respectively. Our IRTF spectrum shows the asteroid to be C-class with no evidence of hydration. Thermal models from the IRTF data provide a size and visual albedo consistent with the radar-derived estimate. We estimate the spin-pole to be within a few tens of degrees of λ=317°, β=−20°. Our radar observations reveal a secondary approximately 0.3 km in diameter, giving this binary one of the largest size differentials of any known NEA. The secondary is in a near-circular orbit with period 15.6±0.1 h and a semi-major axis of 4.7±0.2 km. Estimates of the binary orbital pole and secondary rotation rate are consistent with the secondary being in a spin-locked equatorial orbit. The orbit corresponds to a primary mass of M=1.95±0.25×¹⁰¹³ kg, leading to a primary bulk density of , one of the lowest values yet measured for a main-belt or near-Earth asteroid.     

Photometric study of Centaur (60558) 2000 EC98 and trans-neptunian object (55637) 2002 UX25 at different phase angles

We present photometric observations of Centaur (60558) 2000 EC₉₈ and trans-neptunian object (55637) 2002 UX₂₅ at different phase angles and with different filters (mainly R but also V and B for some data). Results for 2000 EC₉₈ are: (i) a rotation period of 26.802±0.042 h if a double-peaked lightcurve is assumed, (ii) a lightcurve amplitude of 0.24±0.06 for the R band, (iii) a phase curve with H=9.03±0.01 and G=−0.39±0.08 (R filter) and H=9.55±0.04 and G=−0.50±0.35 (V filter) or a slope of (R filter) and 0.22±0.06 (V filter), (iv) the color indices B-V=0.76±0.15 and V-R=0.51±0.09 (for α=0.1-0.5°) and 0.55±0.08 (for α=1.4-1.5°). The rotation period is amongst the longest ever measured for Centaurs and TNOs. We also show that our photometry was not contaminated by any cometary activity down to magnitude ?27/arcsec². For 2002 UX₂₅ the results are: (i) a rotation period of 14.382±0.001 h or 16.782±0.003 h (if a double-peaked lightcurve is assumed) (ii) a lightcurve amplitude of 0.21±0.06 for the R band (and the 16.782 h period), (iii) a phase curve with H=3.32±0.01 and G=+0.16±0.18 or a slope of (R filter), (iv) the color indices B-V=1.12±0.26 and V-R=0.61±0.12. The phase curve reveals also a possible very narrow and bright opposition surge. Because such a narrow surge appears only for one point it needs to be confirmed.     

Estimating Mercury's 88-day libration amplitude from orbit

In this paper it is derived that the libration of Mercury can be described by where Φ₀ is the unknown libration amplitude, M is Mercury's mean anomaly and K=−9.483. Φ₀ can be determined by comparing pairs of images of the same landmarks taken by an orbiter at different positions of Mercury. If the angle between the orbit plane of a polar orbiter and Mercury's line of periapsis is between −60° and 60° and if one landmark at the equator is imaged per day with a relative precision of , then the libration amplitude can be determined in two Mercury years (176 days) with an accuracy of or better, which is sufficient to answer the question whether Mercury has a solid or fluid core.     

The orbit, mass, size, albedo, and density of (65489) Ceto/Phorcys: A tidally-evolved binary Centaur

Hubble Space Telescope observations of Uranus- and Neptune-crossing object (65489) Ceto/Phorcys (provisionally designated 2003 FX₁₂₈) reveal it to be a close binary system. The mutual orbit has a period of 9.554±0.011 days and a semimajor axis of 1840±48 km. These values enable computation of a system mass of (5.41±0.42)×¹⁰¹⁸ kg. Spitzer Space Telescope observations of thermal emission at 24 and 70 μm are combined with visible photometry to constrain the system's effective radius and geometric albedo . We estimate the average bulk density to be , consistent with ice plus rocky and/or carbonaceous materials. This density contrasts with lower densities recently measured with the same technique for three other comparably-sized outer Solar System binaries (617) Patroclus, (26308) 1998 SM₁₆₅, and (47171) 1999 TC₃₆, and is closer to the density of the saturnian irregular satellite Phoebe. The mutual orbit of Ceto and Phorcys is nearly circular, with an eccentricity ?0.015. This observation is consistent with calculations suggesting that the system should tidally evolve on a timescale shorter than the age of the Solar System.     

The Planetary Fourier Spectrometer (PFS) onboard the European Mars Express mission

The Planetary Fourier Spectrometer (PFS) for the Mars Express mission is an infrared spectrometer optimised for atmospheric studies. This instrument has a short wave (SW) channel that covers the spectral range from 1700 to (1.2-) and a long-wave (LW) channel that covers 250- (5.5-). Both channels have a uniform spectral resolution of . The instrument field of view FOV is about 1.6^° (FWHM) for the Short Wavelength channel (SW) and 2.8^° (FWHM) for the Long Wavelength channel (LW) which corresponds to a spatial resolution of 7 and 12 km when Mars is observed from an height of 250  km. PFS can provide unique data necessary to improve our knowledge not only of the atmosphere properties but also about mineralogical composition of the surface and the surface-atmosphere interaction.The SW channel uses a PbSe detector cooled to 200-220 K while the LW channel is based on a pyroelectric (L_iT_aO₃) detector working at room temperature. The intensity of the interferogram is measured every 150 nm of physical mirrors displacement, corresponding to 600 nm optical path difference, by using a laser diode monochromatic light interferogram (a sine wave), whose zero crossings control the double pendulum motion. PFS works primarily around the pericentre of the orbit, only occasionally observing Mars from large distances. Each measurements take 4 s, with a repetition time of 8.5 s. By working roughly 0.6 h around pericentre, a total of 330 measurements per orbit will be acquired 270 looking at Mars and 60 for calibrations. PFS is able to take measurements at all local times, facilitating the retrieval of surface temperatures and atmospheric vertical temperature profiles on both the day and the night side.     

Near infrared imaging spectroscopy of Venus with the Anglo-Australian Telescope

We have obtained full-disk spatially resolved spectra of the Venus nightside at near-infrared wavelengths during July 2007 using the Anglo-Australian Telescope and Infrared Imager and Spectrograph 2 (IRIS2). The data have been used to map the intensity and rotational temperature of the O₂(a¹Δ_g) airglow band at . The temperatures agree with those obtained in earlier IRIS2 observations and are significantly higher than expected from the Venus International Reference Atmosphere (VIRA) profile. We also report the detection of the corresponding ν=0-1O₂ airglow band at with a similar spatial distribution to the ν=0-0 band. Observations in the thermal window have been used to image surface topography using two different methods of cloud correction. We have also obtained images that can be used to study cloud motion.     

Non-constant rotation period of Comet C/2001 K5 (LINEAR)

The article presents results of CCD photometry in R-band of a dynamically new Comet C/2001 K5 (LINEAR), obtained at a heliocentric distance of about 5.6 AU, after the perihelion passage. Being so distant from the Sun, this comet was extremely active (Afρ close to 2000 cm), exhibiting quite well developed dust coma and tail. During the observations, general photometric behavior of the comet with heliocentric distance r was well described by the 2.5nlog(r) function with coefficient n=5. The radial profiles of the coma were found to be undulated, with mean slope of the dependence between cometary magnitude and 2.5log of aperture radius (at comet distance) equal to . The light curve of Comet LINEAR exhibited short-term variability which we attributed to cyclic changes of dust emission, induced by nucleus rotation. Model computations by some authors have revealed that active comets can change their spin status quite substantially even during a single orbital revolution. Thus, attempting to search for a rotation frequency, we have modified the classical PDM approach by including the spin acceleration term. Such DynamicalPDM (DPDM) method revealed the most reliable solution for the frequency f₀=0.019048±0.000013 h⁻¹ and its first time-derivative (index “zero” denotes reference to the mid time of the whole observing run), indicating a rapid spin-down of the nucleus. These parameters are equivalent to the rotation period of 52.499±0.036 h and its relative increment of 0.02729±0.00013. We present the most probable evolution of the rotation frequency of Comet LINEAR, based on the results of periodicity analysis and a simple, almost parameter independent, dynamical model of nucleus rotation. It is also shown that the DPDM may be an effective tool for determination of a nucleus radius, which provided us with the value of 1.53±0.25 km for Comet LINEAR.     

New near-aphelion light curves of Comet 2P/Encke

We present new, near-aphelion, time series of photometry of Comet 2P/Encke in Cousins-R band. With these light curves we find that the dominant, synodic rotational periodicity is either P₀=11.079±0.009 h or 2P₀=22.158±0.012 h. This is in contrast to data from the 1980s published by others that are consistent with 15.08- and 22.6-h periods. Those periods do not satisfy our phased light curves, and also the 1980s data are not easily reconciled with our periods. This could be due to P/Encke having non-principal axis rotation or due to a drift in the rotation period caused by outgassing torques. We observed the comet at five epochs: July, August, September, and October 2001, and September 2002, and the comet was at times intrinsically brighter than expected for a bare nucleus, due to an apparent contribution from an unresolved coma. Three-quarters of the data were obtained in the second and fifth epochs, and we analyzed these two time series using both the phase-dispersion minimization and “WindowCLEAN” techniques. At both epochs and with both techniques strong periodicities were found near frequencies and . By then using visual inspection of the phased light curves to corroborate these frequencies, and by using the data from the other three epochs to properly align light curve features, we were able to derive P₀ and 2P₀ as the only solutions that satisfy all our observations. The periodicity due to f₁ is clearly seen in our data, but we cannot tell from our data alone whether it is a manifestation of the nucleus's shape, non-principal axis rotation, or both.