A&A 447, L21-L24 (2006)
S. A. Levshakov1 - M. Centurión2 - P. Molaro2,3 - M. V. Kostina4
1 - Department of Theoretical Astrophysics, Ioffe Physico-Technical Institute, 194021 St. Petersburg, Russia
2 - Osservatorio Astronomico di Trieste, via G. B. Tiepolo 11, 34131 Trieste, Italy
3 - Observatoire de Paris 61, avenue de l'Observatoire, 75014 Paris, France
4 - Sobolev Astronomical Institute, St. Petersburg State University, 198504 St. Petersburg, Russia
Received 20 September 2005 / Accepted 4 January 2006
Aims. We analyzed the C I lines associated with the damped Ly system observed at = 1.15 in the spectrum of HE 0515-4414 to derive the 12C/13C ratio.
Methods. The spectrum was obtained by means of the UV-Visual Echelle Spectrograph (UVES) at the ESO Very Large Telescope (VLT).
Results. The obtained lower limit 12C/13C > 80 (2 C.L.) shows for the first time that the abundance of 13C in the extragalactic intervening clouds is very low. This rules out a significant contribution from intermediate-mass stars to the chemical evolution of matter sampled by this line of sight. The estimated low amount of 13C is in agreement with low abundances of nitrogen observed in damped Ly systems - the element produced in the same nuclear cycles and from about the same stars as 13C.
Key words: cosmology: observations - line: profiles - stars: nucleosynthesis - quasars: absorption lines - quasars: individual: HE 0515-4414
The evolution of the chemical composition of matter in the Universe is closely related to the history of star formation and destruction, stellar nucleosynthesis and enrichment of the interstellar/intergalactic medium (ISM/IGM) with processed material. The behavior of the stable CNO isotopes during the cosmic time is of particular interest since their ratios trace the production of primary and secondary elements which depends on the principal chemical evolution parameters: the stellar initial mass function (IMF), the rate of the mass loss from evolved stars, and the dredge-up episodes leading to the mixture of the surface layers with deeper layers of the star (Wannier 1980; Renzini & Voli 1981; Marigo 2001; Meynet et al. 2006).
12C is a primary product of stellar nucleosynthesis and is formed in the triple-process during hydrostatic helium burning (van den Hoek & Groenewegen 1997; El Eid 2005). 13C is supposed to have mainly a secondary origin, and is produced in the hydrogen burning shell when the CN cycle converts 12C into 13C (Wannier 1980). However, chemical evolution models considered by Prantzos et al. (1996) match observations better if a mixture of primary+secondary origin in intermediate-mass stars is assumed for 13C. Evolution of massive rotating stars at very low metallicities may also contribute to primary 13C (Meynet et al. 2006). The predicted values for the ratio = 12C/13C in the rotationally enhanced winds diluted with the supernova ejects are between 100 and 4000, whereas the 13C yields of massive non-rotating stars are negligible, . Besides, Meynet et al. show that the rotating AGB and massive stars have about the same effects on the isotope production. The difference between them is only in the isotope composition of the massive star wind material and the AGB star envelopes: the former is characterized by very low values of , while the latter have ranging between 19 and 2500 (the lower values correspond to the most massive AGB stars).
Chemical evolution models predict a decrease of the isotope ratio 12C/13C with time and an increase with galactocentric distance at a fixed time (Audouze et al. 1975; Dearborn et al. 1978; Tosi 1982; Romano & Matteucci 2003).
For instance, a photospheric solar ratio (Asplund et al. 2005), representative of the local ISM 5 billion years ago, is higher than the present value obtained through optical, UV, and IR absorption line observations as well as radio emission line measurements (see, e.g., Hawkins & Jura 1987; Centurión & Vladilo 1991; Centurión et al. 1995; Goto et al. 2003, and references therein). It should be noted, however, that on scales of 100 pc the local ISM is probably chemically inhomogeneous (Casassus et al. 2005) making this difference uncertain.
The high resolution observations ( ) of quasar absorption-line spectra available nowadays at large telescopes allow us to probe the isotope composition of the intervening damped Ly (DLA) systems through the analysis of C I lines, and, hence, to perform an important test of models of stellar nucleosynthesis outside the Milky Way. The most convenient C I transitions are from the 2u ( Å) and 3u ( Å) multiplets (Morton 2003) which cover the redshift range from to (look-back time 7.7-11.8 Gyr) in optical spectra of distant QSOs.
The analysis of the isotopic composition of the DLA systems is also connected to the interpretation of the hypothetical variation of the fine-structure constant ( ) at early cosmological epochs (Ashenfelter et al. 2004a,b; Fenner et al. 2005, hereafter FMG).
In the present letter we report constraints on the carbon isotope abundance in the = 1.15 sub-DLA system toward HE 0515-4414 (Reimers et al. 1998). The metallicity of this sub-DLA is found to be lower than solar: [Zn/H according to de la Varga et al. (2000), or [Zn/H as measured by Reimers et al. (2003) who revised the H I column density. Sections 2 and 3 describe the observations and the analysis. The obtained results and their astrophysical implications are discussed in Sect. 4.
The observations were acquired with the UV-Visual Echelle Spectrograph (UVES) at the VLT 8.2 m telescope at Paranal, Chile, and the spectral data were retrieved from the ESO archive. The seven selected exposures were described in detail in Levshakov et al. (2006, hereafter Paper I).
The C I lines analyzed in the present paper are located near the center of the echelle orders. This minimizes possible distortions of the line profiles caused by the decreasing spectral sensitivity at the edges of echelle orders.
Following Paper I, where all details of the data reduction are given, we work with the reduced spectra with the original pixel sizes in wavelength. The residuals of the calibrations give mÅ. The observed wavelength scale of each spectrum was transformed into vacuum, heliocentric wavelength scale (Edlén 1966). The instrumental profile is dominated in this case by the slit width (0.8 arcsec). The spectral resolution calculated from the narrow lines of the arc spectra is km s-1.
After the normalization to the local continuum,
the spectra from different exposures were rebinned with the step
equal to the mean pixel size. To eliminate instrumental velocity
calibration errors, we used the first exposure as a reference frame
and calculated residual differences in the radial velocity offsets of other
exposures through the cross-correlation analysis.
The aligned C I spectra were
co-added with weights inversely proportional to
The resulting average C I profiles are shown in Fig. 1 by dots with
1 error bars.
|Figure 1: Normalized intensities (dots with error bars) of the C I lines selected from the VLT/UVES spectrum of HE 0515-4414. The zero radial velocity is fixed at z = 1.150789 (QBR). The over-plotted smooth curves show the best fitted synthetic profiles ( , the number of degrees of freedom ). The dashed vertical lines mark positions of the C I sub-components. Note that the absorption feature consists of two sub-components with different oscillator strengths: f1560.68 = 0.0581 and f1560.70 = 0.0193. The 12C I oscillator strengths, laboratory vacuum wavelengths and isotopic shifts ( IS) are taken from Morton (2003).|
|Open with DEXTER|
The key points in the study of the isotopic abundance from the UV C I lines are as follows:
The heterogeneity of the absorbing medium may cause small relative velocity shifts between C I transitions with different J similar to the shifts observed between H2 lines in the direction of Ori A (Jenkins & Peimbert 1977). This would add an additional noise in the measurements of the disparity of the C I line positions. To avoid this putative effect we selected from the observed eight C I lines (QBR) only two J=1 transitions which show most pronounced isotopic shifts. These are just the lines C I and mentioned above. Other C I lines ( and ), as well as the lines with J=0 and 2 are less sensitive to the presence of 13C since their isotopic shifts are smaller than those of and .
Following QBR, we use two-component model to constrain the isotope ratio from the minimization of the objective function defined by Eq. (2) in Paper I. The values (normalized per degree of freedom) are calculated in the interval and shown as a function of in Fig. 2. This function gradually decreases with increasing and tends to a limit at . This global minimum lies within uncertainty range since at the expected mean value of is equal to . The optimized values and formal uncertainties of the model parameters for both C I absorption components (labeled by subscripts "1'' and "2'') at zero 13C abundance are as follows: the column density cm-2, cm-2, the Doppler parameter km s-1, km s-1, and the radial velocity difference km s-1. The corresponding synthetic profiles are shown in Fig. 1 by smooth curves. The model parameters do not change for .
|Figure 2: Confidence intervals in the "-'' plane calculated from the simultaneous fit of the C I lines shown in Fig. 1.|
|Open with DEXTER|
We also used a 3-component model to test robustness of the obtained results. Being applied to the lines of C I , C I , and C I this model yields and with at . The concordance of all available C I profiles is demonstrated in Fig. 3.
We re-analyzed the profiles of the C I lines associated with the sub-DLA system observed at = 1.15 in the spectrum of HE 0515-4414. Our main purpose was to set a constraint to the carbon isotopic abundance outside the Milky Way, in distant intervening clouds at the cosmological epoch corresponding to 8.2 Gyr of the look-back time. Similar tasks were discussed by Carlsson et al. (1995), Labazan et al. (2005), and FMG.
The present analysis gives 12C/13C > 80 (2 C.L.). This low abundance of 13C does not support the enrichment of gas by the wind from rotating massive stars (Meynet et al. 2006). Besides, it can constraint the fraction of intermediate-mass (IM) stars ( ) in the IMF. For instance, the FMG models for the chemical evolution of the Galaxy with normal IMF attain the solar ratio 12C/13C 90 already at [Fe/H for the outer radius and [Fe/H for the solar radius case, respectively. On the contrary, the models with an enhanced population of IM stars produce oversolar abundances of 13C at all metallicities [Fe/H] < -0.3 (cf. Fig. 9 in FMG). Thus, these models can be ruled out by our measurements.
|Figure 3: Same as Fig. 1 but for a 3-component Voigt profile model fitted to all available C I lines. The over-plotted smooth curves show the best fitted synthetic profiles ( , the number of degrees of freedom ). The dashed vertical lines mark positions of the C I sub-components.|
|Open with DEXTER|
Closely related to 12C/13C ratio is the problem of nitrogen, 14N, which is expected to be produced in the same way and from about the same stars as 13C. Measurements in DLAs show very low N abundances which is consistent with the absence of an enhanced population of IM stars to the chemical evolution of the gas (Centurión et al. 2003). Thus, low 13C and low 14N seem to be in agreement. It should be noted, however, that the chemical evolution of nitrogen is still not well understood since the existing models predict abundances higher than observed in DLAs (see, e.g., Fig. 6 in FMG).
The present bound to the abundance of 13C and, as a result, to the contribution of the AGB stars to the chemical evolution of the = 1.15 system may be indirectly related to the claims of a possible time variation of the fine-structure constant, (see Murphy et al. 2004, and references therein). A significant portion of the sample used by Murphy et al. involves the comparison of Mg II and Fe II wavelength shifts. Later on, it was shown that over-solar abundances (0.3 dex) of 25,26Mg isotopes with respect to 24Mg in the absorbing material can imitate an apparent variation of in the redshift range between 0.5 and 1.8 (Ashenfelter et al. 2004a,b). The production of Mg isotopes is believed to occur in about the same stars which produce 13C. Thus, the obtained bound to the amount of 13C, taken as a typical value for the DLA systems, poses a limit to the role of the AGB stars in mimicking the variations. However, the variability of has not been supported by more recent studies (Quast et al. 2004; Chand et al. 2004; Levshakov et al. 2005, 2006). This may suggest the presence of other systematic errors to explain the different results concerning the variability of .
The authors thank Edward Jenkins and Ralf Quast for useful comments. This work is supported by the RFBR Grant No. 03-02-17522 and by the RLSS Grant No. 1115.2003.2.