Issue 
A&A
Volume 561, January 2014



Article Number  A32  
Number of page(s)  5  
Section  Atomic, molecular, and nuclear data  
DOI  https://doi.org/10.1051/00046361/201322751  
Published online  20 December 2013 
log gf values for astrophysically important transitions Fe II
^{1} School of Mathematics and Physics, Queen’s University Belfast, Belfast BT7 1NN, UK
email: a.hibbert@qub.ac.uk
^{2} Indian Maritime University, Kolkata Campus, P19, Taratala Road, 700088 Kolkata, INDIA
Received: 25 September 2013
Accepted: 31 October 2013
Aims. In a recent measurement, Meléndez & Barbuy (2009, A&A, 497, 611) report accurate log gf values for 142 important astrophysical lines with wavelengths in the range 4000 Å to 8000 Å. Their results include both solar and laboratory measurements. In this paper, we describe a theoretical study of these lines.
Methods. The CIV3 structure codes, combined with our “finetuning” extrapolation process, are used to undertake a largescale CI calculation involving the lowest 262 finestructure levels belonging to the 3d^{6}4s, 3d^{7}, 3d^{5}4s^{2}, 3d^{6}4p, and 3d^{5}4s4p configurations.
Results. We find that many of the 142 transitions are very weak intercombination lines. Other transitions are weak because the dominant configurations in the two levels differ by two orbitals.
Conclusions. The comparison between our log gf values and the experimental values generally shows good agreement for most of these transitions, with our theoretical values agreeing slightly more closely with the solar than with the laboratory measurements. A detailed analysis of the small number of transitions for which the agreement between theory and experiment is not as good shows that such disagreements largely arise from severe cancellation due to CI mixing.
Key words: atomic data / relativistic processes / methods: numerical
© ESO, 2013
1. Introduction
The Goddard High Resolution Spectrograph (GHRC) and the more recent Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST) have produced many iron lines, particularly lines of Fe II, over a wide wavelength range. With the groundbased highresolution fibrefed echelle spectrograph mounted on the 3.6m telescope at ESOLa Silla Observatory under the HARPS GTO scheme, new spectroscopic parameters for different ions, including Fe II, are emerging. In one such measurement, Sousa et al. (2008) has reported log gf values for 263 Fe I and 36 Fe II weak lines and were subsequently used to determine iron abundances. Meléndez & Barbuy (2009) have presented log gf values for weak 142 Fe II lines. These solar and laboratory measurements, claimed to be accurate and precise, included the Fe II lines measured by Sousa et al. (2008). For the transitions in common in these two independent measurements, there is generally good agreement. In the present investigation, we have undertaken a largescale calculation in which we represent the wave functions of the levels involved in these (and other) transitions by means of configuration interaction (CI) expansions. We are then able to compare our calculated oscillator strengths with the experimental results given by Sousa et al. (2008) and by Meléndez & Barbuy (2009). A preliminary report was given in the proceedings of the XXVIIth ICPEAC conference (Deb & Hibbert 2012).
2. Method of calculation
Our group has undertaken a number of recent CI calculations of transitions in Fe II. There are many similarities in the approaches we have adopted in these calculations, but the method of optimisation of some of the orbitals has varied from calculation to calculation. Corrégé & Hibbert (2005, 2006) used orbitals up to 7s, 7p, 6d, and 4f to study a number of transitions involving 3d^{6}nl and 3d^{5}4s^{2} or 3d^{5}4s4p states. One of the difficulties that arise is the significant variation in the form of the optimal 3d function between different states. For transitions between levels of these configurations, the 3d orbital that appears in the dipole matrix elements is represented best by the HartreeFock orbital of the 3d^{6}4s ^{6}D state. Subsequently, we (Deb & Hibbert 2010a,b, 2011) studied forbidden transitions, but since many of the transitions involve the 3d^{7} levels, we chose the 3d function to be the HF function for the 3d^{7}^{4}F state. However, in all our calculations, the remaining dfunctions were used to account for the variation in the form of the optimal 3d function in different configurations. In the present work we have reverted to the orbitals given by Corrégé & Hibbert (2005), since the transitions reported by Meléndez & Barbuy (2009) mostly have the form 3d^{6}4s–3d^{6}4p.
The CIV3 program of Hibbert (1975), Hibbert et al. (1991) has been used to carry out the present structure calculation of Fe II. The form of the CI wave functions is, in LSJ coupling, (1)In (1), Φ_{i} represents a configuration state function (CSF) which is constructed from a common set of oneelectron orbitals of the form (2)where Y is a spherical harmonic, χ a spin function (often denoted by α or β for or , respectively), and the radial functions in (2) are expressed in analytic form as linear combinations of normalised Slater orbitals (STOs): (3)and where the STOs take the form (4)Also in (1), α_{i} represents the angular momentum coupling scheme and other necessary labelling, and a_{i} is the corresponding component of the eigenvector, associated with this wave function, of the diagonalised Hamiltonian matrix whose typical element is H_{ij} = ⟨ Φ_{i}  H  Φ_{j} ⟩.
We have sought to represent the 262 finestructure levels belonging to the 3d^{6}4s, 3d^{7}, 3d^{5}4s^{2}, 3d^{6}4p, and 3d^{5}4s4p configurations with energies up to 88 615 cm^{1}. We treat this set of configurations as our reference set. Beyond this energy, levels of the 3d^{6}5p configuration begin to appear in the spectrum, and some of the wave functions require careful treatment of CI mixing, particularly between 3d^{6}5p and 3d^{5}4s4p. The work of Corrégé & Hibbert (2005) was specifically focussed on this type of mixing, but the transitions reported by Meléndez & Barbuy (2009) involve levels lying well below any of the 3d^{6}5p levels. Below energies of 88 615 cm^{1}, there are a few levels that belong to the 3d^{6}4d and 3d^{6}5s configurations, but they do not show any significant coupling with the even levels considered in the present calculations. Similarly, because we have not included 3d^{6}5p in the reference set, the mixing between 3d^{6}5p and 3d^{5}4s4p will not be correctly represented.
Corrégé & Hibbert (2005, 2006) were able to include only a few sextet and quartet states, in order to keep the configuration size manageable, and so no doublets were included in their calculation. Nor did they include the 3d^{7} configuration. In the present calculation, we include all doublets, quartets, sextets, and one octet belonging to the above configurations and below 88 615 cm^{1}. The configuration set in the present calculation is therefore significantly larger that those in Corrégé & Hibbert (2005, 2006).
The set of configurations was chosen in two steps: (a) keeping subshells up to 3p filled, we included single and double replacements of orbitals from the outer shells by any of the available orbitals, for each configuration in the reference set; (b) then we included configurations formed by replacing up to two orbitals (in total) from the 3s and 3p subshells by available orbitals up to n = 6. This produced a large number of separate configuration state functions (CSFs).
For each LSπ symmetry, we first undertook a calculation in LS coupling. We found that a significant number of the CSFs had very small eigenvector components, so we deleted those whose  a_{i}  had a value less than 0.001. In previous calculations, we found that such a cutoff made very little difference to the calculated energies or to the eigenvector components of the remaining CSFs.
Relativistic effects were introduced in the BreitPauli approximation, with the following operators: mass correction and Darwin terms, and a modified spinorbit term H_{so} that allows for both the nuclear spinorbit effect and the main part of the spinotherorbit effect: (5)where ζ(l) is a parameter that only depends on the lvalue of the interacting electrons in the BreitPauli Hamiltonian matrix element. The values of the parameters ζ(l) were chosen so that, for certain key CSFs, the matrix element of H_{so} reproduces the sum of the matrix elements of the full spinorbit and spinotherorbit operators. This process is then independent of the numbers or types of the CSFs included in the CI expansions. The most suitable parameters were found to be ζ(s) = 0.0, ζ(p) = 0.88, ζ(d) = 0.59242, and ζ(f) = 0.52908. The number of CSFs included in our calculation, for each Jvalue, is given in Table 1.
Number of CSFs used.
Finally, we adopted our customary practice of “finetuning”, whereby we make small adjustments to the Hamiltonian matrix elements so as to bring the calculated energy differences (i.e., above the ground state) as closely as possible into line with experimental values.
3. Results and discussions
The oscillator strengths determined by Meléndez & Barbuy (2009) are all weak. Most correspond either to intercombination lines (E1 transitions that in LS coupling would be zero but become nonzero when spinorbit mixing is introduced, resulting in the upper and lower levels of the transitions being described by wave functions containing CSFs with the same J but different L or S) or to transitions in which the dominant CSFs of the two levels differ by more than one orbital. In the latter case, a single configuration approximation would again result in a zero oscillator strength, assuming orthogonal orbitals, but configuration mixing makes a nonzero result possible. Moreover, in the calculation of many of these transitions, configuration mixing results in cancellation of contributions from different CSFs, and in some cases this cancellation is substantial, lowering the oscillator strength by one or occasionally more orders of magnitude.
An interesting and important example is provided by the three transitions of multiplet 42. These are usually described as 3d^{5}4s^{2}^{6}S_{5/2}–3d^{6}4p ^{6}P, with J = 3/2, 5/2, 7/2: a notionally twoelectron change between upper and lower levels. The nonzero nature of the calculated oscillator strength arises because of mixing of the dominant CSF of the lower level with 3d^{5}4p^{2} and the dominant CSF of the upper level with 3d^{5}4s4p CSFs, primarily the latter. These additional CSFs have direct dipole interactions with the dominant CSF of the opposite parity and are the main contributors to the oscillator strengths. The calculated oscillator strengths are therefore very sensitive to the degree of this mixing.
log gf values for the multiplet 42 transitions.
Some comparative results for the three transitions are displayed in Table 2. Our work is in very close agreement with the experimental determinations of Meléndez & Barbuy (2009). The earlier calculations of Raassen & Uylings (1998, 1999) give gf values that are almost a factor of two lower than ours, and are somewhat similar to the values given previously by Kurucz (1988) and in later revisions of his database, though these are even lower. On the other hand, the current (2013) tabulation on the CD of Kurucz & Bell (1995) shows values that have been substantially revised; although these are closer to ours, the ratios of the three oscillator strengths (or differences between log gf values) do not match either those of other authors or of the present calculation. Our work therefore shows a clear preference for the results of Meléndez & Barbuy (2009), and the closeness of the agreement is encouraging, given the sensitivity to CI mixing in our work.
The full comparison between our results and the experimental results of Meléndez & Barbuy (2009) is shown in Table 3. The close agreement between theory and experiment seen for the multiplet 42 continues for many of the 142 transitions displayed in Table 3, even though the calculated values involve substantial CI mixing and often CI cancellation. For some transitions, the agreement between calculated and experimental oscillator strengths is within 10%, though for many transitions the results differ by around 25%, which is still good for such weak transitions. In any case for some transitions, this size of difference occurs between the results of Meléndez & Barbuy (2009) and of Sousa et al. (2008). For some, though, there is a difference of about a factor of two, while in two cases, there are greater discrepancies: at λ6433, the difference is two orders of magnitude; at λ6508, the difference is more than four orders of magnitude.
We briefly examine three specific transitions, to see why such a variation in difference occurs.
3.1. λ6416.91 Å
For this transition, our result is almost identical to that of Meléndez & Barbuy (2009). It corresponds to the transition 3d^{6}(^{3}D)4s ^{4}D_{5/2}–3d^{6}(^{5}D)4p ^{4}P, with the levels labelled by the dominant CSF in each case. With the assumption of orthogonal orbitals, the oscillator strength would be zero, because of the different resultant angular momenta of the two 3d^{6} cores. The actual nonzero result of our calculations arises through the mixing of 3d^{6}(^{5}D)4p ^{4}P with 3d^{6}(^{3}D)4p ^{4}P, even though the for the second CSF is only 0.02. This second CSF in the odd parity state has a direct dipole interaction with the dominant even parity CSF. All other such interactions are minor in comparison. The agreement with the experimental result suggests that our finetuning process is successful in obtaining mixing coefficients that are reasonably accurate.
3.2. λ4515.33 Å
This transition corresponds to 3d^{6}(^{3}F2)4s ^{4}F_{5/2}–3d^{6}(^{5}D)4p ^{4}F. Our oscillator strength is more than double the value of Meléndez & Barbuy (2009). The lower level exhibits strong mixing with 3d^{6}(^{3}F1)4s ^{4}F_{5/2} (the 3d cores can also be written in terms of seniority as F and F). This mixing gives contributions to the dipole matrix element that are additive, and if the strength of the mixing were in error, there would be little effect on the total. However, these contributions arise from the mixing of both 3d^{6}(^{3}F2)4p ^{4}F and 3d^{6}(^{3}F1)4p ^{4}F in the wave function for 3d^{6}(^{5}D)4p ^{4}F. Again, this mixing has an additive effect on the oscillator strength, which makes it less sensitive to errors in the CI mixing coefficients. Other mixing in the lower level with 3d^{6}(^{5}D)4s ^{4}D and 3d^{7}^{4}F has a smaller but again noncancelling effect on the oscillator strength. It would seem that the calculated oscillator strength could only be reduced through the CI mixings being too strong, though changes by a factor of two would be surprising.
3.3. λ6508.12 Å
This is the transition for which our oscillator strength is four orders of magnitude below the experimental value determined by Meléndez & Barbuy (2009). It corresponds to the transition 3d^{6}(^{5}D)4p ^{4}F–3d^{5}4s^{2}^{4}D_{3/2}. There are two orbitals that are different in the dominant CSFs of the two levels so that again the single configuration result for the oscillator strength would be zero: the nonzero result is obtained through CI mixing.
To compare the effects of mixing arising from different CSFs, it is convenient to define where l and u denote the lower and upper levels of the transition, ΔE is the transition energy, and Op is the appropriate dipole operator. Then the oscillator strength is given simply by For this particular transition, the main contributions to the dipole integral are displayed in Table 4.
There is severe cancellation within the final four interactions owing to the interaction between 3d^{6}(^{5}D)4p ^{4}F with 3d^{6}(^{5}D)4p ^{4}D, resulting in a net contribution of 0.009 to the total of 100F_{ij}. A similar degree of cancellation occurs amongst the first three interactions listed, with a combined contribution of − 0.049, resulting in a total of − 0.040 for these seven interactions, which is a factor of more than 20 smaller than the largest of the individual contributions. At this point, other contributions of a similar magnitude to the total of these seven make a significant contribution: with the seven alone, the oscillator strength would be 1.6 × 10^{7}, whereas when all contributions are included, the oscillator strength reduces still further to 2.1 × 10^{9}.
One possible source of error in the calculation is therefore the loss of significant figures arising from these cancellations. Moreover, the relevant mixing coefficients are themselves small. Their values are shown in Table 5.
Main contributions F_{ij} to the dipole integral in line λ6508.12 Å.
Significant CI mixing coefficients relevant to the line λ6508.12 Å.
The magnitude of at least one of the a_{i} in each of the key interactions is itself very small, and is perhaps open to change if more CSFs were to be added to the wave function expansions. In view of the extent of the cancellation effects, the calculated oscillator strength for this particular transition is much less reliable than those of other transitions in Table 3.
4. Conclusions
In the previous section we discussed three exemplars of transitions, all of which (as E1 transitions) might normally be expected to be strong but which in fact are weak, sometimes very weak. The type of accuracy that can be expected of theoretical results depends, not on the wavelength of the transition, but on the makeup of the wave functions of the levels involved in each transition and, in particular, on the angular momenta of the 3d^{n} cores.
The characteristic feature of transitions such as the 6416.91 Å line (Sect. 3.1) is that the 3d^{6} cores, while having different LS symmetries in the 4s and 4p states, are mixed mainly with CSFs containing the other symmetry. That is, the main mixing in the 3d^{6}(L_{1}S_{1})4s level is with 3d^{6}(L_{2}S_{2})4s, while the one for 3d^{6}(L_{2}S_{2})4p is with 3d^{6}(L_{1}S_{1})4p. The close agreement between our oscillator strengths and those of Meléndez & Barbuy (2009) encourages us to believe that, for such lines, we have obtained this mixing accurately and therefore also the oscillator strengths.
The characteristic feature of transitions such as 4515.33 Å (Sect. 3.2) is that, in addition, the core angular momenta L_{1}S_{1} and/or L_{2}S_{2} can occur with more than one seniority, and this introduces a further CI mixing that is usually very strong (in that case ^{3}F2 and ^{3}F1). Our experience is that such mixing is
difficult to achieve accurately. When this occurs, we would be more hesitant about the accuracy of the oscillator strengths.
Lines such as 6508.12 Å (Sect. 3.3) are characterised by a twoelectron change between the main CSFs of the two levels (3d4p → 4s^{2}). Frequently for such lines there is severe cancellation in the dipole matrix elements due to strong CI mixing, giving rise to calculated oscillator strengths that can be abnormally low (by several orders of magnitude). Then, different calculations may give very different results, pointing to a need to undertake more extensive and perhaps more focussed calculations for such lines.
In Table 3 we have not given the usual labels for the levels involved in the transitions. In the near future, we plan to submit for publication our full set of results, comprising oscillator strengths of all transitions between levels of the configurations 3d^{6}4s, 3d^{7}, 3d^{5}4s^{2}, 3d^{6}4p, 3d^{5}4s4p, and also 3d^{6}5p. The labels of the lower and upper levels of each transition will then be specified.
Acknowledgments
The authors would like to thank STFC, UK, for support under the Rolling Grant PP/D00103X/1.
References
 Corrégé, G., & Hibbert, A. 2005, ApJ, 627, L157 [NASA ADS] [CrossRef] [Google Scholar]
 Corrégé, G., & Hibbert, A. 2006, ApJ, 636, 1166 [NASA ADS] [CrossRef] [Google Scholar]
 Deb, N. C., & Hibbert, A. 2010a, ApJ, 711, L104 [NASA ADS] [CrossRef] [Google Scholar]
 Deb, N. C., & Hibbert, A. 2010b, A&A, 524, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deb, N. C., & Hibbert, A. 2011, A&A, 536, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deb, N. C., & Hibbert, A. 2012, J. Phys. Conf. Ser., 388, 152005 [NASA ADS] [CrossRef] [Google Scholar]
 Hibbert, A. 1975, Comp. Phys. Commun., 9, 141 [NASA ADS] [CrossRef] [Google Scholar]
 Hibbert, A., Glass, R., & Fischer, C. F. 1991, Comput. Phys. Commun., 64, 455 [NASA ADS] [CrossRef] [Google Scholar]
 Kroll, S., & Kock, M. 1987, A&AS, 67, 225 [NASA ADS] [Google Scholar]
 Kurucz, R. L. 1988, Trans. IAU, 28, 168 [Google Scholar]
 Kurucz, R. L., & Bell, B. 1995, Kurucz CDROM No. 23 (Cambridge, Mass.: Smithsonian Astrophysical Observatory) [Google Scholar]
 Meléndez, J., & Barbuy, B. 2009, A&A, 497, 611 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Raassen, A. J. J., & Uylings, P. H. M. 1998, J. Phys. B: Atom. Mol. Opt. Phys., 31, 3137 [NASA ADS] [CrossRef] [Google Scholar]
 Raassen, A. J. J., & Uylings, P. H. M. 1999, ftp://ftp.wins.uva.nl/pub/orth/iron/FeII.E1 [Google Scholar]
 Schnabel, R., SchultzJohanning, M., & Kock, M. 2004, A&A, 414, 1169 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sousa, S. G., Santos, N. C., Mayor, M., et al. 2008, A&A, 487, 373 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
All Tables
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.