Next Article in Journal
Shape and Satellite Studies of Highly Charged Ions X-ray Spectra Using Bayesian Methods
Previous Article in Journal
Cyanide Molecular Laser-Induced Breakdown Spectroscopy with Current Databases
Previous Article in Special Issue
Ideas and Tools for Error Detection in Opacity Databases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Atomic Data Assessment with PyNeb: Radiative and Electron Impact Excitation Rates for [Fe ii] and [Fe iii]

by
Claudio Mendoza
1,*,
José E. Méndez-Delgado
2,
Manuel Bautista
3,
Jorge García-Rojas
4,5 and
Christophe Morisset
6
1
Physics Center, Venezuelan Institute for Scientific Research (IVIC), Caracas 1020, Venezuela
2
Astronomisches Rechen-Institut, Zentrum für Astronomie der Universität Heidelberg, D-69120 Heidelberg, Germany
3
Department of Physics, Western Michigan University, Kalamazoo, MI 49008, USA
4
Instituto de Astrofísica de Canarias (IAC), E-38205 La Laguna, Spain
5
Departamento de Astrofísica, Universidad de La Laguna, E-38206 La Laguna, Spain
6
Instituto de Astronomía, Universidad Nacional Autónoma de México, Ensenada 22860, Mexico
*
Author to whom correspondence should be addressed.
Submission received: 6 February 2023 / Revised: 14 March 2023 / Accepted: 23 March 2023 / Published: 1 April 2023
(This article belongs to the Special Issue Development and Perspectives of Atomic and Molecular Databases)

Abstract

:
We use the PyNeb 1.1.16 Python package to evaluate the atomic datasets available for the spectral modeling of [Fe ii] and [Fe iii], which list level energies, A-values, and effective collision strengths. Most datasets are reconstructed from the sources, and new ones are incorporated to be compared with observed and measured benchmarks. For [Fe iii], we arrive at conclusive results that allow us to select the default datasets, while for [Fe ii], the conspicuous temperature dependency on the collisional data becomes a deterrent. This dependency is mainly due to the singularly low critical density of the 3 d 7 a 4 F 9 / 2 metastable level that strongly depends on both the radiative and collisional data, although the level populating by fluorescence pumping from the stellar continuum cannot be ruled out. A new version of PyNeb (1.1.17) is released containing the evaluated datasets.

1. Introduction

PyNeb1 [1,2] is a widely used Python tool for the analysis of emission lines in nebular plasmas to determine the temperature, density, and chemical abundances [3]. It essentially solves the equations of statistical equilibrium by addressing an extensive atomic database of radiative and collisional rates; therefore, PyNeb’s diagnostic potential relies directly on the accuracy of this database. In this respect, its data curation strategy is based on the selection of recommended default datasets and keeping access to a historical archive of radiative and collisional files rather than discarding them. This singular scheme allows modelers to estimate the uncertainties of the plasma diagnostics and chemical abundances due to the scatter of the atomic parameters. It also provides a versatile environment for atomic data assessment, which has been previously used productively to benchmark the accuracy of the radiative rates in N- and P-like ions and effective collision strengths in C-like ions [4] and to determine the atomic-data impact on density diagnostics [5].
As a follow-up to [4], we apply PyNeb’s capabilities to assess the radiative and collisional rates associated with the transition arrays of [Fe ii] and [Fe iii]. These ions give rise to forbidden lines in the optical and infrared spectra of H ii regions, planetary nebulae, and quasars that are used to determine plasma conditions and iron abundance, the latter particularly in the context of depletion factors in the life cycle of dust grains [6,7,8]. These iron transition arrays are also useful for studying the shocked plasma caused by outflows from newborn stars that form bright Herbig–Haro (HH) objects perturbing stellar accretion [9,10,11]. The high-resolution spectra of such HH objects—e.g., HH 202S [9] and HH 204 [11] in the Orion Nebula—taken with the Ultraviolet Visual Echelle Spectrograph (range 3100–10,400 Å) [12] present an observational opportunity to benchmark the Fe ii and Fe iii atomic data. The optical emission spectra of these two HH objects are dominated by photoionization from θ 1 Ori C, the main ionization source of the Orion Nebula, and can therefore be treated as small-scale H ii regions [10].
Due to the astrophysical importance of Fe ii, the computation of radiative and collisional rates for the spectral modeling of this system has received considerable attention since the early semi-empirical computations of A-values (radiative transition probabilities) for the forbidden [13] and dipole-allowed [14] transitions and electron-impact collision strengths with the distorted wave method [15,16]. Our intention here is not to evaluate this vast body of data but to concentrate on complete and accurate datasets to model the spectra of [Fe ii] and [Fe iii] observed in nebular plasmas.
A landmark systematic effort to compute radiative and collisional data for the fine-structure spectral modeling of Fe ions was undertaken by the Iron Project2 using the R-matrix method [17] and relativistic multi-configuration structure codes such as superstructure [18], civ3 [19], and hfr  [20]. For Fe ii and Fe iii, A-values were computed for both allowed and forbidden lines [21,22,23] and electron-impact effective collision strengths (ECS) for infrared, optical, and ultraviolet transitions [24,25,26]. Although these R-matrix calculations were originally performed in L S coupling followed by algebraic transformations to intermediate coupling, more recent scattering calculations for these ionic species [27,28,29,30,31] have relied on R-matrix implementations based on the Breit–Pauli Hamiltonian (bprm [32], bsr [33]) and Dirac–Fock Hamiltonian (darc [34]). In some cases, e.g., the collisional excitation of levels within the Fe + ground configuration, low-temperature ECS have been computed with several R-matrix implementations—intermediate coupling frame transformation (icft, [35]), bprm, and darc—to obtain uncertainty estimates [36].
Extensive multi-configuration Breit–Pauli (mcbp) calculations of g f -values (weighted oscillator strengths) and A-values for allowed and forbidden transitions in Fe ii and Fe iii [37,38,39,40,41,42,43] with the code civ3 have shown that fine-tuning corrections to the Hamiltonian matrix by fitting to the experimental level energies lead to more accurate configuration expansions. However, significant A-value differences (greater than 20–30%) are found for some transitions computed with comparable multi-configuration methods that also implement fine-tuning corrections (e.g., [23,44,45]), which are caused by strong configuration mixing, double-excitation representations, and subtle core–valence correlation involving the closed 3p 6 and open 3d n subshells.
Since the Fe ii radiative data display a large scatter, experimentalists and astronomical observers have been enticed to measure accurate radiative parameters such as lifetimes and branching fractions to determine A- and g f -values; however, most of them address allowed transitions involving odd-parity states. Lifetimes have been measured for 21 Fe ii levels with a time-resolved nonlinear laser-induced fluorescence technique to an accuracy of 2–3%, and by means of branching fractions, absolute A-values were reported to be within 6% and 26% for the strong and weak transitions, respectively [46]. Accurate g f -values for 142 Fe ii lines were determined from laboratory and solar spectroscopic multiplets and theoretical fine-structure ratios, assuring that oscillator-strength error no longer hinders stellar abundance studies [47]. Branching fractions have been measured for 121 ultraviolet Fe ii lines (2250–3280 Å), and through published lifetimes, g f -values have been derived and applied to iron abundance determinations in the Sun and the metal-poor star HD 84937 [48]. Good agreement is obtained with the Fe standard solar abundance. On the other hand, the only lifetime measurements of [Fe ii] even-parity metastable states to derive transition probabilities have been carried out by the Ferrum Project using a laser probing technique on a stored ion beam, the branching fractions being derived from astrophysical spectra [49,50,51]. Consistent differences between the observed intensity ratios and theoretical branching fractions have suggested that the [Fe ii] lines might be used to determine reddening.
In the present data assessment, we revise the Fe ii and Fe iii atomic datasets—namely, energy levels, A-values, and ECS—available in PyNeb 1.1.16. In most cases, we reconstruct the atomic models from the sources by mapping the rates to the NIST level structure. We renormalize the number of levels to a predefined value to enable one-to-one comparisons among the theoretical datasets and benchmarks with the HH 202S and HH 204 spectra and the aforementioned experimental data. We also extend the PyNeb database with additional new datasets, the main objective being the default selection for a new version: PyNeb 1.1.17.
We give a brief description of PyNeb’s file structure and methods in Section 2, and carry out the data assessments and discussions of Fe iii and Fe ii in Section 3 and Section 4 followed by the conclusions in Section 5.

2. PyNeb

A brief description of the PyNeb package is given in [4], in particular of its object-oriented architecture. The ion object is thereby represented by the Atom class, offering a set of methods through the atomicData manager. Regarding energy levels, A-values, and ECS, PyNeb has three file types:
  • x x x _ levels.dat: lists NIST energy levels for ionic species x x x ;
  • x x x _ atom_ref.dat: lists A-values for transitions between energy levels of ion x x x contained in the ref dataset;
  • x x x _ coll_ref.dat: lists ECS for transitions between energy levels of ion x x x contained in the ref dataset.
A printout of the NIST energy levels for a specific ion, Fe iii say, may be obtained with the function getLevelsNIST():
# Call thePyNeb Python module
import pyneb as pn
# List NIST energy levels for Fe III
levels = pn.getLevelsNIST(’Fe3’)
print(levels)
and a list of all the Fe iii atom and coll files is currently available with:
pn.atomicData.getAllAvailableFiles(’Fe3’)
  • [’* fe_iii_atom_Q96_J00.dat’,
  •  ’* fe_iii_coll_Z96.dat’,
  •  ’fe_iii_atom_BBQ10.dat’,
  •  ’fe_iii_atom_NP96.dat’,
  •  ’fe_iii_coll_BB14.dat’,
  •  ’fe_iii_coll_BBQ10.dat’]
where the recommended default files are prefixed with an asterisk. Alternative files may be easily accessed by the user to replace the default files with the function
pn.atomicData.setDataFile(’fe_iii_atom_BBQ10.dat’)
pn.atomicData.setDataFile(’fe_iii_coll_BBQ10.dat’)
We can now instantiate the Fe iii ion with the new atomic datasets
pn.atomicData.getDataFile(’Fe3’)
Fe3=pn.Atom(’Fe’,3,NLevels=34)
and apply a series of methods to determine some of its plasma properties:
# Emissivity for transition 2-1 at T = 1.0e4 K and n_e = 1.0e4 cm-3
Fe3.getEmissivity(tem=1.0e4,den=1.0e4,lev_i=2,lev_j=1)
# Emissivity for transition with wavelength 4658.17 A at T = 1.0e4 K
# and n_e = 1.0e4 cm-3
Fe3.getEmissivity(tem=1.0e4,den=1.0e4,wave=4658.17)
# Level critical densities at T = 1.0e4 K
Fe3.getCritDensity(1.0e4)
# Level populations at T = 1.0e4 K and n_e = 1.0e4 cm-3
Fe3.getPopulations(tem=1.0e4, den=1.0e4)
In the ion instantiating command above, the option NLevels specifies the number of levels of the atomic model, which can be the total (default) or a reduced set.
A complete list of methods is included in the PyNeb manual. The Atom class also addresses the ion atomic properties, e.g., the ECS at a selected temperature
# Collision strength for transition 2-1 at T = 1.0e4 K
Fe3.getOmega(1.0e4,2,1)
This function is useful as it returns the interpolated value at the prescribed input temperature allowing comparisons between different ECS datasets with tabulations on diverse temperature ranges. The original temperature and ECS arrays may be listed with the commands
# Temperature and collision strength original arrays
Fe3.getTemArray()
Fe3.getOmegaArray()

3. Fe iii

3.1. Fe iii Atomic Datasets in PyNeb

Fe iii atomic models in PyNeb 1.1.16 comprise 34 fine-structure levels (energy E < 60,000  cm 1 ) from the 3 d 6 and 3 d 5 4 s electron configurations, whose attributes were obtained from the NIST server (April 2017). Levels 3 d 6 b 1 D 2 and b 1 S 0 at E > 77,000  cm 1 are therefore not taken into account, and from the 3 d 5 4 s configuration, only the lowly lying 7 S 3 and 5 S 2 are included. Transitions among the levels of these two even-parity configurations are electric dipole (E1) forbidden and thus occur via electric quadrupole (E2) and magnetic dipole (M1) operators. Shortcomings in the spectral fits resulting from the number of levels in the atomic models, 34 levels in this case, will be examined by comparing with extended models of 144 levels (see Section 3.2).
Two fundamental interactions in the computation of accurate radiative (A-values) and collisional (ECS) rates for atomic ions are electron correlation (configuration interaction, CI) and relativistic coupling. We give below brief descriptions of the numerical methods and approximations that were used to compute the datasets available in PyNeb 1.1.16.
  • atom_NP96 —Contains A-values for E2 and M1 transitions computed in an extensive CI framework with the multi-configuration Breit–Pauli (mcbp) code superstructure [22].
  • coll_Z96 —ECS were calculated in an 83-term non-relativistic R-matrix calculation including the 3 d 6 , 3 d 5 4 s , and  3 d 5 4 p configurations. ECS for 219 fine-structure levels were then obtained through algebraic recoupling [25].
  • atom_Q96_J00 —The Pauli Hartree–Fock hfr code was used to compute A-values for the radiative transitions within the 3 d 6 configuration [44]. Wave functions were generated with CI expansions adjusting the electrostatic and spin–orbit integrals to fit the spectroscopic level energies. A-values for the 3 d 5 4 s 7 S 3 and 5 S 2 levels were obtained independently with hfr using empirically adjusted Slater parameters [52].
  • atom_BBQ10, coll_BBQ10 —A 36-configuration CI expansion including pseudo-orbitals ( 4 f ¯ , 5 s ¯ , 5 p ¯ , and  5 d ¯ ) was used to compute A-values with the mcbpautostructure atomic structure code [27,53]. Term-energy corrections were introduced to fine-tune the wave functions. ECS were computed with the Dirac–Coulomb R-matrix package (darc) based on a 3-configuration target constructed with the grasp92 multi-configuration Dirac–Hartree–Fock structure code [54] under the extended average level approximation.
  • coll_BB14 —ECS were computed with the icft R-matrix method using a 136-term (322 fine-structure levels) target representation [28]. This mcbp 3-configuration target was generated with autostructure using a Thomas–Fermi–Dirac–Amaldi model potential.

3.2. Fe iii Revised and New Atomic Datasets

The present work intends to revise the datasets currently available in PyNeb, incorporate new ones, and implement procedures for data evaluation. For such purposes, we have collected the data from the original sources, verified consistency, and implemented new atom and coll files. To begin with, we have introduced a new level file, fe_iii_levels.dat, with updated spectroscopic measurements (improved accuracy and new levels such as 3 d 6 b 1 D 2 ) from the NIST database (July 2022) [55]. We list below the revised and new atomic datasets with brief descriptions; the suffix “ _ n ” is now appended to the dataset label, where n is the number of levels in the atomic model. For data comparisons, we have normalized all Fe iii datasets to 34 levels, but to estimate model convergence, we have extended some models to 144 levels.
An important issue is how to deal with missing energy levels in both the theoretical and spectroscopic reports since the PyNeb energy-level files are based on the NIST database, which is incomplete at the higher energies. In spectrum modeling, levels with zero decay rates may lead to numerical problems; we have, therefore, limited the number of levels in the atomic models to those that have been both computed and measured. From the user’s point of view, modeling turnaround time is also an important issue; therefore, the number of levels in the atomic models must take this point into consideration.
  • atom_Q96_34 —The main differences with the atom_Q96_J00 dataset described in Section 3.1 are the A-values for transitions involving the 3 d 5 4 s 7 S 3 and 5 S 2 levels. The radiative rates in atom_Q96_J00 for transitions decaying from the 7 S 3 upper level are set to A = 0.0 , while those from 5 S 2 were calculated with hfr with empirically adjusted Slater parameters [52]. In this dataset, radiative rates for transitions involving these levels have been updated with A-values from atom_BB14_34. This choice is to a certain extent arbitrary as inferred from the wide A-value scatter shown for the transitions 5 S 2 5 D J in Table 1. We include in this comparison A-values computed with extensive CI using both hfr and autostructure (atom_FBQh16 and atom_FBQa16) [45]; however, these datasets lack the completeness required to model collisionally excited nebular lines and, consequently, will not be further considered in the present data assessment.
  • coll_Z96_34 —ECS listed in the supplementary data associated with the original publication [25] were downloaded directly from the CDS to build a new file (see Section 3.1 for details of the calculation).
  • coll_Z96_144 —The coll_Z96_34 dataset was extended to include 144 levels from the 3 d 6 , 3 d 5 4 s , and  3 d 5 4 p configurations with E < 121,500  cm 1 . For higher energies, inconsistencies in the reported list of measured and computed levels begin to appear. Although this model extension is not expected to contribute to the collisionally excited nebular lines, it may be relevant if the populations of the 3 d 5 4 p levels generate E1 arrays.
  • atom_DH09_34A-values for the E2 and M1 transitions between levels of the 3 d 6 configuration have been computed with the mcbp code civ3 [19,56] in a CI scheme spanning single and double excitations [38]. To improve accuracy, the diagonal elements of the Hamiltonian matrix were fine-tuned to fit the experimental energies. Similarly to atom_Q96_34, A-values for transitions involving the 3 d 5 4 s 7 S 3 and 5 S 2 levels are taken from atom_BB14_34.
  • atom_BBQ10_34, coll_BBQ10_34A-values and ECS were downloaded from the xstar database [57], and the PyNeb atom and coll files were rebuilt.
  • atom_BB14_34, coll_BB14_34 —Computed A-values and ECS were obtained from the crlike_nrb13#fe2.dat adf04 file downloaded from the open-adas3 database. Details of the ECS computation [28] are described in Section 3.1. No information is given in the adf04 file on the provenance of the radiative data; thus, we assume they are coproducts of the target calculation.
  • atom_BB14_144, coll_BB14_144 —The BB14_34 model has been extended to include 144 levels from the 3 d 6 , 3 d 5 4 s , and  3 d 5 4 p configurations with E < 121,500  cm 1 .

3.3. Fe iii Observational Benchmarks

3.3.1. A-Value Ratio Benchmark

It is well known that the intensity ratio of two spectral lines arising from a common upper level can be expressed as
I ( j , i ) I ( j , i ) = A ( j , i ) A ( j , i ) × λ ( j , i ) λ ( j , i )
enabling an observational benchmark of the theoretical A-value ratio A ( j , i ) / A ( j , i ) . This measure provides indications of the accuracy of both the observed intensities and computed radiative rates, signposting unreliable lines to be excluded from the spectrum fits.
In Table 2 and Table 3, we compare the observed line-intensity ratios from HH 202S and HH 204  [9,11] with the A-value ratios from the BBQ10_34, BB14_34, DH09_34, and Q96_34 atom files. The relative uncertainties in the observed intensity ratios from HH 202S are, in general, larger than those from HH 204. Although both spectra are of comparable quality, HH 204 was observed under photometric conditions whereby the absolute calibration error is slightly smaller; moreover, the HH 204 integrated emission comprised twice the area of HH 202S with longer exposure times.
For both sources, the observed intensity ratios for transitions arising from levels 3 d 6 a 3 F J or lower and ending up at a level within the ground configuration 3 d 6 a 5 D with wavelengths λ > 4600  Å are generally more accurate (better than 15%). For transitions arising from higher levels, uncertainties as large as 70% have been listed. On the other hand, most theoretical ratios agree to within 10% except for significantly discrepant cases in atom_BBQ10_34: I ( λ 3286 ) / I ( λ 8729 ) ; I ( λ 3366 ) / I ( λ 9960 ) ; and I ( λ 3357 ) / I ( λ 8838 ) . As a whole, the theoretical ratios lie within the observed ratio error bars except for some questionable lines in HH 202S ( λ 8729 , affected by telluric absorptions) and in HH 204 ( λ λ 4080 , 4097 blended with O ii lines and λ λ 7078 , 7088 , 9204 being misidentifications).

3.3.2. [Fe iii] Spectrum Fits

Following [4], we take advantage of PyNeb’s methods for emission-line modeling to implement a further observational benchmark to assess the accuracy of the atomic datasets described in Section 3.1 and Section 3.2. For this purpose, we again use the high-resolution spectra of the bright HH 202S and HH 204 sources.
The electron temperature and density of the HH object are obtained by a least-squares fit of all the reliable spectral lines as proposed by [58]. The theoretical emissivity and observed intensity of each line are respectively normalized with respect to the emissivity and intensity sums of the transition array. The fit measure is given in terms of
χ 2 = i = 1 N ( I i ϵ i ) 2 N σ I i 2 ,
where N is the number of spectral lines, I i and ϵ i are the normalized intensity and emissivity, respectively, and  σ I i 2 is the variance of the normalized intensity [59]. No explicit error is considered in the theoretical emissivities. The least-squares minimization is carried out with the Python function scipy.optimize.minimize finding no degeneracies in the χ 2 convergence.
PyNeb provides two methods to compute the emissivity of a spectral line (e.g., λ 5270.57 between levels i = 6 and j = 2 of Fe iii) at ( T e , n e ) by addressing, say, the datasets atom_Q96_J00 and coll_Z96. The Python commands to run these two methods (e1 and e2) are as follows:
# Call the PyNeb Python module
import pyneb as pn
# Select the radiative and collisional datasets
pn.atomicData.setDataFile(’fe_iii_atom_Q96_J00.dat’)
pn.atomicData.setDataFile(’fe_iii_coll_Z96.dat’)
# Instantiate Fe III
Fe3=pn.Atom(’Fe’,3)
# Check the datafiles
print(Fe3)
# Determine and print the two emissivities
e1=Fe3.getEmissivity(tem=1.e4,den=1.0e4,wave=5270.57)
e2=Fe3.getEmissivity(tem=1.e4,den=1.0e4,lev_i=6,lev_j=2)
print(e1,e2)
In the present work, consistency between e1 and e2 for the whole spectrum was tested to avoid incorrect line identifications, but we found that e2 avoided the problem of transitions with the same wavelength in atomic models with a large number of levels.
In the spectrum fits of HH 202S and HH 204, a handful of lines are excluded due to the following problems:
  • Misidentifications: λ λ 7078 , 7088 , 9204
  • Line blending or telluric contamination [9]: λ λ 4080 , 4097 , 4925 , 4931, 4987, 8729, 8838
  • Small A-values ( 10 3  s 1 ): λ λ 4008 , 4047.
We first examine in Table 4 the spectral fits with the atomic datasets available in PyNeb 1.1.16 (see Section 3.1). For each coll dataset BBQ10, BB14, and Z96, we run the fits in turn with three atom datasets: BBQ10, NP96, and Q96_J00. For HH 202S, we obtain an average temperature and density of T e = ( 7 ± 1 ) × 10 3  K and n e = ( 6 ± 5 ) × 10 4  cm 3 , which indicate large uncertainties due to the atomic data. In particular, coll_BB14 gives, on average, a questionable low temperature ( 5.4 × 10 3  K) and high density ( 1.3 × 10 5  cm 3 ). The quality of the fits is generally poor, the best being atom_BBQ10–coll_Z96 with χ 2 = 1.74 . The fits for HH 204 are as unimpressive: T e = ( 8 ± 2 ) × 10 3  K) and n e = ( 5 ± 4 ) × 10 4  cm 3 , the best ( χ 2 = 16.9 ) being with the file pair atom_BBQ10–coll_Z96.
In a similar fashion, we proceed to benchmark the revised and new datasets described in Section 3.2 by running a grid of fits with different atom–coll file pairs: four radiative and three collisional datasets (see Table 5). The improvement is considerable; the average temperature and density for HH 202S is now T e = ( 7.9 ± 0.3 ) × 10 3  K and n e = ( 2.1 ± 0.3 ) × 10 4  cm 3 and for HH 204 T e = ( 8.9 ± 0.6 ) × 10 3  K and n e = ( 1.5 ± 0.3 ) × 10 4  cm 3 . Outstanding fits for HH 202S are obtained with the atom_DH09_34–coll_Z96_34 and atom_Q96_34–coll_Z96_34 dataset pairs, and the revised coll_BB14_34 ECS dataset now gives more reasonable temperatures ( 8 × 10 3  K) and densities ( 2 × 10 4  cm 3 ). We find questionable ECS differences between coll_BB14 and coll_BB14_34 for transitions involving levels a 3 P J , a 3 F J , and  a 1 G 4 ; furthermore, for levels with E > 43,000  cm 1 , the ECS in coll_BB14 are Υ ( j , i ) = 0.0 . We also found ECS differences for transitions to level 3 d 5 4 s 7 S 3 in coll_BBQ10 and coll_BBQ10_34 due to improvements performed by [27] after publication. On the other hand, the extended datasets atom_BB14_144, coll_BB14_144, and coll_Z96_144, as expected, make no difference in these plasma conditions.
The temperature and density scatters are now within 6% and 20%, respectively, and our average values compare satisfactorily with those adopted in the observational papers [9,11]: T e = ( 8.760 ± 0.180 ) × 10 3  K and n e = ( 1.354 ± 0.121 ) × 10 4  cm 3 for HH 204 and n e = ( 1.743 ± 0.236 ) × 10 4  cm 3 for HH 202S. No adopted temperature is quoted for the latter, but our average value compares well with their T e ( [ S II ] ) = ( 8.25 ± 0.54 ) × 10 3  K, T e ( [ Ar III ] ) = ( 8.26 ± 0.41 ) × 10 3  K, and  T e ( He I ) = ( 7.95 ± 0.20 ) × 10 3  K.

3.3.3. Line-Ratio Diagnostics

Line-ratio diagnostics are widely used in nebular astrophysics to determine the plasma thermodynamic properties (temperature, density, and ionic abundances). They are based on the assumption that the observed line-intensity ratio is equivalent to the theoretical emissivity ratio
I ( j , i ) I ( j , i ) = ϵ ( j , i ) ϵ ( j , i ) .
Using our benchmarked atomic datasets for [Fe iii], we investigate the impact of the atomic data scatter on the line-ratio temperature and density behaviors.
A common [Fe iii] emissivity ratio to determine the plasma density in the range 3 log n e 6 is I ( λ 4658 ) / I ( λ 4702 ) (see, for instance, [10,11]). In Figure 1, we show its dependency on the atomic data. In the left panel, we select dataset atom_DH09_34 and plot I ( λ 4658 ) / I ( λ 4702 ) as a function of density for three coll datasets: BBQ10_34, BB14_34, and Z96_34. Large discrepancies are found in the low-density regime ( log n e < 4 ), particularly in coll_BBQ10_34 (40% lower than coll_Z96_34 at log n e < 3 ). In the right panel, we select dataset coll_Z96_34 and plot the emissivity ratio for four atom datasets: DH09_34, BBQ10_34, BB14_34, and Q96_34. Differences are not as large, atom_BBQ10_34 being the more discordant ( 15 % lower at log n e < 3 ).
For log n e 3 , the emissivity ratio I ( λ 4658 ) / I ( λ 4702 ) is practically density insensitive, being mostly dependent on the ECS of the involved atomic levels. To discern a reliable magnitude in this regime in the light of the discrepancies due to the collisional datasets, we select from the DESIRED database (Méndez-Delgado et al., unpublished) the following low-density sources: 30 Doradus [60]; Sh 2-311 [61]; M20 [62]; M17 [63]; HII-1, HII-2, UV-1 [64]; K932, NGC 5461, NGC 604, NGC 2363, NGC 4861, NGC 1741-C, NGC 5447, VS44, H1013 [65]; NGC 3125, POX 4, TOL 1924-416, TOL 1457-262, NGC 6822 (HV), NGC 5408 [66]; Sh 2-100, Sh 2-288 [67]; NGC 5471, NGC 5455  [68], Sh 2-152 [69]; and N44C, N11B, N66A, NGC 1714, IC 2111, N81 [70]. We then plot the observed line-intensity ratios in Figure 2 as a function of the [S ii] density. In spite of the scatter, there is a definite trend for I ( λ 4658 ) / I ( λ 4702 ) > 3.5 , which would reinforce the validity of the coll_BB14_34 and coll_Z96_34 datasets.
Regarding a useful [Fe iii] temperature diagnostic, we choose I ( λ 4658 ) / I ( λ 3240 ) due to its small variations in density and atomic data. In Figure 3, we plot this emissivity ratio as a function of temperature for the same atomcoll grid used in the density diagnostic. The curves show only a weak dependency on the choice of the collisional and radiative datasets.
A relevant question is whether the density and temperature estimates obtained from the I ( λ 4658 ) / I ( λ 4702 ) and I ( λ 4658 ) / I ( λ 3240 ) line ratios, respectively, compare with those obtained from the spectral fits described in Section 3.3.2. For HH 202S, as an example, the observed I ( λ 4658 ) / I ( λ 4702 ) = 2.81 and I ( λ 4658 ) / I ( λ 3240 ) = 13.6 . In the density diagnostic, we set the temperature at T e = 7.9 × 10 3  K and obtain n e = 2.34 × 10 4  cm 3 for the pair atom_DH09_34–coll_Z96_34 and n e = 1.86 × 10 4  cm 3 for atom_DH09_34–coll_BB14_34, values that are fairly close to those obtained from the fits, respectively n e = 1.97 × 10 4  cm 3 and n e = 1.78 × 10 4  cm 3 . However, if we designate the pair atom_DH09_34–coll_BBQ10_34, the theoretical I ( λ 4658 ) / I ( λ 4702 ) < 2.78 impairs a density reading for HH 202S. Furthermore, if we set the density at n e = 2.1 × 10 4  cm 3 , we obtain from the temperature diagnostic 7.89 × 10 3  K, 7.95 × 10 3  K, and  8.33 × 10 3  K for the atomcoll pairs DH09_34–Z96_34, DH09_34–BBQ10_34, and DH09_34–BB14_34, respectively, which compare well with 7.71 × 10 3  K, 7.94 × 10 3  K, and  8.18 × 10 3  K from the spectral fits.

3.4. Radiative Lifetimes of 3 d 6 and 3 d 5 4 s Levels

As previously mentioned, the [Fe iii] nebular spectrum displays transitions between levels of the 3 d 6 and 3 d 5 4 s configurations. The PyNeb Fe iii atomic models are based on the lower 34 levels with E < 60,000  cm 1 ; that is, the cut is imposed just before the onset of the 3 d 5 4 s level array. As a result, in the present lifetime analysis, we obviate the 3 d 6 b 1 D 2 level at 77,044.67 cm 1 and the unobserved 3 d 6 1 S 0 estimated theoretically at E > 99,000  cm 1 . From  3 d 5 4 s , only the lowly lying 7 S 3 and 5 S 2 levels are taken into account.
Radiative lifetimes derived from the A-values in the atom datasets BB14_34, BBQ10_34, DH09_34, and Q96_34 are compared in Figure 4. For each level, we determine the lifetime differences of these four datasets with respect to the average value. Such differences are, in general, within 0.1 dex except for five levels in atom_BB14_34, nine in atom_BBQ10_34, and two in each atom_DH09_34 and atom_Q96_34. Notably large discrepancies (as large as 0.8  dex) are found for the 3 d 6 3 H J levels, which would affect the physical conditions and Fe 2 + abundance derived from the λ λ 4881, 4987, 4924, 5032, 4986 line intensities. Furthermore, the  3 d 6 b 3 P J levels in atom_BB14_34 and 3 d 6 a 3 P 2 and a 1 G 4 in atom_BBQ10_34 have lifetimes almost 0.2 dex lower. Significant differences also stand out for the 3 d 5 4 s 7 S 3 and 5 S 2 levels between these two datasets. A reassuring outcome is the agreement between the atom_DH09_34 and atom_Q96_34 lifetimes: within 20% except for the problematic 3 H J levels. Most of the conclusions from this comparison are in line with [27].

3.5. Radiative Lifetimes of 3 d 5 4 p Levels

Apart from even-parity levels, dataset atom_BB14_144 contains 44 odd-parity levels belonging to the 3 d 5 4 p configuration. Their radiative lifetimes can be compared with those derived from the A-values computed with extensive CI (single and double excitations up to n 5 , 3 , and 6p as well as selective promotions from the 3s and 3p subshells) with the mcbp code civ3  [39]. As shown in Figure 5, the atom_BB14_144 lifetimes are on average 0.15 dex shorter. The larger differences are found for levels belonging to the 3 d 5 ( 4 P ) 4 p 3 P o , 3 d 5 ( 4 P ) 4 p 5 P o , and  3 d 5 ( 4 D ) 4 p 5 F o spectroscopic terms.

3.6. Effective Collision Strengths

As shown in Table 5, we have three revised coll datasets containing ECS for the 34-level atomic model of [Fe iii]: Z96_34, BBQ10_34, and BB14_34. We compare these data at T = 10 4  K in Figure 6 using the latter as the reference. In general, there is a reasonable agreement for most transitions, but for several others, huge differences are encountered. Transitions involving the 3 d 5 ( 6 S ) 4 s 7 S 3 level are problematic. For instance, in Figure 6 (right panel), we have excluded from coll_Z96_34 transitions of the type 3 d 5 ( 6 S ) 4 s 7 S 3 3 d 6 1 L J as they appear with Υ = 0.0 ; furthermore, coll_Z96_34 transitions with Δ S = 3 / 2 have abnormally small ECS ( log Υ < 4 ). These shortcomings may be due to the non-relativistic method used in [25]. On the other hand, as discussed in [28], the main differences between the ECS in coll_BBQ10_34 and coll_BB14_34 are due to a mixup by the former in the assignments of the 3 d 5 ( 6 S ) 4 s 7 S 3 and 3 d 6 3 G 3 adjacent levels. Large discrepancies are also found in transitions of the type 3 d 6 a 3 P J 3 d 6 3 H J and 3 d 6 b 1 G 4 3 d 6 5 D J .

3.7. Fe iii Discussion

Our initial step in the present atomic data assessment was the comparison of observed line-intensity ratios with the corresponding theoretical A-value ratios for transitions arising from a common upper level (see Section 3.3.1). This procedure readily yields a first exposure of observational line misidentifications, blending, and contamination (e.g., telluric or instrumental) and of inherent computational difficulties such as strong valence–core correlation, double-excitation couplings, and strong CI mixing. Therefore, the agreement between theoretical ratios, on the one hand, and between observed and theoretical ratios, on the other hand, are both revealing.
Although the spectra of HH 202S and HH 204 are of comparable quality, the reported observational errors of the former are on average larger (36%) than the latter (12%). Comparisons among theoretical A-value ratios and with the observed counterparts lead to an overall accuracy level for the theoretical ratios of around 20–30%, except for odd theoretical outliers pinpointed in Section 3.3.1 believed to be caused by the aforementioned computational difficulties.
The spectrum fits performed in Section 3.3.2 bring to the fore a frequent drawback in modeling codes: the error-prone mapping of radiative or collisional datasets onto an atomic model with different level numbering as illustrated with dataset coll_BB14 that led to file deprecation. The four revised atom datasets (BBQ10_34, BB14_34, DH09_34, Q96_34) and three coll datasets (BBQ10_34, BB14_34, Z96_34) appear to give comparable temperatures and densities within error bands of 6% and 20%, respectively, which also match satisfactorily estimates obtained from [Fe iii] line-ratio diagnostics such as I ( λ 4658 ) / I ( λ 4702 ) and I ( λ 4658 ) / I ( λ 3240 ) (see Section 3.3.3) and from alternative ionic diagnostics reported in [9,11]. However, the density-sensitive ratio I ( λ 4658 ) / I ( λ 4702 ) computed with coll_BBQ10_34 does not match the ratios observed in low-density ( log n e 3 ) nebular sources (see Figure 2).
As shown in Section 3.4, radiative lifetimes for levels of the 3 d 6 and 3 d 5 4 s configurations computed with A-values from the BBQ10_34, BB14_34, DH09_34, and Q96_34 atom datasets agree to within 0.1  dex except for 3 d 6 3 H J and 3 d 5 4 s 7 S J , 5 S 2 . However, the best agreement ( 20 % ) is found between A-values from structure calculations, namely atom_DH09_34 and atom_Q96_34, rather than those associated with scattering targets. This may be due to the convergence and accuracy of the CI expansions implemented in each calculation type. In structure calculations of radiative lifetimes for transitions involving levels of the 3 d 6 and 3 d 5 4 s configurations in Fe iii, for instance, there is no need to include odd-parity configurations, while in a Fe iii scattering target, a reduced set of configurations of both parity types must be considered. A similar argument may be applied to the lifetimes of 3 d 5 4 p levels (see Section 3.5), where the atom_BB14_144 lifetimes are on average shorter than those from a lengthy mcbp structure calculation [39].
The main finding of the ECS comparison in Section 3.6 is the large discrepancies displayed by coll_Z96_34 ECS with respect to coll_BB14_34 for some transitions with log Υ < 1 (see Figure 6), which may be due to the non-relativistic R-matrix method used by [25] in the computation of the former dataset. The larger ECS differences between coll_BBQ10_34 and coll_BB14_34 are mainly due to level misassignments.
Analyzing results from these data comparisons and spectral fits, we select as the [Fe iii] default datasets in PyNeb 1.1.17 the atom_BB14_144 and coll_BB14_144 files, while files coll_BBQ10 and coll_BB14 are deprecated. Deprecated files are still available with the command pn.atomicData.includeDeprecatedPath(). Files atom_BBQ10 and coll_Z96 remain active, and files coll_BBQ10_34, atom_BB14_144, atom_DH09_34, and atom_Q96_34 are incorporated.

4. Fe ii

Similarly to the Fe iii atomic data assessment, we implement observational benchmarks to evaluate the atomic datasets for Fe ii available in PyNeb 1.1.16, transcribing them directly from the source and implementing new models. NIST energy levels in this version of PyNeb are the same as those listed in the database of July 2022 [55].

4.1. Fe ii Atomic Datasets in PyNeb

The following datasets are available for Fe ii in PyNeb 1.1.16:
  • atom_VVKFHF99, coll_VVKFHF99 —These datasets list radiative and collisional rates for transitions among 80 levels with energies E < 44,760  cm 1 compiled from different sources [71]. A-values for the dipole allowed and forbidden transitions are from [21] and [23], respectively. ECS are mainly from [24]; however, the source data are not available and A-values and ECS for transitions with spin change Δ S = 2 have been obviated.
  • atom_B15, coll_B15A-values and ECS for transitions among 52 levels belonging to the 3 d 7 , 3 d 6 4 s , and  3 d 5 4 s 2 configurations ( E < 31,500  cm 1 ) [29]. A-values were computed with the mcbp atomic structure code autostructure using different expansions and orbital optimization strategies and a dipole correction to the Thomas–Fermi–Dirac–Amaldi potential. Fine-tuning was carried out by means of term-energy corrections adjusted by fitting to the experimental term energies. ECS were computed with the Dirac–Coulomb R-matrix package darc [34] using a target generated with the fully relativistic atomic structure code grasp0 [54] using a six-configuration expansion comprising 329 levels.
  • atom_SRKFB19, coll_SRKFB19A-values and ECS for 250 levels from the 3 d 7 , 3 d 6 4 s , 3 d 5 4 s 2 , 3 d 6 4 p , and  3 d 5 4 s 4 p configurations [31]. A-values were computed with grasp0 in a 20-configuration expansion and ECS with darc using a 716-level target representation.

4.2. Fe ii Revised and New Atomic Datasets

For comparison purposes, models were normalized to 52 levels ( E < 31,500  cm 1 ), i.e., the number of levels considered in atom_B15 that includes the last even quartet term before the appearance of the 3 d 6 4 p level array. Two datasets have been extended to comprise 173 levels ( E < 66,100  cm 1 ), a cut before the onset of unmeasured levels in the NIST level list, and 225 levels ( E < 73,800  cm 1 ) to discuss some assignment inconsistencies. Although ECS from [24] are included in VVKFHF99, we have not considered their dataset individually as it excludes the doublet terms in their atomic model.
  • atom_QDZa96, atom_QDZh96 —E2 and M1 A-values calculated for transitions among the 63 levels of the configurations 3 d 6 4 s , 3 d 7 , and  3 d 5 4 s 2  [23]. Extensive CI expansion was used with the Breit–Pauli superstructure [18] (QDZa96) and Pauli hfr [20] (QDZh96) atomic structure codes. These datasets are not complete enough to be used in PyNeb spectral modeling as they obviate transitions such as a 4 F a 6 D and a 2 P a 6 D ; therefore, they will only be used in radiative data comparisons.
  • atom_VVKFHF99_52, coll_VVKFHF99_52 —The VVKFHF99 datasets containing A-values and ECS described in Section 4.1 are reduced to 52 levels from even-parity configurations.
  • atom_DH11_52 —E2 and M1 A-values for transitions between levels of the 3 d 6 4 s , 3 d 7 , and  3 d 5 4 s 2 configurations were computed in a large-scale CI framework with the mcbp structure code civ3 [42]. Wave functions were fine-tuned to fit the experimental energies.
  • atom_B15_52, coll_B15_52 —Data (A-values and ECS) from the B15 calculation described in Section 4.1 were downloaded directly from the CDS, and the datasets were reconstructed.
  • atom_TZ18_52, coll_TZ18_52A-values and ECS were computed with the close-coupling B-spline Breit–Pauli R-matrix method for 340 levels of the 3 d 6 4 s , 3 d 7 , 3 d 5 4 s 2 , 3 d 6 4 p , and  3 d 5 4 s 4 p configurations [30]. The CI target representation was constructed with a Hartree–Fock method. Semi-empirical fine-tuning in the structure and scattering calculations was introduced by fitting tothe experimental energies. The final Hamiltonian of the scattering calculation contains the spin–orbit term.
  • atom_SRKFB19_52, coll_SRKFB19_52 —Radiative and collisional data from the SRKFB19 calculation were described in Section 4.1. New files were constructed with data extracted from the adf04_DARC716_2.10.18 file (Cathy Ramsbottom, private communication) and reduced to 52 levels.
  • atom_SRKFB19_173, coll_SRKFB19_173, atom_TZ18_173, coll_TZ18_173 —The SRKFB19 and TZ18 datasets have been extended to 173 levels ( E < 66,100  cm 1 ) to test model convergence.
  • atom_SRKFB19_225, atom_TZ18_225 —The theoretical SRKFB19 and TZ18 datasets have been extended to 225 levels ( E < 73,800  cm 1 ) to bring out term assignment inconsistencies involving levels with total orbital angular momentum quantum number L = 5 (see Table 6). With respect to NIST and other theoretical datasets, the term assignments of levels at 21,430.36 and 21,581.62  cm 1 ( a 4 H ) and 26,170.18 and 26,352.77  cm 1 ( b 2 H ) have been interchanged in atom_SRKFB19_225. Since the total angular momentum quantum number J of each level coincides in the different datasets, the incongruent term labeling would be inconsequential if the A-values for the transitions involving these levels were comparable. As shown in Table 7, this is not the case, as discrepancies as large as an order of magnitude were encountered. A similar mixup was found with the odd-parity z 4 H o and y 2 H o levels, which in this case are even energetically misplaced (see Table 6). Assignment discrepancies of this sort could be due to strong CI admixture in terms with high orbital angular momentum ( L > 4 say) in lowly ionized systems that the grasp0 multi-configuration Dirac–Hartree–Fock structure code has problems resolving.

4.3. Observational and Spectroscopic Benchmarks

Due to the wide astrophysical interest in Fe ii, this ion presents interesting possibilities for implementing observational and spectroscopic benchmarks. There is interest not only in the forbidden-line spectrum but also in the allowed transition arrays involving levels from the odd-parity configurations 3 d 6 4 p and 3 d 5 4 s 4 p , which may be populated by fluorescence pumping from the stellar continuum [72,73,74]. Fe ii ultraviolet emission in active galactic nuclei is not fully understood, impeding quasar classification schemes and Fe abundance estimates [75].

4.3.1. A-Value Ratio Benchmark

In Table 8 and Table 9, we compare the observed line-intensity ratios in the spectra of HH 202S and HH 204 involving transitions from a common upper level with the corresponding A-value ratios (see Equation (1)) from atom datasets B15_52, DH11_52, TZ18_52, SRKFB19_52, QDZa96, and QDZh96. In this comparison, we must point out the worrisome discrepancies displayed by some A-value ratios from atom_SRKFB19_52 with both the observed and other theoretical values; thus, the integrity of this radiative dataset is questionable.
For two cases, the observed line-intensity ratios are very different from the corresponding theoretical values: I ( λ 3968.27 ) / I ( λ 4305.90 ) and I ( λ 4889.71 ) / I ( λ 5412.68 ) , which are caused by the misidentification of the stronger [Fe ii] λ λ 3968.66 , 4889.62 lines with the weaker [Fe ii] λ λ 3968.27 , 4889.71 , which are emitted from different atomic transitions.
Excluding atom_SRKFB19_52, the A-value ratios in Table 8 and Table 9 agree to 20 25 % except for two cases in HH 202S, I ( λ 7637.52 ) / I ( λ 9051.95 ) and I ( λ 7686.93 ) / I ( λ 8891.93 ) , and one in HH 204, I ( λ 7733.13 ) / I ( λ 9033.49 ) , which are caused by lines with A-values of order 10 3 subject to large numerical uncertainties. Most theoretical ratios are within the estimated observational error bars except for those involving A-values strongly affected by level mixing effects that are frequent in Fe ii [40].

4.3.2. [Fe ii] Spectrum Fits

In a similar manner to [Fe iii] (see Section 3.3.2), we fit the reliable [Fe ii] lines of the high-resolution spectra of the bright objects HH 202S and HH 204 with PyNeb theoretical emissivities optimized through a least-square procedure in terms of the electron temperature and density. The theoretical emissivity and observed intensity of each line are, respectively, normalized to the sum of the emissivities and the sum of intensities.
In the spectrum fit of HH 202S and HH 204, a selection of lines were excluded mainly due to observational problems [9,11]: λ λ 4413.78, 4416.27, 4452.10, 4509.61, 4514.90, 4528.38, 4889.71, and 4905.35 in addition to λ λ 4276.84, 4319.62, 4359.33, 4452.10, 7665.28, 9399.04, and 9997.49 in the latter source. The line λ 4382.74 is also questionable due to an incorrect identification.
We show in Table 10 the results of the spectral fits with the atomic datasets available in PyNeb 1.1.16. In contrast to [Fe iii] (see Section 3.3.2), the optimization procedure does not always converge to the desired accuracy (e.g., SRKFB19 for HH 204), and the density variation may not show a minimum or may depend on the initial input value. The temperature scatter is worrisome, and the HH 204 fits, in particular, show large values of χ 2 .
The revised and new Fe ii atomic models comprise four radiative (atom) and four collisional (coll) files normalized to 52 levels. The spectral fits for HH 202S and HH 204 are shown in Table 11, where each collisional dataset appears to lead to a diverse temperature and density pair ( T e , n e ) . For HH 202S in units of 10 3  K and 10 4  cm 3 , respectively,
  • coll_B15_52: ( 12.0 ± 0.1 , 6.3 ± 0.6 )
  • coll_TZ18_52: ( 8.3 ± 0.4 , 3.9 ± 0.5 )
  • coll_SRKFB19_52: ( 10.3 ± 0.6 , 2.4 ± 0.3 )
  • coll_VVKFHF99_52: ( 10.2 ± 0.7 , 1.6 ± 0.2 )
where the scatter in both temperature and density due to the atomic data is surprisingly broad, and only the temperature from coll_TZ18_52, T e = ( 8.3 ± 0.4 ) × 10 3  K, agrees with that determined from the [Fe iii] fit, T e = ( 7.9 ± 0.3 ) × 10 3  K (see Section 3.3.2). On the other hand, the density n e = ( 2.1 ± 0.3 ) × 10 4  cm 3 obtained in the [Fe iii] fit matches the values by coll_SRKFB19_52, ( 2.4 ± 0.3 ) × 10 4  cm 3 , and coll_VVKFHF99_52, ( 1.6 ± 0.2 ) × 10 4  cm 3 .
For HH 204,
  • coll_B15_52: ( 18.4 ± 0.3 , 3.8 ± 0.2 )
  • coll_TZ18_52: ( 11.9 ± 0.1 , 3.3 ± 0.6 )
  • coll_SRKFB19_52: ( 15.7 ± 0.1 , 1.0 ± 0.2 )
  • coll_VVKFHF99_52: ( 13.6 ± 0.2 , 0.85 ± 0.07 )
whereby the temperatures are significantly higher than the [Fe iii] value T e = ( 8.9 ± 0.6 ) × 10 3  K, while the [Fe iii] density, n e = ( 1.5 ± 0.3 ) × 10 4  cm 3 , favors coll_SRKFB19_52.
This confusing situation may lead to different analyses. For instance, in Figure 7, we plot for lines in HH 202S the ratio of the normalized theoretical emissivity to the normalized observed intensity as a function of the upper-level energy of the transition. In the left panel, the spectrum was fitted with datasets on atom_TZ18_52 and coll_TZ18_52, resulting in the low-temperature T e = 7.94 × 10 3  K and density n e = 3.29 × 10 4  cm 3 . For transitions with upper-level energies E < 2 × 10 4  cm 1 , the number of ratios above and below 1.0 is roughly equal, while at the higher energies, the number of ratios above 1.0 progressively diminishes. This behavior is similar to that depicted in Figure 4 of [73], interpreted as a signature of missing fluorescence pumping by the stellar continuum in the plasma model.
In the right panel, the spectral fit was performed with datasets atom_DH11_52 and coll_SRKFB19_52, resulting in a higher temperature T e = 11.1 × 10 3  K and comparable density n e = 2.58 × 10 4  cm 3 ; however, the missing continuum-pumping signature is not present. The diverse, high temperatures of Table 11 might thus be interpreted as attempts to fit the spectra without accounting for all the level-populating mechanisms.
A second reading is to look at the critical densities of the levels. In Figure 8, we plot the critical densities at T e = 10 4  K for the lower 52 levels of [Fe ii] computed with dataset atom_TZ18_52 and three coll files: B15_52, SRKFB19_52, and TZ18_52. For levels with index i > 10 , critical densities by coll_B15_52 are a factor of two higher, and for i > 35 , they are all significantly different from each other. A key feature is the low critical density ( log ρ crit < 2.6 ) of the 3 d 7 a 4 F 9 / 2 metastable level (level i = 6 ). For an electron density n e = 2 × 10 4  cm 3 , the population of this level may be larger than that of the ground level for temperatures as low as T e 7 × 10 3  K; therefore, this level plays a dominant role in the spectral formation of [Fe ii].
Moreover, in Figure 9, we plot the critical density of level 3 d 7 a 4 F 9 / 2 as a function of temperature using different atomic datasets. In the right panel, we show results computed with the dataset coll_TZ18_52 and three atom files (TZ18_52, B15_52, and DH11_52). Although the sensitivity to the radiative data is extreme, for  T e > 5 × 10 3  K, it is approximately temperature independent. In the left panel, the critical density is computed with the atom_TZ18_52 and three coll datasets (B15_52, SRKFB19_52, and TZ18_52). The dependency on the collisional data is not as marked as on the radiative data, but it is temperature dependent.
In Table 11, we also list the fit results using the extended dataset pairs atom_TZ18_173–coll_TZ18_173 and atom_TZ18_173–coll_SRKFB19_173. For the HH 202S spectrum, the resulting temperature (in 10 3  K units), density (in 10 4  cm 3 units), and  χ 2 values ( T e , n e , χ 2 ) = ( 7.81 , 3.29 , 7.91 ) and ( T e , n e , χ 2 ) = ( 9.75 , 2.32 , 2.84 ) , respectively, are, as expected, close to those from the atom_TZ18_52–coll_TZ18_52 and atom_TZ18_52–coll_SRKFB19_52 pairs: ( T e , n e , χ 2 ) = ( 7.94 , 3.29 , 7.90 ) and ( T e , n e , χ 2 ) = ( 9.91 , 2.21 , 2.75 ) . However, for HH 204, the smaller χ 2 values in atom_TZ18_173–coll_TZ18_173 and atom_TZ18_173–SRKFB19_173 indicate improved fits: χ 2 = 40.2 . and 64.2 as compared to χ 2 = 154 . and 84.7 in the pairs atom_TZ18_52–coll_TZ18_52 and atom_TZ18_52–coll_SRKFB19_52, respectively. Moreover, although the temperatures and densities with coll_SRKFB19_52 and coll_SRKFB19_173 are close, the densities with coll_TZ18_52 and coll_TZ18_173 differ substantially: 2.63 and 1.51, respectively.

4.3.3. Line-Ratio Diagnostics

An extensive discussion of [Fe ii] line-ratio diagnostics in the infrared, near-infrared, and optical is given in [29]. Using lines from the HH 202S and HH 204 spectra, we have examined several line ratios to be used as density diagnostics involving “safe” low-energy levels to avoid predicted fluorescence effects [73]. Since temperature diagnostics inherently rely on levels at higher energies, we have refrained from their treatment in the present atomic data comparisons.
We have selected as a showcase in Figure 10 the emissivity ratio I ( λ 8892 ) / I ( λ 9267 ) that is density sensitive in the range 3 log n e 5  cm 3 . In the left panel, we plot the emissivity ratio computed with the atom_TZ18_52 radiative dataset and several coll ECS files highlighting its dependency on the collisional datasets. We do not show its behavior below log n e = 3 as it becomes more complicated, manifesting divergences also due to the radiative data. We also include in Figure 10 the observed intensity ratio in HH 202S, I ( λ 8892 ) / I ( λ 9267 ) = 2.36 ± 0.42 , which, except for the curve computed with coll_TZ18_52 ECS, seems to degrade the diagnostic potential of this ratio for log n e 4.5 since the density variations lie within its error band. On the other hand, the temperature dependency of the ratio above T e 10 4  K is generally weak, as shown in the right panel.
If this observed line-intensity ratio in HH 202S is used to determine the density, we obtain 7.94 × 10 4 , 1.32 × 10 4 , 1.41 × 10 4 , and  3.55 × 10 4  cm 3 for the coll datasets B15_52, TZ18_52, SRKFB19_52, and VVKFHF99_52, respectively, which compare poorly with estimates from the spectrum fits of Section 4.3.2: respectively, 5.98 × 10 4 , 3.29 × 10 4 , 2.21 × 10 4 , and  1.47 × 10 4  cm 3 . Nonetheless, if we take into consideration the gross uncertainties brought about by the collisional datasets, on average, both methods coincide on a similar poor density diagnostic: ( 4 ± 3 ) × 10 4  cm 3 .

4.4. Radiative Lifetimes of the 3 d 6 4 s , 3 d 7 , and  3 d 5 4 s 2 Levels

To compare the lifetimes of levels belonging to the 3 d 6 4 s , 3 d 7 , and  3 d 5 4 s 2 configurations predicted by the atom datasets B15_52, DH11_52, TZ18_52, SRKFB19_52, QDZa96, and QDZh96, we again determine an average lifetime for each level by comparing the respective differences for each dataset (see Figure 11). A problematic and identifiable case in this plot is the lowest level of the first excited configuration, 3 d 7 a 4 F 9 / 2 , with the longest lifetime ( log τ ¯ = 4.22  s) and large discrepant values ( 0.4  dex). This is an important metastable level inasmuch as being dominant in the plasma radiative and collisional equilibrium and in opening the routes to continuum pumping. In general, the scatter is 0.1  dex, except for lifetimes derived from the A-values of the questionable atom_SRKB19_52 dataset.
Lifetimes for a handful of metastable levels of [Fe ii] have been measured with a laser probing technique, and A-values were derived using branching fractions from the spectra of the Eta Carinae ejecta taken with the Hubble Space Telescope [49,50,51]. A comparison between the experimental lifetimes and those derived from the PyNebatom datasets is given in Table 12. Lifetimes from atom_SRKFB19_52 are markedly discrepant from both the measured and other theoretical values and are thus excluded from this table. The theoretical lifetime of the 3 d 6 4 s b 2 H 11 / 2 level shows a large scatter. Although there is good agreement between theory and experiment for the 3 d 5 4 s 2 a 6 S 5 / 2 level, the rest of the lifetimes from B15_52, DH11_52, and TZ18_52 are generally longer (as long as 33%) than experiment. The best overall agreement with the experiment ( 7 % ) is displayed by atom_QDZh96. In Table 13, we compare experimental and theoretical A-values. QDZh96 values are within the experimental error bars except for A ( a 4 G 9 / 2 , a 4 F 9 / 2 ) and A ( b 4 D 7 / 2 , a 4 F 7 / 2 ) for which the theoretical A-values agree to within 12%. The average agreement between the experiment and the rest of the theoretical data is around 20%.

4.5. Radiative Lifetimes of the 3 d 6 ( 5 D ) 4 p Odd-Parity Levels

Radiative lifetimes of the Fe ii odd-parity levels 3 d 6 ( 5 D ) 4 p have been measured with a time-resolved non-linear laser-induced fluorescence technique, whereby an improved signal-to-noise ratio reduces uncertainties to 2–3% [46]. In Table 14, we compare these measurements with lifetimes derived from A-values from atom_TZ18_173 and atom_RKFB19_173. Theoretical lifetimes agree to 30 % but are, on average, around 14% shorter than measurements and well outside the aforementioned error bars.
Experimental E1 A-values were reported using the lifetimes in Table 14 and branching fractions determined with a Fourier transform spectrometer and a high-resolution grating spectrometer to precisions of 6% and 26% for the strong and weak transitions, respectively [46]. We compare these measurements with A-values from atom_TZ18_173 and atom_SRKFB19_173 in Figure 12. A-values from the former dataset with log A > 6 agree with the experiment to within 0.3 dex, while large discrepancies (as large as 3 dex) are displayed by the latter.

4.6. Fe ii  g f -Values

Oscillator strengths ( g f -values) for 142 Fe ii lines ( λ λ 4000 8000 ) have been derived from both laboratory and computed data and benchmarked with accurate spectra from the Sun and metal-poor stars [47]. These g f -values have been compared with data computed with extensive CI using the mcbp code civ3, including fine-tuning [43]. This calculation included 262 fine-structure levels from the 3 d 6 4 s , 3 d 7 , 3 d 5 4 s 2 , 3 d 6 4 p , and  3 d 5 4 s 4 p configurations. In general, good agreement was found except for a few ill-behaved transitions susceptible to effects difficult to constrain computationally such as severe cancellation due to CI mixing.
In Figure 13, we compare the observed g f -values (MB09) [47] with theoretical data from Ref. [43] (DH14), atom_TZ18_173, and atom_SRKFB19_173. If problematic transitions are excluded, the average agreement of DH14 and atom_TZ18_173 with MB09 is within 0.2 dex. Dataset atom_SRKFB19_173, on the other hand, contains several transitions with discrepancies greater than 0.5 dex.
Table 15 shows three transitions for which DH14 shows sizable discrepancies with respect to MB09. Datasets atom_TZ18_173 and atom_SRKFB19_173 do not perform much better; thus, the hindrance from cancellation effects in computational estimates put forward in Ref. [43] is reinforced. This table also lists three transitions we were not able to identify in atom_TZ18_173 or atom_SRKFB19_173 for which the agreement between DH14 and MB09 is reasonable.

4.7. Fe ii Branching Fractions

Branching fractions
B F j i = I ( j , i ) i I ( j , i ) = A ( j , i ) i A ( j , i )
for 121 Fe ii UV lines originating from levels of the 3 d 6 ( 5 D ) 4 p configuration have been measured in low-current spectra from Fe–Ne and Fe–Ar hollow cathode discharge lamps. The spectra were taken with a UW 3 m echelle spectrograph with a resolving power of 250,000 [48]. We compare these branching fractions in Figure 14 with those derived from A-values of the atom_TZ18_173 and atom_SRKFB19_173 atom datasets. For transitions with log ( B F expt ) > 2 , the branching fractions from atom_TZ18_173 agree with the experiment, on average, to better than 30%, while significantly larger discrepancies are displayed by atom_SRKFB19_173. However, most theoretical branching fractions, in general, do not match the accuracy implicit in the experimental error bars.

4.8. [Fe ii] Effective Collision Strengths

We have constructed three revised coll datasets containing ECS for the [Fe ii] 52-level atomic model of PyNeb: B15_52, TZ18_52, and SRKFB19_52. We determine for each transition the average ECS from these three datasets at T = 10 4  K, and as depicted in Figure 15, we compare the logarithmic differences with respect to this average. ECS from coll_B15_52 appear to lie below average, displaying differences as large as Δ log Υ 4 , while those from coll_SRKFB19_52 are above average with a scatter within 1.0 dex. The most interesting result is coll_TZ18_52 with an ECS average well within 0.5 dex. This outcome supports the findings in the spectrum fits (see Table 11), whereby the electron densities predicted by coll_B15_52 are higher, coll_SRKFB19_52 lower, and coll_TZ18_52 in between.
Due to the diagnostic potential of [Fe ii] emission lines in interstellar-medium plasmas, accurate ECS at low temperatures were recently computed for transitions among levels of the 3 d 6 ( 5 D ) 4 s a 6 D J ground term [36]. To obtain uncertainty estimates, different target models and relativistic R-matrix methods (bprm, icft, darc) were used in these calculations. In Figure 16, we compare selected ECS from [36] with those from coll datasets B15_52, TZ18_52, and SRKFB19_52. For this purpose, we use the PyNebgetOmega(T,j,i) function that allows ECS linear interpolation at the input temperature T; if the latter lies outside the temperature range of the tabulated ECS, the function lists the end value. This is the case of coll_B15_52 and coll_SRKFB19_52 for log T < 3 (see Figure 16). As expected, the best agreement is with coll_SRKFB19_52, although there are points outside the error bars. Large discrepancies are also found with coll_B15_52 and coll_TZ18_52 at higher temperatures ( log T > 4 ).

4.9. Fe ii Discussion

In the implementation of new datasets for [Fe ii], we found term misassignments for levels with total orbital angular momentum quantum number L = 5 in the atomic model of atom_SRKFB19_225 (see Table 6). A further comparison in Table 7 of A-values for transitions arising from the 3 d 6 4 s a 4 H J and 3 d 6 4 s b 2 H J ( J = 9 / 2 and 11 / 2 ) levels indicate large discrepancies with respect to other atom datasets (TZ18_225, DH11_52, B15_52, QDZa96, and QDZh96). This is an indication of faulty data that can be corroborated with the discrepancies displayed by dataset atom_SRKFB19_52 in the comparison of observed and computed A-value ratios carried out in Section 4.3.1.
This comparison also brings out observed lines with questionable identifications (e.g., λ λ 3968.27 , 4889.71 ) and unreliable line-intensity ratios—e.g., I ( λ 7637.52 ) / I ( λ 9051.95 ) , I ( λ 7686.93 ) / I ( λ 8891.93 ) , and  I ( λ 7733.13 ) / I ( λ 9033.49 ) —due to A-values with magnitude log A 3 . Theoretical A-value ratios agree to within 20–25% and, with respect to the observed line-intensity ratios, to within the error bars except for some outliers subject to strong level mixing, as pointed out in [40].
We have found atomic data files, namely atom_B15 and coll_B15, with faulty data caused again by the incorrect mapping to atomic models with different level numbering. This problem is exacerbated by data producers who use several multi-configuration methods implementing a variety of relativistic Hamiltonians (Dirac–Fock, Breit–Pauili, Pauli). This multi-code approach may enable accuracy estimates of the radiative and collisional data, but it often leads to level structures with conflicting numbering and identifications that are difficult to remap. The use of fine-tuning in atomic calculations to match the experimental level order does mitigate this problem.
In contrast to [Fe iii], for which the spectrum fits using different atomic datasets resulted in acceptable temperature (6%) and density (20%) uncertainties thus sustaining their reliability, for [Fe ii], the collisional datasets lead to temperature variations of around 15–20% and density discrepancies as large as a factor of four. However, the reasonable temperature and density accord between coll_SRKFB19_52 and coll_VVKFHF99_52 must be noted. These results appear to indicate that, for [Fe ii], the electron density is not a good fitting parameter or, as previously proposed [72,73,74], fluorescence continuum pumping might be a level-populating mechanism we have not taken into account. Furthermore, the extremely low critical density of level 3 d 7 a 4 F 9 / 2 (see Figure 8) and its sensitivity to both the radiative and collisional datasets (Figure 9) give this level a protagonistic role in the [Fe ii] spectrum formation, notwithstanding the poor accuracy of its small A-value.
The potential of line-ratio density diagnostics in [Fe ii] is compromised by the mismatch with the spectrum fits, and the temperature diagnostics are discouraged by the looming role of fluorescence pumping. In Section 4.3.3, we use the I ( λ 8892 ) / I ( λ 9267 ) ratio to illustrate a density diagnostic problem: the density variations in the range of interest ( 3 log n e 5  cm 3 ) lie within the observed line-intensity error band. As performed in Section 3.3.3 for the analysis of the [Fe iii] density diagnostics, an observational sample of low-density ionized nebulae with [Fe ii] detections would be helpful to constrain the predictions of both the radiative and collisional parameters. The lack of such a sample in the literature prevents us from drawing further conclusions.
We have performed extensive comparisons of computed and measured radiative data for both the 3 d 6 4 s , 3 d 7 , and  3 d 5 4 s 2 even-parity levels and the 3 d 6 4 p and 3 d 5 4 s 4 p odd-parity levels. Regarding even-parity levels (see Section 4.4), the computed radiative lifetimes agree to within 0.1 dex except for the singular 3 d 7 a 4 F 9 / 2 metastable level with log τ 4  s and an uncertainty of 0.4 dex. It must be noted that the A-values in atom_SRKB19_52 display questionable lifetimes for some of these levels (see Figure 11). The best overall agreement between theory and experiment ( 7 % ) is with the atom_QDZh96 dataset, while other atom datasets (B15_52, DH11_52, and TZ18_52) in general yield longer (as long as 33%) lifetimes than experiment. As discussed in [40], the large uncertainty of the 3 d 6 4 s b 2 H 11 / 2 lifetime is caused by imprecise level mixing. A-values derived from the experimental lifetimes, on average, agree to around 20% with those in the theoretical datasets.
Accurate (2–3%) radiative lifetimes have been measured for the 3 d 6 ( 5 D ) 4 p odd-parity levels [46], which have been compared in Section 4.5 with those derived from the atom_TZ18_173 and atom_SRKFB19_173 datasets. Theoretical lifetimes agree to within 30%, but on average they are around 14% shorter than measurements. E1 A-values have also been derived from these lifetime measurements, and those with log A > 6 agree with atom_TZ18_173 to within 0.3 dex while larger discrepancies (as large as 3 dex) are found in atom_SRKFB19_173. Moreover, the comparison of the observed g f -values [47] in Section 4.6 with those derived from the theoretical datasets shows good agreement (0.2 dex) with atom_TZ18_173 and discrepancies larger than 0.5 dex with atom_SRKFB19_173. A similar outcome is obtained in the comparison of measured branching ratios [48] in Section 4.7, whereby those derived from A-values in atom_TZ18_173 with log ( B F expt ) > 2 agree to better than 30%, while some from atom_SRKFB19_173 show larger discrepancies. In general, the theoretical radiative datasets do not reach anywhere near the accuracy quoted for the measured lifetimes, A-values, g f -values, and branching fractions.
The comparison of collisional datasets in Figure 15 of Section 4.8 gives a macro-measure of the overall ECS magnitudes that has a direct impact on the critical densities. In this comparison, we determined for each transition an average ECS using three coll datasets (B15_52, SRKFB19_52, and TZ18_52) and then plotted the dispersion for each dataset. Figure 15 shows that the smallest dispersion (within 0.5 dex) is by coll_TZ18_52 while differences in coll_B15_52 and coll_SRKFB19_52 are larger and, respectively, mostly negative and positive. The fairly large ECS discrepancies found in Figure 16 for transitions within the 3 d 6 ( 5 D ) 4 s a 6 D ground term at low temperatures gives a further indication of the patchy quality of the collisional data for Fe ii.
From the extensive data tests hereby carried out on the Fe ii datasets, we find the radiative data in atom_SRKFB19_52 to be faulty, and this could indeed have an impact on the coll_SRKFB19_52 data. We would therefore select as the PyNeb 1.1.17 defaults the atom_TZ18_173 and coll_ TZ18_173 datasets.

5. Conclusions

We have performed extensive data evaluations of the levels, atom, and coll files for [Fe iii] and [Fe ii] in PyNeb 1.1.16 to examine their worthiness in nebular spectral modeling and to select the package defaults. In most cases, we have reconstructed the datasets from the sources and introduced newly published data, thus leading to a more reliable atomic database for these species in the new release PyNeb 1.1.17.
The present data assessment comprises: fits of high-resolution spectra from the HH 202S and HH 204 Herbig–Haro objects taken with the Ultraviolet Visual Echelle Spectrograph [9,11]; comparisons of theoretical and observed line-intensity ratios; and comparisons with measurements of radiative lifetimes, A-values, g f -values, and branching ratios. Observed lines may suffer from misidentifications, line blending, and telluric or instrumental contamination. On the other hand, the computation of accurate radiative data for lowly ionized ions of the iron group with open 3 d n shells is hindered by strong configuration mixing, double-excitation representations, and subtle core–valence correlation involving closed and open subshells; thus, broadly speaking, we would not expect an A-value accuracy of better than 20%. Similarly, the ECS at low temperatures are affected by complex resonance structures sensitive to target models and channel convergence that seriously compromise accuracy ratings.
For Fe iii, we introduced a new levels file from the NIST database with more complete level assignments and accurate energies. The previous version of the levels file is in the deprecated directory, renamed as fe_iii_levels_2023.dat. As described in Section 3.2, we implemented the new atom files Q96_34, DH09_34, BBQ10_34, BB14_34, and BB14_144 and the coll files Z96_34, BBQ10_34, BB14_34, Z96_144, and BB14_144. In PyNeb 1.1.17 atom_BBQ10_34, atom_BB14_34, coll_Z96_34 are not included as they are, respectively, made redundant by atom_BBQ10, atom_BB14_144, and coll_Z96_144 while coll_BB14 and coll_BBQ10 are deprecated.
A comparison of the observed line-intensity ratios for transitions arising from the same upper level in [Fe iii] with the corresponding A-value ratios constrains the overall accuracy of the latter to around 20–30%. Spectrum fits with the new datasets led to temperatures and densities within 6% and 20%, respectively, validating their usefulness. Radiative lifetimes computed with these atom datasets agree to within 0.1 dex. A revision of coll files brings out the incompleteness of the Z96_34 dataset and its poor accuracy for ECS with log Υ < 1 . For [Fe iii], we have selected the atom_BB14_144 and coll_BB14_144 datasets as the defaults in PyNeb 1.1.17.
For Fe ii, we have implemented the atom files VVKFHF99_52, DH11_52, B15_52, TZ18_52, SRKFB19_52, TZ18_173, SRKFB19_173 and coll files VVKFHF99_52, B15_52, TZ18_52, SRKFB19_52, TZ18_173, and SRKFB19_173. After extensive tests and much pondering, the atom_SRKFB19_52 and atom_SRKFB19_173 datasets were found to contain faulty A-values, while atom_VVKFHF99_52 and coll_VVKFHF99_52 have missing data; therefore, these files are not included in PyNeb 1.1.17. Apart from this drawback, the theoretical A-value ratios for transitions arising from a common upper level derived from the rest of the atom files agree to within 20–25%, and most lie within the experimental error bars. The spectrum fits yielded temperature and density variations of around 20% and a factor of 4, respectively, which appear to indicate that the density is not a good fitting parameter and fluorescence continuum pumping might be a level-populating contributor not taken into account. However, the poor theoretical characterization of the key level 3 d 7 a 4 F 9 / 2 , inasmuch as its inaccurate A-value and extremely low critical density sensitive to both the radiative and collisional datasets, deters us from more conclusive remarks. Extensive comparisons of the radiative data for both forbidden and allowed transitions with experiment bring out the computational difficulties to match the measured accuracy (better than 10%); thus, a global theoretical A-value accuracy of 20–30% would be statistically sound. Finally, the examined coll datasets lead to different critical densities around T e = 10 4  K, a situation that is difficult  to analyze.
In PyNeb 1.1.17, the atom files VVKFHF99_52, TZ18_52 and coll files VVKFHF99_52, SRKFB19_52, SRKFB19_173 are made redundant by atom_VVKFHF99, atom_TZ18_173, coll_VVKFHF99, and coll_SRKFB19. Furthermore, atom_B15, coll_B15, atom_VVKFHF99, coll_VVKFHF99 are deprecated. We select atom_TZ18_173 and coll_TZ18_173 as defaults.
The status of all the data files for Fe iii and Fe ii in PyNeb 1.1.16 and 1.1.17 are summarized in Table A1 of Appendix A.
A final recommendation emerging from the present work must be put forward regarding the curation of the atomic databases associated with spectral modeling codes. The potential of the latter in astrophysical analyses relies in great part on the integrity (accuracy and completeness) of their tabulations of radiative and collisional rates, which, as shown here, demand detailed evaluation and upgrading. For complex ionic spectra such as [Fe iii] and [Fe ii], this may imply a lengthy endeavor.

Author Contributions

C.M. (Claudio Mendoza) conceptualized and coordinated the full paper. J.E.M.-D. and J.G.-R. provided the observational data and extensive discussions about their validity and provenance. M.B. provided data, performed calculations, and led discussions about fluorescence continuum pumping in [Fe ii]. C.M. (Christophe Morisset) researched data, performed ancillary calculations, and software support regarding the PyNeb package. Investigation, methodology, analysis, and validation: C.M. (Claudio Mendoza), J.E.M.-D., J.G.-R., M.B., C.M. (Christophe Morisset). Writing and original draft preparation: C.M. (Claudio Mendoza), J.E.M.-D. Review and editing: C.M. (Claudio Mendoza), J.E.M.-D., C.M. (Christophe Morisset), J.G.-R., M.B. Visualization: C.M. (Claudio Mendoza), J.E.M.-D. All authors have read and agreed to the published version of the manuscript.

Funding

J.E.M.-D. acknowledges funding from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) in the form of an Emmy Noether Research Group (grant number KR4598/2-1, PI Kathryn Kreckel). Ch.M. acknowledges the support from UNAM/DGAPA/PAPIIT grant IN101220. J.G.-R. acknowledges support from an Advanced Fellowship under the Severo Ochoa excellence program CEX2019-000920-S; financial support from the Canarian Agency for Research, Innovation, and Information Society (ACIISI) of the Canary Islands Government and the European Regional Development Fund (ERDF) under grant with reference ProID2021010074; and support under grant P/308614 financed by funds transferred from the Spanish Ministry of Science, Innovation, and Universities (MCIU) charged to the General State Budgets and with funds transferred from the General Budgets of the Autonomous Community of the Canary Islands by the MCIU.

Data Availability Statement

The data files revised and compiled for this revision are publicly available in PyNeb 1.1.17.

Acknowledgments

We acknowledge extensive private communications and data exchange with Catherine Ramsbottom (Queen’s University Belfast, UK) regarding the radiative and collisional data for Fe ii. We are indebted to Grażyna Stasińska (Observatoire de Paris, France) for data exchange and useful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

    The following abbreviations are used in this manuscript:
CDSCentre de Données astronomiques de Strasbourg
NISTNational Institute os Standards and Technology

Appendix A

Table A1. Status of the Fe iii and Fe ii atomic data files in the PyNeb library for versions 1.1.16 and 1.1.17. ‘X’ means the file is present; ‘*’ means the file is present and is the default; ‘D’ means the file is deprecated (but may still be accessible with the command pn.atomicData.includeDeprecatedPath()).
Table A1. Status of the Fe iii and Fe ii atomic data files in the PyNeb library for versions 1.1.16 and 1.1.17. ‘X’ means the file is present; ‘*’ means the file is present and is the default; ‘D’ means the file is deprecated (but may still be accessible with the command pn.atomicData.includeDeprecatedPath()).
Fe iii File Namev1.1.16v1.1.17Fe ii File Namev1.1.16v1.1.17
fe_iii_atom_BB14_144.dat *fe_ii_atom_TZ18_173.dat *
fe_iii_atom_BBQ10.datXXfe_ii_atom_B15.dat*D
fe_iii_atom_DH09_34.dat Xfe_ii_atom_B15_52.dat X
fe_iii_atom_NP96.datXXfe_ii_atom_DH11_52.dat X
fe_iii_atom_Q96_34.dat Xfe_ii_atom_SRKFB19.datXD
fe_iii_atom_Q96_J00.dat*Xfe_ii_atom_VVKFHF99.datXD
fe_iii_coll_BB14_144.dat *fe_ii_coll_TZ18_173.dat *
fe_iii_coll_BB14.datXDfe_ii_coll_B15.dat*D
fe_iii_coll_BBQ10.datXDfe_ii_coll_B15_52.dat X
fe_iii_coll_BBQ10_34.dat Xfe_ii_coll_B15_old.datXD
fe_iii_coll_Q96.datXXfe_ii_coll_SRKFB19.datXX
fe_iii_coll_Z96.dat*Xfe_ii_coll_VVKFHF99.datXD
fe_iii_coll_Z96_144.dat X

Notes

1
2
3

References

  1. Luridiana, V.; Morisset, C.; Shaw, R.A. PyNeb: A New Software for the Analysis of Emission Lines. In IAU Symposium; Cambridge University Press: Cambridge, UK, 2012; Volume 283, pp. 422–423. [Google Scholar] [CrossRef] [Green Version]
  2. Luridiana, V.; Morisset, C.; Shaw, R.A. PyNeb: A new tool for analyzing emission lines. I. Code description and validation of results. Astron. Astrophys. 2015, 573, A42. [Google Scholar] [CrossRef] [Green Version]
  3. Peimbert, M.; Peimbert, A.; Delgado-Inglada, G. Nebular Spectroscopy: A Guide on H II Regions and Planetary Nebulae. Publ. Astron. Soc. Pac. 2017, 129, 082001. [Google Scholar] [CrossRef]
  4. Morisset, C.; Luridiana, V.; García-Rojas, J.; Gómez-Llanos, V.; Bautista, M.; Mendoza, C. Atomic Data Assessment with PyNeb. Atoms 2020, 8, 66. [Google Scholar] [CrossRef]
  5. Juan de Dios, L.; Rodríguez, M. Atomic data and the density structures of planetary nebulae. Mon. Not. R. Astron. Soc. 2021, 507, 5331–5339. [Google Scholar] [CrossRef]
  6. Perinotto, M.; Bencini, C.G.; Pasquali, A.; Manchado, A.; Rodriguez Espinosa, J.M.; Stanga, R. The iron abundance in four planetary nebulae. Astron. Astrophys. 1999, 347, 967–974. [Google Scholar]
  7. Delgado Inglada, G.; Rodríguez, M.; Mampaso, A.; Viironen, K. The Iron Abundance in Galactic Planetary Nebulae. Astrophys. J. 2009, 694, 1335–1348. [Google Scholar] [CrossRef]
  8. Delgado-Inglada, G.; Mesa-Delgado, A.; García-Rojas, J.; Rodríguez, M.; Esteban, C. The Fe/Ni ratio in ionized nebulae: Clues on dust depletion patterns. Mon. Not. R. Astron. Soc. 2016, 456, 3855–3865. [Google Scholar] [CrossRef] [Green Version]
  9. Mesa-Delgado, A.; Esteban, C.; García-Rojas, J.; Luridiana, V.; Bautista, M.; Rodríguez, M.; López-Martín, L.; Peimbert, M. Properties of the ionized gas in HH 202—II. Results from echelle spectrophotometry with Ultraviolet Visual Echelle Spectrograph. Mon. Not. R. Astron. Soc. 2009, 395, 855–876. [Google Scholar] [CrossRef]
  10. Méndez-Delgado, J.E.; Esteban, C.; García-Rojas, J.; Henney, W.J.; Mesa-Delgado, A.; Arellano-Córdova, K.Z. Photoionized Herbig-Haro objects in the Orion Nebula through deep high-spectral resolution spectroscopy—I. HH 529 II and III. Mon. Not. R. Astron. Soc. 2021, 502, 1703–1739. [Google Scholar] [CrossRef]
  11. Méndez-Delgado, J.E.; Henney, W.J.; Esteban, C.; García-Rojas, J.; Mesa-Delgado, A.; Arellano-Córdova, K.Z. Photoionized Herbig-Haro Objects in the Orion Nebula through Deep High Spectral Resolution Spectroscopy. II. HH 204. Astrophys J. 2021, 918, 27. [Google Scholar] [CrossRef]
  12. D’Odorico, S.; Cristiani, S.; Dekker, H.; Hill, V.; Kaufer, A.; Kim, T.; Primas, F. Performance of UVES, the echelle spectrograph for the ESO VLT and highlights of the first observations of stars and quasars. In Discoveries and Research Prospects from 8- to 10-Meter-Class Telescopes; Bergeron, J., Ed.; Society of Photo-Optical Instrumentation Engineers (SPIE): Bellingham, WA, USA, 2000; Volume 4005, pp. 121–130. [Google Scholar] [CrossRef]
  13. Garstang, R.H. Transition probabilities for forbidden lines of Fe II. Mon. Not. R. Astron. Soc. 1962, 124, 321. [Google Scholar] [CrossRef]
  14. Kurucz, R.L.; Peytremann, E. A Table of Semiempirical gf Values. Pt 1: Wavelengths: 5.2682 NM to 272.3380 nm; Pt 2: Wavelengths: 272.3395 NM to 599.3892 nm; Pt 3: Wavelengths: 599.4004 NM to 9997.2746 NM; SAO Special Report; Smithsonian Astrophysical Observatory: Cambridge, MA, USA, 1975. [Google Scholar]
  15. Nussbaumer, H.; Storey, P.J. Atomic data for Fe II. Astron. Astrophys. 1980, 89, 308–313. [Google Scholar]
  16. Nussbaumer, H.; Pettini, M.; Storey, P.J. Sextet transitions in Fe II. Astron. Astrophys. 1981, 102, 351–358. [Google Scholar]
  17. Hummer, D.G.; Berrington, K.A.; Eissner, W.; Pradhan, A.K.; Saraph, H.E.; Tully, J.A. Atomic data from the IRON project. I. Goals and methods. Astron. Astrophys 1993, 279, 298–309. [Google Scholar]
  18. Eissner, W.; Jones, M.; Nussbaumer, H. Techniques for the calculation of atomic structures and radiative data including relativistic corrections. Comput. Phys. Commun. 1974, 8, 270–306. [Google Scholar] [CrossRef]
  19. Hibbert, A. CIV3—A general program to calculate configuration interaction wave functions and electric-dipole oscillator strengths. Comput. Phys. Commun. 1975, 9, 141–172. [Google Scholar] [CrossRef]
  20. Cowan, R.D. The Theory of Atomic Structure and Spectra; University of California: Berkeley, CA, USA, 1981. [Google Scholar]
  21. Nahar, S.N. Atomic data from the Iron Project. VII. Radiative dipole transition probabilities for Fe II. Astron. Astrophys. 1995, 293, 967–977. [Google Scholar]
  22. Nahar, S.N.; Pradhan, A.K. Atomic data from the Iron Project. XVII. Radiative transition probabilities for dipole allowed and forbidden transitions in Fe III. Astron. Astrophys. Suppl. Ser. 1996, 119, 509–522. [Google Scholar] [CrossRef]
  23. Quinet, P.; Le Dourneuf, M.; Zeippen, C.J. Atomic data from the IRON Project. XIX. Radiative transition probabilities for forbidden lines in Fe II. Astron. Astrophys. Suppl. Ser. 1996, 120, 361–371. [Google Scholar] [CrossRef]
  24. Zhang, H.L.; Pradhan, A.K. Atomic data from the Iron Project. VI. Collision strengths and rate coefficients for Fe II. Astron. Astrophys. 1995, 293, 953–966. [Google Scholar]
  25. Zhang, H. Atomic data from the Iron Project. XVIII. Electron impact excitation collision strengths and rate coefficients for Fe III. Astron. Astrophys. Suppl. Ser. 1996, 119, 523–528. [Google Scholar] [CrossRef]
  26. Bautista, M.A.; Pradhan, A.K. Atomic data from the Iron project. XIII. Electron excitation rates and emissivity ratios for forbidden transitions in Ni II and Fe II. Astron. Astrophys. Suppl. Ser. 1996, 115, 551–559. [Google Scholar]
  27. Bautista, M.A.; Ballance, C.P.; Quinet, P. Atomic Data and Spectral Model for Fe III. Astrophys. J. 2010, 718, L189–L193. [Google Scholar] [CrossRef] [Green Version]
  28. Badnell, N.R.; Ballance, C.P. Electron-impact Excitation of Fe2+: A Comparison of Intermediate Coupling Frame Transformation, Breit-Pauli and Dirac R-matrix Calculations. Astrophys. J. 2014, 785, 99. [Google Scholar] [CrossRef] [Green Version]
  29. Bautista, M.A.; Fivet, V.; Ballance, C.; Quinet, P.; Ferland, G.; Mendoza, C.; Kallman, T.R. Atomic Data and Spectral Model for Fe II. Astrophys. J. 2015, 808, 174. [Google Scholar] [CrossRef] [Green Version]
  30. Tayal, S.S.; Zatsarinny, O. Electron-impact excitation of forbidden and allowed transitions in Fe II. Phys. Rev. A 2018, 98, 012706. [Google Scholar] [CrossRef] [Green Version]
  31. Smyth, R.T.; Ramsbottom, C.A.; Keenan, F.P.; Ferland, G.J.; Ballance, C.P. Towards converged electron-impact excitation calculations of low-lying transitions in Fe II. Mon. Not. R. Astron. Soc. 2019, 483, 654–663. [Google Scholar] [CrossRef] [Green Version]
  32. Berrington, K.A.; Eissner, W.B.; Norrington, P.H. RMATRX1: Belfast atomic R-matrix codes. Comput. Phys. Commun. 1995, 92, 290–420. [Google Scholar] [CrossRef]
  33. Zatsarinny, O. BSR: B-spline atomic R-matrix codes. Comput. Phys. Commun. 2006, 174, 273–356. [Google Scholar] [CrossRef]
  34. Ait-Tahar, S.; Grant, I.P.; Norrington, P.H. Electron scattering by Fe XXII within the Dirac R-matrix approach. Phys. Rev. A 1996, 54, 3984–3989. [Google Scholar] [CrossRef]
  35. Griffin, D.C.; Badnell, N.R.; Pindzola, M.S. R-matrix electron-impact excitation cross sections in intermediate coupling: An MQDT transformation approach. J. Phys. B At. Mol. Phys. 1998, 31, 3713–3727. [Google Scholar] [CrossRef]
  36. Wan, Y.; Favreau, C.; Loch, S.D.; McLaughlin, B.M.; Qi, Y.; Stancil, P.C. Electron-impact fine-structure excitation of Fe II at low temperature. Mon. Not. R. Astron. Soc. 2019, 485, 2252–2258. [Google Scholar] [CrossRef]
  37. Deb, N.C.; Hibbert, A. Oscillator strengths for transitions among Fe III levels belonging to the three lowest configurations. J. Phys. Conf. Ser. 2008, 130, 012006. [Google Scholar] [CrossRef] [Green Version]
  38. Deb, N.C.; Hibbert, A. Electric quadrupole and magnetic dipole transitions among 3d6 levels of Fe III. J. Phys. B Atom. Mol. Phys. 2009, 42, 065003. [Google Scholar] [CrossRef]
  39. Deb, N.C.; Hibbert, A. Weighted f-values, A-values, and line strengths for the E1 transitions among 3d 6, 3d 54s, and 3d 54p levels of Fe III. Atom. Data Nucl. Data Tables. 2009, 95, 184–303. [Google Scholar] [CrossRef]
  40. Deb, N.C.; Hibbert, A. Importance of level mixing on accurate [Fe II] transition rates. Astron. Astrophys. 2010, 524, A54. [Google Scholar] [CrossRef] [Green Version]
  41. Deb, N.C.; Hibbert, A. Calculation of Intensity Ratios of Observed Infrared [Fe II] Lines. Astrophys. J. 2010, 711, L104–L107. [Google Scholar] [CrossRef]
  42. Deb, N.C.; Hibbert, A. Radiative transition rates for the forbidden lines in Fe II. Astron. Astrophys. 2011, 536, A74. [Google Scholar] [CrossRef] [Green Version]
  43. Deb, N.C.; Hibbert, A. log gf values for astrophysically important transitions Fe II. Astron. Astrophys. 2014, 561, A32. [Google Scholar] [CrossRef] [Green Version]
  44. Quinet, P. Transition probabilities for forbidden lines of Fe III. Astron. Astrophys. Suppl. Ser. 1996, 116, 573–578. [Google Scholar] [CrossRef] [Green Version]
  45. Fivet, V.; Quinet, P.; Bautista, M.A. Radiative rates for forbidden M1 and E2 transitions of astrophysical interest in doubly ionized iron-peak elements. Astron. Astrophys 2016, 585, A121. [Google Scholar] [CrossRef] [Green Version]
  46. Schnabel, R.; Schultz-Johanning, M.; Kock, M. Fe II lifetimes and transition probabilities. Astron. Astrophys. 2004, 414, 1169–1176. [Google Scholar] [CrossRef] [Green Version]
  47. Meléndez, J.; Barbuy, B. Both accurate and precise gf-values for Fe II lines. Astron. Astrophys. 2009, 497, 611–617. [Google Scholar] [CrossRef] [Green Version]
  48. Den Hartog, E.A.; Lawler, J.E.; Sneden, C.; Cowan, J.J.; Brukhovesky, A. Atomic Transition Probabilities for UV and Blue Lines of Fe II and Abundance Determinations in the Photospheres of the Sun and Metal-poor Star HD 84937. Astrophys. J. Suppl, Ser. 2019, 243, 33. [Google Scholar] [CrossRef] [Green Version]
  49. Rostohar, D.; Derkatch, A.; Hartman, H.; Johansson, S.; Lundberg, H.; Mannervik, S.; Norlin, L.O.; Royen, P.; Schmitt, A. Lifetime Measurements of Metastable States in Fe+. Phys. Rev. Lett. 2001, 86, 1466–1469. [Google Scholar] [CrossRef] [Green Version]
  50. Hartman, H.; Derkatch, A.; Donnelly, M.P.; Gull, T.; Hibbert, A.; Johansson, S.; Lundberg, H.; Mannervik, S.; Norlin, L.O.; Rostohar, D.; et al. The FERRUM Project: Experimental transition probabilities of [Fe II] and astrophysical applications. Astron. Astrophys. 2003, 397, 1143–1149. [Google Scholar] [CrossRef] [Green Version]
  51. Gurell, J.; Hartman, H.; Blackwell-Whitehead, R.; Nilsson, H.; Bäckström, E.; Norlin, L.O.; Royen, P.; Mannervik, S. The FERRUM project: Transition probabilities for forbidden lines in [Fe II] and experimental metastable lifetimes. Astron. Astrophys. 2009, 508, 525–529. [Google Scholar] [CrossRef] [Green Version]
  52. Johansson, S.; Zethson, T.; Hartman, H.; Ekberg, J.O.; Ishibashi, K.; Davidson, K.; Gull, T. New forbidden and fluorescent Fe III lines identified in HST spectra of eta Carinae. Astron. Astrophys. 2000, 361, 977–981. [Google Scholar]
  53. Badnell, N.R. A Breit-Pauli distorted wave implementation for AUTOSTRUCTURE. Comput. Phys. Commun. 2011, 182, 1528–1535. [Google Scholar] [CrossRef]
  54. Parpia, F.A.; Fischer, C.F.; Grant, I.P. GRASP92: A package for large-scale relativistic atomic structure calculations. Comput. Phys. Commun. 1996, 94, 249–271. [Google Scholar] [CrossRef]
  55. Kramida, A.; Yu, R.; Reader, J.; NIST ASD Team. NIST Atomic Spectra Database, version 5.9; Reader, J., Ed.; National Institute of Standards and Technology: Gaithersburg, MD, USA, 2021. [Google Scholar]
  56. Hibbert, A.; Glass, R.; Froese Fischer, C. A general program for computing angular integrals of the Breit-Pauli Hamiltonian. Comput. Phys. Commun. 1991, 64, 455–472. [Google Scholar] [CrossRef]
  57. Mendoza, C.; Bautista, M.A.; Deprince, J.; García, J.A.; Gatuzz, E.; Gorczyca, T.W.; Kallman, T.R.; Palmeri, P.; Quinet, P.; Witthoeft, M.C. The XSTAR Atomic Database. Atoms 2021, 9, 12. [Google Scholar] [CrossRef]
  58. Storey, P.J.; Sochi, T. Electron temperatures and free-electron energy distributions of nebulae from C II dielectronic recombination lines. Mon. Not. R. Astron. Soc. 2013, 430, 599–610. [Google Scholar] [CrossRef]
  59. Sochi, T. Atomic and Molecular Aspects of Astronomical Spectra. Ph.D. Thesis, University College, London, UK, 2012. [Google Scholar]
  60. Peimbert, A. The Chemical Composition of the 30 Doradus Nebula Derived from Very Large Telescope Echelle Spectrophotometry. Astrophys. J. 2003, 584, 735–750. [Google Scholar] [CrossRef] [Green Version]
  61. García-Rojas, J.; Esteban, C.; Peimbert, A.; Peimbert, M.; Rodríguez, M.; Ruiz, M.T. Deep echelle spectrophotometry of S 311, a Galactic HII region located outside the solar circle. Mon. Not. R. Astron. Soc. 2005, 362, 301–312. [Google Scholar] [CrossRef] [Green Version]
  62. García-Rojas, J.; Esteban, C.; Peimbert, M.; Costado, M.T.; Rodríguez, M.; Peimbert, A.; Ruiz, M.T. Faint emission lines in the Galactic H II regions M16, M20 and NGC 3603. Mon. Not. R. Astron. Soc. 2006, 368, 253–279. [Google Scholar] [CrossRef] [Green Version]
  63. García-Rojas, J.; Esteban, C.; Peimbert, A.; Rodríguez, M.; Peimbert, M.; Ruiz, M.T. The chemical composition of the Galactic H II regions M8 and M17. A revision based on deep VLT echelle spectrophotometry. Rev. Mex. Astron. Astrofis. 2007, 43, 3–31. [Google Scholar]
  64. López-Sánchez, Á.R.; Esteban, C.; García-Rojas, J.; Peimbert, M.; Rodríguez, M. The Localized Chemical Pollution in NGC 5253 Revisited: Results from Deep Echelle Spectrophotometry. Astrophys. J. 2007, 656, 168–185. [Google Scholar] [CrossRef]
  65. Esteban, C.; Bresolin, F.; Peimbert, M.; García-Rojas, J.; Peimbert, A.; Mesa-Delgado, A. Keck HIRES Spectroscopy of Extragalactic H II Regions: C and O Abundances from Recombination Lines. Astrophys. J. 2009, 700, 654–678. [Google Scholar] [CrossRef]
  66. Esteban, C.; García-Rojas, J.; Carigi, L.; Peimbert, M.; Bresolin, F.; López-Sánchez, A.R.; Mesa-Delgado, A. Carbon and oxygen abundances from recombination lines in low-metallicity star-forming galaxies. Implications for chemical evolution. Mon. Not. R. Astron. Soc. 2014, 443, 624–647. [Google Scholar] [CrossRef] [Green Version]
  67. Esteban, C.; Fang, X.; García-Rojas, J.; Toribio San Cipriano, L. The radial abundance gradient of oxygen towards the Galactic anti-centre. Mon. Not. R. Astron. Soc. 2017, 471, 987–1004. [Google Scholar] [CrossRef] [Green Version]
  68. Esteban, C.; Bresolin, F.; García-Rojas, J.; Toribio San Cipriano, L. Carbon, nitrogen, and oxygen abundance gradients in M101 and M31. Mon. Not. R. Astron. Soc. 2020, 491, 2137–2155. [Google Scholar] [CrossRef]
  69. Esteban, C.; García-Rojas, J. Revisiting the radial abundance gradients of nitrogen and oxygen of the Milky Way. Mon. Not. R. Astron. Soc. 2018, 478, 2315–2336. [Google Scholar] [CrossRef]
  70. Domínguez-Guzmán, G.; Rodríguez, M.; García-Rojas, J.; Esteban, C.; Toribio San Cipriano, L. The homogeneity of chemical abundances in H II regions of the Magellanic Clouds. Mon. Not. R. Astron. Soc. 2022, 517, 4497–4514. [Google Scholar] [CrossRef]
  71. Verner, E.M.; Verner, D.A.; Korista, K.T.; Ferguson, J.W.; Hamann, F.; Ferland, G.J. Numerical Simulations of Fe II Emission Spectra. Astrophys. J. Suppl. Ser. 1999, 120, 101–112. [Google Scholar] [CrossRef] [Green Version]
  72. Rodríguez, M. Fluorescence of [Fe II] in H II regions. Astron. Astrophys. 1999, 348, 222–226. [Google Scholar]
  73. Verner, E.M.; Verner, D.A.; Baldwin, J.A.; Ferland, G.J.; Martin, P.G. Continuum Pumping of [Fe II] in the Orion Nebula. Astrophys. J. 2000, 543, 831–839. [Google Scholar] [CrossRef]
  74. Hartman, H. Fluorescence in Astrophysical Plasmas. In New Trends in Atomic and Molecular Physics; Springer: Berlin/Heidelberg, Germany, 2013; p. 189. [Google Scholar] [CrossRef] [Green Version]
  75. Sarkar, A.; Ferland, G.J.; Chatzikos, M.; Guzmán, F.; van Hoof, P.A.M.; Smyth, R.T.; Ramsbottom, C.A.; Keenan, F.P.; Ballance, C.P. Improved Fe II Emission-line Models for AGNs Using New Atomic Data Sets. Astrophys. J. 2021, 907, 12. [Google Scholar] [CrossRef]
Figure 1. Density behavior of the [Fe iii] I ( λ 4658 ) / I ( λ 4702 ) emissivity ratio at T e = 7.9 × 10 3  K for different atomic datasets. Left panel: atom_DH09_34 dataset with three coll datasets ( Z96_34, BBQ10_34, and BB14_34). Right panel: coll_Z96_34 dataset with four atom datasets (DH09_34, BBQ10_34, BB14_34, and Q96_34).
Figure 1. Density behavior of the [Fe iii] I ( λ 4658 ) / I ( λ 4702 ) emissivity ratio at T e = 7.9 × 10 3  K for different atomic datasets. Left panel: atom_DH09_34 dataset with three coll datasets ( Z96_34, BBQ10_34, and BB14_34). Right panel: coll_Z96_34 dataset with four atom datasets (DH09_34, BBQ10_34, BB14_34, and Q96_34).
Atoms 11 00063 g001
Figure 2. Observed [Fe iii] I ( λ 4658 ) / I ( λ 4702 ) line-intensity ratio from several low-density sources as a function of the [S ii] density. The theoretical emissivity ratios determined by PyNeb with dataset atom_DH09_34 and three coll datasets (BBQ10_34, BB14_34, and Z96_34) are also depicted.
Figure 2. Observed [Fe iii] I ( λ 4658 ) / I ( λ 4702 ) line-intensity ratio from several low-density sources as a function of the [S ii] density. The theoretical emissivity ratios determined by PyNeb with dataset atom_DH09_34 and three coll datasets (BBQ10_34, BB14_34, and Z96_34) are also depicted.
Atoms 11 00063 g002
Figure 3. Temperature behavior of the [Fe iii] I ( λ 4658 ) / I ( λ 3240 ) emissivity ratio at n e = 2.1 × 10 4  cm 3 for different atomic datasets. Left panel: dataset atom_DH09_34 with three coll datasets (Z96_34, BBQ10_34, and BB14_34). Right panel: dataset coll_Z96_34 with four atom datasets (DH09_34, BBQ10_34, BB14_34, and Q96_34).
Figure 3. Temperature behavior of the [Fe iii] I ( λ 4658 ) / I ( λ 3240 ) emissivity ratio at n e = 2.1 × 10 4  cm 3 for different atomic datasets. Left panel: dataset atom_DH09_34 with three coll datasets (Z96_34, BBQ10_34, and BB14_34). Right panel: dataset coll_Z96_34 with four atom datasets (DH09_34, BBQ10_34, BB14_34, and Q96_34).
Atoms 11 00063 g003
Figure 4. Comparison of Fe iii radiative lifetimes computed from the A-values in the atom datasets BB14_34, BBQ10_34, DH09_34, and Q96_34 atom files. For each level, lifetime differences for these datasets with respect to the average value are plotted as a function of the average value.
Figure 4. Comparison of Fe iii radiative lifetimes computed from the A-values in the atom datasets BB14_34, BBQ10_34, DH09_34, and Q96_34 atom files. For each level, lifetime differences for these datasets with respect to the average value are plotted as a function of the average value.
Atoms 11 00063 g004
Figure 5. Comparison of radiative lifetimes for 44 levels of the 3 d 5 4 p configuration in Fe iii derived from the A-values from atom_BB14_144 ( τ 1 ) and Ref. [39] ( τ 2 ).
Figure 5. Comparison of radiative lifetimes for 44 levels of the 3 d 5 4 p configuration in Fe iii derived from the A-values from atom_BB14_144 ( τ 1 ) and Ref. [39] ( τ 2 ).
Atoms 11 00063 g005
Figure 6. Comparison of [Fe iii] ECS at T = 10 4 K in the coll_BBQ10_34 and coll_Z96_34 datasets with respect to those in coll_BB14_34.
Figure 6. Comparison of [Fe iii] ECS at T = 10 4 K in the coll_BBQ10_34 and coll_Z96_34 datasets with respect to those in coll_BB14_34.
Atoms 11 00063 g006
Figure 7. Ratio of the normalized theoretical emissivity and observed intensity for lines in HH 202S as a function of the upper-level energy of the transition. Left panel: spectral fit with datasets atom_TZ18_52 and coll_TZ18_52 at T e = 7.94 × 10 3 K and n e = 3.29 × 10 4 cm 3 . Right panel: spectral fit with datasets atom_DH11_52 and coll_SRKFB19_52 at T e = 11.1 × 10 3 K and n e = 2.58 × 10 4 cm 3 .
Figure 7. Ratio of the normalized theoretical emissivity and observed intensity for lines in HH 202S as a function of the upper-level energy of the transition. Left panel: spectral fit with datasets atom_TZ18_52 and coll_TZ18_52 at T e = 7.94 × 10 3 K and n e = 3.29 × 10 4 cm 3 . Right panel: spectral fit with datasets atom_DH11_52 and coll_SRKFB19_52 at T e = 11.1 × 10 3 K and n e = 2.58 × 10 4 cm 3 .
Atoms 11 00063 g007
Figure 8. Critical densities for the lower 52 levels of [Fe ii] predicted at T = 10 4  K by using ECS from the coll datasets B15_52, SRKFB19_52, and TZ18_52 and A-values from atom_TZ18_52.
Figure 8. Critical densities for the lower 52 levels of [Fe ii] predicted at T = 10 4  K by using ECS from the coll datasets B15_52, SRKFB19_52, and TZ18_52 and A-values from atom_TZ18_52.
Atoms 11 00063 g008
Figure 9. Critical density of level 3 d 7 a 4 F 9 / 2 of [Fe ii] as a function of temperature computed with different atomic datasets. Left panel: atom_TZ18_52 file and three coll files: B15_52, SRKFB19_52, and TZ18_52. Right panel: coll_TZ18_52 file and three atom files: TZ18_52, B15_52, and DH11_52.
Figure 9. Critical density of level 3 d 7 a 4 F 9 / 2 of [Fe ii] as a function of temperature computed with different atomic datasets. Left panel: atom_TZ18_52 file and three coll files: B15_52, SRKFB19_52, and TZ18_52. Right panel: coll_TZ18_52 file and three atom files: TZ18_52, B15_52, and DH11_52.
Atoms 11 00063 g009
Figure 10. Density behavior of the [Fe ii] I ( λ 8892 ) / I ( λ 9267 ) emissivity ratio using different atomic datasets. The observed ratio in HH 202S is 2.36 ± 0.42 , which is depicted as a gray band. Left panel: atom_TZ18_52 radiative dataset with four coll datasets (SRKFB19_52, B15_52, TZ18_52, VVKFHF99_52) at T e = 10 4  K. Right panel: atom_TZ18_52 and coll_TZ18_52 datasets at four electron temperatures.
Figure 10. Density behavior of the [Fe ii] I ( λ 8892 ) / I ( λ 9267 ) emissivity ratio using different atomic datasets. The observed ratio in HH 202S is 2.36 ± 0.42 , which is depicted as a gray band. Left panel: atom_TZ18_52 radiative dataset with four coll datasets (SRKFB19_52, B15_52, TZ18_52, VVKFHF99_52) at T e = 10 4  K. Right panel: atom_TZ18_52 and coll_TZ18_52 datasets at four electron temperatures.
Atoms 11 00063 g010
Figure 11. [Fe ii] radiative lifetimes for levels within the 3 d 6 4 s , 3 d 7 , and  3 d 5 4 s 2 configurations computed from the A-values in the atom datasets B15_52, DH11_52, TZ18_52, SRKFB19_52, QDZa96, and QDZh96. For each level, the lifetime differences for these datasets with respect to the average value are plotted as a function of the average value.
Figure 11. [Fe ii] radiative lifetimes for levels within the 3 d 6 4 s , 3 d 7 , and  3 d 5 4 s 2 configurations computed from the A-values in the atom datasets B15_52, DH11_52, TZ18_52, SRKFB19_52, QDZa96, and QDZh96. For each level, the lifetime differences for these datasets with respect to the average value are plotted as a function of the average value.
Atoms 11 00063 g011
Figure 12. Comparison of experimental A-values (s 1 ) for E1 transitions arising from the 3 d 6 ( 5 D ) 4 p levels of Fe ii [46] with those from the atom datasets TZ18_173 and SRKFB19_173.
Figure 12. Comparison of experimental A-values (s 1 ) for E1 transitions arising from the 3 d 6 ( 5 D ) 4 p levels of Fe ii [46] with those from the atom datasets TZ18_173 and SRKFB19_173.
Atoms 11 00063 g012
Figure 13. Comparison of observed g f -values for Fe ii (MB09) [47] with those calculated from atom_TZ18_173, atom_SRKFB19_173, and Ref. [43] (DH14).
Figure 13. Comparison of observed g f -values for Fe ii (MB09) [47] with those calculated from atom_TZ18_173, atom_SRKFB19_173, and Ref. [43] (DH14).
Atoms 11 00063 g013
Figure 14. Comparison of Fe ii experimental branching fractions ( B F expt ) [48] with theoretical values ( B F th ) derived from atom_TZ18_173 and atom_SRKFB19_173.
Figure 14. Comparison of Fe ii experimental branching fractions ( B F expt ) [48] with theoretical values ( B F th ) derived from atom_TZ18_173 and atom_SRKFB19_173.
Atoms 11 00063 g014
Figure 15. [Fe ii] ECS logarithmic differences at T = 10 4  K for each transition in the coll datasets B15_52, SRKFB19_52, and TZ18_52 relative to the average value.
Figure 15. [Fe ii] ECS logarithmic differences at T = 10 4  K for each transition in the coll datasets B15_52, SRKFB19_52, and TZ18_52 relative to the average value.
Atoms 11 00063 g015
Figure 16. ECS Υ ( J J ) for the transitions 3 d 6 ( 5 D ) 4 s a 6 D J a 6 D J within the ground term of [Fe ii]. Blue circles: Ref. [36]. Yellow crosses: coll_B15_52. Green squares: coll_SRKFB19_52. Red triangles: coll_TZ18_52.
Figure 16. ECS Υ ( J J ) for the transitions 3 d 6 ( 5 D ) 4 s a 6 D J a 6 D J within the ground term of [Fe ii]. Blue circles: Ref. [36]. Yellow crosses: coll_B15_52. Green squares: coll_SRKFB19_52. Red triangles: coll_TZ18_52.
Atoms 11 00063 g016
Table 1. A-value comparison for the 3 d 5 4 s 5 S 2 3 d 6 5 D J transitions in [Fe iii]. A-values in J00 were calculated with the semi-relativistic hfr code using empirical Slater parameters [52] and included in the atom_Q96_J00 dataset (see Section 3.1). The methods used to compute the radiative data in atom_BBQ10_34 and atom_BB14_34 are described below. A-values in atom_FBQh16 and atom_FBQa16 were, respectively, computed with hfr and autostructure with extensive CI [45].
Table 1. A-value comparison for the 3 d 5 4 s 5 S 2 3 d 6 5 D J transitions in [Fe iii]. A-values in J00 were calculated with the semi-relativistic hfr code using empirical Slater parameters [52] and included in the atom_Q96_J00 dataset (see Section 3.1). The methods used to compute the radiative data in atom_BBQ10_34 and atom_BB14_34 are described below. A-values in atom_FBQh16 and atom_FBQa16 were, respectively, computed with hfr and autostructure with extensive CI [45].
J λ air (Å)A-Value (s 1 )
J00BBQ10_34BB14_34FBQh16FBQa16
42438.2822.225.738.432.031.6
32464.4816.018.728.023.323.1
22483.0111.012.819.115.915.8
12495.016.47.4311.19.299.22
02500.932.02.443.66
Table 2. Comparison of observed [Fe iii] line intensity ratios in HH 202S (Obs) for transitions sharing a common upper level with those estimated from the A-values in the BBQ10_34, BB14_34, DH09_34, and Q96_34 atom files. The error of the least significant figure of the observed intensity ratio is indicated in brackets.
Table 2. Comparison of observed [Fe iii] line intensity ratios in HH 202S (Obs) for transitions sharing a common upper level with those estimated from the A-values in the BBQ10_34, BB14_34, DH09_34, and Q96_34 atom files. The error of the least significant figure of the observed intensity ratio is indicated in brackets.
Line 1Line 2ObsBBQ10_34BB14_34DH09_34Q96_34
UpperLower λ air (Å)Lower λ air (Å)Line Intensity Ratio
3 d 6 3 D 3 3 d 6 5 D 4 3239.79 3 d 6 5 D 3 3286.244(1)3.583.593.283.60
3 d 6 3 D 3 3 d 6 5 D 3 3286.24 3 d 6 a 3 P 2 8728.842.0(5)4.883.123.803.29
3 d 6 3 D 1 3 d 6 5 D 1 3355.50 3 d 6 5 D 0 3366.221.6(8)1.131.121.161.14
3 d 6 3 D 1 3 d 6 5 D 0 3366.22 3 d 6 a 3 P 1 9959.855(4)9.075.426.375.54
3 d 6 3 D 2 3 d 6 5 D 2 3334.95 3 d 6 5 D 1 3356.591.2(2)1.141.151.201.18
3 d 6 3 D 2 3 d 6 5 D 1 3356.59 3 d 6 a 3 P 2 8838.144.1(8)6.014.284.394.19
3 d 5 4 s 7 S 3 3 d 6 5 D 4 3322.47 3 d 6 5 D 3 3371.351.4(3)1.381.37
3 d 5 4 s 7 S 3 3 d 6 5 D 3 3371.35 3 d 6 5 D 2 3406.112.2(5)2.332.31
3 d 6 3 G 3 3 d 6 5 D 3 4046.49 3 d 6 5 D 2 4096.682.3(8)2.272.452.562.53
3 d 6 3 G 4 3 d 6 5 D 4 4008.34 3 d 6 5 D 3 4079.693.8(7)3.153.633.623.92
3 d 6 a 3 F 2 3 d 6 5 D 3 4667.11 3 d 6 5 D 2 4734.000.29(4)0.280.280.290.28
3 d 6 a 3 F 2 3 d 6 5 D 2 4734.00 3 d 6 5 D 1 4777.702.0(3)2.062.072.092.08
3 d 6 a 3 F 3 3 d 6 5 D 4 4607.12 3 d 6 5 D 3 4701.640.19(3)0.180.180.190.17
3 d 6 a 3 F 3 3 d 6 5 D 3 4701.64 3 d 6 5 D 2 4769.532.8(3)2.892.902.932.94
3 d 6 a 3 F 4 3 d 6 5 D 4 4658.17 3 d 6 5 D 3 4754.815.3(6)5.285.405.325.49
3 d 6 a 3 P 1 3 d 6 5 D 2 5011.41 3 d 6 5 D 0 5084.855.9(9)5.726.125.975.95
3 d 6 3 H 4 3 d 6 5 D 4 4881.07 3 d 6 5 D 3 4987.295.4(7)4.995.335.505.76
3 d 6 a 3 P 2 3 d 6 5 D 3 5270.57 3 d 6 5 D 1 5412.0611(2)10.211.411.411.0
Table 3. Comparison of observed [Fe iii] line intensity ratios in HH 204 (Obs) for transitions sharing a common upper level with those estimated from A-values in the BBQ10_34, BB14_34, DH09_34, and Q96_34 atom files. The error of the least significant figure of the observed intensity ratio is indicated in brackets.
Table 3. Comparison of observed [Fe iii] line intensity ratios in HH 204 (Obs) for transitions sharing a common upper level with those estimated from A-values in the BBQ10_34, BB14_34, DH09_34, and Q96_34 atom files. The error of the least significant figure of the observed intensity ratio is indicated in brackets.
Line 1Line 2ObsBBQ10_34BB14_34DH09_34Q96_34
UpperLower λ air (Å)Lower λ air (Å)Line Intensity Ratio
3 d 6 3 D 3 3 d 6 5 D 4 3239.79 3 d 6 5 D 3 3286.243.6(9)3.583.593.283.60
3 d 6 3 D 3 3 d 6 5 D 3 3286.24 3 d 6 5 D 2 3319.271.0(4)1.351.361.471.41
3 d 6 3 D 3 3 d 6 5 D 2 3319.27 3 d 6 a 3 P 2 8728.843.1(8)3.622.302.592.35
3 d 6 3 D 1 3 d 6 5 D 1 3355.50 3 d 6 5 D 0 3366.221.5(3)1.131.121.161.14
3 d 6 3 D 2 3 d 6 5 D 2 3334.95 3 d 6 5 D 1 3356.591.2(3)1.141.151.201.18
3 d 6 3 D 2 3 d 6 5 D 1 3356.59 3 d 6 a 3 P 2 8838.145.3(7)6.014.284.394.19
3 d 6 1 I 6 3 d 6 3 H 6 9701.87 3 d 6 3 H 5 9942.381.4(2)1.581.541.561.55
3 d 5 4 s 7 S 3 3 d 6 5 D 4 3322.47 3 d 6 5 D 3 3371.351.5(2)1.381.37
3 d 5 4 s 7 S 3 3 d 6 5 D 3 3371.35 3 d 6 5 D 2 3406.111.7(2)2.332.31
3 d 6 3 G 3 3 d 6 5 D 3 4046.49 3 d 6 5 D 2 4096.683.4(9)2.272.452.562.53
3 d 6 3 G 4 3 d 6 5 D 4 4008.34 3 d 6 5 D 3 4079.694.4(4)3.153.633.623.92
3 d 6 a 3 F 2 3 d 6 5 D 3 4667.11 3 d 6 5 D 2 4734.000.29(1)0.280.280.290.28
3 d 6 a 3 F 2 3 d 6 5 D 2 4734.00 3 d 6 5 D 1 4777.702.1(1)2.062.072.092.08
3 d 6 a 3 F 3 3 d 6 5 D 4 4607.12 3 d 6 5 D 3 4701.640.18(1)0.180.180.190.17
3 d 6 a 3 F 3 3 d 6 5 D 3 4701.64 3 d 6 5 D 2 4769.532.9(1)2.892.902.932.94
3 d 6 a 3 F 4 3 d 6 5 D 4 4658.17 3 d 6 5 D 3 4754.815.3(2)5.285.405.325.49
3 d 6 a 3 P 1 3 d 6 5 D 2 5011.41 3 d 6 5 D 0 5084.855.9(4)5.726.125.975.95
3 d 6 3 H 4 3 d 6 5 D 4 4881.07 3 d 6 5 D 3 4987.296.1(2)4.995.335.505.76
3 d 6 a 3 P 2 3 d 6 5 D 3 5270.57 3 d 6 5 D 1 5412.0610.8(5)10.211.411.411.0
Table 4. Temperature and density fits of the [Fe iii] observed spectra of HH 202S and HH 204 with the radiative (atom) and collisional (coll) datasets available in PyNeb 1.1.16.
Table 4. Temperature and density fits of the [Fe iii] observed spectra of HH 202S and HH 204 with the radiative (atom) and collisional (coll) datasets available in PyNeb 1.1.16.
DatasetsHH 202SHH 204
AtomColl T e ( 10 3 K ) n e ( 10 4 cm 3 ) χ 2 T e ( 10 3 K ) n e ( 10 4 cm 3 ) χ 2
BBQ10BBQ107.991.9212.09.771.5366.9
NP96BBQ108.373.7913.910.52.2457.2
Q96_J00BBQ107.801.9614.49.601.3658.1
BBQ10BB145.2613.114.95.9511.388.7
NP96BB145.5416.815.36.1510.560.2
Q96_J00BB145.268.5415.85.955.5160.4
BBQ10Z967.342.521.748.531.7016.9
NP96Z967.447.0915.29.544.2069.9
Q96_J00Z967.352.5512.08.791.6837.7
Table 5. Temperature and density fits of the [Fe iii] observed spectra in HH 202S and HH 204 with the improved or new radiative (atom) and collisional (coll) datasets.
Table 5. Temperature and density fits of the [Fe iii] observed spectra in HH 202S and HH 204 with the improved or new radiative (atom) and collisional (coll) datasets.
DatasetsHH 202SHH 204
AtomColl T e ( 10 3 K ) n e ( 10 4 cm 3 ) χ 2 T e ( 10 3 K ) n e ( 10 4 cm 3 ) χ 2
BBQ10_34BBQ10_347.681.946.319.491.5147.6
BB14_34BBQ10_348.241.803.849.911.4428.1
DH09_34BBQ10_347.941.312.249.250.8221.8
Q96_34BBQ10_347.971.663.549.391.1935.6
BBQ10_34BB14_347.832.175.579.012.1650.0
BB14_34BB14_348.472.144.569.031.4933.2
DH09_34BB14_348.181.783.239.001.169.95
Q96_34BB14_348.121.913.978.971.4911.9
BBQ10_34Z96_347.342.511.748.531.7116.9
BB14_34Z96_347.942.322.548.231.6128.5
DH09_34Z96_347.711.970.688.341.304.36
Q96_34Z96_347.682.241.008.111.734.62
BB14_144BB14_1448.472.124.579.031.4933.2
BB14_144Z96_1447.942.332.548.241.6128.5
Table 6. Comparison of NIST energies (cm 1 ) for levels in Fe ii displaying a total orbital angular momentum quantum number L = 5 with those listed in different atom datasets (see Section 4.2). Questionable level assignments are found in atom_SRKFB19_225.
Table 6. Comparison of NIST energies (cm 1 ) for levels in Fe ii displaying a total orbital angular momentum quantum number L = 5 with those listed in different atom datasets (see Section 4.2). Questionable level assignments are found in atom_SRKFB19_225.
LevelNISTSRKFB19_225TZ18_225DH11_52B15_52QDZh96
3 d 7 a 2 H 11 / 2 20,340.2520,340.2520,327.9120,340.3020,340.3020,340
3 d 7 a 2 H 9 / 2 20,805.7620,805.7620,824.9920,805.7720,805.7720,806
3 d 6 ( 3 H ) 4 s a 4 H 13 / 2 21,251.5821,251.5821,283.5221,251.6121,251.6121,252
3 d 6 ( 3 H ) 4 s a 4 H 11 / 2 21,430.3626,170.1721,433.7021,430.3621,430.3621,430
3 d 6 ( 3 H ) 4 s a 4 H 9 / 2 21,581.6226,352.7721,560.1721,581.6421,581.6421,582
3 d 6 ( 3 H ) 4 s a 4 H 7 / 2 21,711.9021,711.9021,671.6321,711.9221,711.9221,712
3 d 6 ( 3 H ) 4 s b 2 H 11 / 2 26,170.1821,430.3626,186.4826,170.1826,170.1826,170
3 d 6 ( 3 H ) 4 s b 2 H 9 / 2 26,352.7721,581.6126,331.5026,352.7726,352.7726,353
3 d 6 ( 3 H ) 4 p z 4 H 13 / 2 o 60,837.5660,837.5560,832.43
3 d 6 ( 3 H ) 4 p z 4 H 11 / 2 o 60,887.61 60,880.58
3 d 6 ( 3 H ) 4 p z 4 H 9 / 2 o 60,989.44 61,012.21
3 d 6 ( 3 H ) 4 p z 4 H 7 / 2 o 61,156.8361,156.8261,178.84
3 d 6 ( 3 H ) 4 p z 2 H 11 / 2 o 65,363.6165,363.6065,333.65
3 d 6 ( 3 H ) 4 p z 2 H 9 / 2 o 65,556.2765,556.2665,551.82
3 d 6 ( 3 G ) 4 p y 4 H 13 / 2 o 66,411.7166,411.7066,306.76
3 d 6 ( 3 G ) 4 p y 4 H 11 / 2 o 66,463.5466,463.5466,602.84
3 d 6 ( 3 G ) 4 p y 4 H 9 / 2 o 66,589.0466,589.0366,733.10
3 d 6 ( 3 G ) 4 p y 4 H 7 / 2 o 66,672.3466,672.3266,705.11
3 d 6 ( 3 G ) 4 p y 2 H 11 / 2 o 67,516.3367,516.3267,942.69
3 d 6 ( 3 G ) 4 p y 2 H 11 / 2 o 67,709.96
3 d 6 ( 3 G ) 4 p y 2 H 9 / 2 o 68,000.7968,000.7767,609.99
3 d 6 ( 3 G ) 4 p y 2 H 9 / 2 o 68,201.16
3 d 6 ( 1 G 2 ) 4 p x 2 H 9 / 2 o 72,130.3872,130.3672,407.94
3 d 6 ( 1 G 2 ) 4 p x 2 H 11 / 2 o 72,261.7472,261.7371,962.96
3 d 6 ( 1 I ) 4 p w 2 H 11 / 2 o 73,603.5473,603.5373,435.97
3 d 6 ( 1 I ) 4 p w 2 H 9 / 2 o 73,751.2873,751.2773,868.28
Table 7. A-values for transitions arising from the 3 d 6 4 s a 4 H J and 3 d 6 4 s b 2 H J levels ( J = 9 / 2 and 11 / 2 ) upper levels of [Fe ii] listed in different atom datasets, which show the large discrepancies in atom_SRKFB19_225.
Table 7. A-values for transitions arising from the 3 d 6 4 s a 4 H J and 3 d 6 4 s b 2 H J levels ( J = 9 / 2 and 11 / 2 ) upper levels of [Fe ii] listed in different atom datasets, which show the large discrepancies in atom_SRKFB19_225.
λ air (Å)A-Value (s 1 )
SRKFB19_225TZ18_225DH11_52B15_52QDZa96QDZh96
4114.48 1.52 × 10 2 1.44 × 10 1 1.61 × 10 1 1.05 × 10 1 7.80 × 10 2 1.03 × 10 1
4178.96 7.45 × 10 4 3.26 × 10 2 2.02 × 10 2 1.89 × 10 2 9.95 × 10 3 1.56 × 10 2
4211.11 1.81 × 10 2 5.40 × 10 2 5.93 × 10 2 3.15 × 10 2 3.37 × 10 2 4.44 × 10 2
4251.45 6.12 × 10 3 2.86 × 10 2 1.96 × 10 2 1.33 × 10 2 1.28 × 10 2 1.87 × 10 2
5111.64 4.68 × 10 1 1.01 × 10 1 1.07 × 10 1 1.19 × 10 1 1.20 × 10 1 1.31 × 10 1
5220.08 8.90 × 10 2 1.09 × 10 1 1.11 × 10 1 1.32 × 10 1 1.34 × 10 1 1.44 × 10 1
5261.63 1.62 × 10 + 0 3.25 × 10 1 3.39 × 10 1 3.95 × 10 1 4.01 × 10 1 4.29 × 10 1
5333.66 8.85 × 10 1 2.68 × 10 1 2.78 × 10 1 3.23 × 10 1 3.32 × 10 1 3.51 × 10 1
Table 8. Comparison of [Fe ii] line-intensity ratios for transitions arising from a common upper level in the HH 202S spectrum (Obs) with those estimated from A-values in the following atom datasets: B15_52 (T1); DH11_52 (T2); TZ18_52 (T3); SRKFB19_52 (T4); QDZa96 (T5); QDZh96 (T6). The error of the least significant figure of the observed intensity ratio is indicated in brackets.
Table 8. Comparison of [Fe ii] line-intensity ratios for transitions arising from a common upper level in the HH 202S spectrum (Obs) with those estimated from A-values in the following atom datasets: B15_52 (T1); DH11_52 (T2); TZ18_52 (T3); SRKFB19_52 (T4); QDZa96 (T5); QDZh96 (T6). The error of the least significant figure of the observed intensity ratio is indicated in brackets.
Line 1Line 2ObsT1T2T3T4T5T6
UpperLower λ air (Å)Lower λ air (Å)Line Intensity Ratio
b 2 H 9 / 2 a 4 F 7 / 2 4178.96 a 4 F 5 / 2 4251.451.3(9)1.441.051.160.120.790.85
b 2 H 11 / 2 a 4 F 9 / 2 4114.48 a 4 F 7 / 2 4211.112.4(8)3.422.782.730.862.372.37
a 4 G 7 / 2 a 4 F 5 / 2 4319.62 a 4 F 3 / 2 4372.432.4(8)2.001.951.971.881.981.96
a 4 G 9 / 2 a 4 F 9 / 2 4177.20 a 4 F 7 / 2 4276.840.28(8)0.240.210.230.290.250.24
a 4 G 9 / 2 a 4 F 7 / 2 4276.84 a 4 F 5 / 2 4352.782.1(6)2.162.152.242.072.212.19
a 4 G 11 / 2 a 4 F 9 / 2 4243.97 a 4 F 7 / 2 4346.865(1)4.404.544.674.244.604.59
a 6 S 5 / 2 a 6 D 9 / 2 4287.39 a 6 D 7 / 2 4359.331.3(3)1.371.371.381.351.371.38
a 6 S 5 / 2 a 6 D 7 / 2 4359.33 a 6 D 5 / 2 4413.781.4(3)1.441.431.441.421.421.44
a 6 S 5 / 2 a 6 D 5 / 2 4413.78 a 6 D 3 / 2 4452.101.6(4)1.581.571.581.561.571.58
a 6 S 5 / 2 a 6 D 3 / 2 4452.10 a 6 D 1 / 2 4474.902.2(7)2.072.062.062.052.062.06
b 4 F 3 / 2 a 6 D 3 / 2 4509.61 a 4 F 5 / 2 4950.760.4(2)0.310.400.380.110.280.28
b 4 F 3 / 2 a 4 F 5 / 2 4950.76 a 4 F 3 / 2 5020.240.9(4)0.960.970.980.950.970.97
b 4 F 5 / 2 a 6 D 7 / 2 4432.45 a 6 D 5 / 2 4488.750.6(3)0.360.360.360.360.370.36
b 4 F 5 / 2 a 6 D 5 / 2 4488.75 a 6 D 3 / 2 4528.383(2)3.443.443.433.383.423.45
b 4 F 5 / 2 a 6 D 3 / 2 4528.38 a 4 F 7 / 2 4874.500.3(2)0.230.300.280.090.220.21
b 4 F 5 / 2 a 4 F 7 / 2 4874.50 a 4 F 5 / 2 4973.401.3(5)1.251.241.271.261.261.24
b 4 F 5 / 2 a 4 F 5 / 2 4973.40 a 4 F 3 / 2 5043.531.5(7)2.012.071.991.611.931.97
b 4 F 7 / 2 a 6 D 7 / 2 4457.95 a 6 D 5 / 2 4514.904(2)4.304.354.324.224.304.34
b 4 F 7 / 2 a 6 D 5 / 2 4514.90 a 4 F 9 / 2 4774.730.6(3)0.460.600.560.190.440.42
b 4 F 7 / 2 a 4 F 9 / 2 4774.73 a 4 F 7 / 2 4905.350.6(2)0.580.580.600.720.590.59
b 4 F 9 / 2 a 6 D 9 / 2 4416.27 a 6 D 7 / 2 4492.647(2)7.617.547.707.417.607.74
b 4 F 9 / 2 a 6 D 7 / 2 4492.64 a 4 F 9 / 2 4814.540.15(5)0.140.190.170.070.130.12
b 4 F 9 / 2 a 4 F 9 / 2 4814.54 a 4 F 7 / 2 4947.397(2)7.576.357.000.887.097.17
b 4 F 9 / 2 a 4 F 7 / 2 4947.39 a 4 D 7 / 2 6809.244(2)3.493.223.0770.23.984.16
b 4 P 3 / 2 a 6 D 5 / 2 4728.07 a 4 F 7 / 2 5158.010.9(3)1.231.832.000.571.291.19
b 4 P 3 / 2 a 4 F 7 / 2 5158.01 a 4 F 5 / 2 5268.893(1)1.571.561.571.481.571.56
a 4 H 7 / 2 a 4 F 5 / 2 5296.84 a 4 F 3 / 2 5376.470.3(1)0.320.340.340.130.340.34
a 4 H 9 / 2 a 4 F 7 / 2 5220.08 a 4 F 5 / 2 5333.660.4(1)0.420.410.420.100.410.42
a 4 H 11 / 2 a 4 F 9 / 2 5111.64 a 4 F 7 / 2 5261.630.32(6)0.310.320.320.300.310.31
a 2 D 2 3 / 2 a 6 D 3 / 2 4889.71 a 4 F 5 / 2 5412.683(1)0.000.000.000.000.000.00
a 2 D 2 3 / 2 a 4 F 5 / 2 5412.68 a 4 F 3 / 2 5495.842(1)1.891.911.901.891.891.90
b 4 P 5 / 2 a 4 F 9 / 2 5273.36 a 4 F 7 / 2 5433.153(1)3.322.883.193.273.283.26
a 2 D 2 5 / 2 a 4 F 7 / 2 5527.35 a 4 F 5 / 2 5654.8710(4)9.269.099.199.129.068.98
a 2 G 7 / 2 a 4 F 7 / 2 7172.00 a 4 F 5 / 2 7388.171.4(3)1.351.351.361.321.351.35
a 2 G 9 / 2 a 4 F 9 / 2 7155.17 a 4 F 7 / 2 7452.563.1(5)3.173.223.233.103.193.24
a 4 P 1 / 2 a 4 F 5 / 2 9033.49 a 4 F 3 / 2 9267.550.8(2)0.770.780.780.500.780.78
a 4 P 3 / 2 a 6 D 5 / 2 7686.93 a 4 F 7 / 2 8891.930.27(7)0.260.830.440.400.360.44
a 4 P 3 / 2 a 4 F 7 / 2 8891.93 a 4 F 5 / 2 9226.631.7(3)1.761.791.780.881.781.80
a 4 P 5 / 2 a 6 D 7 / 2 7637.52 a 4 F 7 / 2 9051.950.5(1)0.591.961.080.180.891.15
a 4 P 5 / 2 a 4 F 7 / 2 9051.95 a 4 F 5 / 2 9399.047(2)5.135.325.212.645.465.67
Table 9. Comparison of [Fe ii] line-intensity ratios for transitions arising from a common upper level in the HH 204 spectrum (Obs) with those estimated from A-values in the following atom datasets: B15_52 (T1); DH11_52 (T2); TZ18_52 (T3); SRKFB19_52 (T4); QDZa96 (T5); QDZh96 (T6). The error of the least significant figure of the observed intensity ratio is indicated in brackets.
Table 9. Comparison of [Fe ii] line-intensity ratios for transitions arising from a common upper level in the HH 204 spectrum (Obs) with those estimated from A-values in the following atom datasets: B15_52 (T1); DH11_52 (T2); TZ18_52 (T3); SRKFB19_52 (T4); QDZa96 (T5); QDZh96 (T6). The error of the least significant figure of the observed intensity ratio is indicated in brackets.
Line 1Line 2ObsT1T2T3T4T5T6
UpperLower λ air (Å)Lower λ air (Å)Line Intensity Ratio
a 2 F 7 / 2 a 4 D 7 / 2 5163.96 a 2 G 9 / 2 8715.809(2)7.5010.511.12.816.737.48
b 2 H 9 / 2 a 4 F 7 / 2 4178.96 a 4 F 5 / 2 4251.451.2(2)1.441.051.160.120.790.85
b 2 H 11 / 2 a 4 F 9 / 2 4114.48 a 4 F 7 / 2 4211.112.6(2)3.422.782.730.862.372.37
a 4 G 5 / 2 a 6 D 3 / 2 3968.27 a 4 F 5 / 2 4305.907(1)0.000.000.000.000.000.00
a 4 G 5 / 2 a 4 F 5 / 2 4305.90 a 4 F 3 / 2 4358.370.56(8)0.440.420.440.480.460.45
a 4 G 7 / 2 a 4 F 5 / 2 4319.62 a 4 F 3 / 2 4372.432.0(2)2.001.951.971.881.981.96
a 4 G 9 / 2 a 4 F 9 / 2 4177.20 a 4 F 7 / 2 4276.840.25(2)0.240.210.230.290.250.24
a 4 G 9 / 2 a 4 F 7 / 2 4276.84 a 4 F 5 / 2 4352.782.1(1)2.162.152.242.072.212.19
a 4 G 11 / 2 a 4 F 9 / 2 4243.97 a 4 F 7 / 2 4346.864.9(2)4.404.544.674.244.604.59
a 6 S 5 / 2 a 6 D 9 / 2 4287.39 a 6 D 7 / 2 4359.331.39(6)1.371.371.381.351.371.38
a 6 S 5 / 2 a 6 D 7 / 2 4359.33 a 6 D 5 / 2 4413.781.42(6)1.441.431.441.421.421.44
a 6 S 5 / 2 a 6 D 5 / 2 4413.78 a 6 D 3 / 2 4452.101.55(6)1.581.571.581.571.571.58
a 6 S 5 / 2 a 6 D 3 / 2 4452.10 a 6 D 1 / 2 4474.902.0(1)2.072.062.062.052.062.06
b 4 F 3 / 2 a 6 D 3 / 2 4509.61 a 4 F 5 / 2 4950.760.42(9)0.310.400.380.110.280.28
b 4 F 3 / 2 a 4 F 5 / 2 4950.76 a 4 F 3 / 2 5020.241.4(3)0.960.970.980.950.970.97
b 4 F 5 / 2 a 6 D 7 / 2 4432.45 a 6 D 5 / 2 4488.750.49(7)0.360.360.360.360.370.36
b 4 F 5 / 2 a 6 D 5 / 2 4488.75 a 6 D 3 / 2 4528.383.6(7)3.443.443.433.383.423.45
b 4 F 5 / 2 a 6 D 3 / 2 4528.38 a 4 F 7 / 2 4874.500.25(5)0.230.300.280.090.220.21
b 4 F 5 / 2 a 4 F 7 / 2 4874.50 a 4 F 5 / 2 4973.401.4(1)1.251.241.271.261.261.24
b 4 F 5 / 2 a 4 F 5 / 2 4973.40 a 4 F 3 / 2 5043.531.9(3)2.012.071.991.611.931.97
b 4 F 7 / 2 a 6 D 9 / 2 4382.74 a 6 D 7 / 2 4457.950.19(2)0.200.210.200.200.210.20
b 4 F 7 / 2 a 6 D 7 / 2 4457.95 a 6 D 5 / 2 4514.904.3(4)4.304.354.324.224.304.34
b 4 F 7 / 2 a 6 D 5 / 2 4514.90 a 4 F 9 / 2 4774.730.56(6)0.460.600.560.190.440.42
b 4 F 7 / 2 a 4 F 9 / 2 4774.73 a 4 F 7 / 2 4905.350.57(4)0.580.580.600.720.590.59
b 4 F 9 / 2 a 6 D 9 / 2 4416.27 a 6 D 7 / 2 4492.647.7(5)7.617.547.707.417.607.74
b 4 F 9 / 2 a 6 D 7 / 2 4492.64 a 4 F 9 / 2 4814.540.15(1)0.140.190.170.070.130.12
b 4 F 9 / 2 a 4 F 9 / 2 4814.54 a 4 F 7 / 2 4947.397.8(5)7.576.357.000.887.097.17
b 4 F 9 / 2 a 4 F 7 / 2 4947.39 a 4 D 7 / 2 6809.243.2(4)3.493.223.0770.213.984.16
b 4 P 3 / 2 a 6 D 5 / 2 4728.07 a 4 F 5 / 2 5268.892.3(4)1.932.853.140.842.021.85
a 4 H 7 / 2 a 4 F 5 / 2 5296.84 a 4 F 3 / 2 5376.470.36(5)0.320.340.340.130.340.34
a 4 H 9 / 2 a 4 F 7 / 2 5220.08 a 4 F 5 / 2 5333.660.44(5)0.420.410.420.100.410.42
a 4 H 11 / 2 a 4 F 9 / 2 5111.64 a 4 F 7 / 2 5261.630.34(3)0.310.320.320.300.310.31
a 2 D 2 3 / 2 a 6 D 3 / 2 4889.71 a 4 F 3 / 2 5495.8412(2)0.000.000.000.000.000.00
b 4 P 5 / 2 a 4 F 9 / 2 5273.36 a 4 F 7 / 2 5433.152.7(2)3.322.883.193.273.283.26
b 4 P 5 / 2 a 4 F 7 / 2 5433.15 a 4 D 7 / 2 7764.717(1)7.225.644.8822.637.017.42
a 2 G 7 / 2 a 4 F 9 / 2 6896.17 a 4 F 7 / 2 7172.000.10(1)0.100.090.090.090.100.09
a 2 G 7 / 2 a 4 F 7 / 2 7172.00 a 4 F 5 / 2 7388.171.34(9)1.351.351.361.321.351.35
a 2 G 9 / 2 a 4 F 9 / 2 7155.17 a 4 F 7 / 2 7452.563.2(2)3.173.223.233.103.193.24
a 4 P 1 / 2 a 6 D 3 / 2 7665.28 a 6 D 1 / 2 7733.132.8(7)3.263.313.284.443.263.28
a 4 P 1 / 2 a 6 D 1 / 2 7733.13 a 4 F 5 / 2 9033.490.11(2)0.100.320.170.100.140.17
a 4 P 1 / 2 a 4 F 5 / 2 9033.49 a 4 F 3 / 2 9267.550.75(7)0.770.780.780.500.780.78
a 4 P 3 / 2 a 6 D 5 / 2 7686.93 a 6 D 1 / 2 7874.237(1)7.067.317.0511.24
a 4 P 3 / 2 a 6 D 1 / 2 7874.23 a 4 F 7 / 2 8891.930.03(1)0.040.110.060.04
a 4 P 3 / 2 a 4 F 7 / 2 8891.93 a 4 F 5 / 2 9226.631.7(2)1.761.791.780.881.781.80
a 4 P 5 / 2 a 6 D 3 / 2 7926.88 a 4 F 7 / 2 9051.950.05(1)0.040.140.080.00
a 4 P 5 / 2 a 4 F 7 / 2 9051.95 a 4 F 5 / 2 9399.045.4(6)5.135.325.212.645.465.67
Table 10. Temperature and density fits of the [Fe ii] observed spectra of HH 202S and HH 204 with the radiative (atom) and collisional (coll) datasets currently available in PyNeb 1.1.16.
Table 10. Temperature and density fits of the [Fe ii] observed spectra of HH 202S and HH 204 with the radiative (atom) and collisional (coll) datasets currently available in PyNeb 1.1.16.
DatasetsHH 202SHH 204
AtomColl T e ( 10 3 K ) n e ( 10 4 cm 3 ) χ 2 T e ( 10 3 K ) n e ( 10 4 cm 3 ) χ 2
B15B1513.06.434.4024.63.9290.8
VVKFHF99VVKFHF999.871.527.8012.80.98105.
SRKFB19SRKFB1918.82.6719.6
Table 11. Temperature and density fits of the [Fe ii] observed spectra in HH 202S and HH 204 with the improved and new radiative (atom) and collisional (coll) datasets.
Table 11. Temperature and density fits of the [Fe ii] observed spectra in HH 202S and HH 204 with the improved and new radiative (atom) and collisional (coll) datasets.
DatasetsHH 202SHH 204
AtomColl T e ( 10 3 K ) n e ( 10 4 cm 3 ) χ 2 T e ( 10 3 K ) n e ( 10 4 cm 3 ) χ 2
B15_52B15_5211.65.991.5116.53.9067.0
DH11_52B15_5213.46.004.6022.74.0094.7
TZ18_52B15_5211.45.981.9417.13.6970.1
VVKFHF99_52B15_5211.77.122.1017.43.5275.5
TZ18_52TZ18_527.943.297.9011.22.63154.
B15_52TZ18_528.144.107.3711.53.62145.
DH11_52TZ18_528.923.959.0613.13.02160.
VVKFHF99_52TZ18_528.094.437.3711.83.95147.
B15_52SRKFB19_5210.02.052.2014.50.8985.8
DH11_52SRKFB19_5211.12.584.7717.30.8996.3
TZ18_52SRKFB19_529.912.212.7515.21.1584.7
VVKFHF99_52SRKFB19_5210.12.732.7115.61.2490.4
VVKFHF99_52VVKFHF99_529.921.667.7912.70.89106.
B15_52VVKFHF99_529.851.437.1912.10.8188.3
DH11_52VVKFHF99_5211.21.7512.516.70.93164.
TZ18_52VVKFHF99_529.941.477.7013.00.79101.
TZ18_173TZ18_1737.813.297.9110.41.5140.2
TZ18_173SRKFB19_1739.752.322.8415.51.1964.2
Table 12. Comparison of [Fe ii] experimental lifetimes (in s) with those derived from the revised PyNeb atom datasets. The error of the least significant figure of the observed lifetime is indicated in brackets.
Table 12. Comparison of [Fe ii] experimental lifetimes (in s) with those derived from the revised PyNeb atom datasets. The error of the least significant figure of the observed lifetime is indicated in brackets.
ExperimentTheory
LevelRef. [49]Ref. [50]Ref. [51]B15_52DH11_52TZ18_52QDZa96QDZh96
3 d 5 4 s 2 a 6 S 5 / 2 0.23(3) 0.2410.2330.2250.2620.220
3 d 6 4 s a 4 G 11 / 2 0.75(10)0.9090.9600.9340.7740.706
3 d 6 4 s a 4 G 9 / 2 0.65(2) 0.8490.8620.8760.7550.694
3 d 6 4 s b 2 H 11 / 2 3.8(3) 5.313.714.186.595.20
3 d 6 4 s b 4 D 1 / 2 0.54(3)0.5920.6270.6320.6160.550
3 d 6 4 s b 4 D 7 / 2 0.53(3) 0.5480.5780.5750.5680.501
Table 13. Comparison of experimental A-values (s 1 ) for transitions in [Fe ii] with those listed in the revised PyNebatom datasets. The error of the least significant figure of the observed A-value is indicated in brackets.
Table 13. Comparison of experimental A-values (s 1 ) for transitions in [Fe ii] with those listed in the revised PyNebatom datasets. The error of the least significant figure of the observed A-value is indicated in brackets.
ExperimentTheory
UpperLower λ air  (Å)Ref. [50]Ref. [51]B15_52DH11_52TZ18_52QDZa96QDZh96
a 4 G 9 / 2 a 4 F 9 / 2 4177.200.29(5) 0.160.140.140.180.19
a 4 F 7 / 2 4276.840.83(7) 0.660.650.650.750.82
a 4 F 5 / 2 4352.780.36(4) 0.310.310.290.350.38
a 6 S 5 / 2 a 6 D 9 / 2 4287.391.53(22) 1.501.551.621.371.65
a 6 D 7 / 2 4359.331.19(21) 1.111.151.191.021.22
a 6 D 5 / 2 4413.780.84(13) 0.780.810.840.730.86
a 6 D 3 / 2 4452.100.53(8) 0.500.520.540.470.55
a 6 D 1 / 2 4474.900.26(4) 0.240.250.260.230.27
b 4 D 7 / 2 a 6 D 9 / 2 3175.380.23(3) 0.210.240.210.190.21
a 4 F 9 / 2 3376.200.96(10) 0.880.770.810.850.98
a 4 F 7 / 2 3440.990.23(3) 0.300.270.280.290.33
a 4 P 5 / 2 5551.310.18(4) 0.150.140.140.150.18
a 4 P 3 / 2 5613.270.10(3) 0.0780.0690.0750.0810.093
a 4 G 11 / 2 a 4 F 9 / 2 4243.97 1.05(15)0.850.810.841.021.12
a 4 F 7 / 2 4346.85 0.25(5)0.200.180.180.230.25
Table 14. Comparison of experimental lifetimes (in ns) of the 3 d 6 ( 5 D ) 4 p levels of Fe ii with those derived from A-values listed in the revised PyNeb atom datasets. The error of the least significant figure of the observed lifetime is indicated in brackets.
Table 14. Comparison of experimental lifetimes (in ns) of the 3 d 6 ( 5 D ) 4 p levels of Fe ii with those derived from A-values listed in the revised PyNeb atom datasets. The error of the least significant figure of the observed lifetime is indicated in brackets.
ExperimentTheory
LevelRef. [46]TZ18_173SRKFB19_173
z 6 D 9 / 2 o 3.68(7)3.173.39
z 6 D 7 / 2 o 3.67(9)3.193.43
z 6 D 5 / 2 o 3.69(5)3.203.44
z 6 D 3 / 2 o 3.73(7)3.213.44
z 6 D 1 / 2 o 3.68(11)3.213.45
z 6 F 11 / 2 o 3.20(5)2.582.94
z 6 F 9 / 2 o 3.28(4)2.632.98
z 6 F 7 / 2 o 3.25(6)2.653.09
z 6 F 5 / 2 o 3.30(5)2.663.03
z 6 F 3 / 2 o 3.45(12)2.673.04
z 6 P 7 / 2 o 3.71(4)2.933.53
z 6 P 5 / 2 o 3.75(10)2.913.60
z 6 P 3 / 2 o 3.70(12)2.893.57
z 4 D 7 / 2 o 2.97(4)2.742.10
z 4 D 5 / 2 o 2.90(6)2.762.15
z 4 D 3 / 2 o 2.91(9)2.742.13
z 4 F 9 / 2 o 3.72(10)3.293.14
z 4 F 7 / 2 o 3.59(10)3.143.00
z 4 F 5 / 2 o 3.55(8)3.162.99
z 4 P 5 / 2 o 3.27(6)3.202.74
z 4 P 3 / 2 o 3.23(9)3.212.55
Table 15. Problematic transitions in the g f -value comparison of MB09 [47] with atom_TZ18_173, atom_SRKFB19_173, and DH14 [43] showing large discrepancies or unidentified lines.
Table 15. Problematic transitions in the g f -value comparison of MB09 [47] with atom_TZ18_173, atom_SRKFB19_173, and DH14 [43] showing large discrepancies or unidentified lines.
MB09TZ18_173SRKFB19_173DH14
LowerUpper λ air  (Å) log ( g f )
3 d 6 4 p z 4 F 5 / 2 o 3 d 5 4 s 2 c 4 D 3 / 2 6508.13−3.45−3.93−2.56−7.90
3 d 6 4 s c 4 F 9 / 2 3 d 6 4 p x 4 G 9 / 2 o 6433.81−2.37−3.82−2.91−4.41
3 d 6 4 p z 4 F 7 / 2 o 3 d 5 4 s 2 c 4 D 5 / 2 6371.13−3.13−3.54−2.48−4.19
4635.32−1.42 −1.75
4625.89−2.35 −2.41
4549.19−1.62 −1.86
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mendoza, C.; Méndez-Delgado, J.E.; Bautista, M.; García-Rojas, J.; Morisset, C. Atomic Data Assessment with PyNeb: Radiative and Electron Impact Excitation Rates for [Fe ii] and [Fe iii]. Atoms 2023, 11, 63. https://doi.org/10.3390/atoms11040063

AMA Style

Mendoza C, Méndez-Delgado JE, Bautista M, García-Rojas J, Morisset C. Atomic Data Assessment with PyNeb: Radiative and Electron Impact Excitation Rates for [Fe ii] and [Fe iii]. Atoms. 2023; 11(4):63. https://doi.org/10.3390/atoms11040063

Chicago/Turabian Style

Mendoza, Claudio, José E. Méndez-Delgado, Manuel Bautista, Jorge García-Rojas, and Christophe Morisset. 2023. "Atomic Data Assessment with PyNeb: Radiative and Electron Impact Excitation Rates for [Fe ii] and [Fe iii]" Atoms 11, no. 4: 63. https://doi.org/10.3390/atoms11040063

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop