1. Introduction: NORAD-Atomic-Data
The atomic database, NORAD-Atomic-Data [
1,
2], contains parameters of atomic processes, as described below, that are dominant in astrophysical plasma. A major part of the data were obtained under the two international collaborations of the Opacity Project (OP) [
3,
4] and the Iron Project (IP) [
5], and constitute the atomic data source, either for new data or revised and more expanded data for the OP and the IP.
The author (SNN) is the leading member of the Opacity Project ([
4]) and the Iron Project (IP [
5]) for studying the radiative atomic processes and has been computing extensive sets of radiative data. To meet the need of astronomers accessing atomic data all time from an independent platform the online database NORAD-Atomic-Data [
1,
2] was founded in 2007 at the Ohio State University. NORAD-Atomic-Data provides data which are not available or which are updated with larger wavefunction expansion and higher accuracy than those in the two atomic databases, TOPbase [
7] founded under the OP and TIPbase [
8] founded under the IP. All data are processed for clarity and easy applications before uploading them to NORAD-Atomic-Data as explained in [
1]. This report presents several recently added enhancements and under development in the latest form of the database.
Details of format and data structure, and summary of theoretical background, accuracy estimation of data for the processes have been described in Nahar [
1]. More detailed theoretical background are available from Pradhan and Nahar in textbook "Atomic Astrophysics and Spectroscopy" [
9] which is a culmination of findings in atomic processes revealed in the study under the OP and IP. Here we present summary outlines for experimental methods, and more details of one particular experimental method with Advanced Light Source (ALS) at Lawrence Berkeley National Laboratory (LBNL) synchrotron for measurement of high-resolution photoionization cross sections.
We provide below outlines for general guidance on the atomic processes in plasma and relevant contents of the database.
2. Atomic Processes and Parameters in NORAD-Atomic-Data
There are four atomic process that occur commonly in astrophysical plasma - photo-excitation,photoionization, electron-ion recombination, electron-impact excitation. Study of an atomic process requires wavefunctions and energy levels of the atomic system being studied to calculate various parameters of the processes. These parameters at database NORAD-Atomic-Data have been computed using energy levels and wavefunctions generated either by the R-matrix method using its package of and associated codes [
10,
11,
12,
13] or by central-field approximation using atomic structure code SUPERSTRUCTURE (SS) [
14,
15]. Followings are the atomic processes and associated data files.
-
i)
Energy levels:
NORAD-Atomic-Data provides energy table for each atomic species for which atomic processes have been studied. Although basically the same, there are some differences in representations between the energies obtained from the R-matrix method and from atomic structure calculations.
The R-matrix method calculates absolute values for the energy levels and can show the clear presence of ionization threshold as negative energies for bound levels turn positive for continuum levels when energies reach to the threshold. This option is enabled in the R-matrix codes for nonrelativistic LS coupling approximation [
10]. The latest version of the relativistic Breit-Pauli R-matrix (BPRM) codes [
13] computes only the bound levels, that is, no positive energies are computed. In contrast to R-matrix, SS computes energy values relative to the ground level which is set to zero, and hence show positive values. This representation is the same as that of for energies compiled and available at the website of the National Institute of Standards and Technology (NIST [
16]). Atomic structure calculations, such as of SS, do not determine the ionization threshold and computes all energy levels belonging to the set of configurations provided. Hence the computed energies can continue from bound to beyond the ionization threshold. NIST provides ionization thresholds obtained separately.
Energy levels computed in the R-matrix method [
10,
13] go from ground level to highly excited levels with
10,
9, and of parities even and odd, and hence correspond to a large set of levels. The number of energy levels from SS [
14,
15] is typically less since the excitation goes up to
4-6 although it includes all possible combination of quantum numbers.
R-matrix computations include all quantities, such as, oscillator strengths, photoionization cross sections, electron-ion recombination rate coefficients for all levels going up to n=10. SS computes all transition parameters for levels going up to the highest orbitals with n=4-6.
Formats for the energy files are described in [
1].
-
ii)
Radiative atomic transition of excitation and de-excitation:
is the ion
X of charge
q and asterisk indicates excitation. The parameters of interest are: i) Line strength (
S), ii) Oscillator strength (
f), iii) Radiative decay rate (
A-value), for the transitions and iv) Lifetimes (
) of levels.
Radiative data for bound-bound transitions are of the following types: allowed electric dipole (E1) of the same spin (for
LS-allowed and fine structure) and different (intercombination) spin multiplicities of the final state from the initial one, electric quadrupole (E2), electric octupole (E3), magnetic dipole (M1), and magnetic quadrupole (M2). Details of these are described in [
9] For many ions at the database, E1 transitions have been obtained from using BPRM method while parameters for forbidden transitions have been obtained from Breit-Pauli atomic structure calculations using code SS. When R-matrix method has not been used, E1 transitions obtained using SS are presented. Individual transitions form spectral lines. However data for large number of transitions are needed to determine the opacity in plasma, creation of synthetic spectrum of ions, etc.
Formats for the transition parameters are described in [
1].
-
iii)
Photoionization (PI):
The process proceeds directly by photoabsorption as
which produces smooth photoionization cross section (
) feature, such as the background
. Photoionization can occur indirectly via formation of an intermediate doubly-excited autoionizing state before ionization as
This takes place when the photon energy matches to that of a Rydberg state and introduces a resonance, known as a Rydberg resonance, in the cross section
. A resonance has considerable impact, not only by its magnitude of very high probability, but also on the shape of the background cross section. In addition to Rydberg series of resonances, Seaton resonance also forms when the core ion goes through a dipole allowed transition before breaking down for photoionization.
Central-field approximation, such as, distorted wave, can produce only the background cross sections, not the resonances. Only the close-coupling (CC) approximation using the R-matrix method can produce resonances automatically along with background cross sections. We use relativistic Breit-Pauli R-matrix (BPRM) codes [
13,
17] with CC wavefunction expansion to compute
. The resonances are treated with radiation damping effect [
18] for precise treatment. However, we also use nonrelativistic R-matrix codes [
10] when relativistic fine structure calculations become computationally prohibitive.
The parameters calculated for photoionization are cross sections for all bound levels from ground to 10. They are usually of two types. i) Partial cross section for photoionization leaving the core ion in its ground level and ii) total cross section for photoionization leaving the core ion in the ground and various excited states as the photon energy crosses the excitation thresholds. Total cross section is the sum of all partial cross sections for leaving the core ion in various states. iii) We also calculate the individual partial cross sections leaving the core ion in various excited levels. However, such partials are saved only for the ground level for which these cross sections have more applications.
Formats for photoionization cross sections are described in [
1]
-
iv)
Electron-ion recombination:
It is the inverse process of photoionization where electron-ion can recombine directly, known as radiative recombination (RR), as
This provides the smooth background feature of recombination cross section. Recombination can also occur through an intermediate recombination (DR), as
The intermediate autoionizing state can lead to autoionization (AI) where the electron goes free, or to radiative decay channel of DR by emission of a photon.
For electron-ion recombination, we calculate the atomic parameters using the unified method of Nahar and Pradhan [
17,
19] for: i) Level-specific total recombination (includes both RR and DR) rate coefficients (
) of all bound levels going up to
10, ii) total recombination rate coefficient with respect to temperature
) which is the summed contributions of all levels with
10 and of levels with 10
obtained using extension of Bell and Seaton theory [
17,
20] on DR, iii) total recombination cross sections (
) and iv) total recombination rate coefficient (
) with respect to photo-electron energy
E. Recombination resonances can be seen in emission spectra as dielectronic recombination (DES) lines and total
is used in determination of ionization fractions in photoionized or collisional plasma. Data format is explained in [
1].
-
v)
Electron-impact excitation (EIE):
Electron impact excitation of a target leads to emission of a photon as the excited electron drops down:
The line of the emitted photon carries information of the process. Similar to photoionization, an autoionizing state that introduces a resonance can form if the impact energy equals to that of a Rydberg state. EIE parameters are: i) Collision strength (
), ii) Effective collision strengths (
), iii) Collision rate coefficients (
). Emission lines from levels excited by EIE provides diagnostics of thin plasma. Details can be found in, e.g., [
9].
3. Merged Beam Technique at the Advanced Light Source at the Lawrence Berkeley National Laboratory
For a general concept of measurement of photoionization cross sections, we provide summary description of the experimental set-up that was installed at the Advanced Light Source (ALS) Synchrotron Facility at the LBNL. This is one most accurate experimental set-up with third generation synchrotron radiation for high-resolution photoionization cross section.
For a review about the merged beam technique (MBT), we refer to the article by Phaneuf
et al. [
21]. For more information about synchrotron technology we refer to the article by Margaritondo [
22]. For a detailed description of this experimental set-up, we refer to references [
23,
24]. Here, we present a brief but complete description of the method.
To implement the experimental realization of photoionization Eqs. (
2,
3 ), an ion beam and a photon beam are merged. Also, parameters, such as, their width, intensity and overlap length have to be measured.
3.1. The Photon Beam
The experimental set-up was located at the ALS, a third-generation synchrotron accelerator at the LBNL. The Merged Beam Technique (MBT) was implemented in one of the end-lines of the synchrotron ring (beam line 10.0.1.). The end-station was intended to be a multi-user experiment. This MBT end-station was decommissioned in 2015 and by the end of the MBT operational time, the synchrotron’s electron’s storage-ring operated with 500 mA current at 1.9 GeV in top-off mode, i. e., with a quasi-continuous injection of electrons that results in a constant current and stability improvement.
The photon-beam was generated by a 10 cm period undulator inserted in the electron’s trajectory inside the synchrotron’s ring. As a result, the relativistic electrons oscillate producing a photon-beam with a width less than 1 mm and with a divergence less than 0.05°. This photon-beam was then directed to shine over a grazing spherical-grating monochromator, producing wavelength separation as a function of the grazing angle. In this manner, it was possible to scan the photon energy by rotating the grating. Also, the exit slit of the monochromator and the undulator gap were adjusted to maximize the photon-beam intensity.
A calibrated silicon photodiode [
25] with a 5% uncertainty was used to measure the photon flux. This diode generated a current that was read by a precision current meter that provided a normalization signal to the data acquisition system. Collisions of the
ions with the residual gas in the vacuum system may generate
signal. To separate this background signal from the actual photoionization signal, the photon-beam was mechanically pulsed using a controlled chopper-wheel to produce a signal correlated to photons and a signal correlated only to background
, and in turn subtract it from the signal generated when the photons were on.
s is the order of photoionization. For instance,
s = 1 corresponds to single photoionization.
The energy resolution of the photon-beam was typically between 10 meV and 40 meV and it was to nominally programmed on the control system to later derive a more realistic value from FWHM (full width half maximum) of resonance peaks. The energy of the photon beam was nominally set by the control interface system associated to the beam line. The energy calibration was further improved with the aid of an ionization gas-cell located on a side branch of the experimental set-up, in this cell He [
26] and Kr [
27] gases were injected. Their ionization spectra were used as a reference in the energy interval from 21.218 eV to 63.355 eV. Estimates of the energy error from this technique, varied from about ± 30 meV to ± 10 meV, depending on different author’s implementations. For instance, this calibration was employed for the single photoionization spectrum measurement of Cl
+ [
28] resulting in an agreement with NIST database [
29] up to three decimal digits (in eV).
Although, the photon flux and resolution of this facility were uniquely fitted to measure photoionization, higher-order harmonic components in the photon beam, especially in the lower energy interval of its operational range (17 eV to 80 eV) can have considerable effects. This issue was assessed in references [
23,
25,
30].
3.2. The Ion Beam
To generate the ion beam, an initial stage to create the ions was implemented by installing an ion-source. Different types of ion sources were positioned in the end-station. For instance, an ECR (electron cyclotron resonance) ion source was used to generate multicharged ion beams. For example C
2+, N
3+ and O
4+ [
31]. In the initial stage of the experiment, a hot-filament ion-source was used to generate O
+ [
32] and Ne
+ [
33] ions; a caesium-sputter ion source was used for Na
− [
34] and a plasma emission source to produce Li
+ ions [
35].
The ion source was mounted on a floated voltage circuit, in a way, that the whole ion source was biased to the acceleration voltage. As a consequence, ions generated inside the ion source were repelled or extracted by Coulomb repulsion. The standard repulsion voltage was 5 kV or 6 kV to produce ion beams of energies of the order of keV because at these energies, the ionization by collision with the residual gas is expected to be low.
ions extracted from the ion source were accelerated toward a cylindrical electrostatic lens set that provided focusing to the ions. Following, the ions were selected by mass-to-charged analysis with a 60° sector magnet. It is necessary to have fine control of the ion beam direction, for this purpose, sets of electric fields perpendicular to the ion beam trajectory were implemented by applying small voltages to steering plate electrodes (not shown in
Figure 1) located at the entrance and the exit of a deflector made of two spherical-sector electrodes. This deflector produced a 90° change in the ion beam trajectory for the ion beam to face the counter-propagating photon beam from the ALS.
In this stage, the photon beam and the ion beam overlap, and as a consequence, photons can ionize the ion beam. Ions resulting from the overlapping will be called photo-ions. Photo-ions and the ion beam traveled through a voltage bias (typically, 2 kV) metallic mesh or interaction region. Photo-ions produced in this section, suffer a different electric field from those photo-ions produced outside, thereby defining a fixed interaction distance. In this manner, photo-ions produced in this region had a different energy and it was possible to separate them from the parent ion beam by a second analyzing magnet. The interaction region has entrance and exit apertures to help define an effective interaction length. The experiment operated under ultrahigh vacuum conditions.
Intensity distributions of both beams were measured by rotating-wire beam profilers installed before and after the interaction region. A third beam profiler consisted of a mechanically translating-slit scanner positioned halfway between the interaction region.
A second analyzing magnet separated the ion beam from the photo-ions. The parent ion beam was monitored by an elongated Faraday cup. The second magnetic field was adjusted such that the photo-ions passed through an aperture in the back of the Faraday cup. A spherical 90° electrostatic deflector directed them onto a stainless steel plate biased at -550 V, from which secondary electrons were accelerated and detected by a microsphere-plate electron multiplier used in a pulse-counting mode.
The detector’s planes were perpendicular to the beam trajectory, allowing the photo-ions be scanned across the detector, and verify full collection by measuring a plateau as a function of the pertinent parameters. The efficiency of the photo-ion detector was calibrated in situ using an averaging sub-femtoampere meter to record the photo-ion current, which was then compared to the measured photo-ion count rate. The primary ion beam current was measured by a precision current meter.
The overlap of the photon and ion beams, in addition to the interaction region dimensions are needed to measure the cross-section. For this purpose, the beam intensity distributions measured with the beam profilers,
of the ion beam and
of the photon beam, were measured, and a form factor
was derived according to,
where
z is the axis assigned to the propagation direction of the ion beam.
was sampled at the three positions where profilers were mounted: the entrance, the center, and the exit of the interaction region. These values,
, were used to derive
by interpolation along the total interaction region length and derive, by integration over
z, the spatial overlap of the photon and ion beams the interaction region path [
33].
The photoionization cross-section,
, was derived from the experimental parameters as follows,
where
R is the photo-ion count rate,
q is the charge state of
parent-ion
C,
is the ion beam velocity in
,
is the responsivity of the photo-diode (electrons per photon),
is the ion beam current (A), and
is the photo-diode current (A).
Past Reviews on Measured Photoionization of Ions
We present a list of references for photoionization cross sections for various ions measured at ALS in
Table 1. We discuss cross sections from other set-ups and the sources of information below.
There exists some reviews on the subject earlier measurements of photoionization cross sections. The review by Schmith in 1992 [
36], presents a comprehensive treatment of both theory and experimentation with synchrotron radiation up to that date. There is still no reference to the merged-beam technique. It is a handy reference for general aspects and terminology of first principles theoretical methods. The Sonntag and Zimmermann review in 1992 [
37] concentrates on PI experiments of metallic vapors with synchrotron radiation above 20 eV. It presents an interesting review of available techniques to that date, like He-discharge and lasers. The authors mentioned that R-matrix is a more convenient method to describe PI over extended energy regions.
In 2001 West’s review [
38] it is mentioned that there had been two main methods to measure photoionization spectra and cross sections, plasma discharge and merged beam technique. He brings forward that the first experimental measurements of photoionization cross sections were reported by Lyon and collaborators [
39]. Up to the date of the review, the majority of measurements are concentrated in spectra measured with the plasma discharge method. A description of the plasma ionization technique evolution is presented. In this review, the merged beam technique is described to a good extent. This review emphasizes cases where both methods, plasma ionization and the merged beams technique have been complementary. It raises the fact that the amount of experimental data available to date was not comparable with the data generated by theory, to a certain degree, this statement is still valid. By 2006, Kjeldsen [
40] covers the state of the art focused in MBT experimentation with synchrotron radiation. Finally, the photoionization of metal atoms is also reviewed by [
41] and relevant to astrophysics PI data is reviewed by Schippers and Müller [
42].
4. New Additions to NORAD-Atomic-Data
The database NORAD-Atomic-Data has been enhanced by additions of new data, as described below, of atomic processes.
1.) Among the data added, a signiicant amount of new data is related to spectral formation from photo-excitation cross sections of atoms and ions. There are of three types of such data. i) Atomic data for x-ray spectroscopy, especially for K-
and K-
transitions for resonance fluorescence effect predicted by Pradhan et al [
43] for biomedical applications and was later observed experimentally [
44] as shown in Figure 1. ii) Photoiexcitation data for spectral formation by ions of bio-signature elements, such as phosphorus, e.g. [
46]. iii) Atomic data for spectroscopy of lanthanides of neutral to low ionization stages which may relate to the electromagnetic waves created during a kilonovae event of merger of two neutron stars or black holes, e.g. [
47].
Figure 2.
K
α resonance fluorescence in Al plasma [
44]. Top: Theoretically computed K
α absorption features in Al ions isoelectronic with F to He are shown (ionization states are labeled in the top panel and roman numerals in the bottom panel). All contributing K
α transition strengths are added together and shown in terms of X-ray attenuation coefficients (cm2/g) (convoluted over a Gaussian FWHM of 10 eV). Bottom: Experimental measurements from the LCLS-XFEL [
45] (reproduced with permission). The dashed line on the right shows the resonance transitions with increasing K
α emission intensity corresponding to the theoretical absorption complexes. The observed K
α features appear when the XFEL energy equals the emitted photon energy.
Figure 2.
K
α resonance fluorescence in Al plasma [
44]. Top: Theoretically computed K
α absorption features in Al ions isoelectronic with F to He are shown (ionization states are labeled in the top panel and roman numerals in the bottom panel). All contributing K
α transition strengths are added together and shown in terms of X-ray attenuation coefficients (cm2/g) (convoluted over a Gaussian FWHM of 10 eV). Bottom: Experimental measurements from the LCLS-XFEL [
45] (reproduced with permission). The dashed line on the right shows the resonance transitions with increasing K
α emission intensity corresponding to the theoretical absorption complexes. The observed K
α features appear when the XFEL energy equals the emitted photon energy.
2.) The second addition is inclusion of collisional data for collision strengths and rate coefficients for electron impact excitation (EIE). NORAD-Atomic-Data considered almost no collisional data in the past since these data for various atomic systems are available from TIPbase [
8] and from other investigators. However, for possible plasma diagnostics, NORAD-Atomic-Data has started to include these collisional data, e.g. for Ca IV [
48].
3.) The the third addition is the inclusion of experimental data for energy levels, oscillator strengths. NIST database [
16] is the main source for compiled and evaluated experimental oscillator strengths and radiative decay rates. We have added large sets of experimental energy levels and transition parameters of a number of heavy ions, obtained mainly from the leading spectroscopy laboratory of Aligarh Muslim University in India. These files contain complete sets of transitions that were measured in comparison to the evaluated sets available at NIST. The added data are currently for eight ions, Ag III [
49], Ag-IV [
50], Au V [
55], Bi IV [
52], Hg III [
53], Hg-IV [
54], Hg V [
55], Hg-VI [
55]. For the transition parameters, they used vacuum spectrograph for measurement of spectral lines from exposures created by discharge of electrodes, containing the sample, with a trigger module of high voltage electrical pulse. Details of their experiments, calibration and analysis are given in their papers.
Aligarh lab has produced spectral data for many ions. However, currently available data are largely in pdf-formatted files which were saved with publications. More data will be added to NORAD-Atomic-Data when the authors provide the data files converted to ASCII formatted files.
4.) The fourth addition is the inclusion of experimental data for photoionization cross sections. Most experimental set-ups are typically synchrotron based, such as, Advanced Light Source (ALS) of synchrotron at Lawrence Berkeley National Lab (LBNL), set-ups at University of Paris-SUD, Aahus University, BESSY and PETRA in Germany. Another set-up used is laser-induced excited state photoionization in Pakistan and University of Nebraska. Among the experimental photoionization cross section data for ions available at NORAD-Atomic-Data are Cl III [
56], Cl II [
57], K [
58], etc.
5.) The fifth addition of NORAD-Atomic-Data is a web-based interactive feature to obtain line ratios for plasma diagnostics. NORAD-Atomic-Data has all the basic radiative, except for the collisional data for EIE, to calculate the line ratio diagnostics for many atoms and ions. Our recent study on line ratios of P II [
59], with inclusion of electron-ion recombination cascade effect which was known to have impact not considered in the astrophysical models before, demonstrated that the effect can lower or increase the line ratios.
Figure 3 presents an exemplary impact of electron-ion recombination on line ratios (red) curve compared to that without the effect (dark red) on the low lying transitions in P II [
59].
Figure 3.
Emissivity ratios of [P II] lines at ionospheric and nebular temperatures 500K, 1000K, and 10,000K lying within and potentially observed by by JWST spectrometers. While 500K and 1000 K lines (dark red) appear to be identical, they differ by a few percentage points at all, far less than the difference from the 10,000 K line [
59].
Figure 3.
Emissivity ratios of [P II] lines at ionospheric and nebular temperatures 500K, 1000K, and 10,000K lying within and potentially observed by by JWST spectrometers. While 500K and 1000 K lines (dark red) appear to be identical, they differ by a few percentage points at all, far less than the difference from the 10,000 K line [
59].
We have computed the relevant atomic quantities in large sets that can used to compute the line ratios for diagnostics over a large range temperatures and densities. A preliminary version of the code is in operaturation at the the web-based interactive feature for computing the line ratios at any temperature and density.
5. Conclusions
An overview with general outlines of contents on the new additions to database NORAD-Atomic-Data is presented. There are six additions that have been discussed for the expansion of the database. Significant enhancements of data are due to additions of experimental data, not included before, and data for photo-absorption spectral features.
Testing of data for accuracy is given a high priority. Majority of data have been obtained from ab initio R-matrix method. Results using atomic structure calculations have also been added for cases they provide good accuracy, such as for x-ray or high energy region where perturbations and electron-electron interactions are relatively negligible, and for cases, such as large data need for large atomic systems of lanthanides where numerical extends are prohibitive for R-matrix codes.
While theoretical overview along with description of data, accuracy and benchmarking of the online database were presented in the earlier paper on the database [
1], here we present general overview of experimental background, particularly of photoionization.
Future direction are:
i) include more data, both radiative and collisional, with publications
ii) include computer programs that can read the data files and process the data to calculate the quantities of interest
iii) add more data for the web-based feature to compute spectral line ratios of other atoms and ions for astrophysical applications,
iv) an option for selection of partial data, such as, for a particular wavelength range