Free Access
Issue
A&A
Volume 533, September 2011
Article Number A62
Number of page(s) 15
Section Atomic, molecular, and nuclear data
DOI https://doi.org/10.1051/0004-6361/201117471
Published online 29 August 2011

© ESO, 2011

1. Introduction

Recent measurements of neutron(n)-capture element (atomic number Z > 30) emission lines in ionized astrophysical nebulae have spurred laboratory astrophysics investigations of these species. Trans-iron elements were first detected in planetary nebulae (PNe) 35 years ago (Treffers et al. 1976), but it was not until nearly two decades later that their emission lines were identified (Péquignot & Baluteau 1994).

The detection of these elements marked a new era in the field of PNe, since the progenitor stars of these objects can produce trans-iron elements via slow n-capture nucleosynthesis (the s-process) during the asymptotic giant branch (AGB) phase of evolution. Abundance determinations of PNe complement spectroscopic investigations of AGB stars, enabling previously unexplored details of s-process nucleosynthesis to be revealed (e.g., see Sterling & Dinerstein 2008; Sterling & Witthoeft 2011). Indeed, many of the n-capture elements that are accessible in nebular spectra, such as Se, Kr, and Xe, cannot be detected in cool giants such as AGB stars (and in fact had not previously been studied in any of their sites of origin).

The pioneering work of Péquignot & Baluteau (1994) inspired other groups (Dinerstein 2001; Dinerstein & Geballe 2001; Sharpee et al. 2007) to identify previously unrecognized nebular emission lines as transitions of trans-iron elements. In spite of the low cosmic abundances of n-capture elements (Asplund et al. 2009), deep ultraviolet (Sterling et al. 2002; Sterling & Dinerstein 2003; Williams et al. 2008), optical (e.g., Hyung et al. 2001; Sharpee et al. 2003; Liu et al. 2004; Zhang et al. 2005; Sharpee et al. 2007; Sterling et al. 2009; Otsuka et al. 2010, 2011) and near-infrared spectroscopy (Likkel et al. 2006; Sterling et al. 2007; Sterling & Dinerstein 2008) of several Galactic and Local Group PNe revealed spectral features of these species. Abundances of these elements have now been estimated in approximately 100 PNe, allowing for an analysis of s-process elemental enrichments in a statistically meaningful sample of objects.

These detections are not limited to PNe. Indeed, n-capture element emission lines have been detected in a variety of other environments, including H II regions (Aspin 1994; Lumsden & Puxley 1996; Luhman et al. 1998; Baldwin et al. 2000; Puxley et al. 2000; Okumura et al. 2001; Blum & McGregor 2008; Roman-Lopes et al. 2009) and starburst galaxies (Vanzi et al. 2008). These detections demonstrate the promise of nebular spectroscopy as a tool to study n-capture nucleosynthesis in low-mass stars, the chemical evolution of trans-iron elements in the Universe, and the heavy-element nucleosynthetic histories of other galaxies.

Kr is one of the most readily detected n-capture elements in ionized nebulae, along with Se and Xe, with relatively strong transitions in both the optical (Péquignot & Baluteau 1994) and near-infrared (Dinerstein 2001). It has been detected in dozens of Galactic PNe and H II regions, as well as a handful of extragalactic nebulae (e.g., Wood et al. 2007; Otsuka et al. 2011).

Beyond the relative ease of detection, Kr is a sensitive probe of s-process nucleosynthesis. It can be highly enriched by the s-process, due in part to the 86Kr isotope, which has a magic number of neutrons and acts as a bottleneck in the s-process path. Isotopes with magic numbers of neutrons have very small n-capture cross sections and cause the peaks seen in the element-by-element s-process enrichment distribution near Z = 40, 56, and 82 (e.g., Busso et al. 1999; Burris et al. 2000; Sneden et al. 2008). Typically, s-process nucleosynthesis enriches Kr by a significantly larger factor than Se (Sterling & Dinerstein 2008; Karakas et al. 2009), which does not have a stable isotope with a magic number of neutrons. In addition, because it is a noble gas, Kr is not depleted into dust grains, and hence its gas-phase abundance is representative of the total elemental abundance.

However, nebular spectroscopy requires a foundation of atomic data to determine both ionic abundances and ionization balance solutions to correct for unobserved ions. Transition probabilities and effective collision strengths have been determined for many detected n-capture element ions (Biémont & Hansen 1986b,a, 1987; Biémont et al. 1988, 1995; Schöning 1997; Schöning & Butler 1998), enabling their ionic abundances to be derived. However, one of the primary hurdles to nebular studies of trans-iron elements is that the ionization equilibria of these species cannot be accurately solved due to the lack of photoionization (PI) and recombination data for nearly all n-capture element ions.

Motivated by this clear need for atomic data, we have embarked on a program using both theoretical and experimental methods to determine PI cross sections and rate coefficients for radiative recombination (RR), dielectronic recombination (DR), and charge transfer (CT). Collisional ionization and high-temperature DR are of negligible importance for ionization balance solutions in photoionized nebulae such as PNe and H II regions.

This paper is the second in a series presenting atomic data calculations for n-capture element ions. In the first (Sterling & Witthoeft 2011), multi-configuration Breit-Pauli (MCBP) distorted-wave PI cross sections and RR and DR rate coefficients were given for the first six ions of Se. We furnish CT rate coefficients for the first five ions of Ge, Se, Br, Kr, Rb, and Xe in the third paper (Sterling & Stancil 2011), computed in the Demkov and Landau-Zener approximations. In this study, we present MCBP distorted-wave PI cross sections and RR and DR data for the first six Kr ions, calculated with the atomic structure code AUTOSTRUCTURE (Badnell 1986, 1997, 2011). We do not investigate more highly-charged Kr species, since PN central stars and massive young stars (the ionizing sources of H II regions) are not sufficiently hot to significantly photoionize Kr more than six times (Kr6+ has an ionization potential of 111 eV, and PN central stars generally have effective temperatures less than 200 000 K; Napiwotzki 1999; Phillips 2003; Villaver et al. 2007).

Photoionization and recombination studies of low-charge Kr ions are few and far between, particularly in terms of theoretical calculations. On the other hand, neutral Kr has been the subject of numerous experimental (e.g., Henke et al. 1967; Lang & Watson 1975; Marr & West 1976; West & Marr 1976; Samson & Yin 1989; Samson et al. 1991; Sorokin et al. 2000, and references therein) and theoretical (e.g., Amusia et al. 1971; Miller et al. 1977; Parpia et al. 1984; Tulkki et al. 1992) PI studies. However, neutral Kr is a trace species in PNe and H II regions, and therefore its PI properties are not expected to play a significant role in Kr ionization balance solutions.

In terms of recombination, little is known about the ions of interest in our work. Mattioli et al. (2006) compiled RR and DR rate coefficients for all Kr ions, but the data for low-charge ions are rough estimates. For these species, DR rate coefficients were computed using the Burgess formula (Burgess 1965), as modified by Merts et al. (1976), which is appropriate for collisionally-ionized plasmas (indeed, the work of Mattioli et al. was inspired by the utility of Kr as a coolant for magnetically-confined fusion plasmas), but is known to badly underestimate the rate coefficient at the low temperatures of photoionized plasmas. Their RR rate coefficients for low-charge Kr ions are approximate, with RR into the valence and (hydrogenic) excited shells treated separately.

Experimental measurements of absolute PI cross sections provide the most useful comparisons to our calculations. Concurrent to our theoretical investigations, we have experimentally measured PI cross sections for Se and Xe ions (Esteves et al. 2009, 2010; Sterling et al. 2011; Esteves et al. 2011 and in prep.) at the Advanced Light Source synchrotron radiation facility, for the purpose of constraining our calculations. PI investigations were previously carried out for Kr3+ to Kr5+ (Lu et al. 2006a,b; Lu 2006), and more recently for Kr+ (Bizau et al. 2011). We compare our computed PI cross sections to these measurements to benchmark and gauge the accuracy of our direct PI cross sections.

Because no comprehensive studies of the PI and recombination properties of low-charge Kr ions in the photoionized regime existed prior to our study, we provide estimates for the internal uncertainties of our computed atomic data. These uncertainties can affect ionization equilibrium calculations, and hence elemental abundance determinations. Beyond comparisons to experimental PI cross-section measurements, we computed data using three different configuration-interaction (CI) expansions for each Kr ion to gauge the sensitivity of the results to CI effects. Moreover, we tested the effects of other parameters and assumptions in our calculations (e.g., whether the radial wavefunctions were orthogonalized, the radial scaling parameters for continuum orbitals, etc.) on the resulting data.

In the remainder of this paper, we present the computed PI cross sections and RR and DR rate coefficients for low-charge Kr ions. In Sect. 2, we briefly describe the calculations and methodology for each of these processes, and detail the calculated electronic structure for each ion. The resulting data are presented in Sect. 3, along with estimated uncertainties and comparisons to previous studies. We summarize our investigation in Sect. 4.

Table 1

Configuration-interaction expansions used for each Kr ion.

2. Calculations and methodology

The electronic structure, distorted-wave PI cross sections, and RR and DR rate coefficients have been computed for the first six ions of Kr using the atomic structure code AUTOSTRUCTURE (v23.25; Badnell 1986, 1997, 2011). We refer the reader to Badnell et al. (2003) and Badnell (2006b) for a detailed account of the theoretical background for these calculations.

2.1. Electronic structure

We constructed three separate CI expansions (which we name “small”, “medium”, and “large”) for each ion, in order to test the sensitivity of our results to the adopted CI expansion. Based on comparisons to level energies and ionization potentials in NIST (Ralchenko et al. 2008), the medium configuration sets provide the best compromise between accuracy and computational efficiency, and the results in Sect. 3 were calculated with the medium CI expansions. In Table 1, the configurations and numbers of levels for the three CI expansions of each ion are given.

The electronic structure of each ion was computed with Thomas-Fermi-Dirac-Amaldi potentials, in intermediate coupling with κ-averaged relativistic wavefunctions (Cowan & Griffin 1976). The radial wavefunctions were Schmidt orthogonalized, but calculations for each process were also conducted without orthogonalization (see Sect. 3). We determined orbital radial scaling parameters (Table 2) by optimizing on the average of LS term energies in order to best reproduce experimental energies and ionization potentials tabulated by NIST. We generally found best agreement with experimental values by optimizing the scaling parameter for each orbital on terms in the lowest-energy configuration containing electrons in the that orbital. In a few cases, other optimization strategies were required to improve the accuracies of ground configuration level energies and the ionization potential (e.g., see the discussion of the Kr+ structure in Sect. 3.1.2). The computed energies for the medium and large configuration sets typically agree with experimental determinations to within 3−4% (see Table 3), while the small CI expansion energies are in slightly worse agreement. The energies for some levels are more discrepant, especially 1D2 and 2D3/2,5/2 in the ground configurations of Kr2 + −Kr4 +  (the disagreement is as high as 22% for the medium CI expansion of Kr3+ 4s2 4p4 2D3/2).

Table 2

Radial scaling parameters used for CI expansions of each ion.

Calculated ionization potentials (Table 3) agree with NIST values to within 2% for the medium and large configuration sets, with the exceptions of Kr2+ and Kr3+. The ionization potential of Kr2+ is strongly sensitive to the radial scaling parameter of the 4p orbital. Decreasing the 4p scaling parameter leads to an improved ionization potential, but the accuracy of the ground term energies rapidly worsens. If the 4p scaling parameter optimization is conducted so as to minimize its value, the computed ionization potential would improve from an accuracy of 7.1% to 6.2%, while the errors in the 3P metastable state energies would more than triple. We therefore chose to optimize the 4p scaling parameter to best reproduce the experimental energy level structure, at the expense of a slightly less accurate ionization potential. This choice does not significantly affect our results – the direct Kr2+ PI cross section computed with the two values of the 4p scaling parameter differs by less than 2% near the ground state threshold. The Kr3+ ionization potential is quite sensitive to the 4s and 4p scaling parameters, improving when those values are decreased. However, in this case it was not possible to optimize the 4s and 4p scaling parameters so as to reduce their values below those listed in Table 2, and hence we were unable to further improve the accuracy of the ionization potential.

Table 3

Comparison of selected calculated and experimental energies (in Rydbergs).

In Table 4, we compare the Einstein A-coefficients we computed with previously determined values. For forbidden transitions in the ground configuration, our A-coefficients agree to within 50% or better with previous results, with the exception of transitions involving levels whose energy we were not able to compute to within 3−4% accuracy (e.g., 1D2 for Kr2+ and Kr4+, and 2D3/2,5/2 for Kr3+). For permitted transitions, the disagreement is generally larger (especially for neutral and singly-ionized Kr).

Table 4

Calculated Einstein A-coefficients (in s-1) compared to literature values.

2.2. Photoionization, radiative recombination, and dielectronic recombination

In low-density nebulae such as PNe and H II regions, nearly all atomic ions reside in ground configuration states. Therefore, we computed MCBP distorted-wave PI cross sections for each level of the Kr0 to Kr5+ ground configurations from the ionization threshold up to 100 Rydbergs. Likewise, RR and DR rate coefficients were computed for states within the ground configuration of the target ions (i.e., before recombination) Kr+ to Kr6+ over the temperature range (101−107)z2 K. The RR and DR rate coefficients were determined from the direct and resonant portions of the PI cross sections (respectively), using detailed balance. For DR, only Δn = 0 core excitations were considered, since these dominate Δn ≥ 1 core excitations at photoionized plasma temperatures. However, for DR forming Kr0 and Kr+, whose lowest energy 4s2 4pk 5s and (in the case of Kr0) 4s2 4pk 5p levels (k = 4 and 5 for Kr+ and Kr0, respectively) lie below those of 4s2 4pk 4d, we allowed Δn = 1 excitations into the 5s and 5p orbitals. We ignored the small contributions of core excitations into 5d, 5f, and 5g orbitals to keep the size of the calculations tractable.

The CI expansions listed in Table 1 were used in these calculations, with the (N + 1)-electron ion scaling parameters from Table 2. Our tests using target (N-electron) scaling parameters led to substantial disagreement between our calculated PI cross sections and experimental measurements in all cases. The scaling parameters of s, p, and d continuum orbitals are assumed to be the same as those of the highest principal quantum number bound orbitals with the same orbital angular momentum, and equal to unity otherwise. For RR, the (N + 1)-electron ion’s CI expansion was augmented to include one-electron additions to the target’s ground and mixing configurations. In the case of DR, all relevant core excitations were added to the target configuration set, and the (N + 1)-electron ion’s CI expansion was expanded to include one-electron additions to target states to allow Rydberg electrons to radiate into the core.

AUTOSTRUCTURE automatically switches from length to velocity to acceleration gauge as the photon energy increases. However, because we employed κ-averaged relativistic radial wavefunctions, acceleration gauge could not be used and it was necessary to force velocity gauge at high energies. In Sects. 3.1.1 and 3.2.1, we show that this does not affect our calculations at the energies and temperatures of interest.

RR and DR were treated separately using the independent processes approximation (Pindzola et al. 1992). In our RR calculations, we used a maximum principal quantum number of 1000, and a maximum orbital angular momentum of l = 4 for the first two ions and l = 3 for more highly charged species (hydrogenic radial integrals from the recurrence relations of Burgess (1964) were used for higher values of l). For DR, we used a maximum n = 1000 and l = 200, with l ≥ 9 treated hydrogenically except for the Kr5+ and Kr6+ targets, for which l ≥ 10 and 11 (respectively) were computed hydrogenically in order to allow the calculation to converge. To improve the accuracy of our DR calculations, we utilized experimental target energies whenever available, with theoretical level splittings for terms with incomplete energy information in NIST (Badnell 2006a).

3. Results

3.1. Photoionization cross sections

In Fig. 1, we illustrate the ground state Kr0 to Kr5+ PI cross sections near their ionization thresholds. The AUTOSTRUCTURE output was processed with the codes ADASPI (for direct PI) and ADASPE (for photoexcitation-autoionization), which produce the ADAS adf39 and adf38 output files, respectively (Summers 2005). For ease of presentation, we used the post-processing code xpeppi to add the direct and resonant portions of the PI cross sections.

thumbnail Fig. 1

Ground state PI cross sections for Kr0−Kr5+ near their ionization thresholds.

The PI cross sections are available at the CDS, with two data files for each ion. Files labeled “krq+_XPITOT.dat” (where q is the charge of the (N + 1)-electron ion) tabulate the direct-plus-resonant PI cross sections in five columns. Photon energies relative to the ground and initial state are given in the first two columns (respectively), while the next two columns present the energies of the photo-electron relative to the ground and initial states. All energies are in Rydbergs. The cross section in Mb is given in the final column. For ions with multiple ground configuration levels, the cross sections for different levels are delimited by comment lines (indicated by a “#” in the first column) followed by three numbers. The first is the energy of the initial level relative to the ground state (or the total energy, in the case of the ground state itself), the second is the total energy of the final target state, and the last gives the energy order of the initial level (1 for the ground state, 2 for the first excited state, etc.).

Direct PI cross sections, without photoexcitation-autoionization resonances, are presented in files labeled “krq+_XDPITOT.dat”. These files contain only two columns, the first of which is the photon energy in Rydbergs, and the second the cross section in Mb. Cross sections for different states are delimited as in the krq+_XPITOT.dat files.

3.1.1. Estimation of uncertainties

We estimate uncertainties by testing the sensitivity of our results to the orbital radial scaling parameters, the adopted CI expansion, and to various assumptions in our calculations. We also compare our results to experimental absolute PI cross-section measurements (Sect. 3.1.2).

The PI cross sections are most sensitive to the adopted CI expansion. In Figs. 2 and 3, we show the Kr2+ and Kr3+ cross sections computed with each of the three configuration sets (top three panels), while the bottom panel compares the direct cross sections from each calculation. Most notably, the resonance structure, both in terms of strengths and positions, varies considerably among the different CI expansions. This illustrates the challenge of modeling photoexcitation-autoionization processes in the distorted-wave approximation. It should be noted that any theoretical treatment leads to uncertainties in PI resonance structure, so much so that codes used to numerically simulate astrophysical nebulae often ignore these resonances (Ferland et al. 1998) or use cross sections convolved with a broad Gaussian to smooth out uncertainties in resonance positions (Kallman & Bautista 2001).

thumbnail Fig. 2

Comparison of Kr2+ ground state PI cross sections computed for each of the three CI expansions (as indicated in the top three panels). The direct PI cross sections are compared in the bottom panel.

thumbnail Fig. 3

Comparison of Kr3+ ground state PI cross sections computed for each of the three CI expansions (as indicated in the top three panels). The direct PI cross sections are compared in the bottom panel.

We therefore focus on estimating uncertainties in the direct PI cross sections, particularly near the ground state ionization threshold. The Kr2+ cross section computed with the small configuration set is larger than those from the other CI expansions by approximately 25% near the ground state threshold (Fig. 2), while the Kr3+ direct cross sections agree to within 5% (Fig. 3). Likewise, the direct cross section computed with the medium Kr4+ and Kr5+ CI expansions agree with those from the small and large configuration sets to within 5−10%. The disagreement is largest (up to  ~40%) for neutral and singly ionized Kr.

Other tests revealed smaller uncertainties in the direct PI cross sections. For example, setting the scaling parameters to unity for s, p, and d continuum orbitals (rather than the value of the highest principal quantum number bound orbital with the same orbital angular momentum) led to differences of 10% near the ground state threshold of neutral Kr, and less than 5% for other Kr ions. If the PI cross sections are computed without Schmidt orthogonalizing the radial orbitals, the direct cross sections differ by 5% or less in all cases. To test whether our inability to switch to acceleration gauge at high energies (due to the use of κ-averaged relativistic orbitals) leads to inaccuracies, for each ion we ran a PI calculation solely in the length gauge, and another solely in the velocity gauge. The cross sections showed negligible difference (≪1%), indicating that our PI calculations are not affected by forcing velocity gauge up to 100 Ryd.

Based on these tests and comparisons to experimental measurements, we estimate the uncertainties in our direct PI cross sections to be  ~30−50% for most Kr ions. However, comparison with experimental measurements (Sect. 3.1.2) indicate that the uncertainties in the Kr0 and Kr3+ PI cross sections may be as large as a factor of two.

3.1.2. Comparison to experiment

It is critical to benchmark atomic data calculations against experimental measurements whenever possible. Of the six Kr species considered, experimental absolute PI cross sections have been measured for all but Kr2+.

The measurements were conducted at a variety of synchrotron radiation facilities, and in the case of ionized species utilized the merged beams method. Kr3+ through Kr5+ PI cross sections were measured at the Advanced Light Source synchrotron radiation facility at the Lawrence Berkeley National Laboratory in California (Lu et al. 2006a,b; Lu 2006), while that of Kr+ was recently measured at the ASTRID storage ring at the University of Aarhus in Denmark (Bizau et al. 2011). These measurements were complicated by the population of metastable states in the primary ion beams, and hence the measured cross sections are a linear combination of those from the ground and metastable states. However, Bizau et al. (2011) used an ion trap to significantly limit contributions from the 2P1/2 metastable state, and measured the Kr+ ground state PI cross section over a limited energy range. In the case of neutral Kr, we compare our results to the PI cross sections recommended by Richter et al. (2003), which were determined from measurements of the PI and electron-impact ionization cross section ratio and normalized to the double ion chamber PI measurements of Samson & Yin (1989) and Samson et al. (1991).

thumbnail Fig. 4

Comparison of calculated Kr+ PI cross section to the experimental measurements (dotted line) of Bizau et al. (2011). The dashed line shows our original result before comparison to experiment, while the solid line shows the theoretical cross section after optimization to better reproduce the Bizau et al. (2011) results. In each case, the theoretical curves are statistically-weighted linear combinations of the cross sections from each of the two states in the ground configuration. The calculated cross sections were convolved with a Gaussian of 30 meV to approximate the average resolution of the Bizau et al. (2011) measurements.

We substantially improved our calculated Kr+ PI cross section through comparison with the experimental data of Bizau et al. (2011). In Fig. 4, the experimental measurements are shown as a dotted line, and the dashed line represents our original calculation – statistically weighted and convolved with a Gaussian of 30 meV (the average experimental photon energy resolution near the ground state threshold). The calculated cross section was originally a factor of 2.5 smaller than the experimental one. We added and removed several configurations from the target and (N + 1)-electron ion CI expansions, none of which improved agreement with experiment.

Tests showed the PI cross section to be sensitive to the 4p, 4d, 5p, and continuum d orbital radial scaling parameters, generally increasing (improving) with smaller values for those parameters. We therefore re-optimized the Kr+ np and nd scaling parameters so as to minimize their values without significantly harming the agreement with the experimental level energies and ionization potential. In our original structure calculation, the 4p scaling parameter was optimized on terms in the ground, 4s 4p6, 4s2 4p4 4d, and 4s2 4p4 5s configurations and the 5p scaling parameter was optimized on 4s2 4p3 5p2 terms (this configuration was later removed from the CI expansion) in order to reproduce the energy of the first excited state (2P1/2). The 4d scaling parameter was optimized on 4s2 4p3 4d2 terms, to best reproduce the ionization potential. In our new optimization, we instead optimized the 4p radial scaling parameter only on the ground configuration term, and those for 4d and 5p on terms in the lowest-energy configurations containing an electron in those orbitals (4s2 4p4 4d and 4s2 4p4 5p, respectively). This new optimization reduced the accuracy of the first excited state energy from 0.8% to 11% (higher-energy levels are not appreciably affected), and that of the ionization potential from 0.5% to 1.5%. However, it produced a significantly larger direct PI cross section (solid line in Fig. 4), within 40% of the experimental data. This level of improvement demonstrates how useful absolute experimental PI cross sections can be to theoretical determinations.

We compare our PI cross sections to experimental measurements for other Kr ions in Fig. 5. The depicted theoretical data are linear combinations of cross sections for states in the ground configurations that were visually determined to best reproduce the experimental direct PI cross section and, when possible, families of resonances. The cross sections were weighted by the factors (3P0, 3P1, 3P2, 1D2, 1S0) = (0.32, 0.15, 0.20, 0.3, 0.03) for Kr4+, (2P1/2, 2P3/2) = (0.4, 0.6) for Kr5+, and by statistical weights in the case of Kr3+. The calculated cross sections were convolved with a Gaussian with FWHM equal to the experimental resolution (~1 eV for Kr0, and 0.1 eV for the other three ions).

thumbnail Fig. 5

Comparison of calculated PI cross sections (solid lines) with experimental data (dotted lines) for Kr0 (Richter et al. 2003), Kr3+ (Lu et al. 2006b), Kr4+ (Lu 2006), and Kr5+ (Lu et al. 2006a). The theoretical curves have been convolved with a Gaussian of 1 eV for Kr0 and 100 meV for the other species, to match the experimental resolution. The theoretical data are linear combinations of PI cross sections from each state in the ground configuration, weighted according to statistical weights in the case of Kr3+, and by the factors (3P0, 3P1, 3P2, 1D2, 1S0) = (0.32, 0.15, 0.2, 0.3, 0.03) for Kr4+ and (2P1/2, 2P3/2) = (0.4, 0.6) for Kr5+.

In the cases of Kr4+ and Kr5+, the direct cross sections generally agree within the experimental noise. The theoretical cross section is slightly too large for Kr4+ and perhaps slightly small for Kr5+, but over the full range of energies they agree to within the experimental uncertainties of  ~25%.

The agreement is markedly worse for Kr3+, for which the theoretical cross section is larger than that measured experimentally by approximately a factor of two. Changing the metastable fractions does not significantly improve the agreement with the experimental measurements, since the calculated direct PI cross sections of the ground and metastable states have comparable magnitudes and hence any linear combination of their cross sections will result in poor agreement with the experimental data. As seen in Fig. 3, the direct Kr3+ PI cross section is not strongly sensitive to the CI expansion. We verified this by adding and removing several configurations in the target and (N + 1)-electron ion medium configuration sets, none of which significantly affected the direct PI cross section. Near the ground state threshold, the direct cross section is sensitive to the values of the 4d and continuum s and d scaling parameters, decreasing (coming into better agreement with experiment) as the values of those scaling parameters increase. However, our re-optimization designed to increase the 4d (and hence continuum d) scaling parameter – by optimizing on 4s2 4p 4d2 instead of 4s2 4p2 4d terms – did not effect changes in the direct PI cross section of more than 10−15%. Moreoever, while we were able to reduce the direct PI cross section near threshold, the cross section above 4.5 Ryd increased (worsening agreement with experiment) by 10%. Therefore, we did not alter our calculation, and concede that the computed Kr3+ PI cross section is uncertain by approximately a factor of two.

Likewise, for Kr0 the computed cross section is too large at threshold by  ~70% compared to experiment. The cross section falls below experiment near 1.35 Ryd, and the disagreement increases to a factor of two until  ~2.5 Ryd, where a threshold corresponding to Kr+ 4s2 4p4 5p is seen in the theoretical cross section. The PI cross section at threshold is very sensitive to the value of the 4p and 5p scaling parameters, but altering their values in the sense required to improve the direct cross section badly damages energies of excited levels and the ionization threshold. Like Kr3+, the calculated PI cross section is uncertain by a factor of two, but we do not expect this to harm ionization balance solutions significantly since neutral Kr is a trace species in ionized astrophysical nebulae.

3.2. Radiative recombination rate coefficients

Total and partial final-state resolved RR rate coefficients were derived from the direct portion of the PI cross sections using detailed balance. We used the ADASRR code to process the AUTOSTRUCTURE output and produce ADAS adf48 files (available at the CDS), which are similar in structure to adf09 files for DR (Summers 2005). At the end of the adf48 files, total RR rate coefficients for all ground configuration states of the target ion are printed. Information preceding those rates pertain to partial RR rate coefficients, which are fully J-resolved for n ≤ 8 levels, bundled-nl for n ≤ 10 states, and bundled − n for all higher levels (see Badnell 2006a).

The rate coefficients were fit with the analytical formula (Verner & Ferland 1996; Badnell 2006b) αRR(T)=A×[T/T0(1+T/T0)1B0(1+T/T1)1+B0]-1,\begin{equation} \alpha_{{\rm RR}}(T) = A \times \left[ \sqrt{T/T_0}(1 + \sqrt{T/T_0})^{1-B_0}(1 + \sqrt{T/T_1})^{1+B_0} \right]^{-1}, \label{rr1} \end{equation}(1)where B0=B+Cexp(T2/T).\begin{equation} B_0 = B + C \exp (-T_2/T). \label{rr2} \end{equation}(2)The fit coefficients T and T0,1,2 are in temperature units, A and αRR(T) are in cm3 s-1, and B and C are dimensionless.

As was the case for Se ions (Sterling & Witthoeft 2011), we found that Eqs. (1) and (2) do not always provide satisfactory fits to the Kr RR rate coefficients over the full temperature range (101−107)z2 K. The reason is that the slopes of the rate coefficients can increase significantly at high temperatures. These “bumps” are caused by electron capture into low-energy levels, which are not as effectively screened from the nucleus as higher levels. Equations (1) and (2) can be used to fit the low-temperature behavior of the RR rate coefficients, but a different analytical form is often necessary to fit the high-temperature behavior. For these high-temperature fits, we have chosen to utilize the formula typically used to fit DR rate coefficients (e.g., Zatasrinny et al. 2003): αRR(T)=1T3/2i=1nciexp(Ei/T),\begin{equation} \alpha_{{\rm RR}}(T) = \frac{1}{T^{3/2}} \sum_{i=1}^{n} c_i \exp (-E_i/T), \label{dreq} \end{equation}(3)where T and Ei are in temperature units, ci are in cm3 s-1 K3/2, and αRR(T) is in cm3 s-1. The value for n in Eq. (3) ranges from 5 to 7, depending on the target state.

The fit coefficients for the RR rate coefficients are presented in Table 5 for the low-temperature fits (using Eqs. (1) and (2)) and Table 6 for the high-temperature fits (using Eq. (3)). Table 5 also lists the maximum temperature at which the low-temperature fits are valid for each state; if Tmax is not given, then the low-temperature fits are valid over the entire temperature range (101−107)z2 K. Unless noted otherwise, the fits are accurate to within 5% at low temperatures (Table 5) and 2% at high temperatures (Table 6), and respectively reproduce the correct low-temperature and high-temperature asymptotic behaviors of the rate coefficients.

Table 5

Fit coefficients for radiative recombination rate coefficients at low temperature (see Eqs. (1) and (2)).

Table 6

Fit coefficients for radiative recombination rate coefficients at high temperatures, valid from Tmax (see Table 5) to 107z2 K. See Eq. (3).

We illustrate the RR rate coefficients in Fig. 6. For comparison, we plot the approximate RR rate coefficients of Mattioli et al. (2006) as gray curves. It can be seen that the Mattioli et al. (2006) rate coefficients are at least three orders of magnitude too small over the full range of temperatures. As stated in their paper, their RR rate coefficients are rough estimates, with valence shell states treated separately from excited shells (which were taken to be hydrogenic). In the absence of detailed calculations such an approach is reasonable, but the comparison to our results illustrates the dangers of such approximations for heavy element ions.

thumbnail Fig. 6

Radiative recombination rate coefficients (black lines) as a function of temperature for ground state target ions Kr+ (solid curve), Kr2+ (dotted curve), Kr3+ (short-dashed curve), Kr4+ (dash-dot-dash curve), Kr5+ (long-dashed curve), and Kr6+ (dash-dot-dot-dot-dash curve). For comparison, the approximate rate coefficients of Mattioli et al. (2006) are depicted as gray curves, with the line styles corresponding to the same ions.

3.2.1. Estimation of uncertainties

We have estimated the internal uncertainties of our calculated RR rate coefficients near T = 104 K, the canonical temperature of photoionized astrophysical nebulae. The adopted CI expansion is again the largest source of uncertainty. The small and large configuration sets produce RR rate coefficients different from those of the medium CI expansion by 40−45% for the Kr+ target, and by less than 5−7% for other species.

We also computed the rate coefficients using a larger number of interpolation points per decade (five, instead of three), using scaling parameters for the target ion (rather than those of the (N + 1)-electron ion), and without forcing Schmidt orthogonalization of the radial wavefunctions. These tests resulted in uncertainties no larger than 5%, except for the two lowest-charge species. For the Kr2+ target, using target ion scaling parameters resulted in rate coefficients differing by up to 10% for metastable states. For Kr+ forming Kr0, calculations with target ion scaling parameters led to differences of 20−30%, while the use of non-orthogonal orbitals resulted in differences of 7% for the metastable state rate coefficient.

Finally, the κ-averaged relativistic orbitals used in our calculations prevent AUTOSTRUCTURE from switching to acceleration gauge at high energies. Thus, the calculated RR rate coefficients at high temperatures (which correspond to PI cross sections at high energies) may be inaccurate. To determine the temperatures at which such inaccuracies arise, we computed RR rate coefficients using PI cross sections that were terminated at the maximal velocity gauge energy (above which it is necessary to force velocity gauge in AUTOSTRUCTURE). The rate coefficients calculated in this way exhibited negligible differences with our presented results until T > 105 for the first two ionization stages, or  >106 K for higher ions. This demonstrates that forcing velocity gauge at high energies does not affect our results except at temperatures well above the photoionized regime.

Based on these tests, we estimate the internal uncertainties of the computed RR rate coefficients to be  ≤10% except for the Kr+ target, which is uncertain by 50−60%.

3.3. Dielectronic recombination rate coefficients

We determined total and partial final-state resolved DR rate coefficients from the resonant portion of the PI cross sections, with the aid of the ADASDR post-processing code. The resulting ADAS adf09 output files (see Summers 2005) are available at the CDS.

The rate coefficients were fit with Eq. (3), with n between 3 and 7, and the fit coefficients are listed in Table 7. The non-linear least squares fit algorithm produced fits accurate to within 5% over the temperature range (101−107)z2 K for most target states, with exceptions noted in Table 7. These fits also correctly produce the asymptotic behavior of the rate coefficients outside of this temperature range.

Table 7

Fit coefficients for dielectronic recombination rate coefficients; see Eq. (3).

The DR rate coefficients are plotted in Fig. 7. The rate coefficients for different ions are wildly dissimilar at low temperatures due to differences in the structure of near-threshold autoionizing resonances, but show similar behavior at high temperatures.

Mattioli et al. (2006) also computed DR rate coefficients for low-charge Kr ions (shown in Fig. 7 as gray curves), in order to properly model its cooling rate coefficient in magnetically-confined fusion plasmas. These calculations utilized the Burgess formula (Burgess 1965) as modified by Merts et al. (1976), summed over Δn = 0 and 1 resonances. This formula is well-known to inadequately describe low-temperature DR (e.g., Savin et al. 1999, 2003), since the captured electrons are assumed to have thermal energies comparable to the excitation energies of the first resonance transitions of the target ion (Mattioli et al. 2006). While this approximation is adequate for collisionally-ionized plasmas, it is not sufficient for photoionized plasmas, in which the electrons are captured into low-energy autoionizing states. Given this, it is no surprise that the DR rate coefficients from the Burgess formula (Mattioli et al. 2006) are vanishingly small below 105 K. At high temperatures (≥a few times 106 K), the rate coefficients from the Burgess formula exceed those from our calculations, likely because we ignored Δn = 1 core excitations (except in limited cases for Kr+ and Kr2+).

thumbnail Fig. 7

Dielectronic recombination rate coefficients (black lines) as a function of temperature for ground state target ions Kr+ (solid curve), Kr2+ (dotted curve), Kr3+ (short-dashed curve), Kr4+ (dash-dot-dash curve), Kr5+ (long-dashed curve), and Kr6+ (dash-dot-dot-dot-dash curve). DR rate coefficients computed with the Burgess formula over Δn = 0 and Δn = 1 resonances, from Mattioli et al. (2006), are depicted as gray curves, with the line styles corresponding to the same ions.

Figure 8 compares the DR (dashed lines), RR (dotted lines), and total recombination rate coefficients (solid lines) for each Kr ion we considered. As is true of low-charge Se ions (Sterling & Witthoeft 2011), DR dominates RR for Kr ions near 104 K, with the DR rate coefficient exceeding that of RR by factors of 2−3 up to two orders of magnitude. RR is (occasionally) more important than DR only at low temperatures (≤1−2  ×  103 K). Only in the case of Kr+ is RR the most important recombination mechanism near  ~104 K. These results underline the importance of DR in the ionization balance of trans-iron elements.

thumbnail Fig. 8

Comparison of the radiative (dotted curves), dielectronic (dashed curves), and total recombination rate coefficients (solid curves) for the ground state of the first six Kr ions.

3.3.1. Estimation of uncertainties

DR calculations are inherently more uncertain than those of PI or RR, since the unknown energies of low-lying autoionizing levels usually dominate other sources of uncertainty. The energies of these near-threshold autoionizing states have not been experimentally determined for elements beyond the second row of the Periodic Table, and cannot be determined with sufficient accuracy by even the most sophisticated theoretical methods. To help assuage this uncertainty, we have utilized experimental target energies in all of our DR calculations.

The uncertainty in the DR rate coefficients can be estimated by shifting the continuum relative to near-threshold resonances. We accomplished this with the ADASDR post-processing code, determining the magnitude of the shift to be the maximum disparity between calculated and experimental target energies, and the direction of the shift by whether the calculated energies were too large or too small.

We find that the uncertainties in the rate coefficients near 104 K for the ground (metastable) states are 75% (factor of 20) for Kr+, a factor of 3 (10% to an order of magnitude) for Kr2+, 10% (15−35%) for Kr3+, 70% (50−70%) for Kr4+, 50% (25%) for Kr5+, and a factor of 2 for Kr6+. In contrast, DR calculations with the small CI expansions1 led to rate coefficients different from those computed with the medium configuration sets by a factor of 3 for Kr+, 60% for Kr2+, a factor of 2 for Kr5+, and 10−30% for other Kr ions.

With the exception of Kr3+, for which the uncertainty is probably an underestimate, the DR rate coefficients are uncertain by factors of 2−3 for the ground states of these ions. These uncertainties are considerably larger than those typically found for our PI cross sections (30−50%) or RR (<10%). Given that DR is a much more important recombination mechanism than RR for most Kr ions at photoionized plasma temperatures (Fig. 8), these uncertainties are significant. The best hope for accurate low-temperature DR rate coefficients for low-charge ions lies with experimental measurements made with heavy ion storage rings (Schippers et al. 2010), particularly the forthcoming Cryogenic Storage Ring facility (Wolf et al. 2006). Our planned ionization balance calculations for Se and Kr will take into account uncertainties in the atomic data, demonstrating which Se and Kr ions most critically require experimental DR rate coefficients to enable the most accurate possible nebular abundance determinations.

4. Summary

We have presented MCBP calculations of distorted-wave PI cross sections over the photon energy range 0−100 Ryd, and radiative and dielectronic recombination rate coefficients at temperatures (101−107)z2 K for the first six Kr ions. The results of all of our calculations are available at the CDS, as well as at http://www.pa.msu.edu/astro/atomicdata/kr_data.tar.gz. The RR and DR rate coefficients have been fit to analytical functions, with fit coefficients given in Tables 57.

Kr has been detected in the spectra of planetary nebulae and H II regions, and its abundance is a sensitive tracer of enrichments by n-capture nucleosynthesis. But in order to derive elemental Kr abundances from those of the observed ions, unseen ionization stages must be accounted for via ionization equilibrium solutions. Our atomic data determinations, along with charge transfer rate coefficients presented in the third paper of this series (Sterling & Stancil 2011), enable much more accurate Kr ionization balance computations than previously possible.

This study represents the first comprehensive investigation of the PI and recombination properties of low-charge Kr ions specifically focused on the photoionized plasma regime. As such, we have taken significant efforts to provide realistic internal uncertainties for our results, both through comparison to experimental PI measurements, and by testing the sensitivity of our data to changes in the CI expansion and orbital radial scaling parameter values.

We find typical uncertainties of 30−50% for our computed PI cross sections, though comparison to experimental measurements indicates that the Kr0 and Kr3+ cross sections have larger uncertainties of roughly a factor of two near their ionization thresholds. We note that our Kr+ PI calculations were significantly improved after comparison with experimental measurements, from accuracies of a factor of 2.5 near the ground state threshold to within  ~40%. This demonstrates the utility and importance of experimental absolute PI cross-section measurements for constraining and guiding theoretical calculations.

Our RR rate coefficients are more precisely determined, with uncertainties less than 10% except for the Kr+ target (50−60%). On the other hand, the unknown energies of low-lying autoionizing states lead to significantly larger uncertainties in DR rate coefficients, typically factors of 2−3. These uncertainties are magnified by the fact that DR is the dominant form of recombination for low-charge Kr ions at the temperatures of photoionized astrophysical nebulae, with the rate coefficients exceeding those of RR by as much as two orders of magnitude. The dominance of DR over RR in photoionized nebulae is common for heavy element ions (e.g., Badnell 2006a; Altun et al. 2007; Nikolić et al. 2010; Sterling & Witthoeft 2011).

The estimated uncertainties in these atomic data will play an important role in our forthcoming study of the Se and Kr ionization balance in planetary nebulae. We will numerically simulate the Se and Kr ionization balance for a wide range of nebular parameters to derive new analytical corrections for unobserved ions. Using our estimated uncertainties in the PI cross sections and RR and DR rate coefficients, we will run Monte Carlo simulations to test how sensitive elemental abundance determinations are to these uncertainties, thereby illuminating the atomic processes and ionic systems that require additional theoretical and/or experimental investigation.


1

We did not attempt to compute DR rate coefficients with the large CI expansions due to the prohibitive computional requirements.

Acknowledgments

I acknowledge support from an NSF Astronomy and Astrophysics Postdoctoral Fellowship under award AST-0901432 and from NASA grant 06-APRA206-0049. I am grateful to R. Phaneuf, who provided experimental data for Kr3+, Kr4+, and Kr5+ photoionization cross sections; to J.-M. Bizau, who kindly sent experimental data for Kr+; and to M. Richter and A. Sorokin, who provided an analytical fit to the experimental Kr0 PI cross section. I also thank M. Witthoeft and N. R. Badnell for their careful reading of this manuscript and helpful suggestions that improved its quality.

References

  1. Altun, Z., Yumak, A., Yavuz, I., et al. 2007, A&A, 474, 1051 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Amusia, M. Y., Cherepkov, N. A., & Chernysheva, L. V. 1971, Sov. Phys. J. E. T. P., 33, 90 [NASA ADS] [Google Scholar]
  3. Aspin, C. 1994, A&A, 281, L29 [NASA ADS] [Google Scholar]
  4. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  5. Badnell, N. R. 1986, J. Phys. B: At. Mol. Opt. Phys., 19, 3827 [Google Scholar]
  6. Badnell, N. R. 1997, J. Phys. B: At. Mol. Opt. Phys., 30, 1 [Google Scholar]
  7. Badnell, N. R. 2006a, ApJ, 651, L73 [NASA ADS] [CrossRef] [Google Scholar]
  8. Badnell, N. R. 2006b, ApJS, 406, 1151 [Google Scholar]
  9. Badnell, N. R. 2011, Comp. Phys. Comm., 182, 1528 [Google Scholar]
  10. Badnell, N. R., O’Mullane, M. G., Summers, H. P., et al. 2003, A&A, 167, 334 [Google Scholar]
  11. Baldwin, J. A., Verner, E. M., Verner, D. A., et al. 2000, ApJS, 129, 229 [NASA ADS] [CrossRef] [Google Scholar]
  12. Biémont, E., & Hansen, J. E. 1986a, Phys. Scr, 34, 116 [NASA ADS] [CrossRef] [Google Scholar]
  13. Biémont, E., & Hansen, J. E. 1986b, Phys. Scr, 33, 117 [Google Scholar]
  14. Biémont, E., & Hansen, J. E. 1987, Nucl. Inst. Meth. Phys. Res., 23, 274 [NASA ADS] [CrossRef] [Google Scholar]
  15. Biémont, E., Cowan, R. D., & Hansen, J. E. 1988, Phys. Scr, 37, 850 [Google Scholar]
  16. Biémont, E., Hansen, J. E., Quinet, P., & Zeippen, C. J. 1995, A&AS, 111, 333 [NASA ADS] [Google Scholar]
  17. Bizau, J. M., Blancard, C., Coreno, M., et al. 2011, J. Phys. B: At. Mol. Opt. Phys., 44, 055205 [NASA ADS] [CrossRef] [Google Scholar]
  18. Blum, R. D., & McGregor, P. J. 2008, AJ, 135, 1708 [NASA ADS] [CrossRef] [Google Scholar]
  19. Burgess, A. 1964, ApJ, 139, 776 [NASA ADS] [CrossRef] [Google Scholar]
  20. Burgess, A. 1965, ApJ, 141, 1588 [NASA ADS] [CrossRef] [Google Scholar]
  21. Burris, D. L., Pilachowski, C. A., Armandroff, T. E., et al. 2000, ApJ, 544, 302 [NASA ADS] [CrossRef] [Google Scholar]
  22. Busso, M., Gallino, R., & Wasserburg, G. J. 1999, ARA&A, 37, 239 [NASA ADS] [CrossRef] [Google Scholar]
  23. Cowan, R. D., & Griffin, D. C. 1976, J. Opt. Soc. Am., 66, 1010 [NASA ADS] [CrossRef] [Google Scholar]
  24. Dinerstein, H. L. 2001, ApJ, 550, L223 [Google Scholar]
  25. Dinerstein, H. L., & Geballe, T. R. 2001, ApJ, 562, 515 [NASA ADS] [CrossRef] [Google Scholar]
  26. Esteves, D. A., Sterling, N. C., Kilcoyne, A. L. D., et al. 2009, in XXIX International Conference on Phenomena in Ionized Gases, ed. A. Juarez-Reyes, PA1–7, available at http://www.icpig2009.unam.mx/pdf/PA1-7.pdf [Google Scholar]
  27. Esteves, D. A., Sterling, N. C., Kilcoyne, A. L. D., et al. 2010, in 41st Annual Meeting of the APS Division of Atomic, Mol. Opt. Phys., 55, #OPT.10 [Google Scholar]
  28. Esteves, D. A., Sterling, N. C., Kilcoyne, A. L. D., et al. 2011, Phys. Rev. A, in press [Google Scholar]
  29. Ferland, G. J., Korista, K. T., Verner, D. A., et al. 1998, PASP, 110, 761 [Google Scholar]
  30. Henke, B. L., Elgin, R. L., Lent, R. E., & Ledingham, R. B. 1967, Norelco Reporter, 4314, 112 [Google Scholar]
  31. Hyung, S., Aller, L. H., Feibelman, W. A., & Lee, S.-J. 2001, ApJ, 563, 889 [NASA ADS] [CrossRef] [Google Scholar]
  32. Kallman, T., & Bautista, M. 2001, ApJS, 133, 221 [NASA ADS] [CrossRef] [Google Scholar]
  33. Karakas, A. I., van Raai, M. A., Lugaro, M., Sterling, N. C., & Dinerstein, H. L. 2009, ApJ, 690, 1130 [NASA ADS] [CrossRef] [Google Scholar]
  34. Lang, J., & Watson, W. S. 1975, J. Phys. B: At. Mol. Opt. Phys., 8, L339 [NASA ADS] [CrossRef] [Google Scholar]
  35. Likkel, L., Dinerstein, H. L., Lester, D. F., Kindt, A., & Bartig, K. 2006, AJ, 131, 1515 [NASA ADS] [CrossRef] [Google Scholar]
  36. Liu, Y., Liu, X.-W., Luo, S.-G., & Barlow, M. J. 2004, MNRAS, 353, 1231 [NASA ADS] [CrossRef] [Google Scholar]
  37. Lu, M. 2006, Ph.D. Thesis, University of Nevada, Reno [Google Scholar]
  38. Lu, M., Alna’ Washi, G., Habibi, M., et al. 2006a, Phys. Rev. A, 74, 062701 [NASA ADS] [CrossRef] [Google Scholar]
  39. Lu, M., Gharaibeh, M. F., Alna’ Washi, G., et al. 2006b, Phys. Rev. A, 74, 012703 [NASA ADS] [CrossRef] [Google Scholar]
  40. Luhman, K. L., Engelbracht, C. W., & Luhman, M. L. 1998, ApJ, 499, 799 [NASA ADS] [CrossRef] [Google Scholar]
  41. Lumsden, S. L., & Puxley, P. J. 1996, MNRAS, 281, 493 [NASA ADS] [Google Scholar]
  42. Marr, G. V., & West, J. B. 1976, At. Dat. Nucl. Dat. Tab., 18, 497 [NASA ADS] [CrossRef] [Google Scholar]
  43. Mattioli, M., Mazzitelli, G., Fournier, K. B., Finkenthal, M., & Carraro, L. 2006, J. Phys. B: At. Mol. Opt. Phys., 39, 4457 [NASA ADS] [CrossRef] [Google Scholar]
  44. Merts, A., Cowan, R., & Magee, N. 1976, Report LA-6260-MS, Los Alamos National Laboratory [Google Scholar]
  45. Migdalek, J., & Stanek, M. 1989, J. Quant. Spec. Radiat. Transf., 42, 585 [NASA ADS] [CrossRef] [Google Scholar]
  46. Miller, D. L., Dow, J. D., Houlgate, R. G., Marr, G. V., & West, J. B. 1977, J. Phys. B: At. Mol. Phys., 10, 3205 [NASA ADS] [CrossRef] [Google Scholar]
  47. Morton, D. C. 2000, ApJS, 130, 403 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  48. Napiwotzki, R. 1999, A&A, 350, 101 [Google Scholar]
  49. Nikolić, D., Gorczyca, T. W., Korista, K. T., & Badnell, N. R. 2010, A&A, 516, 97 [Google Scholar]
  50. Okumura, S.-I., Mori, A., Watanabe, E., Nishihara, E., & Yamashita, T. 2001, AJ, 121, 2089 [NASA ADS] [CrossRef] [Google Scholar]
  51. Otsuka, M., Tajitsu, A., Hyung, S., & Izumiura, H. 2010, ApJ, 723, 658 [NASA ADS] [CrossRef] [Google Scholar]
  52. Otsuka, M., Meixner, M., Riebel, D., et al. 2011, ApJ, 729, 39 [NASA ADS] [CrossRef] [Google Scholar]
  53. Parpia, F. A., Johnson, W. R., & Radojevic, V. 1984, Phys. Rev. A, 29, 3173 [NASA ADS] [CrossRef] [Google Scholar]
  54. Péquignot, D., & Baluteau, J. P. 1994, A&A, 283, 593 [NASA ADS] [Google Scholar]
  55. Phillips, J. P. 2003, MNRAS, 344, 501 [NASA ADS] [CrossRef] [Google Scholar]
  56. Pindzola, M. S., Badnell, N. R., & Griffin, D. C. 1992, Phys. Rev. A, 46, 5725 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  57. Puxley, P. J., Ramsay Howat, S. K., & Mountain, C. M. 2000, ApJ, 529, 224 [NASA ADS] [CrossRef] [Google Scholar]
  58. Ralchenko, Y., Kramide, A. E., Reader, J., & Team, N. A. 2008, NIST Atomic Spectra Database (version 3.1.4), National Institute of Standards and Technology, Gaithersburg, MD., http://physics.nist.gov/cgi-bin/AtData/main_asd3 [Google Scholar]
  59. Richter, M., Ulm, G., Gerth, C., et al. 2003, in X-Ray and Inner-Shell Processes, ed. A. Bianconi, A. Marcelli, & N. L. Saini, AIP Conf. Proc., 652, 165 [Google Scholar]
  60. Roman-Lopes, A., Abraham, Z., Ortiz, R., & Rodriguez-Ardila, A. 2009, MNRAS, 394, 467 [NASA ADS] [CrossRef] [Google Scholar]
  61. Samson, J. A. R., & Yin, L. 1989, J. Opt. Soc. Am., 6, 2326 [NASA ADS] [CrossRef] [Google Scholar]
  62. Samson, J. A. R., Lyn, L., Haddad, G. N., & Angel, G. C. 1991, J. Phys. IV, C1, 99 [Google Scholar]
  63. Savin, D. W., Kahn, S. M., Linkemann, J., et al. 1999, ApJS, 123, 687 [NASA ADS] [CrossRef] [Google Scholar]
  64. Savin, D. W., Kahn, S. M., Gwinner, G., et al. 2003, ApJS, 147, 421 [NASA ADS] [CrossRef] [Google Scholar]
  65. Schippers, S., Lestinsky, M., Müller, A., Savin, D. W., & Wolf, A. 2010, Int. Rev. At. Mol. Phys., 1, 109 [Google Scholar]
  66. Schöning, T. 1997, A&AS, 122, 277 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Schöning, T., & Butler, K. 1998, A&AS, 128, 581 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Sharpee, B., Williams, R., Baldwin, J. A., & van Hoof, P. A. M. 2003, ApJS, 149, 157 [NASA ADS] [CrossRef] [Google Scholar]
  69. Sharpee, B., Zhang, Y., Williams, R., et al. 2007, ApJ, 659, 1265 [NASA ADS] [CrossRef] [Google Scholar]
  70. Sneden, C., Cowan, J. J., & Gallino, R. 2008, ARA&A, 46, 241 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Sorokin, A. A., Shmaenok, L. A., Bobashev, S. V., et al. 2000, Phys. Rev. A, 61, 022723 [NASA ADS] [CrossRef] [Google Scholar]
  72. Sterling, N. C., & Dinerstein, H. L. 2003, Rev. Mex. Astron. Ser. Conf., 18, 133 [NASA ADS] [Google Scholar]
  73. Sterling, N. C., & Dinerstein, H. L. 2008, ApJS, 174, 158 [NASA ADS] [CrossRef] [Google Scholar]
  74. Sterling, N. C., & Stancil, P. C. 2011, A&A, in press [arXiv:1008.3838] [Google Scholar]
  75. Sterling, N. C., & Witthoeft, M. C. 2011, A&A, 529, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Sterling, N. C., Dinerstein, H. L., & Bowers, C. W. 2002, ApJ, 578, L55 [NASA ADS] [CrossRef] [Google Scholar]
  77. Sterling, N. C., Dinerstein, H. L., & Kallman, T. R. 2007, ApJS, 169, 37 [NASA ADS] [CrossRef] [Google Scholar]
  78. Sterling, N. C., Dinerstein, H. L., Hwang, S., et al. 2009, PASA, 26, 339 [NASA ADS] [Google Scholar]
  79. Sterling, N. C., Esteves, D. A., Bilodeau, R. C., et al. 2011, J. Phys. B: At. Mol. Opt. Phys., 44, 025701 [NASA ADS] [CrossRef] [Google Scholar]
  80. Summers, H. P. 2005, ADAS Manual v2.7, http://adas.phys.strath.ac.uk [Google Scholar]
  81. Treffers, R. R., Fink, U., Larson, H. P., & Gautier, T. N. I. 1976, ApJ, 209, 793 [NASA ADS] [CrossRef] [Google Scholar]
  82. Tulkki, J., Aksela, S., Aksela, H., et al. 1992, Phys. Rev. A, 45, 4640 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  83. Vanzi, L., Cresci, G., Telles, E., & Melnick, J. 2008, A&A, 486, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Verner, D. A., & Ferland, G. J. 1996, ApJS, 103, 467 [NASA ADS] [CrossRef] [Google Scholar]
  85. Villaver, E., Stanghellini, L., & Shaw, R. A. 2007, ApJ, 656, 831 [NASA ADS] [CrossRef] [Google Scholar]
  86. West, J. B., & Marr, G. V. 1976, Proc. R. Soc. Lond., 349, 397 [Google Scholar]
  87. Williams, R., Jenkins, E. B., Baldwin, J. A., et al. 2008, ApJ, 677, 1100 [NASA ADS] [CrossRef] [Google Scholar]
  88. Wolf, A., Buhr, H., Grieser, M., et al. 2006, Hyperfine Interactions, 172, 111 [Google Scholar]
  89. Wood, J. L., Dinerstein, H. L., Geballe, T. R., & Sterling, N. C. 2007, in BAAS, 38, 1113 [Google Scholar]
  90. Zatasrinny, O., Gorczyca, T. W., Korista, K. T., Badnell, N. R., & Savin, D. W. 2003, A&A, 412, 587 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  91. Zhang, Y., Liu, X.-W., Luo, S.-G., Péquignot, D., & Barlow, M. J. 2005, A&A, 442, 249 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1

Configuration-interaction expansions used for each Kr ion.

Table 2

Radial scaling parameters used for CI expansions of each ion.

Table 3

Comparison of selected calculated and experimental energies (in Rydbergs).

Table 4

Calculated Einstein A-coefficients (in s-1) compared to literature values.

Table 5

Fit coefficients for radiative recombination rate coefficients at low temperature (see Eqs. (1) and (2)).

Table 6

Fit coefficients for radiative recombination rate coefficients at high temperatures, valid from Tmax (see Table 5) to 107z2 K. See Eq. (3).

Table 7

Fit coefficients for dielectronic recombination rate coefficients; see Eq. (3).

All Figures

thumbnail Fig. 1

Ground state PI cross sections for Kr0−Kr5+ near their ionization thresholds.

In the text
thumbnail Fig. 2

Comparison of Kr2+ ground state PI cross sections computed for each of the three CI expansions (as indicated in the top three panels). The direct PI cross sections are compared in the bottom panel.

In the text
thumbnail Fig. 3

Comparison of Kr3+ ground state PI cross sections computed for each of the three CI expansions (as indicated in the top three panels). The direct PI cross sections are compared in the bottom panel.

In the text
thumbnail Fig. 4

Comparison of calculated Kr+ PI cross section to the experimental measurements (dotted line) of Bizau et al. (2011). The dashed line shows our original result before comparison to experiment, while the solid line shows the theoretical cross section after optimization to better reproduce the Bizau et al. (2011) results. In each case, the theoretical curves are statistically-weighted linear combinations of the cross sections from each of the two states in the ground configuration. The calculated cross sections were convolved with a Gaussian of 30 meV to approximate the average resolution of the Bizau et al. (2011) measurements.

In the text
thumbnail Fig. 5

Comparison of calculated PI cross sections (solid lines) with experimental data (dotted lines) for Kr0 (Richter et al. 2003), Kr3+ (Lu et al. 2006b), Kr4+ (Lu 2006), and Kr5+ (Lu et al. 2006a). The theoretical curves have been convolved with a Gaussian of 1 eV for Kr0 and 100 meV for the other species, to match the experimental resolution. The theoretical data are linear combinations of PI cross sections from each state in the ground configuration, weighted according to statistical weights in the case of Kr3+, and by the factors (3P0, 3P1, 3P2, 1D2, 1S0) = (0.32, 0.15, 0.2, 0.3, 0.03) for Kr4+ and (2P1/2, 2P3/2) = (0.4, 0.6) for Kr5+.

In the text
thumbnail Fig. 6

Radiative recombination rate coefficients (black lines) as a function of temperature for ground state target ions Kr+ (solid curve), Kr2+ (dotted curve), Kr3+ (short-dashed curve), Kr4+ (dash-dot-dash curve), Kr5+ (long-dashed curve), and Kr6+ (dash-dot-dot-dot-dash curve). For comparison, the approximate rate coefficients of Mattioli et al. (2006) are depicted as gray curves, with the line styles corresponding to the same ions.

In the text
thumbnail Fig. 7

Dielectronic recombination rate coefficients (black lines) as a function of temperature for ground state target ions Kr+ (solid curve), Kr2+ (dotted curve), Kr3+ (short-dashed curve), Kr4+ (dash-dot-dash curve), Kr5+ (long-dashed curve), and Kr6+ (dash-dot-dot-dot-dash curve). DR rate coefficients computed with the Burgess formula over Δn = 0 and Δn = 1 resonances, from Mattioli et al. (2006), are depicted as gray curves, with the line styles corresponding to the same ions.

In the text
thumbnail Fig. 8

Comparison of the radiative (dotted curves), dielectronic (dashed curves), and total recombination rate coefficients (solid curves) for the ground state of the first six Kr ions.

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.