Free Access
Issue
A&A
Volume 584, December 2015
Article Number A24
Number of page(s) 6
Section Atomic, molecular, and nuclear data
DOI https://doi.org/10.1051/0004-6361/201525738
Published online 13 November 2015

© ESO, 2015

1. Introduction

The present work reports on measurements of transition data involving high-excitation levels and lines of singly ionized iron, Fe ii. The complex and line rich spectra of iron are observed in a wide variety of objects, such as stars, interstellar medium, solar coronae and quasars. In spectroscopy of stellar photospheres the atomic transitions most often appear as absorption lines. To use the spectral line for quantitative analysis of the population distribution and abundance determination in a stellar object, reliable f values must be available.

To meet the demands for accurate atomic data, Johansson et al. (2002) initiated the FERRUM project with the goal of providing evaluated transition data for astrophysical applications. Transitions from levels with different excitation potentials, and with a range of transition probabilities have been measured and calculated for the iron group elements, see e.g. Sikström et al. (1999), Hartman et al. (2005), Gurell et al. (2010), and references therein. The present paper extends the data of highly excited levels, around 10 eV. The amount of data available from experimental studies is limited, and the majority must be provided by calculations. The experimental data is important to benchmark different calculations and to provide the user with reliable uncertainties and priorities between data sets. This is especially true for the complex ions where level mixing is important and hard to predict.

A common technique to derive the oscillator strength is to combine measurements of the radiative lifetime of an upper level with relative intensities for all decay channels. The present paper reports on values using this approach, where the radiative lifetimes are determined with Time-Resolved Laser Induced Fluorescence (TR-LIF) from ions produced in a laser ablation plasma. These are combined with branching fractions obtained from intensity calibrated Fourier transform spectra using a hollow cathode discharge as a light source.

2. Radiative lifetime measurements

Figure 1 presents a partial energy level diagram for Fe ii showing the large number of parent terms. The close similarity in the energy differences between 4s and 4p for all parent terms (the promotion energy), give rise to a high Fe ii line density in the region around 230260 nm. In addition, the 4s4p and 4p5s/4d transition arrays have similar energy differences, and the latter transitions thus also appear in the same wavelength region. The target of the present study is the 3d6(5D)4d e6G term. This term has even parity which is the same as the ground configuration 3d6(5D)4s. This, and the large energy difference makes e6G impossible to excite in a one-photon process. However, the use of an intense laser beam opens the possibility for two-photon excitations from the ground term, which has a high relative population in the ablation plasma.

The TR-LIF set-up at Lund High Power Laser Facility is recently described in detail by Engström et al. (2014), where a study similar to this, but for Cr ii, is reported; here we only emphasize the important aspects. The set-up contains two Nd:YAG lasers operating at 10 Hz. One of them (Continuum Surelite) is frequency doubled and focussed onto a rotating iron target placed inside a vacuum chamber to generate the ablation plasma. The second laser is an injection seeded and Q-switched Continuum NY-82. The 532 nm output from this laser was temporally compressed using stimulated Brillouin scattering in water before pumping a dye laser (Continuum Nd-60) using DCM1. To obtain the 237 nm radiation needed for the two-photon excitation, the output was frequency tripled in KH2PO4 (KDP) and BaB2O4 (BBO) crystals and Raman shifted in a high pressure hydrogen cell. The laser-induced fluorescence was detected by a 1/8 m monochromator with a 280 μm wide entrance slit oriented parallel to the excitation laser beam and perpendicular to the ablation laser. All measurements were performed in the second spectral order giving an instrumental line width of 0.5 nm. The fluorescence signal was recorded with a micro-channel plate photomultiplier (PM) tube (Hamamatsu R3809U) with a rise time of 0.2 ns, and digitized by a Tektronix oscilloscope with 2.5 GHz analogue bandwidth. A second channel on the oscilloscope sampled the excitation laser pulse shape from a fast diode fed by scattered light from an attenuating filter. Each decay curve, as well as the temporal shape of the excitation pulse, was averaged over 1000 laser shots. The lifetimes were extracted using the program DECFIT (Palmeri et al. 2008). The program derives the lifetime from the measured laser pulse and the recorded decay. In the present experiment, the fluorescence data was fitted with a single exponential convoluted with the square of the measured laser pulse (since the upper state is excited through a two-photon absorption) and a constant background. The laser pulse had a temporal width (FWHM) of 2 ns. A typical example is shown in Fig. 2.

thumbnail Fig. 1

A partial level diagram of Fe ii showing the parent term structure, and the similar energies for all 4s4p and 4p5s/4d transitions. The arrow marks the transitions investigated in this study.

The final lifetimes given in Table 1 are the averages of around 20 measurements, performed during different days. The quoted uncertainties, between 10 and 17%, include the statistical contributions from the fitting process as well as variations in the results between the different measurements. The high probability of spectral coincidences between the two-photon pumping channel and a single photon excitation of the 4p levels had to be thoroughly considered. Another consequence of the similarity between the 4s4p and 4p4d wavelengths is that the secondary cascades (4s4p) coincide with the measured primary decay (4p4d). The chosen pump and detection channels are given in Table1, and the most difficult cases will be discussed in detail below. All measurements are checked for blending lines, both in the excitation and decay channels, and three of the levels are discussed in more detail below.

thumbnail Fig. 2

The first 20 ns of the decay of the 3d64d e6G11/2 level in Fe ii following two-photon excitation from the ground state (3d64s a6D9/2). The evaluated lifetime is 2.0 ± 0.2 ns. Background subtracted data points (+) are plotted together with a fitted single exponential decay convoluted by the square of the measured laser pulse (solid line). The dashed curve shows the squared laser pulse.

The high J value of the e6G13/2 level, and the ΔJ = 0, ± 1 restriction for E1 transitions, effectively limits the decay channels to only the (5D)4p level with J = 11/2, z6F11/2. Transitions to 4p levels belonging to other parent terms are several orders of magnitudes weaker. The pumping transition is located at 238.0 nm, and the primary and secondary decays fall at 237.64 nm and 238.28 nm, respectively. In this case the detection had to be shifted towards shorter wavelengths to avoid the influence from scattered laser light and secondary decay. However, shifting the detection wavelength too much to the blue brings in the fluorescence at 234.8 nm from 4p z4D7/2 with a lifetime of 3.02 ns (Guo et al. 1992), which is excited by a single 238.0 nm photon. This effect was clearly observed in the decay curve as a weak multi-exponential shape which resulted in a longer lifetime in our single exponential analysis compared to detection at longer wavelengths. By stepping the detection wavelength from this decay towards the primary decay at 237.64 nm, the contribution from z4D7/2 could be monitored, both in the derived lifetime and the fitted curve. From this behaviour and the known instrumental line width of 0.5 nm we estimated the maximum contribution from other decays to be 10%. This is reflected in the larger total uncertainty quoted for this level, 14%, compared to 11% for the other levels.

Table 1

Excitation and detection scheme for the measured levels, and results from this work and comparisons with previous calculations.

The level e6G11/2 was pumped at 237.2 nm from the ground state a6D9/2. The primary decay is at 237.0 nm whereas the secondary decay, a6D-z6F, falls at 238.2 nm. The monochromator was tuned to the blue side of the decay channel to minimize scattered laser light in the detection. In addition, the contribution from secondary decays was avoided.

The level e6G7/2 can be populated by two-photon excitation from the lower levels a6D3/2,5/2,7/2. The excitation from J = 3/2 was the best channel. The transitions from J = 7/2 is close to the excitation route to e6G13/2, but still possible. The transition from J = 5/2 is close to the laser frequency and could not be used.

3. Branching fractions

The BF for a transition from an upper level i to a lower level k can be expressed as: BFik=Aikj=1NAij=Iikj=1NIij·\begin{equation} BF_{ik} = \frac{A_{ik}}{\sum\limits_{j\,=\,1}^{N} A_{ij}} = \frac{I_{ik}}{\sum\limits_{j\,=\,1}^{N} I_{ij}}\cdot \label{eq:bf} \end{equation}(1)Here N is the total number of transitions from level i and j is summed over all lower levels. The first equality is the definition of BFik, i.e. the fraction of the atoms in level i that decay through the channel ik. The second equality assumes the case of an optically thin plasma, where the observed photon intensity of a line, Iij measured as photons per second, is proportional to Aij (Thorne et al. 1999). When all intensities are measured on a common calibrated scale, BFik can be derived from the measured intensities.

The lifetime of the upper level, τi, puts the relative intensities on an absolute scale. With

τi=1j=1NAij$$ \tau_i = \frac{1}{\sum\limits_{j\,=\,1}^{N} A_{ij}} $$and Eq. (1)we derive A from Aik=BFikτi·\begin{equation} A_{ik} = \frac{BF_{ik}}{\tau_i} \cdot \label{eq:A} \end{equation}(2)This requires all lines from the upper level to be measured. In practice it is usually impossible to include all transitions, either because they occur outside the detector range or they are too weak to be observed. Thus, we rewrite Eq. (1)as: BFik=11+1Iikj=1,jknIij+1Aikj=n+1NAij·\begin{eqnarray} BF_{ik} &= & \frac{I_{ik}}{\sum\limits_{j\,=\,1}^{N} I_{ij}} = \frac{1}{1+ \sum\limits_{j=1, j \neq k}^{n} \frac{I_{ij}}{I_{ik}} + \sum\limits_{j=n+1}^{N} \frac{I_{ij}}{I_{ik}}} \label{eq:residual} \\ &\approx& \frac{1}{1+ \frac{1}{I_{ik}} \sum\limits_{j=1, j \neq k}^{n} I_{ij} + \frac{1}{A_{ik}}\sum\limits_{j=n+1}^{N} A_{ij}}\notag\cdot \end{eqnarray}(3)Here n is the number of observed lines and the sum over the intensity ratios for the unobserved lines is approximated using the theoretical transition probabilities from Kurucz (2013).

The iron plasma was produced in a hollow cathode discharge (HCD), operated with neon as carrier gas and using currents between 0.6 and 0.8 A and pressures between 1 and 2 torr. The light emitted was recorded with a Fourier transform spectrometer (Chelsea Instruments FT500) in the range 27 64052 290 cm-1 and detected by a solar blind PM tube (Hamamatsu R166). In total, five spectra were measured. While a higher current in the HCD produced a more intense spectrum it also led to slightly wider lines due to the temperature dependent Doppler effect. Line widths between 0.18 and 0.25 cm-1 were observed with most of the spectra having 0.21 cm-1. To calibrate the intensity scale, three spectra with lower resolution were recorded using a deuterium lamp with known relative spectral radiance, measured at the Physicalisch-Technische Bundesanstalt, Berlin, Germany.

Despite the high spectral resolution of the spectrometer (resolving power about 2 × 105) the very line rich spectrum of Fe+ in this region produced a number of blends that must be taken into account to derive accurate experimental intensities. The first case is the z6F11/2e6G11/2 channel at 42 328.78 cm-1, where there are two additional Fe ii lines at 0.185 and 0.352 cm-1 higher wavenumbers, respectively (Johansson 1978; Nave & Johansson 2013). This blending was handled by fitting a sum of three Gaussian shaped line components, constrained to have the same FWHM and their separations fixed to the accurately known values, using the computer code GFit (Engström 2014). One of these fits is shown in Fig. 3. Consistent line widths and relative intensities were obtained in all five measurements.

thumbnail Fig. 3

The 4p z6F11/24d e6G11/2 line at 42 328.78 cm-1 blended by two other Fe ii lines at 0.185 and 0.352 cm-1 higher wavenumber, respectively. These separations have been kept fixed in the fitting of the three Gaussian shaped lines.

thumbnail Fig. 4

Experimental lifetimes are plotted with error bars (+), whereas the theoretical lifetimes (Kurucz 2013) are marked with crosses (×).

A more severe blending is found in the z6F5/2e6G7/2 channel at 42 375.86 cm-1 where another Fe ii line is only 0.11 cm-1 higher. Furthermore, contrary to the previous case the blending line is also significantly more intense. The experimental spectra were analysed using the same technique of constrained fitting as before. Also, in this case, we obtained the same FWHM as for resolved lines but the relative intensities of the two components varied somewhat between the five measurements. In addition to the three lines included in our work, Kurucz (2013) predicts an A value larger than 107 s-1 for four additional transitions (41 472, 40 478, 40 263 and 39 925 cm-1) from the e6G7/2. Unfortunately all these lines are too severely blended in our spectra to be included in our analysis. This causes the so-called residual in Table 2, i.e. the sum of all theoretical BFs from Kurucz (2013) for the unobserved lines from this level, to be as high as 34%. Since the experimental BFs, according to Eq. (3), depends on this residual, our values for the e6G7/2 level must be considered less reliable than for the other levels which is reflected in the higher uncertainties.

The last example of blending is the close coincidence between z6F7/2e6G5/2 at 42 607.80 cm-1 and a Fe ii transition at 42 607.79 cm-1. In addition, the lines can be expected to have similar intensities. Thus the 7/25/2 combination are omitted from the BF analysis.

Table 2

Resulting transition probabilities and oscillator strengths for the lines studied.

We present the final average BFs and transition probabilities obtained in Table 2, where we include the theoretical results by Kurucz (2013) for comparison. The error estimate in the BFs includes the uncertainty in the intensity calibration, in the area determination of the peaks, and a 50% uncertainty in the contribution from the unobserved branches obtained from Kurucz (2013), as described in Eq. (3). Since all but one of the observed decay branches are within 735 cm-1 of each other and, in addition, occur in a wavenumber region where the deuterium lamp calibration spectrum is quite flat, the uncertainty contribution from the calibration is small, about 3%. However, in the present case the total uncertainty of the transition probabilities for the strong lines is dominated by the contribution from the lifetime measurements. For the weak lines, the uncertainty in the BF is the dominating source. A more general discussion of uncertainties in experimental transition probabilities can be found in Sikström et al. (2002).

4. Discussion

Figure 5 shows a Hubble Space Telescope spectrum of the binary system chi Lupi (HD 141556), which has spectral types B9.5p HgMn and A2 Vm with effective temperatures of Teff = 10 650 K and 9200 K, respectively (Wahlgren et al. 1994). The investigated transitions 4p4d e6G are observed as prominent absorption features around 235 nm.

In Table 1 we compare our experimental lifetimes with the available theoretical data from Kurucz (2013) and Raassen & Uylings (1998a). Kurucz (2013) applied a superposition of configurations method using a modified version of the Cowan codes (Cowan 1981), whereas Raassen & Uylings (1998b) applied the orthogonal operator formalism. Both calculations make use of experimental level energies to improve the results. As seen in Table 1, the two theoretical results are within 7%, but they are about 25% shorter and do not overlap with the experimental values, even taking the rather large uncertainties into account. A graphical comparison is shown in Fig. 4, where a weak tendency of lower experimental lifetimes for lower J values is seen. This trend is within the error bars.

The discrepancy between measured and calculated BFs is relatively small for most lines. An interesting exception are the transitions from e6G7/2, where there is a large difference between the two theoretical studies. Although both give similar lifetimes they differ significantly in the transition probabilities for individual lines. For two of the transitions our values support those of Kurucz (2013), whereas the z4P5/2e6G7/2 line at 37 743 cm-1, which is predicted to be quite intense is barely observable in our spectra. In addition, the large residual for this level, as discussed above, leads to a large uncertainty in the experimental A values. The fact that the values by Kurucz (2013) differ significantly from his earlier 2010 calculations indicate theoretical difficulties for this level2. However, for all levels the resulting experimental A values are a combination of lifetimes and BFs, and the discrepancy between the experimental and theoretical results is in most cases a consequence of the lifetime differences.

The general agreement, except for the mentioned e6G7/2, of the experimental f values with the calculations by Kurucz (2013) and Raassen & Uylings (1998a) provides us with confidence to use either set for abundance determinations and also for high-excitation lines represented by the transitions and levels investigated in this work. We have a preference for the calculations by Kurucz (2013) when compared with the experimental values, however, the discrepancy between experiment and theory has to be considered when estimating the uncertainty in the abundance determination or evaluation of the stellar atmosphere models. We recommend that the experimental data is used where available, and complemented with theoretical calculations for the other lines.

thumbnail Fig. 5

Partial absorption spectrum of the chemically peculiar star chi Lupi, with the investigated Fe ii lines marked with tickmarks. The spectrum is observed by the Hubble Space Telescope and its Goddard High Resolution Spectrograph, and is obtained from the MAST data archive at the Space Telescope Science Institute. See Brandt et al. (1999) for a complete identification list and a comparison with a synthetic spectrum.


1

4-(Dicyanomethylene)-2-methyl-6-(4-dimethylaminostyryl)-4H-pyran.

2

The calculations from 2010 (referenced as K10 in Kurucz notation) are not available online any longer, but replaced by the new data Kurucz (2013).

Acknowledgments

We are grateful to the anonymous referee for valuable comments significantly improving the manuscript. This work was supported by the Swedish Research Council through the Linnaeus grant to the Lund Laser Centre and the Knut and Alice Wallenberg Foundation. H.H. gratefully acknowledges the grant No. 621-2011-4206 from the Swedish Research Council and support from the The Gyllenstierna Krapperup’s Foundation. Figure 5 is based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the MAST data archive at the Space Telescope Science Institute. STScI is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract NAS 5-26555.

References

  1. Brandt, J. C., Heap, S. R., Beaver, E. A., et al. 1999, AJ, 117, 1505 [NASA ADS] [CrossRef] [Google Scholar]
  2. Cowan, R. 1981, The Theory of Atomic Structure and Spectra (California: University of California Press, Berkeley) [Google Scholar]
  3. Engström, L. 2014, GFit, http://kurslab-atom.fysik.lth.se/Lars/GFit/Html/index.html [Google Scholar]
  4. Engström, L., Lundberg, H., Nilsson, H., Hartman, H., & Bäckström, E. 2014, A&A, 570, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Guo, B., Ansbacher, W., Pinnington, E. H., Ji, Q., & Berends, R. W. 1992, Phys. Rev. A, 46, 641 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  6. Gurell, J., Nilsson, H., Engström, L., et al. 2010, A&A, 511, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Hartman, H., Schef, P., Lundin, P., et al. 2005, MNRAS, 361, 206 [NASA ADS] [CrossRef] [Google Scholar]
  8. Johansson, S. 1978, Phys. Scr., 18, 217 [NASA ADS] [CrossRef] [Google Scholar]
  9. Johansson, S., Derkatch, A., Donnelly, M. P., et al. 2002, Phys. Scr. T, 100, 71 [NASA ADS] [CrossRef] [Google Scholar]
  10. Kurucz, R. 2013, http://kurucz.harvard.edu/atoms.html, accessed: 2015-03-03 [Google Scholar]
  11. Nave, G., & Johansson, S. 2013, ApJS, 204, 1 [NASA ADS] [CrossRef] [Google Scholar]
  12. Palmeri, P., Quinet, P., Fivet, V., et al. 2008, Phys. Scr., 78, 015304 [NASA ADS] [CrossRef] [Google Scholar]
  13. Raassen, A. J. J., & Uylings, P. H. M. 1998a, http://www.science.uva.nl/pub/orth/iron/FeII.E1, accessed: 2006-02-14 [Google Scholar]
  14. Raassen, A. J. J., & Uylings, P. H. M. 1998b, A&A, 340, 300 [NASA ADS] [Google Scholar]
  15. Sikström, C. M., Schultz-Johanning, M., Kock, M., et al. 1999, J. Phys. B Atom. Mol. Phys., 32, 5687 [NASA ADS] [CrossRef] [Google Scholar]
  16. Sikström, C. M., Nilsson, H., Litzen, U., Blom, A., & Lundberg, H. 2002, J. Quant. Spectry. Radiat. Transf., 74, 355 [Google Scholar]
  17. Thorne, A., Litzén, U., & Johansson, S. 1999, Spectrophysics: principles and applications (Springer Verlag) [Google Scholar]
  18. Wahlgren, G. M., Adelman, S. J., & Robinson, R. D. 1994, ApJ, 434, 349 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Excitation and detection scheme for the measured levels, and results from this work and comparisons with previous calculations.

Table 2

Resulting transition probabilities and oscillator strengths for the lines studied.

All Figures

thumbnail Fig. 1

A partial level diagram of Fe ii showing the parent term structure, and the similar energies for all 4s4p and 4p5s/4d transitions. The arrow marks the transitions investigated in this study.

In the text
thumbnail Fig. 2

The first 20 ns of the decay of the 3d64d e6G11/2 level in Fe ii following two-photon excitation from the ground state (3d64s a6D9/2). The evaluated lifetime is 2.0 ± 0.2 ns. Background subtracted data points (+) are plotted together with a fitted single exponential decay convoluted by the square of the measured laser pulse (solid line). The dashed curve shows the squared laser pulse.

In the text
thumbnail Fig. 3

The 4p z6F11/24d e6G11/2 line at 42 328.78 cm-1 blended by two other Fe ii lines at 0.185 and 0.352 cm-1 higher wavenumber, respectively. These separations have been kept fixed in the fitting of the three Gaussian shaped lines.

In the text
thumbnail Fig. 4

Experimental lifetimes are plotted with error bars (+), whereas the theoretical lifetimes (Kurucz 2013) are marked with crosses (×).

In the text
thumbnail Fig. 5

Partial absorption spectrum of the chemically peculiar star chi Lupi, with the investigated Fe ii lines marked with tickmarks. The spectrum is observed by the Hubble Space Telescope and its Goddard High Resolution Spectrograph, and is obtained from the MAST data archive at the Space Telescope Science Institute. See Brandt et al. (1999) for a complete identification list and a comparison with a synthetic spectrum.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.