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/00046361/201117471  
Published online  29 August 2011 
Atomic data for neutroncapture elements
II. Photoionization and recombination properties of lowcharge krypton ions^{⋆}
Michigan State UniversityDepartment of Physics and Astronomy,
3248 Biomedical Physical Sciences,
East Lansing
MI
488242320,
USA
email: sterling@pa.msu.edu
Received: 13 June 2011
Accepted: 11 July 2011
We present multiconfiguration BreitPauli distortedwave photoionization (PI) cross sections and radiative recombination (RR) and dielectronic recombination (DR) rate coefficients for the first six krypton ions. These were calculated with the AUTOSTRUCTURE code, using semirelativistic radial wavefunctions in intermediate coupling. Kr has been detected in several planetary nebulae (PNe) and H II regions, and is a useful tracer of neutroncapture nucleosynthesis. PI, RR, and DR data are required to accurately correct for unobserved Kr ions in photoionized nebulae, and hence to determine elemental Kr abundances. PI cross sections have been determined for ground configuration states of Kr^{0}−Kr^{5+} up to 100 Rydbergs. Our Kr^{+} PI calculations were significantly improved through comparison with experimental measurements. RR and DR rate coefficients were respectively determined from the direct and resonant PI cross sections at temperatures (10^{1}−10^{7})z^{2} K, where z is the charge. We account for Δn = 0 DR core excitations, and find that DR is the dominant recombination mechanism for all but Kr^{+} at photoionized plasma temperatures. Internal uncertainties are estimated by comparing results computed with three different configurationinteraction expansions for each ion, and by testing the sensitivity to variations in the orbital radial scaling parameters. The PI cross sections are generally uncertain by 30−50% near the ground state thresholds. Near 10^{4} K, the RR rate coefficients are typically uncertain by less than 10%, while those of DR exhibit uncertainties of factors of 2 to 3, due to the unknown energies of nearthreshold autoionizing resonances. With the charge transfer rate coefficients presented in the third paper of this series, these data enable robust Kr abundance determinations in photoionized nebulae for the first time, providing a new tool for studying heavy element enrichments in PNe and for investigating the chemical evolution of transiron elements.
Key words: atomic data / atomic processes / HII regions / nuclear reactions / nucleosynthesis, abundances / stars: evolution / planetary nebulae: general
The calculation results are only available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/533/A62
© 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. Transiron 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 transiron elements via slow ncapture nucleosynthesis (the sprocess) during the asymptotic giant branch (AGB) phase of evolution. Abundance determinations of PNe complement spectroscopic investigations of AGB stars, enabling previously unexplored details of sprocess nucleosynthesis to be revealed (e.g., see Sterling & Dinerstein 2008; Sterling & Witthoeft 2011). Indeed, many of the ncapture 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 transiron elements. In spite of the low cosmic abundances of ncapture 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 nearinfrared 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 sprocess elemental enrichments in a statistically meaningful sample of objects.
These detections are not limited to PNe. Indeed, ncapture 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; RomanLopes et al. 2009) and starburst galaxies (Vanzi et al. 2008). These detections demonstrate the promise of nebular spectroscopy as a tool to study ncapture nucleosynthesis in lowmass stars, the chemical evolution of transiron elements in the Universe, and the heavyelement nucleosynthetic histories of other galaxies.
Kr is one of the most readily detected ncapture elements in ionized nebulae, along with Se and Xe, with relatively strong transitions in both the optical (Péquignot & Baluteau 1994) and nearinfrared (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 sprocess nucleosynthesis. It can be highly enriched by the sprocess, due in part to the ^{86}Kr isotope, which has a magic number of neutrons and acts as a bottleneck in the sprocess path. Isotopes with magic numbers of neutrons have very small ncapture cross sections and cause the peaks seen in the elementbyelement sprocess enrichment distribution near Z = 40, 56, and 82 (e.g., Busso et al. 1999; Burris et al. 2000; Sneden et al. 2008). Typically, sprocess 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 gasphase 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 ncapture 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 transiron 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 ncapture 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 hightemperature 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 ncapture element ions. In the first (Sterling & Witthoeft 2011), multiconfiguration BreitPauli (MCBP) distortedwave 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 LandauZener approximations. In this study, we present MCBP distortedwave 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 highlycharged 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 (Kr^{6+} 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 lowcharge 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 lowcharge 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 collisionallyionized plasmas (indeed, the work of Mattioli et al. was inspired by the utility of Kr as a coolant for magneticallyconfined fusion plasmas), but is known to badly underestimate the rate coefficient at the low temperatures of photoionized plasmas. Their RR rate coefficients for lowcharge 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 Kr^{3+} to Kr^{5+} (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 lowcharge 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 crosssection measurements, we computed data using three different configurationinteraction (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 lowcharge 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.
Configurationinteraction expansions used for each Kr ion.
2. Calculations and methodology
The electronic structure, distortedwave 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 ThomasFermiDiracAmaldi 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 lowestenergy 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 ^{1}D_{2} and ^{2}D_{3/2,5/2} in the ground configurations of Kr^{2 + }−Kr^{4 + } (the disagreement is as high as 22% for the medium CI expansion of Kr^{3+} 4s^{2} 4p^{4} ^{2}D_{3/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 Kr^{2+} and Kr^{3+}. The ionization potential of Kr^{2+} 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 ^{3}P 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 Kr^{2+} PI cross section computed with the two values of the 4p scaling parameter differs by less than 2% near the ground state threshold. The Kr^{3+} 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.
Comparison of selected calculated and experimental energies (in Rydbergs).
In Table 4, we compare the Einstein Acoefficients we computed with previously determined values. For forbidden transitions in the ground configuration, our Acoefficients 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., ^{1}D_{2} for Kr^{2+} and Kr^{4+}, and ^{2}D_{3/2,5/2} for Kr^{3+}). For permitted transitions, the disagreement is generally larger (especially for neutral and singlyionized Kr).
Calculated Einstein Acoefficients (in s^{1}) compared to literature values.
2.2. Photoionization, radiative recombination, and dielectronic recombination
In lowdensity nebulae such as PNe and H II regions, nearly all atomic ions reside in ground configuration states. Therefore, we computed MCBP distortedwave PI cross sections for each level of the Kr^{0} to Kr^{5+} 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 Kr^{6+} over the temperature range (10^{1}−10^{7})z^{2} 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 Kr^{0} and Kr^{+}, whose lowest energy 4s^{2} 4p^{k} 5s and (in the case of Kr^{0}) 4s^{2} 4p^{k} 5p levels (k = 4 and 5 for Kr^{+} and Kr^{0}, respectively) lie below those of 4s^{2} 4p^{k} 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 (Nelectron) 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 oneelectron 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 oneelectron 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 Kr^{5+} and Kr^{6+} 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 Kr^{0} to Kr^{5+} PI cross sections near their ionization thresholds. The AUTOSTRUCTURE output was processed with the codes ADASPI (for direct PI) and ADASPE (for photoexcitationautoionization), which produce the ADAS adf39 and adf38 output files, respectively (Summers 2005). For ease of presentation, we used the postprocessing code xpeppi to add the direct and resonant portions of the PI cross sections.
Fig. 1
Ground state PI cross sections for Kr^{0}−Kr^{5+} near their ionization thresholds. 

Open with DEXTER 
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 directplusresonant 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 photoelectron 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 photoexcitationautoionization 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 crosssection 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 Kr^{2+} and Kr^{3+} 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 photoexcitationautoionization processes in the distortedwave 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).
Fig. 2
Comparison of Kr^{2+} 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. 

Open with DEXTER 
Fig. 3
Comparison of Kr^{3+} 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. 

Open with DEXTER 
We therefore focus on estimating uncertainties in the direct PI cross sections, particularly near the ground state ionization threshold. The Kr^{2+} 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 Kr^{3+} direct cross sections agree to within 5% (Fig. 3). Likewise, the direct cross section computed with the medium Kr^{4+} and Kr^{5+} 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 Kr^{0} and Kr^{3+} 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 Kr^{2+}.
The measurements were conducted at a variety of synchrotron radiation facilities, and in the case of ionized species utilized the merged beams method. Kr^{3+} through Kr^{5+} 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 ^{2}P_{1/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 electronimpact ionization cross section ratio and normalized to the double ion chamber PI measurements of Samson & Yin (1989) and Samson et al. (1991).
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 statisticallyweighted 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. 

Open with DEXTER 
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 reoptimized 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 4p^{6}, 4s^{2} 4p^{4} 4d, and 4s^{2} 4p^{4} 5s configurations and the 5p scaling parameter was optimized on 4s^{2} 4p^{3} 5p^{2} terms (this configuration was later removed from the CI expansion) in order to reproduce the energy of the first excited state (^{2}P_{1/2}). The 4d scaling parameter was optimized on 4s^{2} 4p^{3} 4d^{2} 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 lowestenergy configurations containing an electron in those orbitals (4s^{2} 4p^{4} 4d and 4s^{2} 4p^{4} 5p, respectively). This new optimization reduced the accuracy of the first excited state energy from 0.8% to 11% (higherenergy 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 (^{3}P_{0}, ^{3}P_{1}, ^{3}P_{2}, ^{1}D_{2}, ^{1}S_{0}) = (0.32, 0.15, 0.20, 0.3, 0.03) for Kr^{4+}, (^{2}P_{1/2}, ^{2}P_{3/2}) = (0.4, 0.6) for Kr^{5+}, and by statistical weights in the case of Kr^{3+}. The calculated cross sections were convolved with a Gaussian with FWHM equal to the experimental resolution (~1 eV for Kr^{0}, and 0.1 eV for the other three ions).
Fig. 5
Comparison of calculated PI cross sections (solid lines) with experimental data (dotted lines) for Kr^{0} (Richter et al. 2003), Kr^{3+} (Lu et al. 2006b), Kr^{4+} (Lu 2006), and Kr^{5+} (Lu et al. 2006a). The theoretical curves have been convolved with a Gaussian of 1 eV for Kr^{0} 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 Kr^{3+}, and by the factors (^{3}P_{0}, ^{3}P_{1}, ^{3}P_{2}, ^{1}D_{2}, ^{1}S_{0}) = (0.32, 0.15, 0.2, 0.3, 0.03) for Kr^{4+} and (^{2}P_{1/2}, ^{2}P_{3/2}) = (0.4, 0.6) for Kr^{5+}. 

Open with DEXTER 
In the cases of Kr^{4+} and Kr^{5+}, the direct cross sections generally agree within the experimental noise. The theoretical cross section is slightly too large for Kr^{4+} and perhaps slightly small for Kr^{5+}, but over the full range of energies they agree to within the experimental uncertainties of ~25%.
The agreement is markedly worse for Kr^{3+}, 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 Kr^{3+} 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 reoptimization designed to increase the 4d (and hence continuum d) scaling parameter – by optimizing on 4s^{2} 4p 4d^{2} instead of 4s^{2} 4p^{2} 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 Kr^{3+} PI cross section is uncertain by approximately a factor of two.
Likewise, for Kr^{0} 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^{+} 4s^{2} 4p^{4} 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 Kr^{3+}, 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 finalstate 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 Jresolved for n ≤ 8 levels, bundlednl 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) (1)where (2)The fit coefficients T and T_{0,1,2} are in temperature units, A and α_{RR}(T) are in cm^{3} 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 (10^{1}−10^{7})z^{2} 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 lowenergy levels, which are not as effectively screened from the nucleus as higher levels. Equations (1) and (2) can be used to fit the lowtemperature behavior of the RR rate coefficients, but a different analytical form is often necessary to fit the hightemperature behavior. For these hightemperature fits, we have chosen to utilize the formula typically used to fit DR rate coefficients (e.g., Zatasrinny et al. 2003): (3)where T and E_{i} are in temperature units, c_{i} are in cm^{3} s^{1} K^{3/2}, and α_{RR}(T) is in cm^{3} 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 lowtemperature fits (using Eqs. (1) and (2)) and Table 6 for the hightemperature fits (using Eq. (3)). Table 5 also lists the maximum temperature at which the lowtemperature fits are valid for each state; if T_{max} is not given, then the lowtemperature fits are valid over the entire temperature range (10^{1}−10^{7})z^{2} 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 lowtemperature and hightemperature asymptotic behaviors of the rate coefficients.
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.
Fig. 6
Radiative recombination rate coefficients (black lines) as a function of temperature for ground state target ions Kr^{+} (solid curve), Kr^{2+} (dotted curve), Kr^{3+} (shortdashed curve), Kr^{4+} (dashdotdash curve), Kr^{5+} (longdashed curve), and Kr^{6+} (dashdotdotdotdash 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. 

Open with DEXTER 
3.2.1. Estimation of uncertainties
We have estimated the internal uncertainties of our calculated RR rate coefficients near T = 10^{4} 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 lowestcharge species. For the Kr^{2+} target, using target ion scaling parameters resulted in rate coefficients differing by up to 10% for metastable states. For Kr^{+} forming Kr^{0}, calculations with target ion scaling parameters led to differences of 20−30%, while the use of nonorthogonal 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 > 10^{5} for the first two ionization stages, or >10^{6} 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 finalstate resolved DR rate coefficients from the resonant portion of the PI cross sections, with the aid of the ADASDR postprocessing 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 nonlinear least squares fit algorithm produced fits accurate to within 5% over the temperature range (10^{1}−10^{7})z^{2} 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.
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 nearthreshold autoionizing resonances, but show similar behavior at high temperatures.
Mattioli et al. (2006) also computed DR rate coefficients for lowcharge Kr ions (shown in Fig. 7 as gray curves), in order to properly model its cooling rate coefficient in magneticallyconfined 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 wellknown to inadequately describe lowtemperature 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 collisionallyionized plasmas, it is not sufficient for photoionized plasmas, in which the electrons are captured into lowenergy 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 10^{5} K. At high temperatures (≥a few times 10^{6} 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 Kr^{2+}).
Fig. 7
Dielectronic recombination rate coefficients (black lines) as a function of temperature for ground state target ions Kr^{+} (solid curve), Kr^{2+} (dotted curve), Kr^{3+} (shortdashed curve), Kr^{4+} (dashdotdash curve), Kr^{5+} (longdashed curve), and Kr^{6+} (dashdotdotdotdash 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. 

Open with DEXTER 
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 lowcharge Se ions (Sterling & Witthoeft 2011), DR dominates RR for Kr ions near 10^{4} 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 × 10^{3} K). Only in the case of Kr^{+} is RR the most important recombination mechanism near ~10^{4} K. These results underline the importance of DR in the ionization balance of transiron elements.
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. 

Open with DEXTER 
3.3.1. Estimation of uncertainties
DR calculations are inherently more uncertain than those of PI or RR, since the unknown energies of lowlying autoionizing levels usually dominate other sources of uncertainty. The energies of these nearthreshold 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 nearthreshold resonances. We accomplished this with the ADASDR postprocessing 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 10^{4} K for the ground (metastable) states are 75% (factor of 20) for Kr^{+}, a factor of 3 (10% to an order of magnitude) for Kr^{2+}, 10% (15−35%) for Kr^{3+}, 70% (50−70%) for Kr^{4+}, 50% (25%) for Kr^{5+}, and a factor of 2 for Kr^{6+}. In contrast, DR calculations with the small CI expansions^{1} led to rate coefficients different from those computed with the medium configuration sets by a factor of 3 for Kr^{+}, 60% for Kr^{2+}, a factor of 2 for Kr^{5+}, and 10−30% for other Kr ions.
With the exception of Kr^{3+}, 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 lowtemperature DR rate coefficients for lowcharge 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 distortedwave PI cross sections over the photon energy range 0−100 Ryd, and radiative and dielectronic recombination rate coefficients at temperatures (10^{1}−10^{7})z^{2} 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 5−7.
Kr has been detected in the spectra of planetary nebulae and H II regions, and its abundance is a sensitive tracer of enrichments by ncapture 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 lowcharge 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 Kr^{0} and Kr^{3+} 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 crosssection 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 lowlying 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 lowcharge 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.
Acknowledgments
I acknowledge support from an NSF Astronomy and Astrophysics Postdoctoral Fellowship under award AST0901432 and from NASA grant 06APRA2060049. I am grateful to R. Phaneuf, who provided experimental data for Kr^{3+}, Kr^{4+}, and Kr^{5+} 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 Kr^{0} 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
 Altun, Z., Yumak, A., Yavuz, I., et al. 2007, A&A, 474, 1051 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Amusia, M. Y., Cherepkov, N. A., & Chernysheva, L. V. 1971, Sov. Phys. J. E. T. P., 33, 90 [NASA ADS] [Google Scholar]
 Aspin, C. 1994, A&A, 281, L29 [NASA ADS] [Google Scholar]
 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Badnell, N. R. 1986, J. Phys. B: At. Mol. Opt. Phys., 19, 3827 [NASA ADS] [CrossRef] [Google Scholar]
 Badnell, N. R. 1997, J. Phys. B: At. Mol. Opt. Phys., 30, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Badnell, N. R. 2006a, ApJ, 651, L73 [NASA ADS] [CrossRef] [Google Scholar]
 Badnell, N. R. 2006b, ApJS, 406, 1151 [Google Scholar]
 Badnell, N. R. 2011, Comp. Phys. Comm., 182, 1528 [NASA ADS] [CrossRef] [Google Scholar]
 Badnell, N. R., O’Mullane, M. G., Summers, H. P., et al. 2003, A&A, 167, 334 [Google Scholar]
 Baldwin, J. A., Verner, E. M., Verner, D. A., et al. 2000, ApJS, 129, 229 [NASA ADS] [CrossRef] [Google Scholar]
 Biémont, E., & Hansen, J. E. 1986a, Phys. Scr, 34, 116 [NASA ADS] [CrossRef] [Google Scholar]
 Biémont, E., & Hansen, J. E. 1986b, Phys. Scr, 33, 117 [NASA ADS] [CrossRef] [Google Scholar]
 Biémont, E., & Hansen, J. E. 1987, Nucl. Inst. Meth. Phys. Res., 23, 274 [NASA ADS] [CrossRef] [Google Scholar]
 Biémont, E., Cowan, R. D., & Hansen, J. E. 1988, Phys. Scr, 37, 850 [NASA ADS] [CrossRef] [Google Scholar]
 Biémont, E., Hansen, J. E., Quinet, P., & Zeippen, C. J. 1995, A&AS, 111, 333 [NASA ADS] [Google Scholar]
 Bizau, J. M., Blancard, C., Coreno, M., et al. 2011, J. Phys. B: At. Mol. Opt. Phys., 44, 055205 [NASA ADS] [CrossRef] [Google Scholar]
 Blum, R. D., & McGregor, P. J. 2008, AJ, 135, 1708 [NASA ADS] [CrossRef] [Google Scholar]
 Burgess, A. 1964, ApJ, 139, 776 [NASA ADS] [CrossRef] [Google Scholar]
 Burgess, A. 1965, ApJ, 141, 1588 [NASA ADS] [CrossRef] [Google Scholar]
 Burris, D. L., Pilachowski, C. A., Armandroff, T. E., et al. 2000, ApJ, 544, 302 [NASA ADS] [CrossRef] [Google Scholar]
 Busso, M., Gallino, R., & Wasserburg, G. J. 1999, ARA&A, 37, 239 [NASA ADS] [CrossRef] [Google Scholar]
 Cowan, R. D., & Griffin, D. C. 1976, J. Opt. Soc. Am., 66, 1010 [NASA ADS] [CrossRef] [Google Scholar]
 Dinerstein, H. L. 2001, ApJ, 550, L223 [NASA ADS] [CrossRef] [Google Scholar]
 Dinerstein, H. L., & Geballe, T. R. 2001, ApJ, 562, 515 [NASA ADS] [CrossRef] [Google Scholar]
 Esteves, D. A., Sterling, N. C., Kilcoyne, A. L. D., et al. 2009, in XXIX International Conference on Phenomena in Ionized Gases, ed. A. JuarezReyes, PA1–7, available at http://www.icpig2009.unam.mx/pdf/PA17.pdf [Google Scholar]
 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]
 Esteves, D. A., Sterling, N. C., Kilcoyne, A. L. D., et al. 2011, Phys. Rev. A, in press [Google Scholar]
 Ferland, G. J., Korista, K. T., Verner, D. A., et al. 1998, PASP, 110, 761 [NASA ADS] [CrossRef] [Google Scholar]
 Henke, B. L., Elgin, R. L., Lent, R. E., & Ledingham, R. B. 1967, Norelco Reporter, 4314, 112 [Google Scholar]
 Hyung, S., Aller, L. H., Feibelman, W. A., & Lee, S.J. 2001, ApJ, 563, 889 [NASA ADS] [CrossRef] [Google Scholar]
 Kallman, T., & Bautista, M. 2001, ApJS, 133, 221 [NASA ADS] [CrossRef] [Google Scholar]
 Karakas, A. I., van Raai, M. A., Lugaro, M., Sterling, N. C., & Dinerstein, H. L. 2009, ApJ, 690, 1130 [NASA ADS] [CrossRef] [Google Scholar]
 Lang, J., & Watson, W. S. 1975, J. Phys. B: At. Mol. Opt. Phys., 8, L339 [NASA ADS] [CrossRef] [Google Scholar]
 Likkel, L., Dinerstein, H. L., Lester, D. F., Kindt, A., & Bartig, K. 2006, AJ, 131, 1515 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, Y., Liu, X.W., Luo, S.G., & Barlow, M. J. 2004, MNRAS, 353, 1231 [NASA ADS] [CrossRef] [Google Scholar]
 Lu, M. 2006, Ph.D. Thesis, University of Nevada, Reno [Google Scholar]
 Lu, M., Alna’ Washi, G., Habibi, M., et al. 2006a, Phys. Rev. A, 74, 062701 [NASA ADS] [CrossRef] [Google Scholar]
 Lu, M., Gharaibeh, M. F., Alna’ Washi, G., et al. 2006b, Phys. Rev. A, 74, 012703 [NASA ADS] [CrossRef] [Google Scholar]
 Luhman, K. L., Engelbracht, C. W., & Luhman, M. L. 1998, ApJ, 499, 799 [NASA ADS] [CrossRef] [Google Scholar]
 Lumsden, S. L., & Puxley, P. J. 1996, MNRAS, 281, 493 [NASA ADS] [Google Scholar]
 Marr, G. V., & West, J. B. 1976, At. Dat. Nucl. Dat. Tab., 18, 497 [NASA ADS] [CrossRef] [Google Scholar]
 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]
 Merts, A., Cowan, R., & Magee, N. 1976, Report LA6260MS, Los Alamos National Laboratory [Google Scholar]
 Migdalek, J., & Stanek, M. 1989, J. Quant. Spec. Radiat. Transf., 42, 585 [NASA ADS] [CrossRef] [Google Scholar]
 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]
 Morton, D. C. 2000, ApJS, 130, 403 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Napiwotzki, R. 1999, A&A, 350, 101 [Google Scholar]
 Nikolić, D., Gorczyca, T. W., Korista, K. T., & Badnell, N. R. 2010, A&A, 516, 97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Okumura, S.I., Mori, A., Watanabe, E., Nishihara, E., & Yamashita, T. 2001, AJ, 121, 2089 [NASA ADS] [CrossRef] [Google Scholar]
 Otsuka, M., Tajitsu, A., Hyung, S., & Izumiura, H. 2010, ApJ, 723, 658 [NASA ADS] [CrossRef] [Google Scholar]
 Otsuka, M., Meixner, M., Riebel, D., et al. 2011, ApJ, 729, 39 [NASA ADS] [CrossRef] [Google Scholar]
 Parpia, F. A., Johnson, W. R., & Radojevic, V. 1984, Phys. Rev. A, 29, 3173 [NASA ADS] [CrossRef] [Google Scholar]
 Péquignot, D., & Baluteau, J. P. 1994, A&A, 283, 593 [NASA ADS] [Google Scholar]
 Phillips, J. P. 2003, MNRAS, 344, 501 [NASA ADS] [CrossRef] [Google Scholar]
 Pindzola, M. S., Badnell, N. R., & Griffin, D. C. 1992, Phys. Rev. A, 46, 5725 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Puxley, P. J., Ramsay Howat, S. K., & Mountain, C. M. 2000, ApJ, 529, 224 [NASA ADS] [CrossRef] [Google Scholar]
 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/cgibin/AtData/main_asd3 [Google Scholar]
 Richter, M., Ulm, G., Gerth, C., et al. 2003, in XRay and InnerShell Processes, ed. A. Bianconi, A. Marcelli, & N. L. Saini, AIP Conf. Proc., 652, 165 [Google Scholar]
 RomanLopes, A., Abraham, Z., Ortiz, R., & RodriguezArdila, A. 2009, MNRAS, 394, 467 [NASA ADS] [CrossRef] [Google Scholar]
 Samson, J. A. R., & Yin, L. 1989, J. Opt. Soc. Am., 6, 2326 [NASA ADS] [CrossRef] [Google Scholar]
 Samson, J. A. R., Lyn, L., Haddad, G. N., & Angel, G. C. 1991, J. Phys. IV, C1, 99 [Google Scholar]
 Savin, D. W., Kahn, S. M., Linkemann, J., et al. 1999, ApJS, 123, 687 [NASA ADS] [CrossRef] [Google Scholar]
 Savin, D. W., Kahn, S. M., Gwinner, G., et al. 2003, ApJS, 147, 421 [NASA ADS] [CrossRef] [Google Scholar]
 Schippers, S., Lestinsky, M., Müller, A., Savin, D. W., & Wolf, A. 2010, Int. Rev. At. Mol. Phys., 1, 109 [Google Scholar]
 Schöning, T. 1997, A&AS, 122, 277 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Schöning, T., & Butler, K. 1998, A&AS, 128, 581 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Sharpee, B., Williams, R., Baldwin, J. A., & van Hoof, P. A. M. 2003, ApJS, 149, 157 [NASA ADS] [CrossRef] [Google Scholar]
 Sharpee, B., Zhang, Y., Williams, R., et al. 2007, ApJ, 659, 1265 [NASA ADS] [CrossRef] [Google Scholar]
 Sneden, C., Cowan, J. J., & Gallino, R. 2008, ARA&A, 46, 241 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sorokin, A. A., Shmaenok, L. A., Bobashev, S. V., et al. 2000, Phys. Rev. A, 61, 022723 [NASA ADS] [CrossRef] [Google Scholar]
 Sterling, N. C., & Dinerstein, H. L. 2003, Rev. Mex. Astron. Ser. Conf., 18, 133 [NASA ADS] [Google Scholar]
 Sterling, N. C., & Dinerstein, H. L. 2008, ApJS, 174, 158 [NASA ADS] [CrossRef] [Google Scholar]
 Sterling, N. C., & Stancil, P. C. 2011, A&A, in press [arXiv:1008.3838] [Google Scholar]
 Sterling, N. C., & Witthoeft, M. C. 2011, A&A, 529, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sterling, N. C., Dinerstein, H. L., & Bowers, C. W. 2002, ApJ, 578, L55 [NASA ADS] [CrossRef] [Google Scholar]
 Sterling, N. C., Dinerstein, H. L., & Kallman, T. R. 2007, ApJS, 169, 37 [NASA ADS] [CrossRef] [Google Scholar]
 Sterling, N. C., Dinerstein, H. L., Hwang, S., et al. 2009, PASA, 26, 339 [NASA ADS] [Google Scholar]
 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]
 Summers, H. P. 2005, ADAS Manual v2.7, http://adas.phys.strath.ac.uk [Google Scholar]
 Treffers, R. R., Fink, U., Larson, H. P., & Gautier, T. N. I. 1976, ApJ, 209, 793 [NASA ADS] [CrossRef] [Google Scholar]
 Tulkki, J., Aksela, S., Aksela, H., et al. 1992, Phys. Rev. A, 45, 4640 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Vanzi, L., Cresci, G., Telles, E., & Melnick, J. 2008, A&A, 486, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Verner, D. A., & Ferland, G. J. 1996, ApJS, 103, 467 [NASA ADS] [CrossRef] [Google Scholar]
 Villaver, E., Stanghellini, L., & Shaw, R. A. 2007, ApJ, 656, 831 [NASA ADS] [CrossRef] [Google Scholar]
 West, J. B., & Marr, G. V. 1976, Proc. R. Soc. Lond., 349, 397 [NASA ADS] [CrossRef] [Google Scholar]
 Williams, R., Jenkins, E. B., Baldwin, J. A., et al. 2008, ApJ, 677, 1100 [NASA ADS] [CrossRef] [Google Scholar]
 Wolf, A., Buhr, H., Grieser, M., et al. 2006, Hyperfine Interactions, 172, 111 [NASA ADS] [CrossRef] [Google Scholar]
 Wood, J. L., Dinerstein, H. L., Geballe, T. R., & Sterling, N. C. 2007, in BAAS, 38, 1113 [Google Scholar]
 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]
 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
All Figures
Fig. 1
Ground state PI cross sections for Kr^{0}−Kr^{5+} near their ionization thresholds. 

Open with DEXTER  
In the text 
Fig. 2
Comparison of Kr^{2+} 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. 

Open with DEXTER  
In the text 
Fig. 3
Comparison of Kr^{3+} 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. 

Open with DEXTER  
In the text 
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 statisticallyweighted 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. 

Open with DEXTER  
In the text 
Fig. 5
Comparison of calculated PI cross sections (solid lines) with experimental data (dotted lines) for Kr^{0} (Richter et al. 2003), Kr^{3+} (Lu et al. 2006b), Kr^{4+} (Lu 2006), and Kr^{5+} (Lu et al. 2006a). The theoretical curves have been convolved with a Gaussian of 1 eV for Kr^{0} 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 Kr^{3+}, and by the factors (^{3}P_{0}, ^{3}P_{1}, ^{3}P_{2}, ^{1}D_{2}, ^{1}S_{0}) = (0.32, 0.15, 0.2, 0.3, 0.03) for Kr^{4+} and (^{2}P_{1/2}, ^{2}P_{3/2}) = (0.4, 0.6) for Kr^{5+}. 

Open with DEXTER  
In the text 
Fig. 6
Radiative recombination rate coefficients (black lines) as a function of temperature for ground state target ions Kr^{+} (solid curve), Kr^{2+} (dotted curve), Kr^{3+} (shortdashed curve), Kr^{4+} (dashdotdash curve), Kr^{5+} (longdashed curve), and Kr^{6+} (dashdotdotdotdash 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. 

Open with DEXTER  
In the text 
Fig. 7
Dielectronic recombination rate coefficients (black lines) as a function of temperature for ground state target ions Kr^{+} (solid curve), Kr^{2+} (dotted curve), Kr^{3+} (shortdashed curve), Kr^{4+} (dashdotdash curve), Kr^{5+} (longdashed curve), and Kr^{6+} (dashdotdotdotdash 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. 

Open with DEXTER  
In the text 
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. 

Open with DEXTER  
In the text 
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.