A&A 459, 257-263 (2006)
DOI: 10.1051/0004-6361:20065512
H. Rauer1,2 - M. Weiler1 - C. Sterken3 - E. Jehin4 - J. Knollenberg1 - O. Hainaut4
1 - Institut für Planetenforschung, DLR, Rutherfordstr. 2, 12489 Berlin, Germany
2 -
Zentrum für Astronomie und Astrophysik, TU Berlin, Germany
3 -
Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium
4 -
European Southern Observatory, Alonso de Cordova 3107, Santiago de Chile
Received 26 April 2006 / Accepted 4 August 2006
Abstract
Aims. We present observations of CN emission and the scattered solar light on cometary dust particles around the impact time of the Deep Space spacecraft (NASA) into the nucleus of comet 9P/Tempel 1. The purpose of the observations was to compare post-impact activity to the conditions pre-impact to search for new spectral emission lines after impact, to quantify the increase in gas activity due to the impact and to study the long-term activity changes.
Methods. We performed long-slit spectroscopy observations of comet 9P/Tempel 1 at the VLT, ESO, using the FORS instruments from July 2 to July 12, 2005. A wavelengths range of 370-920 nm was covered using two grisms. Four different position angle settings of the slit were applied each night with the projected Sun-comet line as standard setting, for which we report results here.
Results. The optical spectra of comet 9P/Tempel 1 showed the usual emission bands in the optical wavelengths range of the radicals: CN, C3, C2 and NH2. No new emission bands were detected after impact. The ejecta cloud of gas and dust caused by the impacting spacecraft into the cometary nucleus could be followed over the observing period. The projected expansion velocities have been determined. The night after impact we observed about
molecules of the CN parent in the ejected cloud. However, after five days the appearance of the gas and dust coma was back to pre-impact conditions.
Key words: comets: individual: Tempel 1 - comets: general
Comets are studied to investigate bodies which have been modified
only little since their formation during the early phases of the
solar system. However, although comets have been observed intensely
for a long time, we have little information
about their interior. Observations of the gas and dust coma provide only indirect evidence for the nucleus
composition. The Deep Impact mission (NASA) offered the opportunity to
study fresh sub-surface material for the first time. On July 4,
2005, a projectile of 370 kg mass impacted on comet P/Tempel 1 with
a relative velocity of 10.3 km s-1. A crater was produced on
the surface and a cloud of ejecta material expanded from the nucleus
(A'Hearn et al. 2005). This event was observed
not only from the impactor and fly-by spacecraft, but also in a
world-wide ground-based observing campaign (Meech et al. 2005).
![]() |
Figure 1:
Continuum subtracted spectra taken 1![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
At the European Southern Observatory (ESO) two observing campaigns were performed around impact time to investigate the gas and dust components of the comet. All ESO telescopes were involved in these campaigns to provide information on the target comet in all available wavelength band passes from optical to thermal infrared. Here, we report on the first results from optical long-slit spectroscopic observations of the gas and dust component. The scientific goals of these observations were to investigate whether new emission lines would appear in the spectrum after impact, quantify the gas sublimated during and after impact, search for new active surface areas on the comet and investigate for how long the additional activity is sustained. A comparison with pre-impact conditions is made to investigate possible differences in abundance ratios of the ices in the nucleus. In this paper, we concentrate on the CN and dust activity evolution of the comet. A comparison to other radicals emitting at optical wavelengths and present in our spectra is made in a subsequent publication.
The observations started on the night of July 02/03, two days before impact, and
lasted until July 10 at UT1 of the VLT, ESO, using FORS 2. Two
additional nights of measurements were then performed using FORS1 at
UT2. During the observations, 9P/Tempel 1 was at a heliocentric distance of 1.51 AU and a geocentric distance of 0.88-0.94 AU (see Table 1). Two grisms were used at UT1, covering in
total 370-920 nm. However, the red part of the spectral range, 610-920 nm, was covered only once per night, whereas the blue range
was the standard setting. The resolving power is 780 for the
370-620 nm range and 660 for the 610-920 nm range, using FORS 2.
For the spectra taken at FORS 1 the resolving power is 780. The slit length was 6.8
and a slit width of 1
was used to observe the comet.
The pixel scale is 0.252
pixel-1 (after a
binning)
in the spatial and 1.5 Å pixel-1 in the wavelength direction
for FORS 2. For FORS 1, the corresponding values are 0.20
pixel-1 and 1.2 Å pixel-1. These values correspond to a pixel scale of
162.3-129.2 km pixel-1 for FORS 2 and
135.2-135.9 km pixel-1 for FORS 1, respectively.
The detector of the FORS 2 instrument consists of two individual
pixel CCDs which are separated by a gap of 480
m,
corresponding to 4.03
.
The gap is parallel to the wavelength
direction. The FORS 1 instrument uses a single CCD with
pixels.
The position angle of the projected
solar direction ranges from 291.7
on the evening of July 2,
2005, to 290.3
on the morning of July 12. The slit was oriented at four different position angles, the
reference position being along the Sun-comet line. The additional spectra
were taken perpendicular to the Sun-comet line and at 45
angle in between. Here, we only report on the observations made
along the Sun-comet direction and the perpendicular direction. In addition to
the spectra, images were made at the beginning of each night in
broad band filters to study the dust coma of the comet. Table 1 provides an overview of the spectroscopic observing sequence for each
night.
All frames taken with FORS 2 showed a coherent noise pattern which had to be removed. This noise pattern had a low frequency contribution parallel to the wavelength axis and a high frequency component parallel to the spatial direction. The low-frequency component was basically constant over each night and could be well determined from bias frames, where it was clearly visible. For the high-frequency component parallel to the spatial direction the investigations showed, however, a slight variation in frequency during the night. This noise component was determined using unexposed edges of the CCDs of FORS 2 in each single exposure. A Fourier filter method was applied for both noise components to identify and remove the noise. Then, the standard CCD reduction including bias subtraction and flatfielding was performed. Some straylight contribution was found in frames with high exposure level. Again, the "non-exposed'' edges of each frame could be used to determine the level of the straylight in each spectrum, which was then subtracted. Since no "unexposed edges'' are present in data taken with FORS 1, a straylight correction could not be done for spectra taken on July 10/11 and 11/12.
The spectra were wavelength calibrated, sky background subtracted and corrected for extinction. For the sky background subtraction, separate sky frames taken through the night were used since the cometary coma fills the entire slit length. Absolute flux calibration was performed using spectroscopic standard stars observed prior to and after the comet observations. Observations during all nights, except during the night of July 06/07, were obtained under photometric conditions. To subtract the continuum emission caused by scattered solar light on the cometary dust particles in the coma, a solar catalog spectrum (Kurucz et al. 1984) was convolved with a Gaussian instrumental profile to the resolution of our observations and then fitted to the data in each row by a polynomial fit. An example of calibrated and continuum subtracted spectra taken before and after impact is shown in Fig. 1.
To investigate whether new emission lines appeared after the impact of the spacecraft into the comet, we compared spectra taken before and after the event. To enhance the signal-to-noise ratio in the reference before impact, all spectra taken on July 3/4 were co-added. After impact, the ejecta cloud moves through the field-of-view and only individual spectra were used for comparison. All gas emission lines have increased in intensity (Fig. 1) after impact. However, no new emissions were present, and thus, no new molecules have been detected.
Figure 2 shows the intensity distribution of the
continuum, caused by solar light scattered by dust particles, along
the slit parallel to the Sun-comet line. Two profiles are shown for the night of July 04/05 and one for all other nights. In Fig. 3, the radial continuum intensity profiles along the perpendicular direction are shown. The ejected dust cloud can be seen in the slit
after impact in the sunward direction and 270
direction measured with respect to the projected solar direction. The cloud moves and expands along both directions. After five days the cloud vanished and the coma is back to an intensity profile
similar to pre-impact. No change in the intensity profile is seen in the 90
and 180
directions.
![]() |
Figure 2: Continuum flux along the Sun-comet line, measured at 5540 Å to 5575 Å. The different colors give the profiles at different nights. Isolated peaks in the profiles are caused by stars within the slit. Positive distance is counted towards the projected solar direction, negative in the antisolar direction. |
Open with DEXTER |
![]() |
Figure 3:
Continuum flux perpendicular to the Sun-comet
line, measured at 5540 Å to 5575 Å. The different colors give the profiles at different nights. Isolated peaks in the profiles are caused by stars within the slit. Positive distance is counted towards the 270![]() ![]() |
Open with DEXTER |
The corresponding spatial profiles in the light of CN emission (violet system at 3875 Å) are shown in Figs. 4 and 5. Like in the continuum profiles, a cloud of ejected
material moves in the sunward direction after impact. However, the
CN cloud is visible in all directions (0,
90
,
180
,
and 270
)
already 18 h after
impact, again expanding and diluting in the subsequent nights. The fast lateral
expansion of gas molecules in the coma causes the gas cloud to spread around the nucleus
quickly. The long-slit spectra indicate that then a shell-like cloud of gas
molecules expands around the nucleus. Four days after impact, the
cloud has diluted and is visible just by a higher intensity profile in
comparison to pre-impact. Again, at the end of our observing period
the CN profile is back to pre-impact conditions.
![]() |
Figure 4: CN emission flux along the Sun-comet line for different nights. Positive distance values are counted in the projected solar direction, negative values in the antisolar direction. |
Open with DEXTER |
![]() |
Figure 5:
CN emission flux perpendicular to the Sun-comet line for different nights. Positive distance values are counted in the 270![]() ![]() |
Open with DEXTER |
We notice another difference between the continuum and gas intensity
profiles. After the impact cloud has been ejected, the intensity in
the continuum remains constant in the innermost coma near the
nucleus. In the light of CN, however, it rises again on July 6
(Fig. 4). Since the sky conditions of that night were
non-photometric, it remains unknown whether the observed increase is
indeed caused by an increased cometary activity. Then, the CN intensity at the nucleus position is
going back to the pre-impact value. In order to study short-term
variations, we investigated the brightness in the inner coma. Figure 6 shows the integrated flux over the inner pixels
around the nucleus, as measured in the spectra of the 0-180
direction. Because the pixels at and near the nucleus
position are most affected by any possible mismatch between the
fitted underlying continuum and the spectra, the innermost 7 pixels
have been excluded. Then the adjacent 20 pixels (5
)
on each
side have been summed. This approach was chosen to minimize the influence of the impact cloud. The resulting integrated CN flux around the
nucleus is also shown in Fig. 7, phased with the
rotation period of the comet and normalized to the mean flux value.
The rotation period was taken as (
) h (A'Hearn
et al. 2005).
The data points in the night following the impact (July 04/05) are
higher than intensities obtained on other nights at similar rotational phases. The most likely
explanation is the increase in CN activity due to the impact cloud,
which is also seen in Fig. 4 and which fades during our observing
period.
At small rotation phases, the CN flux is increased by about 30% compared to the levels at larger phase (>0.6). Unfortunately, the data points on July 05/06 have a large error because this night was not photometric. The data points of July 07/08 are also enhanced and show a clear decrease in intensity over the night with increasing rotational phase angle. The observed intensity variation fits to the proposed rotational modulation found from small aperture high-spectral resolution observations obtained at VLT UT2 using UVES (Jehin et al. 2006). This variation, therefore, may be caused by an active area on the nucleus surface rotating in and out of sunlight, which was already present before impact. However, again it is unfortunate that data at small phase angles prior to the impact are not covered by our observations. We therefore can not exclude that we still see the effect of the impact cloud. Such an interpretation would be supported by the low flux level observed at small phase angles on July 10/11, a time period when the spatial profiles show that the impact cloud has already been diluted (Fig. 4). However, since the aperture used for the determination of the data points in Fig. 7 is very small, the ejecta cloud is providing a substantial contribution to the flux in the aperture for at most two days after impact. Therefore, small aperture observations like performed with UVES or in our small integration window should be less sensitive to the impact cloud. By comparison, the wider angle observations of Schleicher et al. (2006) show a stronger contribution from the impact cloud. In such data, therefore, rotational flux variations are diluted. To definitely conclude on the effect of rotational modulation, data at different phase angles and observing times from different sites must be combined in future.
![]() |
Figure 6:
CN line flux in the inner coma (![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 7:
CN line flux over rotational phase. The values displayed are taken at the 0![]() ![]() |
Open with DEXTER |
We measured the projected expansion velocities of the dust and gas
ejecta cloud. This was done by determining the projected
nucleocentric distance where the intensity profiles of the cloud
equals the background coma (see Figs. 2-5). The projected mean expansion velocities are given in Table 2. The velocity of the outermost part of the impact gas cloud can be
determined only in the night after impact, before the cloud expands beyond the edges
of the slit and dilutes. The lower dust velocity can also be determined in following
nights. In the night of July 05/06, the tip of the dust cloud was
within the gap between the two CCDs of FORS 2 and an accurate estimate of
the dust expansion velocity was therefore not possible. In the
following night, the dust velocity could be estimated again. It shows a clear
decrease from night to night, which is caused by radiation pressure starting to act
efficiently on the dust particles, accelerating them anti-sunward. On
July 07/08 and later the dust cloud is too diluted to allow for a clear velocity determination. As expected from the influence of radiation pressure, the deceleration along the sunward direction is more efficient than along the 270 direction.
Table 2:
Projected mean expansion velocities of
the outermost part of the dust and gas cloud.
gives the time since impact and
p the projected distance from the nucleus of the outermost
detectable part of the impact cloud. PA lists the position angle at which the measurement was done, with respect to the projected solar direction.
Table 3:
Projected mean expansion velocities of
the center of the CN gas cloud. gives the time since impact and
p the projected distance from the nucleus. PA lists the position angle at which the measurement was done, with respect to the projected solar direction.
In addition, the velocity of the projected center of the CN cloud was determined on
July 04/05 and July 05/06. When determing the position of the center of the impact cloud, the geometric dilution of the cloud while moving to larger nucleocentric distances has to be taken into account. For a steady-state coma, one would expect a decrease of the column density, and thus the observed line brightness, according to ,
when
is the projected nucleocentric distance. The expansion of the gas cloud resulting from the Deep Impact event is a transient phenomenon, and therefore probably not following a
dependency. Since the distribution of the cloud material with time and nucleocentric distance is not known, in this work a weighting of the measured line flux with
is applied as a first guess. Each line flux value along the radial intensity profiles was multiplied by the square of its distance to the nucleus position. If a larger power of
is used for correction, the derived velocities of the center of the cloud would be lower. Then, a Gaussian profile was fitted to the resulting corrected cloud and the position of the center of the Gaussian was assumed to be the projected center of the CN cloud. The derived radial positions and corresponding velocities of the center of the CN cloud are given in Table 3.
The next step of investigation is to quantify the increase in gas
emission due to the impact. First, we concentrated on the longer
term effects and used a Haser model (Haser 1957) to compute the gas
production rate of CN before and 5-8 days after impact. At this
time, the ejecta cloud had already left the field of view and a Haser model
was a good approximation of the coma column density. Scale lengths
of CN were determined from the spatial coma intensity profiles. Because of the asymmetry of the coma, we treated the radial intensity profiles along each position angle of the slit separately. Because of the rather long
daughter scale length for CN, it could be determined only from one
spectrum where the comet was close to the edge of the slit. The
destruction scale length derived was then assumed to be the same
for both sides and all spectra. The resulting scale lengths for
parents, ,
and daughters,
,
for the comet at a heliocentric distance of
1.51 AU are:
![]() ![]() ![]()
|
Table 4: CN gas production rates determined using a Haser model. Columns give the values derived from spatial profiles along the different position angles with respect to the projected solar direction. All values are in [1024 s-1].
During the days immediately after impact we see the ejecta cloud travelling outwards. Obviously, it is not possible to use a steady-state Haser model to derive the production rates of the comet at this time. Since the details of the observed transient coma phenomenon, such as the distribution of the CN parents and their velocities, are unknown, a simplistic approach is chosen. As long as the complete CN cloud lies within the slit used for observations, as it was the case in the night of July 04/05, the total number of CN molecules in the cloud can be estimated assuming isotropic cloud expansion. It is further assumed that all CN parent molecules in the ejecta cloud were released at the same time, during impact. The approximation of the same age for all CN parent molecules in the impact cloud is justified as long as the time of release of the parent molecules is short compared to the observational time scale.
If a simple two step chemical reaction network is
assumed, with one reaction for the formation of CN and reaction rate coeffient ,
and one for its
destruction with
,
the total number of CN molecules in the impact cloud as
a function of time is obtained from the solution of the differential
equations
![]() |
= | ![]() |
(1) |
![]() |
= | ![]() |
(2) |
![]() |
(4) |
Table 5:
Number of CN radicals within the impact
cloud and the derived number of parent molecules released by the
impact. The values derived from spectra in different directions are
given, assuming a spherically symmetric coma, respectively. t-t0 gives
the time since impact.
corresponds to the total number CN
molecules in the cloud,
is the number of parent
molecules, most likely HCN, derived from Eq. (3). All values are given in [1029 molecules].
The impact of the Deep Impact spacecraft into comet P/Tempel 1 led to the formation of a crater on its surface and an ejecta cloud of gas and dust expanding into the cometary coma (A'Hearn et al. 2005). Our observations of this event performed by long-slit spectroscopy at the ESO VLT prior to and post impact showed the ejecta cloud in the light of the radicals seen in optical spectra of comets and in the continuum light caused by solar radiation scattered on dust particles. It was discussed whether new molecules would become visible from fresh sub-surface ice layers of the nucleus or formed by the high density and temperature conditions during the impact. Our observations could, however, not show any new emission lines after impact.
The ejecta cloud was covered in the long-slit spatial profiles and
was followed until 8 days after the impact event. The cometary coma,
however, showed conditions similar to pre-impact already 5-6 days
after the event. This is supported by the spatial CN and dust
distribution (Figs. 2-5) and a quantitative
comparison of CN production rates (Table 4) of 5-7 days
after the impact event with their values before impact. The projected expansion velocity of the CN
gas cloud is around 1.2 km s-1 for the outermost part (see Table 2) on
the day after impact. The center of the cloud moves with a mean projected velocity of about 0.6 km.
The slower dust cloud moved with about 0.21 km s-1 in the projected solar direction within one night after impact, and decreased to a mean velocity of
0.10 km s-1 until July 07, showing the effect of the solar radiation pressure. In the 270
direction, the projected dust velocities are somewhat higher, decreasing from about 0.37 km s
to 0.21 km s
in the same time.
We estimated the total number of CN parent molecules ejected by the
impact event assuming a two-step formation and destruction model of
CN to take into account the chemical processes in the cloud from
impact time up to our observations. We assume all parent molecules
released at the same time during the impact and isotropic conditions
in the coma. With these assumptions we estimate a total of about
parent molecules ejected by Deep Impact. Assuming HCN
as the main parent species, this would correspond to a mass of (
) tonnes.
The impact event provided a total kinetic energy of
J. With a sublimation enthalpy of 35.6 kJ/mol for HCN
(Stephenson & Malanowski 1987), a total of
molecules of HCN could be sublimated due to the impact event if all
energy would be used for the sublimation of HCN only. This value is of
the same order of magnitude as the number of CN parents produced by
the impact. Since HCN is only a minor component of cometary nuclei
(for 9P/Tempel 1, an abundance ratio HCN/
of (
)% was determined (Mumma et al. 2005)), the kinetic energy provided by the impact spacecraft is not sufficient to explain the observed amount of volatiles released due to the impact. Sublimating icy grains ejected into the coma by the impact event are a likely explanation for the observed amount of volatiles in the impact cloud.
The assumption of an instantaneous release of all parent species may nevertheless be justified since studies of the dust cloud (Schleicher et al. 2006) indicate typical grain sizes below about 2.5
m. If icy particles also have such low sizes and contain dark material, the complete sublimation of the grains takes less than 30 min (Beer et al. 2006) at heliocentric distances between 1-2 AU. Therefore, the approach chosen may provide an approximation of the total amount of sublimated HCN.
Acknowledgements
The authors warmly thank Jean Manfoid for his support during the observations with FORS1.
C. Sterken acknowledges financial support from the Belgian Fund for Scientific Research (FWO).
Table 1:
Observing log including all long-slit spectra.
The Table provides the date and
time of each observation, the exposure time, exp, and the
wavelengths range,
,
covered in each spectrum. The
position angle of the slit, PA, is measured from the projected solar
direction towards North. Observations marked with
were done using FORS 1, others were performed with FORS 2.
Night not photometric.