Benchmarking Atomic Data for Astrophysics: Be-like Ions between B II and Ne VII

Benchmarking Atomic Data for Astrophysics: Be-like Ions between B II and Ne VII


Large-scale self-consistent multiconfiguration Dirac–Hartree–Fock and relativistic configuration interaction calculations are reported for the levels in Be-like ions from B ii to Ne vii. Effects from electron correlation are taken into account by means of large expansions in terms of a basis of configuration state functions, and a complete and accurate data set of excitation energies, lifetimes, wavelengths, and electric dipole, magnetic dipole, electric quadrupole, and magnetic quadrupole line strengths, transition rates, and oscillator strengths for these levels is provided for each ion. Comparisons are made with available experimental and theoretical results. The uncertainty of excitation energies is assessed to be 0.01% on average, which makes it possible to find and rule out misidentifications and aid new line identifications involving high-lying levels in astrophysical spectra. The complete data set is also useful for modeling and diagnosing astrophysical plasmas.

atomic data - atomic processes



College of Science, Hunan University of Technology, Zhuzhou 412000, China \affiliationGroup for Materials Science and Applied Mathematics, Malmö University, SE-20506, Malmö, Sweden;
\affiliationSpace Science Laboratory, University of California, Berkeley, CA 94720, USA

1 Introduction

A wealth of astrophysical spectra has been obtained from different missions, such as Chandra, X-ray Multi-Mirror Mission (XMM-Newton), Hinode, and Hard X-ray Modulation Telescope (HXMT). The analysis of these expensively acquired spectra is often not limited by the capabilities of the spectrometers themselves, but by the lack of accurate reference data in the atomic line database (Bigot & Thévenin, 2006; Kallman & Palmeri, 2007; Ruffoni et al., 2013). To solve this dilemma, we have computed highly accurate energy and radiative transition data for L-shell atomic ions (Jönsson et al., 2013; Ekman et al., 2014; Wang et al., 2014, 2015, 2016a, 2016b, 2017c; Radžiūtė et al., 2015; Si et al., 2016). This work reports accurate atomic data resulting from our systematic calculations for beryllium-like ions from B ii to Ne vii.

These ions have been observed in spectra from different astrophysical objects, such as the sun (Curdt et al., 1997, 2001, 2004; Feldman et al., 1997; Del Zanna & Andretta, 2011; Ko et al., 2002; Parenti et al., 2005; Tian et al., 2009; Thomas & Neupert, 1994), white dwarfs (Raassen et al., 2002; Werner et al., 2004) , as well as nebular regions in RR Tel (Penston & Lago, 1983). The atomic lines do not only tell what elements are present in astrophysical objects and what are their relative abundances, but also reveal the physical conditions in the plasmas, such as density, temperature and radiation fields, along with important excitation mechanisms and photo-processes. For example, the O v lines identified by the Extreme Ultraviolet (EUV) Imaging Spectrometer of the Hinode are density sensitive and can be used for density diagnostics of the solar corona (Landi & Young, 2009). C iii lines have been observed in NGC 2440 and provided carbon abundance information for this planetary nebula (Rubin et al., 2004). Due to the high cosmic abundance from these light elements, the lines often appear strong, increasing their value as probes of physical conditions.

Atomic data for low-lying levels and with and (the complexes) of the ions from B ii to Ne vii have previously been obtained from different calculations (Safronova et al., 1996, 1997, 1999; Tachiev & Froese Fischer, 1999; Komasa et al., 2002; Gálvez et al., 2003; Froese Fischer & Tachiev, 2004; Gu, 2005a; Ho et al., 2006; Cheng et al., 2008; Marques et al., 2012; Verdebout et al., 2014). For many objects, however, emission lines from higher-excitation lines are observed. In the spectrum of the active Seyfert Galaxy NGC 1068 (Brinkman) and the supergiant Orionis  (Raassen & Pollock, 2013), lines from levels with in the light elements such as C, N, O are identified.

There is thus a need for corresponding atomic data for the higher-lying levels to be provided. Our present study is improving existing data, and extending the included levels to as is described below. The new data will extend the probes of physical conditions in different astrophysical objects, as well as of analyzing new observations from different space missions and laboratory experiments (Träbert et al., 2014a, b).

The energies for the () states and transition rates for the transition were calculated by Zhang & Gou (2005) for the beryllium isoelectronic sequence ( = 4–10) using a configuration interaction method. Excitation energies of the and Rydberg states ( = 3–60) in B ii were semi-empirically determined by Sakho et al. (2013) in the framework of the screening-constant-by-unit-nuclear-charge method.  Savukov & Markhotok (2006) performed many-body perturbation theory (MBPT) calculations for a large number of energy levels in Be-like Ne vii. The AUTOSTRUCTURE calculations of Fernández-Menchero et al. (2014) were performed for the levels of Be-like ions with . Data for the states of N iv were provided by  Aggarwal et al. (2016) using the GRASP0 code (Grant et al., 1980), and by Fernández-Menchero et al. (2017) using the B-spline box-based close-coupling method (Zatsarinny & Froese Fischer, 2009).

Among the different studies involving the levels the one by  Savukov & Markhotok (2006) for Ne vii provides energy data of high accuracy. In comparison, although complete sets of data including transition rates are provided, the studies by Fernández-Menchero et al. (2014), for the ions with , and by  Aggarwal et al. (2016) and  Fernández-Menchero et al. (2017), for N iv, are quite inaccurate because of limited configuration interaction effects included in the calculations. For example, the excitation energies of  Fernández-Menchero et al. (2014),  Aggarwal et al. (2016), and  Fernández-Menchero et al. (2017) for N iv depart from the experimental energies compiled in the Atomic Spectra Database (ASD) of the National Institute of Standards and Technology (NIST) (Kramida et al., 2016) by up to 5.8%, 5.7%, and 2.4%, respectively. These uncertainties are too large for identification and deblending of new observations from different space missions.

To fill the gaps of lacking accurate atomic data for high-lying states of the ions from B ii to Ne vii, in particular accurate transition rates, we performed large-scale multiconfiguration Dirac–Hartree–Fock (MCDHF) and relativistic configuration interaction (RCI) calculations using the latest version of the GRASP2K code (Jönsson et al., 2013). The calculations provide a consistent and highly accurate data set of energy structure and radiative transition parameters involving also the levels for these low charged Be-like ions. This work is a continuation of our previous calculations (Wang et al., 2015), in which the corresponding accurate results are reported for the ions in the range of nuclear charges . The excitation energies and lifetimes for the levels of the (with ) and (with and ) configurations, which are below the levels of the configuration, and electric dipole (E1), magnetic dipole (M1), electric quadrupole (E2), and magnetic quadrupole (M2) transition rates among these states are calculated for each ion. To assess the accuracy of the present MCDHF/RCI data, MBPT calculations for Ne vii are carried out using the Flexible Atomic Code (FAC) (Gu, 2008). The two sets of results are in excellent agreement. Compared with previous calculations of Be-like ions, the present work results in a significant extension of accurate energy and transition data for higher-lying states of the configurations, which will greatly aid the analysis of new spectra from astrophysical sources, as well as improving the assessment of blending for diagnostic lines of interest.

2 Theory and Calculations

2.1 Mcdhf

The MCDHF method has been described in detail by Grant (2007) and Froese Fischer et al. (2016). The method is also outlined in our previous calculations (Jönsson et al., 2013; Ekman et al., 2014; Wang et al., 2016a, b, 2017c, 2017a, 2017b; Radžiūtė et al., 2015; Si et al., 2016; Chen & Wang, 2017; Chen et al., 2017). Based on the active space (AS) approach (Olsen et al., 1988; Sturesson et al., 2007) for the generation of the configuration state function (CSF) expansions, separate calculations are done for the even and odd parity states. For the even parity states, the CSF expansions are obtained by allowing single and double (SD) excitations from the multi-reference (MR) configurations , , , , , , , , , , , , , , , , , , , and to an AS of orbitals. For the odd parity states, the CSF expansions are obtained by allowing SD excitations from the MR configurations , , , , , , , , , , , , , , , , , , and to an AS of orbitals. In the first step of the calculations, the AS is


Then, the AS is increased in the following way:




SD excitations from the subshell to active sets with principal quantum numbers are allowed, whereas at most one excitation from the subshell is allowed to active sets with principal quantum numbers . The final model using the AS active set contains 777 325 even and 800 410 odd parity CSFs. The Breit interaction and leading QED effects (vacuum polarization and self-energy) are included in subsequent RCI calculations. All calculations were performed using the GRASP2K code (Jönsson et al., 2007, 2013).

2.2 Mbpt

The MBPT method is explained in Lindgren (1974); Safronova et al. (1996); Vilkas et al. (1999); Gu (2005b, 2007), and it has been implemented in the FAC package (Gu, 2008). The key feature of the MBPT method is the partitioning of the Hilbert space of the system into two subspaces, the model space and the orthogonal space . The configuration interaction effects in the space are exactly considered, while the interaction between the space and is taken into account with the second-order perturbation method. For the MBPT calculation, the model space M contains the even and odd multi-reference configurations of the MCDHF/RCI method, while the space O contains all the possible configurations that are generated by SD virtual excitations of the O space. For single and double excitations, the maximum n values are 125 and 65, respectively, with a maximum value of 25. Just as for the MCDHF/RCI calculations, QED effects are also included.

3 Evaluation of Data

3.1 Energy Levels

The excitation energies for the lowest 138 states of the and (with , with and ) configurations with from our MCDHF/RCI calculations are listed in Table 1, along with the results for the lowest 100 states of these configurations with . All these levels are below the levels of the configurations. In relativistic calculations the wave functions for the states are given as expansions over -coupled CSFs. This labeling system is rarely used in databases or by experimentalists. To overcome this labeling problem, the method developed  Gaigalas et al. (2004, 2017) is used to transform wave functions from the -coupling to the -coupling scheme. In Table 1, for each level numbered by a key , the designation, the total angular momentum and parity , and the radiative lifetime estimated from theoretical transition rates are also included.

Among the six Be-like ions considered here, both experimental and theoretical energy data of the levels up to the configurations are relatively complete for Ne vii. In what follows we will firstly access the accuracy of the MCDHF/RCI excitation energies for Ne vii, by comparing available results. In Table 2, we present the MCDHF/RCI excitation energies of the 138 levels for Ne vii as a function of the increasing active set (AS). For comparison, the compiled values from the NIST ASD (Kramida et al., 2016) are given as well. The mean differences between the MCDHF/RCI and NIST excitation energies and the corresponding standard deviations are cm, cm, cm, , and cm for the calculations based on , , , , and , respectively. This implies that the MCDHF/RCI calculations are well converged with respect to an increasing size of the AS and that the accuracy cannot be further improved by extending the orbital set. Remaining energy differences are due to higher order correlation effects not captured within the framework of SD excitation from the MR. For AS the standard deviation given above corresponds to an average relative difference of .

In Table 3, computed excitation energies of Ne vii from the present MCDHF/RCI and MBPT calculations and from the calculations by Savukov & Markhotok (2006) [hereafter referred to as MBPT2],
 Tachiev & Froese Fischer (1999); Froese Fischer & Tachiev (2004) [MCHF/BP],  Safronova et al. (1996, 1997) [MBPT3], and by  Fernández-Menchero et al. (2014) [AUTOSTRUCTURE] are compared with the NIST compiled values. The present MBPT calculations used the same method and code as in Wang et al. (2015) and  Gu (2005a). Due to these similarities, the latter results are not shown in Table 3.

Compared with the MBPT3 and MCHF/BP calculations, the present MCDHF/RCI and MBPT calculations, as well as the MBPT2 calculations, report data for higher states. Moreover, the latter three calculations show better agreement with the NIST experimental values. The average differences with the standard deviations between the computed excitation energies and the NIST values for the levels are cm for MCDHF/RCI, cm for MBPT, and cm for MBPT2. The corresponding values for MBPT3 and MCHF/BP are cm and cm, respectively.

Inspecting all the results of the MCDHF/RCI, MBPT, and MBPT2 calculations for Ne vii more carefully, one can see that the accuracies of three different calculations, for which the relative average energy differences with the NIST values are (MCDHF/RCI), (MBPT), and (MBPT2), respectively, are generally at the same level. The largest deviations with the NIST values are -776 cm for of MCDHF/RCI, -1111 cm for of MBPT, and 2102 cm of MBPT2. In comparison, although providing complete sets of energy data along with collision strengths, the AUTOSTRUCTURE calculations (Fernández-Menchero et al., 2014), show a much larger average difference with the NIST values, cm, which is far greater than the requirement of analyses of new astrophysical spectra.

Experimental and theoretical energy data are available for many levels along the isoelectronic sequence from B ii to Ne vii. To further assess the accuracy of our computed energy values, a comparison with the NIST experimental values as well as with values from the MBPT3 and MCHF/BP calculations, is given in Table 4. From the table it is clear that the present MCDHF/RCI calculations provide more theoretical data compared with the two previous calculations. The differences with the NIST experimental values are also much smaller for our calculations than those from the MBPT3 and MCHF/BP calculations. More specifically, excluding the levels with , with , and with , for which the differences of theory from experiment are greater than 1000 cm (they will be further discussed below), the average differences with standard deviations, , with the NIST values are cm for MCDHF/RCI, cm for MBPT3, and cm for MCHF/BP. It should be noted that the differences of the MBPT3 and NIST values decrease with increasing for the same level and that the large uncertainties for this method mainly come from the low-charge ions.

Among the levels at , at , and at , where the differences of the MCDHF/RCI and NIST values are greater than 1000 cm, the largest difference is -10451.5 cm for the level at . The NIST energy 498045.5 cm for this level is in good agreement with our MCDHF/RCI value 498010 cm for the level . From Table 1 we can see that the weights for the first two eigenvector components of the level are 70 % () and 27 % (), whereas those of the level are 72 % () and 27 % (). The two levels with relatively pure coupling are mixed with each other, but the mixing is not strong. This implies that the NIST value 498045.5 cm should be designated as the level. Apart from at , we cannot find obvious alternate designations for the other three levels in the present MCDHF/RCI calculations. As an example, the energy differences for the and levels as a function of are shown in Figure 1. Two anomalies exist for the (the difference is 1030.5 cm) and (1989.4 cm ) levels at . Because the same scheme is used for each ion in the present MCDHF/RCI calculations, ensuring that the accuracy of our values for the same level is systematic and consistent, the large differences for these states are caused by relatively large uncertainty of the NIST values.

Lastly we compare the MCDHF/RCI excitation energies for the levels from B ii to Ne vii, which are not computed by Tachiev & Froese Fischer (1999, MCHF/BP); Froese Fischer & Tachiev (2004, MCHF/BP) nor by  Safronova et al. (1996, 1997, MBPT3 ), with the experimental values from the NIST ASD. Removing 3 obvious outliers, = 10 395.6, 2 659, and 6 258 cm for the levels with , and and with , respectively, we find that for the remaining levels in Table 1 cm, which is highly satisfactory considering the range in excitation energy (exceeding 1 500 000 cm).

3.2 Transition rates

Table 5 lists transition rates of the AS active set for the E1, M1, E2, and M2 transitions connecting all the levels included in Table 1. To reduce the amount of data, only values greater than times the sum of values for all transitions from the upper level, i.e., radiative branching ratio (BR) greater than . Also included in this table are wavelengths , line strengths , weighted oscillator strengths , and BR. The general rule is that transition data for electric multipoles (E1,E2) calculated in the Babushkin (length) gauge are preferable over transition data calculated in the Coulomb (velocity) gauge. Whereas generally true, recent work (Wang & Jönsson, 2017; Pehlivan Rhodin et al., 2017) has shown that computed transition data from high states in neutral or near-neutral systems are more stable and show better convergence properties with respect to the increasing active set of orbitals in the Coulomb gauge. For this reason, in Table 5, transition data from the states in B II and C III are computed in the Coulomb gauge. Transitions from all other states are computed in the Babushkin gauge.

Since line strengths for the transitions among the lowest 20 levels up to the configuration are provided by Tachiev & Froese Fischer (1999); Froese Fischer & Tachiev (2004) [MCHF/BP], available in the NIST ASD (Kramida et al., 2016), in Table 6 we compare these data with the values of the present MCDHF/RCI calculations using the AS active set in Ne vii. Our MCDHF/RCI values of the AS active set are also included in the table for comparison. Our two data sets (AS and AS) agree well within 3%, which indicates that the present MCDHF/RCI calculations using the AS active set are well converged. The MCHF/BP values also show good agreement with our results except for the weak and intercombination transitions, which differ by more than 20%. It is also seen that transition data listed in the NIST ASD agree well with our results for most transitions. Differences larger that 20% occur for transitions of relatively small radiative branching ratios only, the transition being an exception. Although the value for this transition is about , which is relatively small, the radiative branching ratio of this transition is 40.3%. This means that the above transition gives an important contribution to the lifetime of the level . An inspection of Table 6 reveals that the accuracy of the NIST values is overestimated for some transitions, especially for the transitions with the estimated accuracy A (), and the present MCDHF/RCI data should be more accurate than those listed in the NIST ASD.

To further access the accuracy of the present MCDHF/RCI transition data, in Figure 2, we compare the line strengths from different resources (MCDHF/RCI (AS), MCHF/BP, and NIST) with the results of the present MCDHF/RCI (AS) calculations for the strong transitions with among the lowest 20 levels up to the configuration for the ions from B ii to Ne vii. Our two data sets (AS and AS) agree very well within 5% for all the transitions of the ions from B ii to Ne vii, which again indicates that the present MCDHF/RCI calculations using the AS active set are well converged. The MCHF/BP values also show good agreement with our results, within 10%, whereas the NIST values differ from the present calculations and the MCHF results by 10% - 60% for a few transitions. For example, the MCDHF/RCI (AS), MCDHF/RCI (AS), and MCHF/BP values for the transition are , , and , respectively, while the corresponding NIST value is about 50% greater (). As a further check of the accuracy of the present transition data, we compare the values of our MCDHF/RCI (AS) and MCDHF/RCI (AS) calculations for transitions involving highly- excited states above the configuration in Ne vii. Our two data sets agree within 10% for most of the transitions. However, large deviations sometimes occur for weak one-photon-two-electrons transitions 1, which are generally very sensitive to electron correlation effects and sometimes to high-order relativistic effects.

3.3 Lifetimes

Lifetimes in Ne vii obtained in the Babushkin gauge in the current work (MCDHF/RCI (AS) and reported by the MCHF/BP calculations  (Tachiev & Froese Fischer, 1999; Froese Fischer & Tachiev, 2004) are compared with the MCDHF/RCI (AS)) calculations in Figure 3. To further check the accuracy of the MCDHF/RCI lifetimes, we also calculate the lifetimes in the Coulomb gauge using the present AS model. The results are also included in Figure 3 for comparison. The agreement between the current lifetimes and the values reported by  Tachiev & Froese Fischer (1999); Froese Fischer & Tachiev (2004) is very good, well within 5%. The differences of the present MCDHF/RCI (AS) and MCDHF/RCI (AS) calculations in the Babushkin gauge are within 5% for almost all levels, with the largest deviation of 6% for the level . The MCDHF/RCI (AS) values in the Coulomb gauge also show good agreement (within 5%) with the Babushkin-gauge values from the same calculation, with the two exceptions: and . These two exceptions, with large deviations, occur for states with lifetimes determined by transitions with relatively small values. For example, the lifetime of the level is determined by the weak intercombination transition with the value of . The contributions from the negative energy continuum, which are important for the transition rates in very weak intercombination transitions, are not included in the present MCDHF/RCI calculations in the Coulomb gauge. Therefore, the accuracy of the values in the Coulomb gauge, but not in the Babushkin gauge, are affected for the weak intercombination transition (Cheng et al., 2008).

Lastly we compare our MCDHF/RCI lifetimes in the Babushkin gauge with the experimental values from different sources and the MCDHF/RCI (Jönsson et al., 1998) and RCI (Cheng et al., 2008) calculations for the levels of the ions from B ii to Ne vii. The agreement between the current lifetimes and lifetimes reported by different experimental and theoretical sources is surprisingly good for both the resonance transition and the intercombination transition . This is a strong indicator of the high accuracy of the present calculations.

3.4 Summary

Self-consistent MCDHF calculations and subsequent RCI calculations taking into account the Breit interaction and leading QED effects have been performed for the states of Be-like ions from B ii to Ne vii. Level energies, transition energies, transition rates and lifetimes have been computed. Labeling of the states has been facilitated by converting from the -coupling to the -coupling scheme. Previous experimental and theoretical data involving low-lying states along the isoelectronic sequence, as well as both low-lying and high-lying states for Ne vii, are used to validate the present computations. An excellent agreement with experimental data is found, which is manifested by an average energy difference with the standard deviation for the states of only cm along the isoelectronic sequence. For all the 138 levels of Ne vii the average relative energy difference is . Lifetimes in the Babushkin and Coulomb gauge agree to within 5% for most states of Ne vii. We have presented calculations for energy and transition data that are more extensive and accurate than previous work.

To summarize, the present work has significantly increased the amount of accurate data for the Be-like isoelectronic sequence, extending our previous calculations (Wang et al., 2015). The complete data set including both energy and transition results provided by the present work and by Wang et al. (2015) can be used for line identification and modeling purposes involving the high-lying states, which fills the gap for lacking accurate atomic data on Be-like ions. Our study can also be considered as a benchmark for other calculations.

Scientific software packages

Scientific software packages including \softwareGRASP2K (Jönsson et al., 2007, 2013) and FAC (Gu, 2008) are used in the present work.

We thank the anonymous referee for many constructive comments on this manuscript. We acknowledge the support from the National Key Research and Development Program of China under Grant No. 2017YFA0403200, the National Natural Science Foundation of China (Grant Grant No. 11703004, No. 11674066, No. 11504421, and No. 11474034) and the Natural Science Foundation of Hebei Province, China (A2017201165). This work is also supported by the Swedish research council under contracts 2015-04842 and 2016-04185. The authors (K.W. and M.F.G.) express their gratitude to the support from the visiting researcher program at the Fudan University.
Figure 1: Energy differences in cm of the MCDHF/RCI and NIST excitation energies as a function of for the and levels. The data come from Table 4.
Figure 2: The differences () in % of line strengths (in length form) for the transitions among the lowest 20 levels up to the configuration for the ions from B ii to Ne vii from different resources, MCDHF/RCI (AS), MCHF/BP, and NIST, relative to the values of the present MCDHF/RCI (AS) calculations.
Figure 3: The differences () in % of lifetimes for all the 138 levels in Ne vii from different resources, MCDHF/RCI (AS) in Babushkin gauge, MCDHF/RCI (AS) in Coulomb gauge, and MCHF/BP, relative to the values of the present MCDHF/RCI (AS) calculations in Babushkin gauge.
Key Level -composition
10 1 0 0.000E+00 0.94 + 0.06 
10 2 111197 0.000E+00 1.00
10 3 111649 5.277E-05 1.00
10 4 112643 1.419E+01 1.00
10 5 215065 2.508E-10 0.99
10 6 289259 3.249E-10 0.99
10 7 289771 3.235E-10 0.99
10 8 290653 3.213E-10 0.99
10 9 317716 1.811E-09 0.99
10 10 393321 1.769E-10 0.92 + 0.06 
10 11 978197 1.788E-11 0.97 + 0.02 
10 12 998117 5.560E-11 0.96 + 0.04 
10 13 1025546 8.302E-12 0.90 + 0.08 
10 14 1028265 1.758E-09 0.98
10 15 1028397 9.983E-10 0.98
10 16 1028650 1.710E-09 0.98
10 17 1054223 4.144E-12 0.98
10 18 1054264 4.149E-12 0.98
10 19 1054337 4.157E-12 0.98
10 20 1071866 6.468E-12 0.98
10 21 1119967 2.374E-11 0.98
10 22 1120462 2.365E-11 0.98
10 23 1121548 2.348E-11 0.98
10 24 1138711 2.016E-11 0.90 + 0.08 
10 25 1148548 1.547E-11 0.99
10 26 1154407 3.165E-11 0.97
10 27 1154989 3.192E-11 0.98
10 28 1155998 3.180E-11 0.98
10 29 1165319 1.790E-11 0.96 + 0.02 
10 30 1171685 1.827E-11 0.99
10 31 1172082 1.827E-11 0.99
10 32 1172665 1.829E-11 0.99
10 33 1177587 1.542E-10 0.93 + 0.06 
10 34 1178262 1.990E-09 0.99
10 35 1179037 7.291E-09 0.99
10 36 1179839 1.013E-11 0.94 + 0.06 
10 37 1184428 1.249E-11 0.98
10 38 1193771 3.455E-12 0.99
10 39 1193989 3.467E-12 0.99
10 40 1194354 3.450E-12 1.00
10 41 1199678 6.441E-12 0.96
10 42 1200102 6.450E-12 0.96
10 43 1200339 6.472E-12 0.97
10 44 1205819 2.755E-11 0.92 + 0.04  + 0.02 
10 45 1212072 3.000E-12 0.97
10 46 1218082 5.155E-12 0.98
10 47 1298906 4.119E-11 0.99
10 48 1306994 4.656E-11 0.96
10 49 1317970 1.274E-10 0.99
10 50 1318016 1.259E-10 0.99
10 51 1318126 1.283E-10 0.99
10 52 1319731 1.633E-11 0.97
10 53 1328043 1.056E-11 1.00
10 54 1328057 1.057E-11 1.00
10 55 1328082 1.058E-11 1.00
10 56 1333874 1.356E-11 0.99
10 57 1335111 2.730E-11 0.99
10 58 1335119 2.730E-11 0.99
10 59 1335132 2.730E-11 0.99
10 60 1336905 2.652E-11 0.98
10 61 1430746 4.405E-11 0.99
10 62 1431042 4.180E-11 0.89 + 0.09 
10 63 1432355 4.336E-11 0.99
10 64 1433369 2.986E-11 0.74 + 0.15  + 0.10 
10 65 1436099 3.333E-11 0.85 + 0.14 
10 66 1438826 1.136E-10 0.89 + 0.10 
10 67 1443735 2.718E-11 0.86 + 0.13 
10 68 1445277 2.697E-11 0.84 + 0.13 
10 69 1445610 2.669E-11 0.97 + 0.02 
10 70 1446554 2.583E-11 0.97 + 0.03 
10 71 1446973 9.947E-11 0.94 + 0.05 
10 72 1446984 9.799E-11 0.94 + 0.06 
10 73 1447007 9.524E-11 0.94 + 0.06 
10 74 1450833 3.157E-11 0.81 + 0.16  + 0.03 
10 75 1450874 2.948E-11 0.99
10 76 1450992 3.201E-11 0.79 + 0.17  + 0.03 
10 77 1451781 2.927E-11 0.98
10 78 1452229 4.105E-11 0.68 + 0.21  + 0.11 
10 79 1452978 2.490E-11 0.98
10 80 1452995 2.532E-11 0.97 + 0.02 
10 81 1453024 5.486E-11 0.60 + 0.37  + 0.02 
10 82 1453025 2.602E-11 0.97 + 0.03 
10 83 1453199 6.466E-11 0.75 + 0.24 
10 84 1453257 6.656E-11 0.58 + 0.18  + 0.18 
10 85 1453375 5.527E-11 0.55 + 0.43 
10 86 1453496 6.837E-11 0.79 + 0.21 
10 87 1453566 7.312E-11 0.61 + 0.19  + 0.15 
10 88 1453809 5.744E-11 0.53 + 0.47 
10 89 1454855 3.425E-11 0.86 + 0.12 
10 90 1454971 8.321E-11 0.62 + 0.37 
10 91 1455742 1.957E-11 0.95 + 0.04 
10 92 1456542 1.456E-11 0.60 + 0.38 
10 93 1457897 4.559E-11 0.58 + 0.38 
10 94 1458067 4.667E-11 0.55 + 0.43 
10 95 1458479 5.025E-11 0.52 + 0.47 
10 96 1460317 8.417E-12 0.97 + 0.02 
10 97 1460552 8.648E-12 0.93 + 0.05 
10 98 1460970 8.447E-12 0.98
10 99 1461766 2.734E-11 0.68 + 0.28  + 0.03 
10 100 1461890 2.823E-11 0.97
10 101 1462168 2.801E-11 0.68 + 0.29 
10 102 1462286 2.844E-11 0.98
10 103 1463091 1.423E-11 0.87 + 0.05  + 0.05 
10 104 1463451 1.444E-11 0.91 + 0.05  + 0.02 
10 105 1463646 1.466E-11 0.93 + 0.05 
10 106 1466327 2.865E-11 0.96
10 107 1466383 3.283E-11 0.74 + 0.23 
10 108 1466386 3.720E-11 0.85 + 0.09 
10 109 1466595 3.349E-11 0.64 + 0.20  + 0.11 
10 110 1466598 2.884E-11 0.80 + 0.16  + 0.02 
10 111 1467077 2.866E-11 0.99
10 112 1467252 3.225E-11 0.78 + 0.21 
10 113 1467610 2.965E-11 0.80 + 0.18 
10 114 1467806 3.740E-11 0.68 + 0.17  + 0.12 
10 115 1469623 7.537E-12 0.86 + 0.12 
10 116 1470147 1.158E-11 0.95 + 0.03 
10 117 1512005 1.107E-10 0.99
10 118 1514042 1.156E-10 0.98
10 119 1517319 2.179E-10 1.00
10 120 1517332 2.167E-10 0.99
10 121 1517363 2.194E-10 1.00
10 122 1517824 4.591E-11 0.98
10 123 1520185 3.602E-11 1.00
10 124 1520189 3.605E-11 1.00
10 125 1520196 3.610E-11 1.00
10 126 1521680 4.097E-11 1.00
10 127 1521981 8.813E-11 1.00
10 128 1521984 8.813E-11 1.00
10 129 1521987 8.813E-11 1.00
10 130 1522131 2.538E-10 1.00
10 131 1522133 2.538E-10 0.71 + 0.29 
10 132 1522134 2.539E-10 0.71 + 0.29 
10 133 1522135 2.539E-10 1.00
10 134 1522380 1.586E-10 0.99
10 135 1522384 1.585E-10 0.99
10 136