Issue |
A&A
Volume 589, May 2016
|
|
---|---|---|
Article Number | A137 | |
Number of page(s) | 4 | |
Section | Atomic, molecular, and nuclear data | |
DOI | https://doi.org/10.1051/0004-6361/201628457 | |
Published online | 26 April 2016 |
Research Note
K-shell photoabsorption and photoionisation of trace elements
I. Isoelectronic sequences with electron number 3 ≤N ≤ 11
1 Physique Atomique et Astrophysique, Université de Mons – UMONS, 20 place du Parc, 7000 Mons, Belgium
e-mail:
patrick.palmeri@umons.ac.be
2 IPNAS, Université de Liège, Campus du Sart Tilman, Bât. B15, 4000 Liège, Belgium
e-mail:
pascal.quinet@umons.ac.be
3 Centro de Física, Instituto Venezolano de Investigación Científicas (IVIC), A.P. 21827, Caracas 1020, Venezuela
e-mail:
claudio@ivic.gov.ve
4 Department of Physics, Western Michigan University, 1903 W Michigan Ave., Kalamazoo, MI 49008, USA
5 ADNET Systems, Inc., Bethesda, MD 20817, USA
6 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA
Received: 8 March 2016
Accepted: 25 March 2016
Context. With the recent launching of the Hitomi X-ray space observatory, K lines and edges of chemical elements with low cosmic abundances, namely F, Na, P, Cl, K, Sc, Ti, V, Cr, Mn, Co, Cu and Zn, can be resolved and used to determine important properties of supernova remnants, galaxy clusters and accreting black holes and neutron stars.
Aims. The second stage of the present ongoing project involves the computation of the accurate photoabsorption and photoionisation cross sections required to interpret the X-ray spectra of such trace elements.
Methods. Depending on target complexity and computer tractability, ground-state cross sections are computed either with the close-coupling Breit-Pauli R-matrix method or with the autostructure atomic structure code in the isolated-resonance approximation. The intermediate-coupling scheme is used whenever possible. In order to determine a realistic K-edge behaviour for each species, both radiative and Auger dampings are taken into account, the latter being included in the R-matrix formalism by means of an optical potential.
Results. Photoabsorption and total and partial photoionisation cross sections are reported for isoelectronic sequences with electron numbers 3 ≤ N ≤ 11. The Na sequence (N = 11) is used to estimate the contributions from configurations with a 2s hole (i.e. [2s]μ) and those containing 3d orbitals, which will be crucial when considering sequences with N > 11.
Conclusions. It is found that the [2s]μ configurations must be included in the target representations of species with N ≥ 11 as they contribute significantly to the monotonic background of the cross section between the L and K edges. Configurations with 3d orbitals are important in rendering an accurate L edge, but they can be practically neglected in the K-edge region.
Key words: atomic data / X-rays: general
© ESO, 2016
1. Introduction
K lines from chemical elements with low cosmic abundances (trace elements) have been observed in the 5–9 keV energy range in spectra taken by the Chandra, XMM-Newton and Suzaku space observatories. Such surprising findings, in spite of trace abundances as low as a factor of 50 relative to Fe, speak for the high sensitivity of the satellite-borne instrumentation but also for puzzling abundance enhancements in some X-ray sources.
K emission lines from Cr and Mn, for instance, were first observed in supernova remnants with Chandra (Hwang et al. 2000) and then with Suzaku (Park et al. 2013) showing ionisation levels consistent with those of Fe, thus implying that these lines were emitted in the same spatial region. The Cr/Fe abundance ratio has been used by Badenes et al. (2008) as an indicator of progenitor metallicity and of the explosion mechanism. Trace-element lines have also been reported in the spectra of galaxy clusters (Tamura et al. 2009) and the Galactic centre (Nobukawa et al. 2010).
The spectra of compact X-ray sources, e.g. accreting black holes and neutron stars and their related warm absorbers, display absorption lines and edges – some of which are blue-shifted by gas outflows – that constrain the elemental abundances, ionisation structure, progenitor mass and velocity dispersion. An interesting example is the galactic black hole GRO J1655-40 where, due to line saturation, K absorption lines from trace elements (P, K, Cr, Mn and Co) impose more constrains on the flow properties than those from abundant species such as Fe, S, Si and Mg (Miller et al. 2006; Kallman et al. 2009). Also in the Chandra HETG spectrum of the galactic black hole GRS1915+105 by Ueda et al. (2009), an unusually broad range of K lines from Cu ions (with electron numbers 2 ≤ N ≤ 6) have been identified in what is believed to be a single zone of ionisation equilibrium.
The Soft X-ray Spectrometer (SXS) aboard the new Hitomi telescope1 will offer high sensitivity and spectral resolution (5 eV) in the 5–10 keV spectral range. It will enable the astrophysical exploitation of K lines and edge structures from trace elements in moderately bright sources with realistic observing times. The SXS energy resolution is expected to be adequate to distinguish between adjacent ionic stages of a K line or edge.
The quantitative interpretation of X-ray spectra requires the understanding of and reliable data for a whole variety of inner-shell processes such as K-shell photoabsorption and electron impact excitation, radiative and Auger decay of K-vacancy levels, K-edge resonance structures and Kα fluorescence emission. Although photoionisation cross sections based on central field approximations are currently available for all subshells of most ions of interest (e.g. Verner et al. 1996), they are not accurate enough as they do not take into account resonances and damping (radiative and Auger) effects that dominate the K-edge behaviour (edge smearing) and location.
As the first stage of a project to compute the required inner-shell data for the spectral modelling of trace elements in X-ray sources, Palmeri et al. (2012, hereafter Paper I) have recently reported – for the complete isonuclear sequences of F, Na, P, Cl, K, Sc, Ti, V, Cr, Mn, Co, Cu and Zn – energies, radiative and Auger widths and fluorescence yields for K-vacancy levels; and wavelengths and radiative A-values for their transitions to valence levels. These data were mainly calculated with the Hartree-Fock method with relativistic corrections (hfr) of Cowan (1981), but data accuracy was assessed with comparisons with the multiconfiguration Breit-Pauli (autostructure, Eissner et al. 1974; Badnell 2011) and the multiconfiguration Dirac-Fock (grasp, Grant et al. 1980; McKenzie et al. 1980) methods.
As a followup, we are currently involved in the extensive computations of the ground-state photoabsorption and photoionisation cross sections of the all the charge states of the aforementioned isonuclear sequences. In this paper we delineate the strategy and methods to accomplish this task, and report cross sections for the simpler systems with electron numbers 3 ≤ N ≤ 11. We conveniently draw from our earlier work on K photoabsorption of cosmic abundant species (Palmeri et al. 2002, 2005; Bautista et al. 2004; García et al. 2005, 2009; Witthoeft et al. 2009, 2011a,b, 2013; Gorczyca et al. 2013; Hasoğlu et al. 2014). The final data sets of the present project are destined to improve the atomic database of the xstar photoionisation modelling code (Bautista & Kallman 2001).
2. Numerical methods
The main problem in rendering the inner-shell photoabsorption and photoionisation cross sections of the ionic systems considered in the present project is computer tractability, as target representations with electron number N> 11 quickly lead to unmanageably large close-coupling calculations. Consequently, some compromises have to be made within the context of the multichannel approach (e.g. LS coupling vs. intermediate coupling), or altogether downsize, in spite of its shortcomings (Badnell et al. 1993), to the isolated-resonance approximation as extensively implemented to treat inner-shell opacity contributions in the Opacity Project (Badnell et al. 2005). These issues have been previously discussed by Witthoeft et al. (2009, 2011a,b), and therefore they will only be briefly reviewed here.
In the close-coupling framework, photoabsorption and photoionisation cross sections are computed with the relativistic (Breit-Pauli) R-matrix method (Berrington et al. 1995), using the rmatrx1 package in the inner region and stgbf0damp7 in the asymptotic region. The latter allows through an optical potential (Gorczyca & Badnell 1996, 2000) for the treatment of spectator-Auger damping effects, which are dominant in the near K-edge region. In the isolated-resonance approximation, the cross sections are calculated with the Breit-Pauli atomic structure code autostructure (Eissner et al. 1974; Badnell 2011).
For each isonuclear series (identified with the atomic number Z), photoabsorption and total and partial photoionisation cross sections are computed for members with 3 ≤ N ≤ Z (the difficult near-neutral members with Z−1 ≤ N ≤ Z will be considered in a more advanced project). We adopt the following computational strategy along isoelectronic sequences.
-
3 ≤ N ≤ 11. Calculations are performed in intermediate coupling with rmatrx1 using the target representations (level energies and Auger widths) listed in Table 9 of Paper I. Target wave functions are computed with the autostructure code.
-
12 ≤ N ≤ 20. Depending on computer tractability, calculations are performed in intermediate coupling with rmatrx1; otherwise LS coupling is implemented but including the non-fine-structure mass–velocity and Darwin interactions. Regarding target representations, two issues are hereby examined with respect to the atomic models tabulated in Paper I: the impact on the background cross section and edge structures due the exclusion of configurations with a 2s hole, i.e. [2s]μ where μ stands for all occupied valence subshells, and those containing 3d orbitals.
-
21 ≤ N ≤ Z−2. Whenever possible, rmatrx1 will be used otherwise we will recur to the isolated-resonance scheme with autostructure.
3. Results
Photoabsorption and total and partial (leaving the daughter ion in the ground and excited levels) photoionisation cross sections have been computed for the ground state and ground-configuration metastable states of ions of F, Na, P, Cl, K, Sc, Ti, V, Cr, Mn, Co, Cu and Zn with electron numbers 3 ≤ N ≤ 11. Target representations (level energies and Auger widths) are taken from Table 9 of Paper I. Although we are mainly interested in the K-edge region, due to the general absence of previous intermediate-coupling multichannel cross sections for these species, we have made an effort to also resolve the L resonances; thus, the cross-section energy tabulations start right from the reaction threshold with a fine mesh step of E/z2 = 0.0001 Ryd where z = Z−N + 1.
We plot in Fig. 1 the K-edge-region photoabsorption cross section of the 2p61S0 ground state of different members of the Ne isoelectronic sequence. Cross sections have been convolved with a Gaussian of width ΔE/E = 0.001 comparable to the SXS spectral resolution of 5 eV (0.4 Ryd). It may be seen that the resonance energy intervals in the K-edge structure increase with Z, but even at low Z and this resolution, the Kβ line can be resolved from the edge structure. This point is further illustrated in Fig. 2 where the closely spaced K lines for phosphorus ions (i.e. systems with a fairly low Z) with electron numbers 7 ≤ N ≤ 10 are shown, the cross sections being again convolved with a Gaussian of the same width. It may be appreciated that the energy intervals in between the Kα lines from the adjacent N-, O- and F-like ions and the Kβ line from the Ne-like ion are larger than 0.4 Ryd, and can therefore be resolved by the SXS.
![]() |
Fig. 1 Photoabsorption cross section in the K-edge region of the 2p61S0 ground state of members of the Ne isoelectronic sequence. The cross sections have been convolved with a Gaussian of width ΔE/E = 0.001. |
![]() |
Fig. 2 Photoabsorption cross section in the K-edge region of the ground state of four P ions. Black curve: P ix (N-like). Red curve: P viii (O-like). Green curve: P vii (F-like). Blue curve: P vi (Ne-like). The cross sections have been convolved with a Gaussian of width ΔE/E = 0.001. |
In the ion representations listed in Table 9 of Paper I for isoelectronic sequences with N ≥ 10, the [2s]μ configurations and those containing 3d orbitals have been excluded. If these target representations are adequate in the present photoionisation computations, it would open the possibility of using the Breit-Pauli R-matrix method in intermediate coupling up to at least N ≤ 20. Therefore, the contributions from these configurations must be evaluated before they are discarded altogether, if at all, and the relatively simple Na isoelectronic sequence presents a convenient way to carry out these tests. We thus consider the following three Ne-like target approximations.
-
Target A:
it is represented by the configurations2p6, [2p] 3s, [2p] 3p, [1s] 3s and [1s] 3p.
-
Target B:
it includes the configurations of Target A plus [2s] 3s and [2s] 3p.
-
Target C:
it includes the configurations of Target B plus [2p] 3d, [2s] 3d and [1s] 3d.
Target A is the one listed in Table 9 of Paper I while Target B additionally includes the [2s] 3s and [2s] 3p configurations with a 2s hole. Target C takes all the possible inner-shell single promotions to the 3ℓ orbitals with ℓ ≤ 2. In Fig. 3 the photoabsorption cross sections of the 2p63s 2S1/2 ground state of P v (Na-like) computed with Target A and Target B are depicted. In order to bring out the main features, the cross sections are convolved with a low resolution (ΔE/E = 0.01) Gaussian profile. It may be appreciated in Fig. 3a that, while the structure of the K edge is similar between the two approximations, Target B gives a significantly lower gradient in the monotonic decrease of the cross section between the L (~15 Ryd) and the K (~160 Ryd) edges. As shown in Fig. 3b, there are also some noticeable differences in the L-edge resonance structure around the spectral head at ~14−15 Ryd. On the other hand, it is found that Target C leads to conspicuous differences with respect to Target B only in the L-edge resonance structure as shown in Fig. 4, where it may be seen that there are significant shifts in the resonance positions and disimilar resonance profile characteristics (width and height); however, we find that an acceptable level of convergence has been attained for this system with Target C (single promotions).
![]() |
Fig. 3 Photoabsorption cross section of the 2p63s 2S1/2 ground state of Na-like P v computed with Target B (black continuous curve) and Target A (red dashed curve). a) Energy region between the L and K edges. b) L-edge region. The cross sections have been convolved with a Gaussian of width ΔE/E = 0.01. |
![]() |
Fig. 4 Photoabsorption cross sections of the 2p63s 2S1/2 ground state of Na-like P v in the L-edge region computed with Target B (black continuous curve) and Target C (blue dashed curve). The cross sections have been convolved with a Gaussian of width ΔE/E = 0.01. |
The photoabsorption and photoionizaton cross sections of the 2p63s 2S1/2 ground state of Na-like Zn xx are plotted in Fig. 5 to illustrate the importance of damping effects. In Fig. 5a damping is excluded and, therefore, the photoabsorption and photoionisation cross sections are alike. Damping is taken into account in Fig. 5b, and as a consequence, the resonances in the photoabsorption cross section show symmetric profiles of constant widths that causes the smearing of the K edge; however, the K edge remains sharp in the photoionisation cross section although most of the resonances are distinctively damped. As shown in Fig. 5c, the position and smearing of the K edge in the photoabsorption cross section depends on the convolution width, in this case a Gaussian of width ΔE/E = 0.001; moreover, the K edge in the photoionisation cross section is now slightly slanted.
![]() |
Fig. 5 Photoabsorption (black continuous curve) and photoionisation (red dashed curve) of the 2p63s 2S1/2 ground state of Na-like Zn xx. a) Undamped cross section. b) Damped cross sections. c) Damped cross sections convolved with a Gaussian of width ΔE/E = 0.001. |
4. Conclusions
-
1.
Photoabsorption and photoionisation (total and partial) crosssections for the ions of F, Na, P, Cl, K, Sc, Ti, V, Cr, Mn, Co, Cu andZn with electron numbers3 ≤ N ≤ 11 have been computed with the Breit-Pauli R-matrix method in intermediate coupling. Both the L and K resonance structures are well resolved with a fine energy mesh. Damping effects are conspicuous in the K-edge resonance structure.
-
2.
It has been shown that, with the SXS nominal resolution of ΔE = 5 eV (0.4 Ryd), it is possible to resolve the Kα and Kβ of these trace elements, even those with low Z (e.g. Z ~ 15).
-
3.
Even though the ion representations listed in Paper I are adequate targets for isoelectronic sequences with N ≤ 10, it has been shown here that for N> 10 the [2s]μ configurations, which were not considered in Paper I, must be additionally included in the close-coupling expansion as they lead to significant contributions to the background cross section in the energy region between the L and K edges.
-
4.
The inclusion of 3d configurations in the target representations of Na-like ions (N = 11) mainly leads to noticeable contributions in the resonance structure of the L edge, but the K edge and the background cross sections are hardly modified. However, they certainly lead to strong configuration interaction among n = 3 valence states, but calculations are expected to grow very rapidly with N making intermediate-coupling R-matrix calculations very large. This implies that the photoionisation cross sections in the L- and K-edge regions of ions with N> 11 will most surely have to be computed separately with different target approximations.
-
5.
Although cross sections for particular species can be obtained directly from the authors, they will only be published as a whole when the data sets are eventually completed. This will avoid data duplication in case there is need to rerun any isoelectronic sequence.
Acknowledgments
This project is sponsored by the NASA grant 12-APRA12-0070 through the Astrophysics Research and Analysis Program.
References
- Badenes, C., Hughes, J. P., Cassam-Chenaï, G., & Bravo, E. 2008, ApJ, 680, 1149 [NASA ADS] [CrossRef] [Google Scholar]
- Badnell, N. R. 2011, Comp. Phys. Comm., 182, 1528 [Google Scholar]
- Badnell, N. R., Griffin, D. C., Gorczyca, T. W., & Pindzola, M. S. 1993, Phys. Rev. A, 48, 2519 [NASA ADS] [CrossRef] [Google Scholar]
- Badnell, N. R., Bautista, M. A., Butler, K., et al. 2005, MNRAS, 360, 458 [NASA ADS] [CrossRef] [Google Scholar]
- Bautista, M. A., & Kallman, T. R. 2001, ApJS, 134, 139 [NASA ADS] [CrossRef] [Google Scholar]
- Bautista, M. A., Mendoza, C., Kallman, T. R., & Palmeri, P. 2004, A&A, 418, 1171 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Berrington, K. A., Eissner, W. B., & Norrington, P. H. 1995, Comp. Phys. Comm., 92, 290 [Google Scholar]
- Cowan, R. D. 1981, The theory of atomic structure and spectra (Berkley: University of California Press) [Google Scholar]
- Eissner, W., Jones, M., & Nussbaumer, H. 1974, Comp. Phys. Comm., 8, 270 [Google Scholar]
- García, J., Mendoza, C., Bautista, M. A., et al. 2005, ApJS, 158, 68 [NASA ADS] [CrossRef] [Google Scholar]
- García, J., Kallman, T. R., Witthoeft, M., et al. 2009, ApJS, 185, 477 [NASA ADS] [CrossRef] [Google Scholar]
- Gorczyca, T. W., & Badnell, N. R. 1996, J. Phys. B Atom. Mol. Phys., 29, L283 [Google Scholar]
- Gorczyca, T. W., & Badnell, N. R. 2000, J. Phys. B Atom. Mol. Phys., 33, 2511 [Google Scholar]
- Gorczyca, T. W., Bautista, M. A., Hasoglu, M. F., et al. 2013, ApJ, 779, 78 [NASA ADS] [CrossRef] [Google Scholar]
- Grant, I. P., McKenzie, B. J., Norrington, P. H., Mayers, D. F., & Pyper, N. C. 1980, Comp. Phys. Comm., 21, 207 [Google Scholar]
- Hasoğlu, M. F., Abdel-Naby, S. A., Gatuzz, E., et al. 2014, ApJS, 214, 8 [NASA ADS] [CrossRef] [Google Scholar]
- Hwang, U., Petre, R., & Hughes, J. P. 2000, ApJ, 532, 970 [NASA ADS] [CrossRef] [Google Scholar]
- Kallman, T. R., Bautista, M. A., Goriely, S., et al. 2009, ApJ, 701, 865 [NASA ADS] [CrossRef] [Google Scholar]
- McKenzie, B. J., Grant, I. P., & Norrington, P. H. 1980, Comp. Phys. Comm., 21, 233 [Google Scholar]
- Miller, J. M., Raymond, J., Fabian, A., et al. 2006, Nature, 441, 953 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Nobukawa, M., Koyama, K., Tsuru, T. G., Ryu, S. G., & Tatischeff, V. 2010, PASJ, 62, 423 [NASA ADS] [Google Scholar]
- Palmeri, P., Mendoza, C., Kallman, T. R., & Bautista, M. A. 2002, ApJ, 577, L119 [NASA ADS] [CrossRef] [Google Scholar]
- Palmeri, P., Kallman, T. R., Mendoza, C., Bautista, M. A., & Krolik, J. H. 2005, in X-ray Diagnostics of Astrophysical Plasmas: Theory, Experiment, and Observation, ed. R. Smith, AIP Conf. Ser., 774, 75 [Google Scholar]
- Palmeri, P., Quinet, P., Mendoza, C., et al. 2012, A&A, 543, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Park, S., Badenes, C., Mori, K., et al. 2013, ApJ, 767, L10 [NASA ADS] [CrossRef] [Google Scholar]
- Tamura, T., Maeda, Y., Mitsuda, K., et al. 2009, ApJ, 705, L62 [NASA ADS] [CrossRef] [Google Scholar]
- Ueda, Y., Yamaoka, K., & Remillard, R. 2009, ApJ, 695, 888 [NASA ADS] [CrossRef] [Google Scholar]
- Verner, D. A., Ferland, G. J., Korista, K. T., & Yakovlev, D. G. 1996, ApJ, 465, 487 [NASA ADS] [CrossRef] [Google Scholar]
- Witthoeft, M. C., Bautista, M. A., Mendoza, C., et al. 2009, ApJS, 182, 127 [NASA ADS] [CrossRef] [Google Scholar]
- Witthoeft, M. C., Bautista, M. A., García, J., et al. 2011a, ApJS, 196, 7 [NASA ADS] [CrossRef] [Google Scholar]
- Witthoeft, M. C., García, J., Kallman, T. R., et al. 2011b, ApJS, 192, 7 [NASA ADS] [CrossRef] [Google Scholar]
- Witthoeft, M. C., Bautista, M. A., García, J., et al. 2013, Atom. Data Nucl. Data, 99, 53 [Google Scholar]
All Figures
![]() |
Fig. 1 Photoabsorption cross section in the K-edge region of the 2p61S0 ground state of members of the Ne isoelectronic sequence. The cross sections have been convolved with a Gaussian of width ΔE/E = 0.001. |
In the text |
![]() |
Fig. 2 Photoabsorption cross section in the K-edge region of the ground state of four P ions. Black curve: P ix (N-like). Red curve: P viii (O-like). Green curve: P vii (F-like). Blue curve: P vi (Ne-like). The cross sections have been convolved with a Gaussian of width ΔE/E = 0.001. |
In the text |
![]() |
Fig. 3 Photoabsorption cross section of the 2p63s 2S1/2 ground state of Na-like P v computed with Target B (black continuous curve) and Target A (red dashed curve). a) Energy region between the L and K edges. b) L-edge region. The cross sections have been convolved with a Gaussian of width ΔE/E = 0.01. |
In the text |
![]() |
Fig. 4 Photoabsorption cross sections of the 2p63s 2S1/2 ground state of Na-like P v in the L-edge region computed with Target B (black continuous curve) and Target C (blue dashed curve). The cross sections have been convolved with a Gaussian of width ΔE/E = 0.01. |
In the text |
![]() |
Fig. 5 Photoabsorption (black continuous curve) and photoionisation (red dashed curve) of the 2p63s 2S1/2 ground state of Na-like Zn xx. a) Undamped cross section. b) Damped cross sections. c) Damped cross sections convolved with a Gaussian of width ΔE/E = 0.001. |
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.