A&A 447, 589-595 (2006)
DOI: 10.1051/0004-6361:20053410
G. A. Galazutdinov1,5 - P. Gnacinski4 - Inwoo Han1 -
Byeong-Cheol
Lee2 - Kang-Min Kim2 - J. Kreowski3
1 - Korea Astronomy and Space Science Institute, Optical
Astronomy Division, 61-1, Whaam-Dong, Yuseong-Gu, Daejeon
305-348, Korea
2 -
Bohyunsan Optical Astronomy Observatory (BOAO) Jacheon
P.O.B.#1, YoungChun, KyungPook 770-820, Korea
3 -
Nicolaus Copernicus University, Center for Astronomy,
Gagarina 11, 87-100 Torun, Poland
4 -
Institute of Theoretical Physics and Astrophysics,
University of Gdansk,
ul. Wita Stwosza 57, 80-952 Gdansk, Poland
5 -
Special astrophysical observatory,
Nizhnij Arkhyz 369167, Russia
Received 12 May 2005 / Accepted 4 October 2005
Abstract
The recently published idea that intensities of some
weak diffuse interstellar bands (DIBs) are related to the C2molecule column density have been examined. We use a set of high
quality echelle spectra of heavily reddened stars, acquired at the
Bohyunsan Optical Astronomical Observatory (South Korea), with a
resolution R=30 000. The high quality (high S/N ratio) of our
spectra is proved by the fact that despite using the most widely used
Phillips (2, 0) band of the C2 molecule (near
8760 Å), we can trace the (3, 0) Phillips band (near 7725 Å) as
well. Equivalent widths of four (5176, 5542, 5546 and 5769
Å) out of 16 examined DIBs demonstrate relatively good
correlation with C2 column density. However, a majority of the
studied DIBs, already reported as "C2'' ones, most likely are
not related to this simplest carbon molecule. A removal of peculiar
objects like HD 34078 from the analyzed sample does not
substantially change the level of correlations.
Key words: ISM: lines and bands - ISM: molecules
Any observationally established relation between any of the
unidentified diffuse interstellar bands (DIBs) and any of the
well-known interstellar molecular absorptions may help to
identify DIBs, revealing the chemical nature of their still unknown
carriers. Thorburn et al. (2003) (hereafter,
Th03) proposed a set of DIBs, all of them very weak, to be closely
related to the C2 molecule. Interstellar space contains much
carbon or hydrocarbon chains discovered by
radioastronomical observations; C2 molecules are very likely
building blocks for such more complex species. Thus such a
relation, if confirmed, may give us some hints as to the chemical
processes that lead to formation of these DIB carriers.
Nevertheless, the absence of diffuse interstellar bands in the
circumstellar shell of the heavily reddened object IRC +10216,
where C2 was found abundant (Kendall 2002, 2004), may cast
doubts on whether any of DIB carriers are
related to C2, despite the lack of quantitative estimates of
C2 abundance.
![]() |
Figure 1: Phillips (2, 0) and (3, 0) bands of C2 molecule observed in spectra of some of our target stars (in top and bottom panels respectively). The telluric lines are divided out. |
Open with DEXTER |
In the context of relations between C2 molecules and at least certain DIBs, Douglas (1977) proposed bare carbon chains as possible DIB carriers. However, the laboratory gas-phase spectra of these species are known only for reasonably short chains i.e. until C5. For longer ones the projects of determining their spectra were conducted in noble gas matrices which makes the experimental spectra difficult to be directly compared to any observational material, but allows one to apply their common characteristics, e.g. intensity ratios. As one may expect abundances of carbon chains to be mutually correlated, a constant strength ratio of some DIB and the C2 column density may help to identify this DIB carrier.
The above considerations motivate a search for correlations between C2 features and strengths of certain DIBs, likely formed in (hydro)carbon chains. However, the existing body of C2measurements is very scarce. Moreover these measurements were performed (by Th03) toward stars that are not heavily reddened; this makes the weak DIBs very weak in their spectra and their intensity measurements uncertain. It is thus important to collect a sample of spectra of high quality, heavily reddened stars; one may expect relatively strong C2 bands as well as reasonably easily detectable weak DIBs. Another important factor should be a broad range of E(B-V)'s which should lead to a broad ranges of intensities of all other interstellar features too. The latter usually correlate more or less tightly to this popular parameter, despite that, as demonstrated in Fig. 3, this correlation is of just rather general sense. A correlation analysis is more reliable if performed on variables spanning relatively broad ranges of values.
The "C2'' DIBs as proposed by Th03 are listed in their Table 2 (18 entries). The weakness of these features makes precise measurements difficult. Apparently this is why Th03 claim that the 4979.58 Å DIB is "new'' (see their Table 2). In fact the feature has been described in two earlier papers: Galazutdinov et al. (2001) and Galazutdinov et al. (2002). The data presented in the two latter papers were acquired both at Terskol Observatory and at ESO which makes the detection certain. However, such a hard-to-detect interstellar feature must be weak. Figure 2 of Th03 shows the weakness of all considered DIBs, demonstrating the difficulties in measuring their equivalent widths. The S/N ratio that one can achieve inside weak DIB profiles must be low, allowing relatively large errors of EWs.
Usually the only sufficiently easily available C2 feature to ground-based observations is the Phillips 2-0 band situated in the near infrared (see Fig. 1). The details of such a band of a homonuclear molecule, full of weak and very narrow features, require high resolution and very high S/N ratio data. Measurements of all spectral features of C2 band(s) may very likely suffer relatively large errors. The BOES spectra, used in this project, allow us to completely remove the fringing (by flatfielding) which is normally a nuisance while highly sensitive, back-illuminated CCDs.
Heavily reddened stars are necessarily quite faint. This is why it is a very time-consuming task to collect a statistically meaningful sample of high S/N spectra of such objects. The sample collected for this project is an intermediate one (20 objects). However, the broad range of reddenings (see Table 1), which should result in a broad range of intensities of both C2 features and DIBs, should reveal the relation suggested by Th03 between the known molecule and the unidentified carriers of DIBs.
Table 1: Observed stars. Columns headings are: Star - object name, V - V magnitude, E(B-V) - reddening, Sp - spectral and luminosity class.
We have collected echelle spectra of 20 heavily reddened stars
(see Table 1) using the fiber fed echelle spectrograph attached to
1.8-m telescope of Bohyunsan Optical Astronomy Observatory (Korea)
during 2004 observing runs. As seen from Table 1 all the
targets are reddened, some of them very heavily. In total
the E(B-V)s of the sample stars range from 0.45 to 1.99 which
should make weak DIBs and molecular features strong enough to
allow detections and reasonably precise measurements but also
covering a broad range of values which facilitates a correlation
analysis.
The spectrograph has 3 observational modes with different
resolving power: 30 000; 45 000 and 90 000. Our spectra are in the
lowest resolution which is caused by the relative faintness of the
heavily reddened objects. The spectrograph allows us to cover in a
single exposure the wavelength range 3700 Å-
10 000 Å, divided into 75-76 spectral orders with CCD
camera equipped with a
pixels matrix (pixel size
15
m
15
m). The achieved resolution allows precise
measurements of wavelengths and intensities of interstellar
features.
Table 2: The equivalent width of DIBs (mÅ) measured in spectra of program stars.
![]() |
Figure 2: The fit of the calculated profile of the Phillips (2-0) C2 band to the observed one in the spectrum of HD 204827. |
Open with DEXTER |
![]() |
Figure 3: A comparison of two spectral ranges observed in HD 204827 and BD+40 4220. Note that intensity ratios of several DIBs do not resemble those of color excesses (E(B-V)s are 1.08 and 1.99 respectively). |
Open with DEXTER |
Our reduction of the echelle spectra have been performed using both the IRAF (Tody 1986) and the DECH codes (Galazutdinov 1992). The programs allow us to perform all standard procedures of CCD spectra processing and analyzing. The wavelength scale was constructed on the basis of a global polynomial of the form described in detail in Galazutdinov et al. (2000).
Our spectra contain the features of the C2 Phillips (2-0) band -
see the examples given in Fig. 1. It is evident that the sequence
of increasing C2 band intensity is not the same as the
sequence of increasing colour excesses. As a demonstration of the
quality of our spectra we present also the Phillips (3-0) band
observed towards our selected targets (Fig. 1). This band was
discovered by van Dishoeck & Black (1986) in the spectrum of
Oph. However, to achieve high enough S/N ratio they
observed this bright object for 7 h and present the C2features between telluric lines. Recently this band was reported in
spectra of heavily reddened stars by Gredel et al. (2001) and by
Gredel (2004). Our Fig. 1 demonstrates the band in the spectra of
stars listed in Table 1 without the telluric lines which have been
removed using the proper standard (Spica). Although the (3-0) band
is weaker than the (2-0) one it behaves in the same fashion, i.e.
the spectra form the same sequence of intensities - see Fig. 1.
The column densities of the C2 rotational levels were derived
from the (2-0) band of the C
Phillips System using the profile fitting technique (Fig. 2). The
molecular cloud velocity (v) and column densities of rotational
levels
to
were simultaneously fitted to the observed
spectra. The Doppler broadening parameter (b) was fitted together
with v and column densities only for sightlines with strong C2lines. For stars with shallow C2 lines the Doppler broadening
parameter was set manually, and an average C2 column density
for different b values was computed. The (2-0) band oscillator
strength
was adopted from Erman
& Iwamae (1995). The wavelengths and oscillator strengths for
individual lines were adopted from van Dishoeck & de Zeeuw
(1984).
The column densities of the higher (unobserved) rotational levels (J
)
of the C2 molecule were estimated from the Bolzmann function:
For three stars the rotational levels
or
were not observed. We have calculated
column densities for these levels from Eq. (1), but the rotational temperature was
obtained for each star independently. For HD 46149
K
(from
), for the
star BD +58 2580 the rotational temperature
K
(from
), and for HD 235825
K
(from
).
The total C2 column density (
)
presented in Table 3 is the sum of rotational levels up to
.
Th03 placed the spectra in sequence after being normalized to EB-V. In Fig. 3 we compare two fragments of the spectra of our heavily reddened targets. Despite the EB-Vs which are different by a factor of 2, several of the DIBs are of almost identical intensity in both spectra (5797, 5828, 5850, 6367, 6379 and 6425 Å). At the same time some rather broad DIBs are evidently stronger towards BD+40 4220 (5780 and 5844) while some narrow features are stronger in the spectrum of HD 204827 (6410 and 6439). This demonstrates the lack of a tight correlation between EB-V and column densities of many DIB carriers. The same concerns known molecular species; e.g., in Fig. 1 it is clearly seen that intensities of C2 features correlate quite poorly with reddening. We thus consider any normalization to EB-V = 1 as risky and we do not apply it in our further considerations.
Evidently all the so-called "C2 DIBs'' are very weak in the
spectrum of the star HD 34078 (AE Aur) where the C2 band is
relatively strong. In the same spectrum another
molecule, CH, is also revealed by very strong spectral features
(Kreowski et al. 1999). The weakness of the selected DIBs in
the target where C2 features are evidently strong casts
serious doubts on whether the set of DIBs is related to
the simple carbon molecule.
![]() |
Figure 4: The variable strength ratio of the strong 4726 DIB and its weak neighbour - 4734. The upper spectrum is the average of BD +59 2735, HD 34078, HD 204827, BD +60 493, HD 226868 and HD 186745 while the lower - of BD +59 456, HD 46149, BD +58 2580, HD 219287, HD 228779 and HD 235825. |
Open with DEXTER |
The Phillips (2-0) band of C2 is usually quite weak and thus its features are difficult to measure. This makes the sample of objects in which it was observed well above the detection level necessarily small. Moreover, all the so-called "C2 DIBs'' are also very weak. The S/N ratio inside their shallow profiles is low and thus any correlation between such weak features must be considered as only a first step toward selecting DIBs that may behave in a similar fashion.
To measure the EWs of weak DIBs we have constructed average spectra involving all our targets. In such very high S/N spectra one can easily trace all the weak DIBs discussed in this paper. The average was then scaled up or down, until it was well-fitted according to a visual inspection, to the observed features in every individual spectrum. The factor established while fitting the average to any of the spectra was then used to multiply the EWs measured on average. This method allows us to avoid uncertainties which otherwise can make measurements uncertain such as: borders of any feature, continuum level, possible strong noise in the feature bottom etc. It makes the sample of measurements more homogeneous than if EWs are determined by integrating very weak and noisy DIB profiles. We are aware of the fact that DIB profiles may vary from object to object, but, unless very large Doppler splitting, such differences should not be visible in our R=30 000 spectra (Galazutdinov et al. 2002). Our attempts to measure the 5004 feature mentioned by Th03 failed. The feature is either very weak or blended with a stellar one which makes possible measurements unreliable.
Table 3:
The column densities (
)
of the C
rotational levels (
).
The values written in boldface are derived from the (2-0)
band of the C
Phillips System. The values written in plain text were calculated from the Boltzmann distribution.
![]() |
Figure 5: The C2 column density versus equivalent width of the selected DIBs (mÅ). The correlation coefficient is given inside each plot. |
Open with DEXTER |
Assuming that the "C2'' DIBs are all related to the C2column densities one may also expect their mutual correlations. If any pair of the DIBs is really related then their strength ratio should be constant inside the measurement scatter. To check this in the case of the strongest "C2'' DIB and the weak new DIB 4734 we have constructed two averages. One of them involves six objects in which the second band seems present: BD +59 2735, HD 34078, HD 204827, BD +60 493, HD 226868 and HD 186745; we also averaged another 6 objects: BD +59 456, HD 46149, BD +58 2580, HD 219287, HD 228779 and HD 235825 in which we could barely trace the weak DIB. The result of the comparison is shown in Fig. 4; the strong 4726 DIB is of very similar strength in both averages while the weak 4734 feature is quite strong in one average and barely seen in the other. This shows that their intensity ratio is very variable and thus their common origin is doubtful. The average colour excesses, E(B-V), are for the two samples 0.99 and 0.93 respectively (see Table 1), so the weak 4734 feature does not relate to E(B-V) as well.
Our measurements of the "C2'' DIB EWs are compared to the C2 column densities in Fig. 5. We present here 16 out of 18 Th03 "C2 DIBs'' i.e. the features which we consider as measured with a reasonably good precision. It is evident that the level of correlation substantially differs from feature to feature. The correlation apparently does not depend on DIB strength. This is important as otherwise one could expect that stronger DIBs, generally measured with higher precision, demonstrate more tight relations. The variable degree of correlation tells us that this is not the case; the observed correlations show that different weak DIBs relate to a very different degree to C2 column densities.
Four weak features: 5176, 5542, 5546 and 5769 seem to be well-correlated to the C2 column density. They are the most likely "C2 DIBs''; however, a much more extensive set of data would be necessary to prove this. On the other hand, DIBs such as 4726, 4969 and 5763 seem not to correlate with the carbon molecule at all. Their EWs are sometimes even higher than those of the well-correlated DIBs. Apparently we observe a real effect. All other observed DIBs lie between good and poor correlation.
Statistically meaningful samples of spectra where absorption
features revealing atomic and molecular gas are detectable must
be observed to test hypotheses concerning chemical reactions that
take place in ISM and the preservation of some more or less complex
species. If the latter are homonuclear (as C2) their features
are numerous but very sharp and shallow. A high signal-to-noise
ratio is very important here. In the case of diffuse bands,
blending of certain features is a serious problem. Members of
blends like 5797+5795 or 6203+6205 (Porceddu et al. 1991; Kreowski et al. 1997) should be measured separately as their
origin is clearly different. Moreover, DIB profiles are not
exactly the same even if free of Doppler splitting (Galazutdinov
et al. 2002) which makes the EW measurements more complicated.
This is why some additional considerations are necessary to make
sure that any two DIBs can be considered as belonging to the
spectrum of a single species (Weselak et al. 2001).
The creation of the C2 "family'' of diffuse bands done by Th03 is certainly premature. Most of the selected features either evidently or likely do not correlate with the C2 column density. Moreover their mutual strength ratios also seem to be variable which makes their membership of any single DIB family not conclusive. Any realistic relation between any two features requires a statistically meaningful sample of objects to be established. On the other hand a single object contradicting any proposed relation, if realistic (like HD 34078), may disprove it. Diffuse bands demonstrate once again the complexity of their identification.
Acknowledgements
G.G. is grateful to the Korea MOST (Ministry of Science and Technology) (grant M1-0222-00-0005) and the KOFST for providing an opportunity to work at KASI through the Brain Poolprogram. J.K. acknowledge the financial support of the Polish State Committee for Scientific Research (grant 2P03D 019 23). J.K. and G.G. acknowledge also the NATO linkage grant PST.CLG.980313.