comets are frozen reservoirs of primitive solar dust grains and ices . analysis of the composition and size distribution of cometary dust grains from infrared imaging and spectroscopic observations expedites an appraisal of the physical characteristics of the solid materials that constituted the primitive solar nebula @xcite .
the study of comets is an indirect probe of the origin of the constituents of the primitive solar system , their subsequent evolution into planetesimals , and their relationship to materials in other astrophysical environments @xcite .
although comets of all types have undergone some amount of post - formation processing , they remain the best preserved sources of material extant during our solar system s epoch of planet formation . in the current paradigm , the nearly isotropic comets ( i.e. , oort cloud and halley - type comets ) formed amongst the giant planets and were scattered into large orbits , in a spherically symmetric manner @xcite .
the ecliptic comets ( including jupiter - family comets ) originate from the kuiper - belt and scattered disk populations and likely formed _ in situ _ or in the transneptunian region @xcite .
comparisons between the nearly isotropic and ecliptic comets may reveal the differences in their post - formation processing or the structure and mineralogy of the proto - planetary disk .
both nearly isotropic and ecliptic comets have been exposed to bombardment by uv photons and cosmic rays , although to varying extents ( the jupiter - family comets have been exposed to @xmath1 times more uv and 100  kev solar cosmic rays than the oort cloud comets ) .
the ecliptic comets have suffered frequent collisions during their residence in the kuiper - belt and are likely to be fragments of larger kuiper - belt bodies @xcite .
the number of comets studied by mid - infrared spectroscopic methods necessary to determine their detailed mineralogies is increasing ( e.g. , see * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ?
* ; * ? ? ?
* ) and we may soon be able to compare comets to each other as groups , rather than individually .
we present _ spitzer space telescope _
@xcite images and spectra of comets 2p / encke , a 3.3  yr period ( @xmath2 ) , jupiter - family comet with a perihelion distance , @xmath3  au , known for an abundance of large dust particles @xcite and an extensive debris trail @xcite ; 67p / churyumov - gerasimenko ( 67p ) , @xmath4  yr , @xmath5  au , a jupiter - family comet and the primary mission target of the european space agency s _ rosetta _ spacecraft ; and c/2001  ht50 ( linear - neat ) ( ht50 ) , a long period , oort cloud comet , @xmath6  yr , @xmath7  au .    comet encke frequently approaches the earth on it s 3.3  yr orbit and is one of the most studied of all comets @xcite .
comet encke was also one of the first comets discovered to have a dust trail @xcite .
dust trails are composed of large ( @xmath8  mm ) , slow moving particles and are precursors to meteor streams [ encke is associated with the taurid meteor stream @xcite ] .
the comet also exhibits weak or non - existent 10   silicate emission @xcite .
the existence of a dust trail , association with a meteor stream , and the lack of a strong silicate feature has led investigators to conclude that encke s dust production is dominated by large particles .
we present _ spitzer space telescope
_ observations of comet encke in  [ spectra - obs - text ] , and derive the comet s dust coma mineralogy at 2.4  au in  [ encke - coma ] . the mineralogy of comet encke is discussed in
[ mineral - discuss ] .
we derive the temperature and effective size of the nucleus of comet encke in  [ encke - results ] and discuss the results in  [ eta - discuss ] .
comet 67p is the primary target of the _ rosetta _ mission .
the spacecraft is designed to intercept and orbit the comet at @xmath9  au ( pre - perihelion ) to study the development of coma activity as the comet approaches the sun .
information on the comet s dust environment is crucial to mission planning , which motivated our _ spitzer _ observations of the comet at 5.0  au ( post - perihelion ) .
the observation is presented in
[ spectra - obs - text ] and the results presented in  [ cg - results ] .
comet ht50 is an oort cloud comet with an orbital period that suggests it has orbited the sun many times ( @xmath6  yr ) .
it was discovered to be cometary at the large heliocentric distance of 7.5  au @xcite .
we observed comet ht50 twice with _
both observations are presented in
[ spectra - obs - text ] and we derive dust mineralogies at both epochs in  [ ht50-results ] .
we discuss ht50 s mineral content and compare it to other oort cloud comets in  [ mineral - discuss ] .
spectra of comets encke , 67p , and ht50 were obtained with the infrared spectrograph ( irs ) instrument @xcite on _
spitzer_. all comets were observed in both low- ( @xmath10128 ) and high- ( @xmath11 ) resolution modes , although not at all wavelengths . the slit widths , aperture extraction sizes ( described below ) , and slit orientations are listed in table  [ apertures ] .
note that the slit parameters force us to sample different areas of each comet s coma ( even though each slit extraction is centered on the nucleus ) and remain as a source of unknown error in our resultant spectra .
a summary of our observations is presented in table  [ spec - log ] .
there were two observations of comet encke , one from our _ spitzer _ guaranteed time observation ( gto ) program , program identification ( pid )  119 , and one from pid  210 , one of comet 67p from our general observer proposal ( pid  2316 ) , and two temporally distinct observations of ht50 from our gto program ( pid  131 ) .
post - pipeline spectral reductions were performed on irs pipeline s12.0.2 basic calibrated data ( bcd ) frames .
our reduction method proceeded as follows : 1 ) the bcd frames were two - dimensionally background subtracted ( when possible , as described below ) and fatally bad pixels were removed by nearest - neighbor interpolation .
2 ) spectral extraction and initial calibration were performed by the _ spitzer
_ irs custom extraction ( spice ) tool , available from the ssc .
if a two - dimensional background subtraction was not possible , we subtracted a background model from the extracted spectra at this step .
3 ) we removed the nucleus contribution , if known , from the extracted spectra .
nuclei were modeled with the near - earth asteroid thermal model @xcite as described in  [ encke - results ] .
4 ) we corrected the combined , nucleus subtracted spectra with our derived extended source calibration ( described below ) .
5 ) we scaled each module to produce one continuous spectrum . the scaling removed module - to - module mis - calibrations and corrected for differences in module slit orientations with respect to the comet comae .
6 ) finally , we scaled the entire spectrum to agree with aperture photometry derived from _ spitzer _ imaging (  [ imaging ] ) , if possible .
all spectra were extracted with a constant width aperture to produce a final spectrum representative of a given aperture size around the comet ( table  [ apertures ] ) .
our sl and ll apertures were 17.5 and 51.5 on the sky ( @xmath12 of the slit ) , i.e. , they encompass the first diffraction ring of a point source at the longest wavelength of each module .
the default point source optimized aperture used by the irs calibration pipeline has an aperture size that varies with wavelength , i.e. , the default aperture width in arc seconds is proportional to wavelength .
such an aperture encompasses a diffraction limited point source ( but a varying amount of sky ) at every extracted wavelength .
therefore , the point source optimized aperture is not optimized for extended source extractions .
we compared spectra of the irs calibration star hd  173511 extracted with our constant width apertures to spectra extracted with the default point source tuned apertures .
the differences in the point - to - point photometry ( aperture averaged at a given wavelength ) between the two spectra were @xmath13% .
we note the irs pipeline uses the entire available slit by default in the high - resolution modules .
an automated spectral calibration of extended sources was not available from the ssc , yet the extent of the comae of encke and ht50 ( described below ) requires us to calibrate these spectra for extended sources .
currently , point sources are used for irs calibration targets .
we estimate the narrow entrance slits of the irs modules block up to 40% of the flux of a point source at the longest wavelengths .
although this does not introduce additional absolute photometric uncertainty into the final data products when calibrating point sources ( as this same fraction is always removed from all point sources ) , the slit - loss does affect the photometric precision for sources larger than the point spread function of the instrument . the emission from comets
is comprised of emission from the nucleus ( generally spatially unresolved ) plus a contribution from the coma , which is a region of radially and azimuthally varying surface brightness . to correctly flux calibrate the comet spectra ,
the nucleus spectrum must be subtracted if the nucleus flux is a significant fraction of the total emission in the slit ( step  3 , above ) .
inspection of the irs peak - up images is used to confirm the validity of this assumption for any individual comet ( cf ,  [ imaging ] ) .
comets encke and 67p required nucleus subtraction , ht50 did not .
the encke nucleus is discussed in  [ encke - results ] and the 67p nucleus is discussed in  [ cg - results ] . in all cases , we assume the coma surface brightness is a slowly varying function over the slit width ( again confirmed by the peak - up images ) .
next , we generate an image of a point source at each extracted wavelength for an obscured primary with _ spitzer _
s optical parameters @xcite .
a slit mask is applied to the image and the flux in the slit mask ( @xmath14 ) is compared to the total flux of the simulated point source ( @xmath15 ) .
the slit - loss correction factor is @xmath16 ( fig .  [ slcf ] ) .
the coefficients of third - order polynomial fits to the slit - loss correction factors ( @xmath17 in units of ) for the four irs modules are given in table  [ slcf - table ] .
this slit - loss correction factor is multiplicatively applied at each wavelength to the nuclear flux subtracted spectra .
the five observations ( table  [ spec - log ] ) were designed with various observing strategies .
the first observation of encke was designed to obtain high - resolution spectra of the comet from 1038   and low - resolution spectra from 514  .
it was executed on 2004 june 25 at 05:28  ut when encke was at a heliocentric distance , @xmath18 , of 2.573  au , a _
spitzer_-comet distance , @xmath19 , of 1.985  au , and a phase angle , @xmath20 , of 21 .
the observation used an irs spectral map astronomical observing request ( aor ) .
we describe spectral map aors with @xmath21 notation , where @xmath22 is the number of steps parallel to the slit s long ( spatial ) direction , and @xmath23 is the number of steps perpendicular to the long direction .
the first comet encke observation used a @xmath24 map of 60 steps perpendicular to the slit direction for all three modules .
the two 60 steps off - source provide background measurements for the on - source position .
a red peak - up image was used to acquire comet encke and center the irs slits .
the second observation of encke only used low - resolution modules to cover the entire 538   spectral region and was executed on 2004 june 25 at 08:27  ut ( @xmath25  au , @xmath26  au , @xmath27 ) .
a blue peak - up image was used to acquire comet encke .
the aor for the 514   low - resolution module ( sl ) executed a @xmath28 spectral map with 20 steps parallel to the slit and the 1438   low - resolution module ( ll ) executed a @xmath29 map with 40 steps parallel to the slit .
the parallel motion was small enough such that the comet was located in the slit for each map position .
each low - resolution module consist of two slits , one of which is centered on the science target and the other offset @xmath30 for sl and @xmath31 for ll from the comet ( fig .
[ bonus - background ] ) .
we use the extra slit position to background subtract our spectra . for example , when the comet is in the 1st order sl slit ( sl1 , 714  ) , sky background is measured in the 2nd order sl slit ( sl2 , 57  ) .
the spacecraft is nodded to place the comet in the sl2 slit , providing a measurement of the sky background in sl1 ( see fig .
[ bonus - background ] ) . for the second observation of comet encke
, the off - source order is used as a background measurement and two - dimensionally subtracted from the science order .
the observation of comet 67p used an irs spectral stare aor in both long wavelength modules and resulted in a spectrum spanning 1438  .
this aor was executed on 2004 july 15 at 08:13  ut ( @xmath32  au , @xmath33  au , @xmath34 ) .
the red peak - up array shows a point source at the location of the comet . to remove the background in ll , we subtracted nod pairs and extracted the resulting spectra .
the high - resolution module did not detect the comet and will no longer be discussed .
two observations of comet ht50 were executed .
the first observation , obtained on 2003 december 17 at 15:39  ut ( @xmath35  au , @xmath36  au , @xmath37 ) , was an irs spectral map aor of size @xmath24 with 7.2 steps ( @xmath38 ) sl , a @xmath39 map of @xmath40 ( @xmath41 ) for the short wavelength , high - resolution module ( sh ) , and a @xmath42 map of @xmath43 ( @xmath44 ) for the long wavelength , high - resolution module ( lh ) .
we extracted spectra from the map positions with the greatest amount of signal .
the second observation was obtained on 2004 july 18 at 10:36  ut ( @xmath45  au , @xmath46  au , @xmath47 ) with an irs spectral stare aor for both high - resolution modules and the first order ( 714  ) of the sl module .
there were no measurements of the background for the sh and lh modules .
instead , the background flux density was estimated using the background model in the _ spitzer planning observations tool _ ( spot ) provided by the _ spitzer _
science center ( ssc ) .
the model background was subtracted after spectral extraction .
restricting the second epoch sl observation of ht50 to a one order spectral stare does not permit us to use the same two - dimensional background subtraction techniques employed for comet encke .
furthermore , the extent of the coma in ht50 does not permit us to directly nod subtract the images to remove the background as described above for comet 67p ( 67p was a point source within the slit ) .
however , the sl module is long enough in the spatial dimension ( @xmath48 long ) such that the coma of ht50 did not entirely fill the slit . to determine the background contribution in sl , a gaussian plus linear background term
was fitted to the spatial dimension for every @xmath49 .
the linear background term was subtracted from the extracted spectrum .
all comets were imaged at the time of spectra acquisition by the peak - up mode of the irs .
comet encke was also observed with the infrared array camera ( irac ) in all band - passes during _
spitzer s _ in - orbit checkout and science verification phase ( w.  t. reach et al .  2006 , in preparation ) on 2003 november 11 at 17:33  ut ( @xmath50  au , @xmath51  au , @xmath52 ) .
the 8.0   images saturated on the nucleus , but the remaining short exposure ( 0.6  s ) 3.65.8   irac images did not saturate and were used in our analysis ( fig .
[ nov - encke - images ] ) .
encke was again imaged with irac during normal science operations on 2004 june 29 at 22:05  ut ( @xmath53  au , @xmath54  au , @xmath55 ) and with the multiband imaging photometer for _ spitzer _ ( mips ) at 24 and 70   on 2004 june 23 at 05:04  ut ( @xmath56  au , @xmath57  au , @xmath58 ; fig .
[ jun - encke - images ] ) .
detailed information regarding the performance and use of irac and mips is provided by @xcite and @xcite , respectively .
the imaging observations are summarized in table  [ image - log ] .
the 2003 november data were processed with irac pipeline s11.0.2 and the 2004 june data were processed with irac and mips pipelines s11.4.0 .
the 2004 june comet encke images were mosaicked in the rest frame of the comet with the mopex software @xcite at the native irac and mips scales ( 1.22  pixel@xmath59 for irac , 2.5  pixel@xmath59 for mips 24 , and 5.1  pixel@xmath59 for mips 70 ) and are presented in fig .
[ jun - encke - images ] .
the nucleus clearly dominates the 5.8 , 8 , and 24   images based on the tenuous coma emission and distinct diffraction rings .
careful inspection also reveals that the nucleus is detected at 3.6 , 4.5 , and 70  .
comet encke s debris trail is seen as a diagonal line of emission extending across the 8 and 24   images @xcite .
also evident are two streaks of material that comprise the remains of encke s activity during this perihelion passage .
the irs peak - up images are presented in fig .
[ peakups ] .
all images use the red peak - up array ( @xmath60  ) except the 2004 june 25 08:27  ut encke peak - up , which uses the blue peak - up array ( @xmath61  ) .
the images are not flux calibrated but are useful in describing the comet morphologies in the following sections .
the spectrum of the encke nucleus must be estimated and subtracted from the raw comet spectrum before extended source calibration , as discussed in  [ spectra - obs - text ] .
the light from a comet nucleus is a combination of reflected light and thermal emission at the shortest wavelengths ( @xmath62  ) and solely comprised of thermal emission at longer wavelengths ( @xmath63  ) .
we approximate the reflected solar spectrum with a 5770  k black body . the thermal emission for comet nuclei is typically modeled with asteroid thermal models @xcite , which have consistently modeled the thermal emission from those comets visited by spacecraft , for example 9p / tempel  1 @xcite and 81p / wild  2 @xcite .
originally , the standard thermal model ( stm ) was successful in modeling the thermal emission from asteroids at small ( @xmath64 ) , phase angles @xcite .
@xcite proposed the near - earth asteroid thermal model ( neatm ) that extends the stm to higher phase angles and introduces a variable infrared ( ir ) beaming parameter , @xmath65 , rather than holding it constant at 0.756 ; a value which reproduces the diameters of asteroids ceres and pallas @xcite .
the ir beaming parameter in models either raises or lowers the sub - solar temperature , @xmath66 , as @xmath67^{1/4}\   \mbox{(k)},\ ] ] where @xmath68 is the bolometric bond albedo , @xmath69 is the incident solar flux , @xmath70 is the infrared emissivity , and @xmath71 is the stefan - boltzmann constant . in our treatment ,
a scaled 5770  k black body spectrum was added to account for any reflected light at the shortest wavelengths ( @xmath72  ) , @xmath73 where @xmath74 is freely variable , unit less scale parameter , @xmath75 is the planck function in w  @xmath76  @xmath59 , @xmath77 is the radius of the comet in km , @xmath78 is the radius of the sun in km , @xmath79 is the comet - observer distance in km , and @xmath18 is the comet - sun distance in km .    the near diffraction limited _ spitzer _ images of comet encke can be used to estimate the size of the comet nucleus if the coma contribution to the observed surface brightness is negligible or can be removed . in the two sets of irac images and the mips 70   image ( fig .
[ nov - encke - images ] and [ jun - encke - images ] ) , the coma was faint or undetected and the nucleus was sufficiently measured through aperture photometry .
aperture corrections from the irac data handbook were applied to the irac photometry to account for the chosen aperture radii and background annuli [ irac is calibrated with a 10  pixel radius aperture and a 1020  pixel annulus , @xcite ] .
aperture corrections applied to the mips photometry were derived from the 500  k mips point source response functions provided by the ssc .
all apertures were 510  pixels in radius except for the 2004 june irac 3.6 and 4.5   apertures , which were set to 2 pixels in radius to avoid neighboring stars .
the mips 24   image was fitted with a point source derived from stars in the field and the irac 8   image was saturated by the nucleus .
the measured fluxes are provided in table  [ encke - nuke ] .
the current estimate of the flux calibration errors is @xmath80% for irac @xcite and @xmath81% and @xmath82% for mips 24 and 70 @xcite .
the errors in table  [ encke - nuke ] were produced from the quadrature addition of the formal photometric errors and the flux calibration errors .
the irac color corrections to the photometry are very sensitive to the underlying spectral energy distribution ( sed ) when dealing with black bodies in the 200400  k range , e.g. , the 3.6   color correction varies from 651% for @xmath83200  k @xcite .
additionally , the 3.6   bandpass likely includes an equal combination of thermal emission and reflected solar radiation from the nucleus ( @xmath84% from each ) [ e.g. , the sed of comet c/1995  o1 ( hale - bopp)@xcite ] .
we fit the sed with an initial neatm model to approximate the true sed and calculate color corrections as prescribed in the irac data handbook @xcite .
the color corrected seds were modeled with a non - linear reduced chi - squared ( @xmath85 ) fitting routine using the neatm to derive the the effective radius of comet encke s nucleus .
the 2004 june epoch had the best spectral coverage , ranging from 3.670  .
three parameters affect the color - temperature of the nucleus , @xmath70 , @xmath68 , and @xmath65 .
for encke s nucleus , we varied @xmath65 as a free parameter , adopted a geometric albedo , @xmath86 , which is related to the bond albedo @xcite , following the discussion by @xcite , and set @xmath87 .
the best - fit effective radius was @xmath88  km ( @xmath89 , @xmath90 ) .
the entire list of fitted parameters is presented in table  [ encke - fit - table ] and the model sed is presented in fig .
[ encke - fit - fig1 ] .
the 2003 november sed was limited in spectral coverage to 3.65.7   and the number of photometry points ( 3 ) permitted us to only fit two parameters at a time .
therefore , we fixed the value of @xmath65 to the best - fit value derived from the 2004 june epoch and fitted @xmath77 and the scale parameter of the reflected light ( @xmath74 ) .
the result was an effective nucleus radius of @xmath91  km with @xmath92 .
the error was derived by varying @xmath65 over the 2004 june range of @xmath93 .
the fit is presented in table  [ encke - fit - table ] and the sed in fig .
[ encke - fit - fig2 ] .
the 2003 november and 2004 june best - fit effective nucleus radii are different by 0.6  km .
if this is evidence for a non - spherical nucleus then the axial ratio is at least @xmath94 .
this value is consistent with the @xcite axial ratio of @xmath95 derived from the visual light curve of the nucleus .
our _ spitzer _ observations were separated by seven months , however , we can not reliably re - phase our data as the rotation period of encke s nucleus is not precisely known @xcite .
previous investigations of comet encke s nucleus with thermal observations provide effective radii consistent with our derived value of @xmath96  km . @xcite
derive a value of @xmath97  km , @xcite posit @xmath986.4  km , and @xcite conclude @xmath99  km .
the 0.6  km effective radius disagreement in our _ spitzer _ estimates can be minimized by varying @xmath65 in the 2003 november fits and forcing the radius to be constant at 2.34  km .
this method results in @xmath100 and @xmath101 .
the new fit is presented in table  [ encke - fit - table ] and in fig .
[ encke - fit - fig2 ] .
the low @xmath85 indicates an under - constrained fit , but the fit appears improved as the scattered light contribution in the model sed is similar to the 2004 june best - fit ( fig .
[ encke - fit - fig1 ] ) .
recent work on near - earth asteroids suggests @xmath65 may have a phase angle dependence @xcite .
we discuss the radius of comet encke and the role of the ir beaming parameter further in  [ eta - discuss ] .      to properly model and interpret the coma dust emission in comet encke requires a useful photometric calibration of the narrow slit spectra .
first , we must take into account the varying slit widths amongst the irs modules ( ranging from 3.7 to 11.1 ) , normalizing to an aperture size of radius 4.5 .
our adopted aperture size is roughly equivalent in area to the module with the smallest slit and extraction aperture ( sl ) , which is @xmath102 as described above (  [ spectra - obs - text ] ) . note that the aperture normalization can not account for the different regions sampled by the different slit sizes and orientations (  [ spectra - obs - text ] ) .
small pointing offsets will cause mis - calibrations of the spectral flux , especially with the sl module , which has a slit width of 3.7 .
this effect is evident when we consider that the inner peak of the _ spitzer _ point spread function ( psf ) has a full - width at half - maximum of about 2.4 at 10  .
the coma flux can be estimated from the mips 24   image where we have subtracted the nucleus .
the resultant flux can then be used to photometrically calibrate all spectra through module - to - module wavelength overlaps .
we measure the coma flux using two methods .
aperture photometry ( method  1 ) on the residual coma of the point source subtracted image yields a coma flux density of @xmath103  w  @xmath76  @xmath59 in a 4.5 radius ( 1.8 pixels ) beam .
the error includes the photometric uncertainty of the nucleus subtraction assuming it is distributed evenly over the psf . alternatively
( method  2 ) , we fit the azimuthally averaged aperture flux density profile from 312  pixels in radius with a power - law , @xmath104 , where @xmath105 is a scale factor to account for the total brightness , @xmath106 is the aperture radius in pixels , and @xmath107 is the logarithmic slope .
our fit yields the values @xmath108  w  @xmath76  @xmath59 and @xmath109 ( @xmath110 ) . the best - fit slope is different than the nominal `` steady - state '' profile ( @xmath111 , * ? ? ?
* ) because of the highly structured dust morphology of comet encke at this epoch and these aperture radii ( see fig .
[ jun - encke - images ] ) . computing the flux enclosed in a 4.5 radius aperture with this surface brightness profile
yields a coma flux density of @xmath112  w  @xmath76  @xmath59 , again including the psf fit error .
the two methods are in good agreement .
we elect to adopt method  2 for scaling the irs spectra as it is less likely influenced by the pixel - to - pixel errors of the psf fitting at the nucleus removal step .
a cut across the psf subtracted image and the aperture photometry profile is presented in figs .
[ encke - psf1 ] and [ encke - psf2 ] .
the best - fit neatm nucleus sed ( @xmath113  km ) was subtracted from the extracted comet encke spectra and then scaled to match the mips coma photometry .
first , we scaled the separate modules by area to the fiducial 4.5 aperture .
next , we scaled the spectra to match the spectra overlap regions . finally , we integrated the spectrum under the mips 24   bandpass to determine the final scale factors ( table  [ scales ] ) .
the sl spectrum from the 2004 june 25 05:28  ut epoch was dominated by the errors after nucleus subtraction ( i.e. , the coma was at low signal to noise ) , hence the large scale factor of 1.53 .
the final , combined 2004 june 25 05:28  ut epoch spectrum is presented in fig .
[ encke - spec ] .
we repeat this procedure for the 2004 june 25 08:27  ut epoch , high - resolution spectrum of encke .
the neatm nucleus predicts a flux larger than observed in the sl spectrum data points by more than one standard deviation ( overall , @xmath114 ) .
an improperly centered ( within the slit ) comet could cause such a photometric error , or the nucleus may be presenting a larger cross section to the telescope than the @xmath115  km model provides .
we can not estimate the coma flux in this module and therefore drop sl from the analysis .
the sh and lh scale factors ( table  [ scales ] ) account for the asymmetric dust morphology of comet encke and the shape and orientation of the slits . the final spectrum is presented in fig .  [ encke - spec ] .
the comet encke spectra are modeled using a thermal grain model developed by @xcite which self - consistently calculates the temperature and thermal emission from a cometary dust mixture using laboratory optical constants .
the minerals include amorphous carbon , amorphous pyroxene , amorphous olivine , and crystalline olivine .
the amorphous carbon is used to represent warm , featureless continuum from deeply absorbing grains .
the amorphous pyroxene and amorphous olivine grains are a 50 - 50 mixture of magnesium and iron which have an increased temperature as compared to pure - magnesium grains .
the increased temperature was required to adequately model the silicate features in the thermal spectrum of comet hale - bopp @xcite .
the crystalline olivine grains are magnesium - pure ( forsterite ) .
the thermal model invokes a hanner grain size distribution ( hgsd ; * ? ? ? * ) to calculate the emission from a population of grains , @xmath116 with radii , @xmath117 , varying from 0.1100  .
the grain population is described by a modified power law , @xmath118 where @xmath117 is the grain radius , @xmath119 is the minimum grain radius ( assumed to be 0.1  ) , @xmath120 is the slope of the distribution at large @xmath117 , and @xmath121 is related to the radius of the peak of the grain size distribution , @xmath122 , by @xmath123 this choice of grain size distribution is known to reproduce the observed seds of many comets in the 3.520   wavelength range .
model grains can be fractally porous , with a density described by @xmath124 where @xmath125 is the bulk density and @xmath126 is the fractal dimension of the dust ( solid spheres have @xmath126 equal to 3 , porous spheres have @xmath127 ) . to keep coma fitting tractable , we chose discrete values of @xmath122 ( 0.1   steps ) and @xmath126 [ one of 3.0 , 2.857 , 2.727 , 2.609 , and 2.5 ( see * ? ? ? * ) ] , and all minerals in a given observation are assumed to have the same grain size distribution . after a best - fit model is derived , we attempt to estimate the error in peak grain size by exploring the curvature of @xmath128 space with respect to @xmath122 . for the encke and ht50 best - fits below ,
the derived @xmath122 errors are @xmath129  .
given the small formal errors in @xmath122 , our choice of 0.1   steps ensures our derived peak grain sizes are within 0.05   from the best values .
the coma of comet encke was fitted ( @xmath1300.9 ) with a porous ( @xmath131 ) , amorphous carbon mineralogy and a peak grain size of 0.4   ( @xmath132 , @xmath133 ) .
the model spectra ( solid curve ) are shown in fig .
[ encke - spec ] .
the agreement in best - fits between the low - resolution and high - resolution data suggest our reduction methods are consistent , regardless of the irs module . the number of peak grains and upper - limits to olivine and pyroxene minerals are presented in table  [ encke - minerals ] .
our best - fit model constrains the sub - micron silicate fraction to @xmath0% by mass .
the small peak grain size suggests that comet encke has not exhausted its reservoir of small particles , but the large particle slope parameter , @xmath134 , does not preclude the importance of large particles . for our derived psd ,
the ratio of the total mass of particles with @xmath135   to the total mass of particles with @xmath136   is 0.076 , i.e. , the sub - micron particles are a minor component of the total coma mass .
discussion of comet encke s mineralogy is presented in
[ mineral - discuss ] .
the peak up image of 67p ( fig .  [ peakups ] ) showed a point source with a full - width at half - maximum of 3 pixels , or 18,000  km at the distance of the comet . although the dust production of jupiter - family comets at 5  au is generally assumed to be minimal or non - existent , a point source does not indicate a bare nucleus .
a point source could also be a combination of nucleus and 1 ) a coma from recent ( hours to weeks ) activity , 2 ) a coma of slowly moving micrometer sized or larger particles ejected this perihelion passage , or 3 ) very large particles ( @xmath137 ) entrained in 67p s debris trail . assuming only a bare nucleus , we fitted the spectrum with the neatm ( @xmath138 , @xmath139 , and @xmath87 ) .
the fit yielded a nucleus radius of @xmath140  km ( @xmath141 ) and is presented along with the spectrum in fig .
[ cg - nuke ] .
our derived best - fit effective radius for this comet is inconsistent with recent estimates of 67p s nucleus size @xcite .
for example , @xcite derive a value of @xmath142  km using the _ hubble space telescope _ , and @xcite estimate a value of @xmath143  km using _ spitzer_/mips .
it is clear our _ spitzer _ spectrum of comet 67p is not that of a bare nucleus .
likely , some amount of dust still enshrouds the nucleus at 5  au .
the nucleus subtracted spectrum is not of high enough quality or spectral range for a detailed coma fit .
when the @xcite neatm fit , calculated for the epoch of the spectrum , is subtracted , the resulting spectrum constrains the dust contribution .
the observed coma flux is @xmath144  w  @xmath76  @xmath59 at 27.9   ( weighted average from 2135  )
.    one component of 67p s dust emission at this large heliocentric distance is the comet s debris trail . to determine the fraction of the flux originating from the trail , we use the trail parameters determined by @xcite from _ iras _ measurements . in the _ spitzer _ observations
, the trail contributes @xmath145  w  @xmath76  @xmath59 , or @xmath146% , to the irs spectrum .
we can further refine this estimate by taking into account the effect of the orbital motion at 5  au in an eccentric orbit . here
, particles in similar orbits will move slower , and hence be closer together , than at @xmath147  au ( the epoch of the _ iras _ observations ) . to estimate the optical depth enhancement caused by this effect
, we differentiate the true anomaly with respect to heliocentric distance @xmath148 where @xmath149 is the true anomaly at the time of observation ( @xmath150 at the iras epoch , @xmath151 at the _ spitzer _ epoch ) and @xmath152 is the eccentricity of the orbit of comet 67p .
the ratio of optical depths from the _ iras _ epoch to the _ spitzer _ epoch is 1.25 , increasing the trail contribution to 3% . the majority of the dust emission detected by _ spitzer
_ likely arises from 1 ) recently ejected dust ( age of order hours to weeks ) , 2 ) large , slowly moving particles from the 2003 perihelion passage [ evidenced by the shallow post - perihelion @xmath18 dependence of dust from optical observations @xcite ] , or 3 ) some combination of the two .
the sed of 67p increases in flux density from 1420  , then flattens ( in slope ) beyond 20   ( fig .
[ cg - nuke ] ) . in principle
, the color - temperature of the coma can constrain the mineral composition of the emitting dust .
we attempted to account for the spectral shape by fitting the spectrum with our thermal model , constrained to one mineral ( amorphous olivine or amorphous carbon ) and constrained to specific peak grain sizes ( @xmath153 , 5 , 15 , or 30  ) .
for similarly sized particles , amorphous carbon dust has a higher color - temperature than amorphous olivine , and thus produces a poorer fit to the data . also , smaller sized particles are warmer and produce poorer fits .
the best - fit model consisted of large ( @xmath154  ) amorphous olivine grains . for models of amorphous olivine dust
, @xmath85 ranged from 0.57 ( @xmath155  , maximum grain size of 30  ) to 0.73 ( @xmath156  , maximum grain size of 100  ) . the increased maximum grain size in the @xmath153   model balances the warmer , @xmath157   particles so that the resultant model spectrum fits the observed color - temperature .
overall , the low signal - to - noise ratio of the sed and the small range of grain temperatures for 1100   grains at 5.0  au ( @xmath158150  k ) prohibits a detailed discussion of 67p s dust coma mineralogy .
inspection of the comet ht50 peak - up images ( fig .
[ peakups ] ) reveals a bright coma at both epochs ( @xmath159 and 4.6  au ) .
we inspected the peak - up images for the signature of a point source .
the peak - up image at @xmath159  au is saturated over the center five pixels and does not yield any information on the nucleus .
the profile of the peak - up image at @xmath160  au does not exhibit a point source within the coma profile .
we therefore proceed with the assumption that the nucleus contribution at each epoch is negligible and do not subtract a model nucleus from the ht50 spectra .
all irs spectra were scaled to the nominal 4.5 aperture to mitigate module photometric mis - matches .
similar to the comet encke spectra , the extracted ht50 spectra must be further scaled to account for varying slit orientations , coma asymmetries ( fig .
[ peakups ] ) , and the overall coma profile ( table  [ scales ] ) .
the spectra and our best - fit thermal models are presented in figs .
[ ht50 - 3au ] and [ ht50 - 4au ] , and in table  [ ht50-minerals ] .
the best - fit thermal models have a large peak grain size of 1.2   and a mixed amorphous carbon and silicate mineralogy ( @xmath161 , @xmath162 , @xmath163 at 3.2  au and @xmath134 , @xmath164 , @xmath165 at 4.6  au ) .
the total number of amorphous carbon grains from 0.110   varies from @xmath166 ( @xmath159  au ) to @xmath167 ( @xmath160  au ) .
a number of factors contribute to the number of observed grains , including the changing grain size distribution , the aperture size at the distance of the comet , the shape of the coma , and the dependence of coma activity on heliocentric distance ( including phenomena such as jet activity and gas / dust outbursts ) .
we can not account for the heliocentric dependence of the dust production without a more rigorous temporal sampling of the coma at these epochs .
our best - fit model indicates a weak detection of crystalline olivine in ht50 at 3.2  au with a signal - to - noise ratio of 4 .
figure  [ ht50-cryst ] presents a closer analysis of the crystalline fit . here ,
the spectrum was smoothed by a seven - point statistically weighted , moving average to increase the signal - to - noise of the spectrum .
the spectrum was then normalized by the best - fit model , excluding the crystalline component .
excess emission ( above unity in the normalized spectrum ) is potentially due to crystalline olivine dust emission ( represented in the figure by the solid line ) .
the chi - squared fitting of the 1-sigma correlated errors indicates that a better fit ( @xmath168 ) is obtained with crystalline olivine present rather than absent ( to 4-@xmath71 or @xmath169 confidence ) .
inspection of fig .
[ ht50-cryst ] suggests the detection of crystals is driven by the shape of the spectrum at 2224  .
altogether , the _ spitzer _
observations suggest the presence of crystalline olivine in comet ht50 .
the mineralogy of comet ht50 is discussed and compared to other oort cloud comets in  [ mineral - discuss ] .
recent observations of near - earth objects ( neo ) suggest a possible dependence of @xmath65 on phase angle @xcite . @xcite
derived an empirical @xmath65-phase angle correlation from observations of neos at phase angles @xmath17060 .
@xcite added more data to the discussion and derived the trend @xmath171 .
however , extensive data to assess whether an @xmath65-phase angle correlation exists for _ single _ objects does not exist .
only one asteroid ( 2002  ny40 ) has been observed in detail at disparate phase angles and no @xmath65-phase dependence was found @xcite .
comet encke is a near - earth object , and observations of its nucleus can be used to assess the validity of the neatm and the @xmath65-phase angle correlation .    to examine the potential correlation of the ir beaming parameter with phase angle
, we will assume the nucleus has a radius equal to the radius derived from our 2004 june sed , @xmath113  km .
the ir beaming parameter is also best constrained by the 2004 june sed to a value of @xmath172 . with these parameters ,
the neatm predicts a flux of 3.9  jy for the @xcite observation on 1997 july 19 at a phase angle of 44 .
however , @xcite observed a flux density of @xmath173  jy .
the 2004 june derived radius and ir beaming parameter do not account for the observed flux in 1997 july , just as they did not account for the 2003 november sed (  [ encke - results ] ) .
the @xcite 1997 july 8.5 , 10.7 , and 11.6   observations were fitted with the 2004 june radius and we derived a best - fit @xmath65 value of @xmath174 .
the 2003 november and 1997 july data sets are limited , but suggest a larger @xmath65 value is required for higher phase angles .
an @xmath65-phase angle dependence can be a direct result of certain physical properties of the nucleus surface .
there are two possible mechanisms that may produce the a change of @xmath65 with phase angle ( see * ? ? ?
* ; * ? ? ?
either , the surface roughness of the nucleus is exacerbated at high phase angles ( i.e. , the observer views an increasing amount of shadowing ) , requiring a lower temperature to model the nucleus ; or alternatively , the night hemisphere of the nucleus does not completely cool to the background temperature and contributes appreciable flux at high phase angles , requiring a lower temperature to model the nucleus .
we consider these possibilities in order .
in the former case , surface roughness causes severe shadowing across the diurnal hemisphere . at high phase angles
, an observer could be viewing the shadowed side of a scarp or other surface topology as observed on the surfaces of comets 1p / halley , tempel  1 , 19p / borrelly , and 81p / wild  2 @xcite .
this geometry presents an overall cooler surface when compared to a smooth sphere at the same phase angle ( see * ? ? ?
* fig .  5 ) .
this requires a lower temperature in the neatm and , therefore , a higher ir beaming parameter . at
low phase angles the contrary is true and an observer may be viewing more surfaces normal to the sun direction , thus requiring a hotter surface and lower ir beaming parameter in the neatm , consistent with the @xcite trend .
the temperature map of the tempel  1 nucleus by @xcite clearly shows how local surface temperature relies on the topology of the nucleus .
surface roughness on smaller scales ( e.g. , in the regolith of small bodies ) is also important in thermal observations @xcite and necessary to explain polarization measurements of atmophereless bodies @xcite .
if surface elements on the night hemisphere contribute significant flux to an observer , they can lower the observed color - temperature . in the formalism of the neatm
, all radiation is assumed to arise from the sun - lit hemisphere and therefore a lower observed color - temperature will raise the derived @xmath65 value . the temperature of the night hemisphere will depend upon the angular rotation rate of the nucleus and the ability of the surface to hold heat , i.e. , its thermal inertia . following @xcite and @xcite we can test whether or not the encke nucleus could be considered a fast or slow rotator .
a slow rotator model uses the same assumptions as the stm and neatm , i.e. , the sub - solar point is the hottest point and the night side emits no light .
a fast rotator model assumes the cooling time scale is longer than the rotation period and , therefore , the object will be isothermal with respect to latitude .
the unit less parameter @xmath175 is used to determine the applicability of the two models , @xmath176 where @xmath177 is the thermal inertia in j  k@xmath59  m@xmath178  s@xmath179 and @xmath180 is the angular rotation rate of the object in s@xmath59 .
slow rotators have @xmath181 and fast rotators have @xmath182 .
we chose 10 and 320  j  k@xmath59  m@xmath178  s@xmath179 as the thermal inertia extrema for short period comets , as suggested by the observations of the sub - solar temperature of tempel  1 @xcite and ground - based observations of ( 3200 )  phaethon @xcite .
@xcite derived two possible periods for the encke nucleus , 11.1 or 22.2  hr .
the 11.1  hr rotation period has been verified by radar observations @xcite . to determine the sub - solar temperature we use our best neatm fits , @xmath183 at 1.1  au and @xmath184 at 2.6  au , which produces @xmath185  k at 1.1  au and 269  k at 2.6  au .
together , @xmath175 ranges 0.020.6 at 1.1  au and 0.051.6 at 2.6  au . with these new measurements we conclude , as did @xcite , that at best we can use a slow rotator model ( if @xmath186  j  k@xmath59  m@xmath178  s@xmath179 ) and at worst the encke nucleus is in an intermediate state between slow and fast rotation ( if @xmath187  j  k@xmath59  m@xmath178  s@xmath179 ) .
additionally , there is no clear transition to a fast rotator as encke recedes from the sun .
we also note the recent mid - infrared measurements of asteroid 2002  ny40 by @xcite at phase angles of 22 and 59 show no @xmath65-phase angle dependence .
although the majority of the observations suggest the dependence @xcite , the relation may not hold true for any particular object .
the case of comet encke is not resolved .
the observed seds are fit by an @xmath65 parameter that varies with phase angle but may also be fit by changing the observed effective radius within the constraints of the derived nucleus shape .
an sed of the encke nucleus at high phase angle that constrains both effective radius and color temperature is required to support either possibility .
the mineralogies of comet ht50 s dust coma at both epochs are consistent with each other given the errors and upper - limits ( table  [ rel - mass ] ) .
also presented in table  [ rel - mass ] are the mineralogies of comets hale - bopp and tempel  1 , pre- and post-_deep impact _ encounter @xcite .
the low mass ratio between the silicate minerals and the deeply absorbing grains ( represented by amorphous carbon ) appears to be intermediate between the aging jupiter - family comet encke and the presumably pristine materials ejected from comet hale - bopp and comet tempel  1 ( post-_deep impact _ ) .
does comet ht50 have a processed surface similar to what may be on the surfaces of jupiter - family comets ?
even comet c/2001  q4 ( neat ) exhibited a large range in silicate - to - carbon ratios , varying from 2.75.7 in @xmath157  hour @xcite .
comet q4 , like hale - bopp , showed strong jets @xcite and q4 s high silicate fractions may originate in the ( local ) jet activity . in comet jets , violent sublimation of ices and volatile gases may excavate and entrain pristine grain materials from sub - surface reservoirs within the nucleus , similar to the _ deep impact _ event . visual inspection of the _ spitzer _ images of ht50 ( fig .
[ peakups ] ) show no distinct jet - features in the coma and there were no reports of any outbursts by this comet .
therefore , it may be that comet ht50 is dominated by `` whole surface '' ( global ) sublimation and that jupiter - family and long period / oort cloud comets are somewhat similar in the structure of their immediate surfaces .
indeed , it has been proposed that oort cloud comets can form a cohesive crust from normal re - surfacing processes , including galactic cosmic rays @xcite , supernovae @xcite , and interstellar grain impacts @xcite .
it is apparent that more work determining the mineralogy of comet comae is needed to understand the extent of the silicate - to - carbon ratio variations in comets , its correlation with nucleus activity , and possible connections to nucleus surface structure .
previous narrow - band photometry of comet encke indicated a weak silicate feature at small heliocentric distances @xcite but no feature at 1.2  au @xcite .
the coma is usually considered to be dominated by large particles , based upon a low coma color - temperature and a lack of a strong silicate feature @xcite . instead , our best - fit models indicate silicates are a minor constituent ( @xmath0% by mass for @xmath188  ) to the coma . to compare our results to previous investigations near perihelion
, we computed a `` silicate upper - limit '' model derived from our best - fit encke mineralogy in table  [ encke - minerals ] .
the upper - limit model sets the number of peak grains for the silicate minerals to their three - sigma upper - limits and decreases @xmath189 for the amorphous carbon component by three - sigma .
the mineral ratio by sub - micron mass becomes 100:1:18:11 ( amorphous carbon : amorphous olivine : amorphous pyroxene : crystalline olivine ) . the upper - limit model and photometry of comet encke at perihelion in 1987 july from @xcite is presented in fig .
[ encke - jul87 ] . also included is the encke nucleus as derived from table  [ encke - fit - table ] ( @xmath115  km , @xmath52 ) .
the coma and nucleus models were computed for the same geometry as the photometry data ( @xmath190  au , @xmath191  au , @xmath52 ) .
the coma model was scaled to account for encke s difference in dust production between @xmath192  au and 2.4  au .
the upper - limit model produces a silicate feature at 10   that approximates the weak silicate excess exhibited by the 913   photometry points .
the general shape of the sed is also approximated by the model , except for the short wavelength flux points , which are underestimated by factors of 1.21.8 .
the model likely requires a contribution from scattered sunlight to account for the 2.2   and 3.6   data points .
only the 5   data point remains unexplained ( the error bar may not reflect variations in atmospheric transparency through the optically thick 5   band - pass ) .
modifying the relative contributions of the silicates and carbonaceous minerals does not produce a satisfactory fit , although the possibility remains that the psd could be varied to account for the spectral shape .
we conclude that our _ spitzer_-derived mineralogy is a robust estimate of 2p / encke s coma composition and that our best - fit model s psd , derived at @xmath193  au , may not be representative of the psd observed throughout encke s entire orbit . indeed ,
the brightness behavior near perihelion is asymmetric and is best explained by different active regions on the nucleus , therefore , a varying psd about perihelion could be expected @xcite .
jupiter - family comets are thought to form in the transneptunian region and beyond @xcite .
we have derived mineralogies for two jupiter - family comets : comet encke , likely dominated by deeply absorbing grains , and comet tempel  1 dominated by silicate grains @xcite . these are striking differences for two comets assumed to arise from the same region of the proto - solar disk .
the difference in mineralogy may be due to each comet s original compositions or may be a reflection of different compositional processing histories .
ultraviolet and ion radiation transforms relatively optically inactive organic material into a dark , carbonized solid @xcite .
outgasing from surfaces or sub - surface layers with different radiation exposures may explain the differences in the coma mineralogy of comets encke and tempel  1 .
alternatively , the outgasing surfaces or sub - surface layers may have original and different compositions .
such layering in jupiter - family comet nuclei has been proposed by @xcite to account for topographical features on comets tempel  1 , borrelly , and wild  2 .
we present _ spitzer_/irs spectra of comets 2p / encke , 67p / churyumov - gerasimenko , and c/2001  ht50 ( linear - neat ) , and _ spitzer_/irac and mips images of comet 2p / encke .
comet encke exhibited a smooth continuum , best modeled by carbonaceous grains with a small peak grain size ( @xmath194  ) .
previous investigations into comet encke s dust coma revealed a weak silicate feature at perihelion ( @xmath195  au ) . we conclude the weak silicate feature is due to the paucity of silicate grains and the preponderance of carbonaceous grains ( or some other warm , deeply absorbing material ) .
we constrain the sub - micron silicate fraction to @xmath0% by mass .
the nucleus of comet encke is fit by the near - earth asteroid thermal model with an effective radius @xmath88 .
the nucleus was observed at phase angles 20 and 63 and may be exhibiting a variation of the infrared beaming parameter with phase angle , which indicates of a rough nucleus surface or appreciable night side temperature .
comet 67p exhibited a significant coma at a heliocentric distance of 5  au , @xmath196  w  @xmath76  @xmath59 at 27.9  .
67p s known dust trail comprises approximately 3% of the measured dust flux density .
the remaining coma flux was due to 1 ) recently ejected dust ( age of order hours to weeks ) , 2 ) large , slowly moving particles from the 2003 perihelion passage , or 3 ) some combination of the two .
comet ht50 displayed a significant silicate mineralogy with a silicate - to - carbon sub - micron mass ratio of 0.6 .
the derived ratio of 0.6 is an order of magnitude lower than the silicate - to - carbon ratios of post-_deep impact _
comet 9p / tempel  1 and other oort cloud comets , c/1995  o1 ( hale - bopp ) and c/2001  q4 ( neat ) .    the differences in silicate - to - carbon mass ratios in comet comae may be linked to strong jet activity in comets .
s derived silicate - to - carbon sub - micron mass ratio is 0.6 , but analysis of comet hale - bopp , which exhibited strong jet activity , derived a ratio of 8.1 . at this time , the wide diversity in comet comae mineralogy likely has not been probed .
the authors thank william reach for useful discussions and comments .
in addition , the authors wish to thank the referee for their insightful comments and suggestions that improved the manuscript .
this work is based on observations made with the spitzer space telescope , which is operated by the jet propulsion laboratory , california institute of technology under a contract with nasa .
support for this work was provided by nasa through contracts 1256406 , 1263741 , 1275835 , issued by jpl / caltech to the university of minnesota .
c.e.w . and m.s.k .
acknowledge support from the national science foundation grant ast-037446 .
m.s.k acknowledges support from the university of minnesota doctoral dissertation fellowship .
hanner , m.  s.  1983 , cometary exploration ; proceedings of the international conference , budapest , hungary , november 15 - 19 , 1982 .
volume 2 ( a84 - 47701 23 - 90 ) .
budapest , akademiai kiado , 1983 , p.  1 - 22 .
nasa - supported research .
, 2 , 1                                                                                + 2004 jun 25.23 & 210 & red , sl , sh , lh & post & 2.573 & 1.985 & 21 + 2004 jun 25.76 & 119 & blue , sl , ll & post & 2.577 & 1.982 & 21
+   + 2004 jul 15.34 & 2316 & red , ll , lh & post & 4.978 & 4.743 & 12 +   + 2003 dec 17.65 & 131 & red , sl , sh , lh & post & 3.238 & 2.652 & 16 + 2004 jul 18.44 & 131 & red , sl , sh , lh & post & 4.598 & 4.368 & 13 +      sl2 & 1&8591233 & -0&3691497 & 0&0425228 & -0&0017076 + sl1 & 0&8424725 & -0&0064947 & -0&0006418 & 0&0000008 + sh & 0&6111256 & 0&0466900 & -0&0040512 & 0&0000857 + ll2 & 1&9293931 & -0&1338932 & 0&0052684 & -0&0000686 + ll1 & 0&8016602 & 0&0051615 & -0&0002820 & 0&0000017 + lh & 1&0151328 & -0&0165096 & 0&0004606 & -0&0000064 +          2003 nov 17.75 & 3.55 & 1.986 & 0.061 & @xmath199 + & 4.49 & 5.06 & 0.11 & @xmath199 + & 5.73 & 1.007 & 0.021 & @xmath200 + 2004 jun 23.21 & 23.7 & 2.30 & 0.24 & @xmath201 + & 71.4 & 1.27 & 0.35 & @xmath202 + 2004 jun 29.92 & 3.55 & 7.0 & 2.8 & @xmath203 + & 4.49 & 3.95 & 0.42 & @xmath202 + & 5.73 & 1.07 & 0.11 & @xmath201 + & 7.87 & 3.99 & 0.12 & @xmath201 +      radius , @xmath77 ( km ) & 2.34 & 0.14 & 1.72 & 0.10 & + geometric albedo , @xmath204 & & & + ir beaming parameter , @xmath65 & 0.735 & 0.046 & & 1.026 & 0.061 + ir emissivity , @xmath70 & & & + reflected light scale , @xmath74 & 0.047 & 0.027 & 0.031 & 0.022 & 0.127 & 0.016 + degrees of freedom , @xmath205 & & & + @xmath85 & & & +              2003 dec 17.65 & 3.2 & 4.2 & 46.2 & 1.2 & 3.000 & 3.12 & @xmath208 & 0.12 & @xmath209 & & 1.27 & @xmath210 & 1686 & 4.2 + 2004 jul 18.44 & 4.6 & 3.7 & 29.6 & 1.2 & 2.857 & 2.98 & @xmath211 & 0.49 & @xmath212 & & & 1690 & 23.8 +      2p / encke & 2&6 & 1&00 & @xmath213 & @xmath2140 & & @xmath2140 & & @xmath2140 & & @xmath2140 & + c/2001  ht50 ( linear - neat ) & 3&2 & 0&63 & @xmath215 & 0&032 & @xmath216 & @xmath2140 & & 0&34 & @xmath217 & 0&59 & @xmath218 + c/2001  ht50 ( linear - neat ) & 4&6 & 0&82 & @xmath219 & 0&18 & @xmath220 & @xmath2140 & & @xmath2140 & & @xmath2142 & + c/1995  o1 ( hale - bopp ) & 2&8 & 0 & & 0 & & 0 & & 0 & & 8 & + 9p / tempel  1 ( pre-_di _ ) & 1&5 & 0 & & 1 & & 0 & & 0 & & + 9p / tempel  1 ( _ _ di__+1.0 ) & 1&5 & 0 & & 0 & & 0 & & 0 & & 3 & +