Next Article in Journal
Cross Sections and Rate Coefficients for Vibrational Excitation of H2O by Electron Impact
Next Article in Special Issue
Revealing the Target Electronic Structure with Under-Threshold RABBITT
Previous Article in Journal
Evaluation of Fe XIV Intensity Ratio for Electron Density Diagnostics by Laboratory Measurements
Previous Article in Special Issue
A Tribute to Oleg Zatsarinny (1953–2021): His Life in Science
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Benchmark Angle-Differential Cross-Section Ratios for Excitation of the 4p5s Configuration in Krypton

1
Department of Physics, California State University, Fullerton, CA 92831, USA
2
Department of Physics and Astronomy, Drake University, Des Moines, IA 50311, USA
*
Author to whom correspondence should be addressed.
Deceased.
Submission received: 10 August 2021 / Revised: 27 August 2021 / Accepted: 29 August 2021 / Published: 2 September 2021

Abstract

:
Benchmark intensity ratio measurements of the energy loss lines of krypton for excitation of the 4 p 6 1 S 0 4 p 5 5 s [ 3 / 2 ] 2 , 4 p 5 5 s [ 3 / 2 ] 1 , 4 p 5 5 s [ 1 / 2 ] 0 , and 4 p 5 5 s [ 1 / 2 ] 1 transitions are reported, these being the lowest electronic excitations for krypton. The importance of these ratios as stringent tests of theoretical electron scattering models for the noble gases is discussed, as well as the role of spin-exchange and direct processes regarding the angular dependence of these ratios. The experimental data are compared with predictions from fully-relativistic B-spline R-matrix (close-coupling) calculations.

1. Introduction

The heavy rare gases are of great interest in both industry and collision physics. There is a wide array of applications involving the use of species such as krypton and xenon, with the most popular being that they act as excellent buffers in plasmas. Krypton is found primarily in the Earth’s atmosphere. It is widely used in electron impact-initiated plasmas for fluorescent lighting [1], flash lamps for use in high-speed photography [2], and high-power gas lasers [3]. In fundamental physics, species such as krypton are excellent candidates for studies in spectroscopy, primarily because of their rich target structure involving spin-exchange and relativistic effects such as spin–orbit coupling, which are important in both the target structure and the scattering continuum [4].
A comprehensive theoretical study of electron scattering for the noble gases Ne, Ar, Kr, and Xe based on the semi-relativistic distorted-wave method was reported by Bartschat and Madison [5]. Exchange effects were approximated by the semi-empirical local model potential suggested by Furness and McCarthy [6], while relativistic effects were accounted for through the one-electron spin–orbit term of the Breit–Pauli Hamiltonian. These potentials were added to the static Hartree potential generated from the final-state electronic configuration of the target, and the distorted waves were then calculated in that potential. A few years later, Zuo et al. [7] developed a fully relativistic distorted-wave approximation (RDWA), yielding results of comparable accuracy relative to those obtained by Bartschat and Madison. In both models, the  ( n p 6 ) 1 S 0 ground state and the n p 5 ( n + 1 ) s excited states were described by single-configuration states. Note that the two excited states with total electronic angular momenta J = 0.2 are actually well L S -coupled, while the two J = 1 states are a linear combination of two states, with the details depending on the coupling scheme used. Bartschat and Madison, for example, adopted the so-called “intermediate coupling scheme” in which a linear combination of L S -coupled 1 P 1 and 3 P 1 state is used.
Table 1 shows the result of a multi-configuration expansion of the states of interest for the present work. The Cowan code [8] only uses one term-averaged 5s orbital. One can clearly observe that a single configuration dominates in all cases, and also that the two J = 1 states are essentially a 50/50 mix of states with singlet and triplet spin character.
A major breakthrough in the treatment of electron–atom collision came in the 1990s with the convergent close-coupling (CCC) [11,12] and R-matrix (close-coupling) with pseudostates (RMPS) [13,14] methods. The general idea is the same, namely to extend the traditional low-energy close-coupling method to intermediate and even high energies by introducing so-called pseudo states. These have the mathematical properties of bound states (most importantly, they are easily normalizable), but due to their confinement to either a hard (R-matrix) or soft (CCC with a Laguerre basis) box, they provide a discretization of the target continuum. While these methods have been very successful for light (quasi-)one-electron and (quasi-)two-electron atomic and ionic targets, even the H 2 molecule [15], the situation is much less satisfactory for more complex, open-shell targets, especially if these targets are also heavy; hence, both electron correlation and relativistic effects need to be described properly by the theoretical approaches.
Concerning the heavy noble gases in particular, this is a situation where the targets are convenient for experimental investigations but difficult for theory. Numerous high-quality experimental datasets have existed already for a while (see, for example, refs. [16,17,18,19] from our own group as just a small selection), but comparisons with those data have been rare due to the challenges faced by theory. The currently available version of the CCC code remains limited to one-electron and two-electron systems outside of a structureless 1 S 0 core. The well-known R-matrix code of the Belfast group [20] has been extended to the RMPS framework, but applications to the heavy noble gases using the full power of the approach have also been limited—likely due to the fact that the target description with strongly term-dependent one-electron valence orbitals remains a challenge, although sufficiently large configuration interaction (CI) expansions should be able to solve this problem even if all one-electron orbitals are forced to be orthogonal to each other.
A very promising method available for these targets is the B-spline R-matrix (BSR) approach that allows the use of non-orthogonal sets of one-electron orbitals and employs B-splines to represent them. Technically, this has the advantages of compact CI expansions yielding sufficient accuracy as well as accurate and efficient integration schemes. A general computer code that can be used for calculations of atomic structure, photo-ionization, and electron collisions in a non-relativistic and semi-relativistic (Breit–Pauli) framework was published by Zatsarinny [21] and used for numerous calculations of the above processes. A review of the method and its applications until 2013 was published by Zatsarinny and Bartschat [22]. Similarly to the Belfast code, the available suite of computer codes is general, i.e., it is applicable to targets such as the heavy noble gases of interest for the present work. A fully-relativistic version described by Zatsarinny and Bartschat [23] also exists. While a comprehensive write-up is not available due to the untimely death of Dr. Zatsarinny in March 2021, the code with instructions is freely available on his GitHub site [24].

2. Cross-Section Ratios

Disagreements observed between the theory and experiment, even for angle-integrated cross sections, are often accredited to the challenges faced in the absolute normalization of the experimental data using conventional electrostatic electron spectrometers [25]. This problem has been difficult to overcome systematically when normalizing inelastic scattering features, with non-zero energy loss E L , to the elastic scattering feature ( E L = 0 ) as a standard. This is due to the required accurate characterization of the spectrometer analyzer’s detection efficiency for different E L values in the spectrum. Therefore, in order to obtain absolute inelastic differential cross sections (DCS), the scattering intensity of inelastic features in the same energy loss spectrum is normalised to the elastic scattering intensity by using DCSs for elastic scattering at the same incident electron energy ( E 0 ) and scattering angle θ , assuming that the detection efficiency across the electron energy loss spectrum is accurately known. Elastic scattering DCSs can be accurately determined by using the standard relative flow method [26] and can be used as a calibration in the observed electron energy loss spectra. However, the added problem to determine accurately normalized inelastic DCSs suffers, as aforesaid, from the accurate determination of the detection efficiency of the scattered electron detector across the observed energy loss spectrum, i.e., as a function of the residual energy of the scattered electron energy E R ( = E 0 E L ) . Knowing this, the absolute inelastic DCSs may be determined from the spectrum if the elastic DCS is known in tandem with the analyzer detection efficiency. This is indeed the case for time-of-flight electron spectrometers, where the detection efficiency is independent of E R of the scattered electron.
At CSUF, we recently built an electron time-of-flight (TOF) spectrometer to determine accurate elastic-to-inelastic differential electron scattering ratios for the excitation of the X 1 Σ g + → b 3 Σ u + transition in H 2 . These benchmark ratios were found to be in excellent agreement with the CCC predictions of the Curtin University group [27]. After normalizing our TOF spectrum’s elastic feature to corresponding well-established elastic electron scattering DCSs from our group [28], these resulted, in turn, in benchmark theoretical and experimental DCSs for elastic and inelastic processes for the H 2 molecule. Since it is much easier and more reliable to measure relative ratios in scattering experiments, we decided to revisit inelastic scattering ratios for the lowest transitions in Kr to provide further experimental benchmark data for testing existing and future collision models. Such ratio measurements were earlier carried out by us for Ne [29], Kr [18,19], and Xe [16]. DWBA calculations for these ratios, particularly regarding the observed deviations from statistical ratios, were also reported by Bartschat and Madison [30]. Referring to the excellent agreement between the experimental ratios of [29] and the predictions from their semi-relativistic BSR calculations for Ne, Zatsarinny and Bartschat [31] recommended such relative inelastic scattering ratios in the heavy noble gases as representing “a very sensitive test to the quality of the theoretical model”.
The present work with a Kr target focuses on inelastic scattering ratios only in order to draw attention to such ratios as benchmarks to aid theory in a critical area of experimental and theoretical electron collision physics. DCS ratios have been determined for the first four excitation levels. We define these ratios as follows.
r ( E 0 , θ ) = I 1 ( E 0 , θ ) I 3 ( E 0 , θ ) σ ( 5 s [ 3 / 2 ] 2 ) σ ( 5 s [ 1 / 2 ] 0 ) ,
r ( E 0 , θ ) = I 2 ( E 0 , θ ) I 4 ( E 0 , θ ) σ ( 5 s [ 3 / 2 ] 1 ) σ ( 5 s [ 1 / 2 ] 1 ) ,
r ( E 0 , θ ) = I 1 ( E 0 , θ ) I 2 ( E 0 , θ ) σ ( 5 s [ 3 / 2 ] 2 ) σ ( 5 s [ 3 / 2 ] 1 ) ,
r ( E 0 , θ ) = I 3 ( E 0 , θ ) I 4 ( E 0 , θ ) σ ( 5 s [ 1 / 2 ] 0 ) σ ( 5 s [ 1 / 2 ] 1 ) .
Here, I m ( E 0 , θ ) ( m = 1 4 , see labels in Figure 1) are the intensities of the electron energy loss features for excitation of the 4 p 6 4 p 5 5 s , 5 s [ K ] J transition, which are proportional to the differential cross section σ . In our notation and discussion below, we use the intermediate-coupling scheme [8,29] that is appropriate for the heavy noble gases in which K is the total angular momentum of the parent ion 4p 5 core and J is the total angular momentum of the core plus valence 5s, 5s’ electron. The ratio r defined in Equation (1) addresses the ratio between the pure triplet 5 s [ 3 / 2 ] 2 ( 3 P 2 ) and 5 s [ 1 / 2 ] 0 ( 3 P 0 ) states (with its approximate L S -coupled component in brackets), respectively. Since the statistical weights are 2 J + 1 , one would expect the DCS ratio to be close to five for a well L S -coupled system, provided one can neglect the energy difference between the two states as well as the potential term dependence of the valence orbital. In order to indicate those possible effects, we follow the standard notation of using 5 s and 5 s for the outer orbital when the parent ion is either the ( 4 p 5 ) 2 P 3 / 2 (for 5 s ) or the ( 4 p 5 ) 2 P 3 / 2 (for 5 s ) state.
On the other hand, the observable r defined in Equation (2) addresses the ratio for the 5 s [ 3 / 2 ] 1 and 5 s [ 1 / 2 ] 1 states. While the appropriate intermediate-coupling coefficients can be deduced from the respective line strengths if these two J = 1 states are written as a linear combination of L S -coupled 1 P 1 and 3 P 1 states [5], the ratio  r can be used as a test of how well a collision model accounts for the respective excitation portions as well as the interference of the singlet and triplet contributions [29]. At sufficiently high energies and small angles, one would expect r a 2 / b 2 , where a and b are the coefficients of the 1 P 1 term in the intermediate-coupling description of the two states [30]. In the present case (see Table 1), this would suggest an asymptotic ratio of  r 1.2 .
Next, r defined in Equation (3) represents the ratio of the L S -coupled 5 s [ 3 / 2 ] 2 ( 3 P 2 ) state and the mixed-coupled 5 s [ 3 / 2 ] 1 state, while r defined in Equation (4), in a similar fashion as  r , is the ratio between the pure σ ( 5 s [ 1 / 2 ] 0 ( 3 P 0 ) state and the other mixed state, ( 5 s [ 1 / 2 ] 1 ) . From an experimental point of view, it is important to note that r and r are less affected by the instrumental transmission efficiency than r and r , since the features in r and r are very closely spaced in energy loss. This is especially important at low E 0 values close to threshold and will be further addressed in the following experimental section.
In this paper, the recently obtained inelastic DCS ratios for Kr are presented in comparison with the respective theoretical predictions, as well as the earlier experimental ratios obtained in our group by Guo et al. [4]. A brief summary of the fully relativistic D(irac)BSR electron scattering model for Kr is given in the next section. More details, also including the Xe target, can be found in [32]. This is followed by a discussion of the experimental apparatus and the measurement procedures for obtaining the above-mentioned DCS. The experimental and theoretical results are presented and discussed in detail in Section 5 before we finish with a summary and our conclusions in Section 6.

3. Theory

Since the details of the calculation were given by Zatsarinny and Bartschat [32], we only summarize the main points here. The calculations reported in this paper were performed using the R-matrix (close-coupling) approach, as implemented in the DBSR suites of computer codes. The initial structure calculation for Kr + was carried out with the GRASP2K relativistic atomic-structure package [33]. After that, the valence orbitals were generated in a frozen-core calculation for Kr + by using the average-term approximation. All these states of Kr + were then used as target states in B-spline bound-state close-coupling calculations in order to generate the low-lying states of atomic Kr (with N = 36 electrons) employing non-orthogonal, term-dependent orbitals for each Kr state.
In the scattering calculations, we included the lowest 31 physical states of Kr, i.e., the fine-structure levels with configurations 4 p 5 5 s , 4 p 5 5 p , 4 p 5 4 d , and 4 p 5 6 s , respectively. This model will be referred to as DBSR-31 below. Since this is a fully-relativistic approach, the mixing coefficients in the multi-configuration expansions are mathematically not the same as in the intermediate-coupling scheme. This is a similar situation as in the RDW model [7]. In a truly complete expansion, the specifics of the numerical basis would not matter. While this ideal situation is not achievable in practice, it still makes some sense to interpret the final results in the intermediate-coupling scheme, which appears to capture the essential physics regarding the orbital and spin angular momentum character of the states in question.
We then used the DBSR version [23] to solve the ( N + 1 ) electron collision problem. We calculated partial-wave contributions up to J = 61 / 2 numerically. No extrapolation scheme to account for even higher partial waves was necessary for all observables presented in this paper.
At this time, the only results from the 31-state model are available. In the original paper [32], 47 state results were also presented, but both calculations were limited to energies below 13.5 eV. Note that the general features seen in the near-threshold region were reproduced in both models (except very close to threshold, where there are no experimental data for the ratios), although the absolute values differed. In the comparison with experiment, the available data were visually normalized to the 47 state results, which made the latter calculation appear more superior. Since we compare ratios in the present paper, the 31 state model should be reasonable. Ideally, one would extend these calculations much further to include a large number of pseudo-states to also couple to the ionization continuum. Due to the passing of Dr. Zatsarinny, however, this is not realistic in the timeframe for this Special Issue dedicated to him.

4. Experiment

The present California State University (CSUF) energy loss system is a moderate-current, high-resolution electrostatic energy loss spectrometer, which has been well tested and described in detail before [34]. The system consists of an electron monochromator and an electron energy analyzer, with both employing hemispherical energy selectors. The collimated gas beam is delivered to the collision region via a moveable gas aperture source, which is aligned and placed about 5 mm below the collision region. The entire spectrometer is housed in a vacuum chamber that is pumped with a 12 inch diffusion pump down to a base pressure of ≈ 1 × 10 7 Torr. In order to create and maintain an environment suitable for measuring stable low-energy electron energy loss spectra over long periods of time, the electron gun and energy analyzer are both baked to about ≈120 C during the experiment, and the vacuum chamber is properly oil-baffled.
The experimental apparatus is computer controlled via a LabVIEW program that is run in the multi-channel scaling mode and monitors the input pressure of the target gas, steers the moveable gas source in and out of the interaction region, drives a stepper motor that sets the scattering angle position of the electron energy analyzer, and controls the energy ramp of the electron energy analyzer. The experimental energy loss spectra acquired by the multi-channel scaling program are analyzed by a separate software that employs a multi-Gaussian instrumental line-profile unfolding technique to fit the energy loss spectra [34].
Krypton energy loss spectra were acquired at E 0 values of 11.5 eV, 12.0 eV, 13.5 eV, 15.0 eV, and 20.0 eV for θ ranging from 15 to 130 with a typical energy resolution of around 35–45 meV for an incident current of about 13 nA to 22 nA. E 0 was calibrated using the He 2 2 S resonance at 19.366 eV at θ = 90 [35,36] to obtain an E 0 with an accuracy of 50 meV or better during the entire run at that E 0 value. The CSUF spectrometer and the earlier spectrometer in [18] differ somewhat in that the former used real apertures in the analyzer as opposed to virtual apertures in the latter. Additionally, the CSUF spectrometer employs an aperture gas target collimation system rather than the earlier hypodermic needle gas collimator. In this experiment, the moveable gas source was kept fixed because the system was not observing elastic scattering as it did in previous work [37] in which elastic background scattering from surfaces was found to be significant. This was also conducted to improve the acquired scattering counts (more acquisition time for inelastic measurements), and thus minimize statistical errors.
An important consideration in this experiment was to control the transmission efficiency for the different features (features one to four in Figure 1), which had respective energy loss values of 9.915 eV, 10.033 eV, 10.563 eV, and 10.644 eV. At small E 0 values, the residual energies of the scattered electrons become significantly different across the energy loss spectrum. Table 2 shows the ratio Δ E R / E R ¯ ( = α ) at different E 0 values for the four ratios, where, e.g., in the case of r, features one and three, Δ E R = E R ( 1 ) E R ( 3 ) is the difference between the residual energies of scattered electrons having excited features one and three, and E R ¯ = ( [ E R ( 1 ) + E R ( 3 ) ] / 2 ) is the mean residual energy of these features. The parameter α indicates whether the transmission efficiency is approximately a linear function of E R , i.e., the approximate value of the fractional difference in the transmission of electrons for the features involved. At E 0 = 11.5 eV, for example, r and r have relative fractional transmission differences of nearly 50% while that of r and r is only ≈2%. This means that the measured ratios r and r are significantly more precise as benchmarks than r and r . In order to circumvent this transmission problem, the scattered electron detector analyzer was first tuned to the average mid-residual energy of the features one and three or two and four, which are similar for both r and r at any E 0 value. This employed tuning the gun and analyzer to the elastic scattering at 90 at E 0 = E R ¯ before setting up the electron gun to deliver a focused beam at E 0 . At a working E 0 of 13.5 eV, for instance, the analyzer was tuned to elastically scattered electrons at E 0 = ( 3.314 + 3.325 ) / 2 eV = 3.320 eV (see Table 2). As a result, the transmission peaked at the mid- E L of the states one and three or two and four, while it fell equally at E L values of the features on either side of this maximum, thus giving consistently reproducible and accurate r and r values by breaking away from a skewed transmission function for the r and r features. This procedure does not affect r and r , which were always well reproduced for any tuning conditions of the analyzer. At E 0 15 eV, this problem is diminished, and the r and r ratios were found to be reproducible in the ≈10–12% region.

5. Results and Discussion

Table 3, Table 4, Table 5, Table 6 and Table 7 list the ratios r, r , r , and r at different E 0 values. Figure 2, Figure 3, Figure 4 and Figure 5 exhibit the ratios from the current measurements compared to our past ratios of Guo et al. [18] and to the DBSR-31 results. These figures will be discussed in detail below.
Figure 2 shows our present r values at several E 0 values compared to our previous measurements of this parameter and theory. We observe excellent agreement for all E 0 values with our earlier r values [18], as well as with the DBSR-31 predictions. We emphasize that the present data have significantly smaller uncertainties than our earlier measurements [18]. This is due to the better statistical counts in the present work that focused only on the lowest four excited states of Kr rather than an extended spectrum obtained before [18,19]. At our lowest E 0 of 12.0 eV, ≈ 1.5 eV above threshold, we observed large deviations from the statistical ratio value of five, with values exceeding 12 (the theoretical prediction is ≈20) at θ 130 . This effect is most likely due to the markedly different residual energies of the scattered electron after exciting the 5 s [ 3 / 2 ] 2 state ( E R = 2.124 eV) compared to the 5 s [ 1 / 2 ] 0 excitation with a relatively much smaller E R of 1.505 eV; see Table 2. The other (likely less important) factor is the presence of spin–orbit interaction in the continuum electron channel. These triplet states are well L S -coupled (cf. Table 1). Even in multi-configuration expansions, both L S -coupling and the 4 p 5 6 s configuration dominate, and the coupling coefficients are very close to unity for these triplet states. Our data for r deviate from five for θ > 0 in an oscillatory manner but converge to r 5 when θ 0 , except at our lowest E 0 of 12.0 eV. This effect, which is observed at all E 0 values of 13.5 eV and above, is at the 20% level for the present E 0 values.
Figure 3 exhibits our present r values at several E 0 values. Once again, we note generally very good agreement with the earlier experimental and the theoretical predictions, except at E 0 = 12.0 eV where we see a somewhat different angular distribution for the present r data, which are ≈50% smaller than the earlier values of [18] and ≈30–50% smaller than the DBSR-31 predictions for θ 30 . This is not due to transmission effects, as can be argued on the basis of our r values, which have essentially the same instrumental transmission dependence (compared with Table 2) of ≈0.48. We did not observe the oscillatory behavior predicted by theory, and the earlier values [18] are scattered at small and large θ . At E 0 = 13.5 eV, the agreement improves considerably. It is excellent at large θ but poorer at θ 50 . The oscillatory behavior of r at this energy is reproduced by all measured and theoretical values. At E 0 = 15.0 eV, the agreement between experiment and theory is excellent, and the oscillatory angular behavior of r is observed in all three data sets. At E 0 = 20.0 eV, the agreement between experiments is excellent, but the measured values are significantly (by ≈20%) below theory. Here, one should expect transmission effects to be small (see Table 2), in the region of a few percent. Nevertheless, the agreement between experiment and theory for such a difficult target is very good and certainly encouraging.
As expected, we observe a trend of r 1.2 = ( 0.74 / 0.67 ) 2 (see Table 1), which is based on the fact that excitation of the singlet part of the wave function for the two states will ultimately dominate that of the triplet part. However, given that this limit is approached better at large than at small angles suggests that spin-dependent effects are still important. While the trend is best for 20 eV, this is not a sufficiently high energy to predict the ratio using such a simple argument.
Figure 4 presents our results for the r ratio, which is our most extensive ratio spanning from threshold of ≈1.5 eV exit energy up to 10 eV (see Table 2). Since transmission for both states is the same, the reproducible characteristics of r is excellent and not influenced by different spectrometer detector tunings. This makes this accurate ratio a very useful parameter for testing theory. This ratio compares the L S -coupled triplet state 5 s [ 3 / 2 ] 2 to the 5 s [ 3 / 2 ] 1 mixed singlet-triplet state. At E 0 = 11.5 eV, we see excellent agreement between the present work and theory, except that the angular dependence is shifted between theory and experiment by a small amount of ≈ 2 . Our r values are also quantitatively smaller by ≈7%, which is within the uncertainty of the measurements. We note (not shown) that the theoretical DCS for the well L S -coupled 5 s [ 3 / 2 ] 2 state, which is only excitable by spin-exchange, is essentially constant in the θ range from 60 to 120 . Thus, it is the decrease in the excitation of the 5 s [ 3 / 2 ] 1 state that results in the rise of r . Furthermore, we see that at small θ 30 rises. As discussed for r above, this result supports the remaining importance of spin-exchange processes at small θ for the energies investigated in this study.
Our findings might also help to interpret the reversal of the angular momentum perpendicular to the scattering plane, ( L ), imparted to the target. This was found at similarly small θ in the excitation of the resonance transition in Ne [38], i.e., of the 2 p 5 3 s [ 3 / 2 ] 1 L S -coupled 1 P 1 component, which is again mixed (but to a lesser extent) with the L S -coupled 3 P 1 component. It is the latter’s excitation via spin-exchange that could enable a physical explanation for the reversal of L compared to the case for a direct singlet → singlet excitation.
At E 0 = 12.0 eV, the theory shows excellent agreement with experiment, but at E 0 of 13.5 eV and 15.0 eV it predicts a clear rise at θ 80 . This is again likely due to a drop in the direct excitation process of the 5 s [ 3 / 2 ] 1 state, which in turn is due to an effective overestimate of the singlet L S -contribution as aforementioned for r . The feature in r that peaks around θ = 30 for E 0 = 20.0 eV is found to be in excellent agreement between both experiments and theory. At this energy, however, the theory predicts oscillations that are not visible to the same extent in the measurements. Overall, we again judge the agreement between experiment and theory as encouragingly good.
For the r ratio, exhibited in Figure 5, we observe a similar behavior as for r . At E 0 = 12.0 eV, the agreement with theory is excellent. The earlier measurements show much more scatter because the DCS (thus the scattering intensity) of the 5 s [ 1 / 2 ] 0 state is the smallest of the four features (see also Figure 2). At small θ , we observe a similar rise in r as for r , which is an indication that spin exchange processes are prevalent at these angles. At E 0 = 13.5 eV, we again note the increased r values at the mid-angles around 80 as well as non-zero r values for θ 0 , similarly as for 12.0 eV. This indicates that small-angle spin-exchange processes exist in this system. At E 0 = 15.0 eV, the overall agreement is excellent, except for the “bump” in r at mid-angles around 80 . At 20.0 eV, a similar oscillatory behavior in r as for r is observed, with some disagreement between experiment and theory for θ 50 . We note again that the feature peaking at θ = 30 (as in r at this E 0 ) is in excellent agreement between both experimental datasets and the DBSR-31 predictions, although somewhat less pronounced compared to r . Clearly, the overall agreement with theory is encouragingly good for this complex target.

6. Conclusions

New and improved measurements of the intensity ratios r, r , r , and an additional useful ratio r for the first four excited states of Kr were presented and compared to earlier ratios taken in our laboratory [18] and compared to benchmark DBSR calculations. The agreement between the measurements is very good for r and r and excellent for r and r . This provides confidence in their overall reliability. The overall agreement with the theoretical predictions is also very good to excellent, but it suggests places where the DBSR model may need to be improved even further.
Similar measurements of ratios will soon be carried out for Xe to further extend the tests of the DBSR results published in [32]. Depending on the outcome, new calculations with a further improved target structure and even more coupled states in the collision part of the problem may be necessary. We hope that the present work on Kr and the planned experiment on Xe will provide an impetus for the theory to further push benchmark calculations for electron scattering from multi-electron targets with more than two electrons outside of a structureless core.

Author Contributions

The project was conceived by M.A.K. who also supervised the experimental part of the study carried out by A.S. and J.-B.F. The calculations were carried out by O.I.Z. and K.B. The manuscript was mostly prepared by M.A.K. and K.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the US National Science Foundation under Grant No. NSF-RUI AMO 1911702 (MAK), which included a post-doctoral fellowship for A. Sakaamini. The work of OZ and KB was supported by the NSF under grant Nos. PHY-1803844, OAC-1834740, PHY-2110023, and the XSEDE allocation TG-PHY-090031.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript or in the decision to publish the results.

Sample Availability

Tabulated results of the cross sections in digital form are available from the authors upon request.

Abbreviations

The following abbreviations are used in this manuscript:
DCSAngle-differential cross section;
DBSRDirac B-Spline R-matrix.

References

  1. Murphy, A.B.; Tam, E. Thermodynamic properties and transport coefficients of arc lamp plasmas: Argon, krypton and xenon. J. Phys. D Appl. Phys. 2014, 47, 295202. [Google Scholar] [CrossRef]
  2. Kawanaka, J.; Shirai, T.; Kubodera, S.; Sasaki, W. 1.5 kW high-peak-power vacuum ultraviolet flash lamp using a pulsed silent discharge of krypton gas. Appl. Phys. Lett. 2001, 79, 3752. [Google Scholar] [CrossRef]
  3. Obenschain, S.; Lehmberg, R.; Kehne, D.; Hegeler, F.; Wolford, M.; Sethian, J.; Weaver, J.; Karasik, M. High-energy krypton fluoride lasers for inertial fusion. Appl. Opt. 2001, 54, F103. [Google Scholar] [CrossRef]
  4. Guo, X.; Trajmar, S.; Zeman, V.; Bartschat, K.; Fontes, C.J. Low-energy electron impact excitation cross sections for the 4p55s levels of krypton—Supplementary test of collisions models for the heavy rare gases. J. Phys. B At. Mol. Opt. Phys. 1999, 32, L155. [Google Scholar] [CrossRef]
  5. Bartschat, K.; Madison, D.H. Electron impact excitation of rare gases: Differential cross sections and angular correlation parameters for neon, argon, krypton and xenon. J. Phys. B At. Mol. Phys. 1987, 20, 5839. [Google Scholar] [CrossRef]
  6. Furness, J.B.; McCarthy, I.E. Semiphenomenological optical model for electron scattering on atoms. J. Phys. B At. Mol. Phys. 1973, 6, 2280. [Google Scholar] [CrossRef]
  7. Zuo, T.; McEachran, R.P.; Stauffer, A.D. Relativistic distorted-wave calculation of electron excitation of heavy noble gases. J. Phys. B At. Mol. Opt. Phys. 1992, 25, 3393. [Google Scholar] [CrossRef]
  8. Cowan, R.D. The Theory of Atomic Structure and Spectra; University of California Academic Press: Riverside, CA, USA, 1981. [Google Scholar]
  9. Fontes, C.J.; Los Alamos National Laboratory, Los Alamos, NM, USA. Personal communication, 1996.
  10. Moore, C.E. Atomic Energy Levels; NSRDS-NBS Publication; US Govt Publishing: Washington, DC, USA, 1957; Volume 35. [Google Scholar]
  11. Bray, I.; Stelbovics, A.T. Convergent close-coupling calculations of electron-hydrogen scattering. Phys. Rev. A 1992, 46, 6995. [Google Scholar] [CrossRef] [Green Version]
  12. Fursa, D.V.; Bray, I. Calculation of electron-helium scattering. Phys. Rev. A 1995, 52, 1279. [Google Scholar] [CrossRef]
  13. Bartschat, K.; Hudson, E.T.; Scott, M.P.; Burke, P.G.; Burke, V.M. Electron—Atom scattering at low and intermediate energies using a pseudo-state/R-matrix basis. J. Phys. B At. Mol. Phys. 1996, 29, 115. [Google Scholar] [CrossRef]
  14. Hudson, E.T.; Bartschat, K.; Scott, M.P.; Burke, P.G.; Burke, V.M. Electron scattering from helium atoms. Phase shifts, resonance parameters and total cross sections. J. Phys. B At. Mol. Opt. Phys. 1996, 29, 5513. [Google Scholar] [CrossRef]
  15. Zammit, M.C.; Savage, J.S.; Fursa, D.V.; Bray, I. Complete solution of electronic excitation and ionization in electron-hydrogen molecule scattering. Phys. Rev. Lett. 2016, 116, 233201. [Google Scholar] [CrossRef] [Green Version]
  16. Khakoo, M.A.; Trajmar, S.; LeClair, L.R.; Kanik, I.; Csanak, G.; Fontes, C.J. Differential cross sections for electron impact excitation of Xe: I. Excitation of the five lowest levels; experiment and theory. J. Phys. B At. Mol. Opt. Phys. 1996, 29, 3455. [Google Scholar] [CrossRef]
  17. Khakoo, M.A.; Trajmar, S.; Wang, S.; Kanik, I.; Aguirre, A.; Fontes, C.J.; Clark, R.E.H.; Abdallah, J. Differential cross sections for electron impact excitation of Xe: II. Excitation of the sixth to twentieth lowest levels; experiment and theory. J. Phys. B At. Mol. Opt. Phys. 1996, 29, 3477. [Google Scholar] [CrossRef]
  18. Guo, X.; Mathews, D.F.; Mikaelian, G.; Khakoo, M.A.; Crowe, A.; Kanik, I.; Trajmar, S.; Zeman, V.; Bartschat, K.; Fontes, C.J. Journal of Physics B: Atomic, Molecular and Optical Physics Differential cross sections for electron-impact excitation of krypton at low incident energies: I. Excitation of the 4p55s configuration. J. Phys. B At. Mol. Opt. Phys. 2000, 33, 1895. [Google Scholar] [CrossRef]
  19. Guo, X.; Mathews, D.F.; Mikaelian, G.; Khakoo, M.A.; Crowe, A.; Kanik, I.; Trajmar, S.; Zeman, V.; Bartschat, K.; Fontes, C.J. Journal of Physics B: Atomic, Molecular and Optical Physics Differential cross sections for electron-impact excitation of krypton at low incident energies: II. Excitation of the 4p55p, 4p54d and 4p556s configurations. J. Phys. B At. Mol. Opt. Phys. 2000, 33, 1921. [Google Scholar] [CrossRef]
  20. Berrington, K.A.; Eissner, W.; Norrington, P.H. RMATRX1: Belfast atomic R-matrix codes. Comp. Phys. Commun. 1995, 92, 290. [Google Scholar] [CrossRef]
  21. Zatsarinny, O. BSR: B-Spline Atomic R-Matrix Codes. Comp. Phys. Commun. 2006, 174, 273. [Google Scholar] [CrossRef]
  22. Zatsarinny, O.; Bartschat, K. XXVIII International Conference on Photonic, Electronic and Atomic Collisions (ICPEAC 2013) IOP Publishing. J. Phys. Conf. Ser. 2014, 488, 012044. [Google Scholar] [CrossRef] [Green Version]
  23. Zatsarinny, O.; Bartschat, K. Relativistic B-spline R-matrix method for electron collisions with atoms and ions: Application to low energy electron scattering from Cs. Phys. Rev. A 2008, 77, 062701. [Google Scholar] [CrossRef]
  24. Oleg Zatsarinny’s GitHub Repository. Available online: https://github.com/zatsaroi (accessed on 1 August 2021).
  25. Zatsarinny, O.; Wang, Y.; Bartschat, K. Electron-impact excitation of argon at intermediate energies. Phys. Rev. A 2014, 89, 022706. [Google Scholar] [CrossRef]
  26. Brunger, M.J.; Buckman, S.J. Electron–molecule scattering cross-sections. I. Experimental techniques and data for diatomic molecules. Phys. Rep. 2002, 357, 215. [Google Scholar] [CrossRef]
  27. Zawadzki, M.; Wright, R.; Dolmat, G.; Martin, M.F.; Hargreaves, L.; Fursa, D.V.; Zammit, M.C.; Scarlett, L.H.; Tapley, J.K.; Savage, J.S.; et al. Time-of-flight electron scattering from molecular hydrogen: Benchmark cross sections for excitation of the X 1 Σ g + b 3 Σ u + transition. Phys. Rev. A 2018, 97, 050702(R). [Google Scholar] [CrossRef] [Green Version]
  28. Muse, J.; Silva, H.; Lopes, M.C.A.; Khakoo, M.A. Low energy elastic scattering of electrons from H2 and N2. J. Phys. B At. Mol. Opt. Phys. 2008, 41, 095203. [Google Scholar] [CrossRef]
  29. Khakoo, M.A.; Wrkich, J.; Larsen, M.; Kleiban, G.; Kanik, I.; Trajmar, S.; Brunger, M.J.; Teubner, P.J.O.; Crowe, A.; Fontes, C.J.; et al. Differential cross sections and cross-section ratios for electron-impact excitation of the neon 2p53s configuration. Phys. Rev. A 2002, 65, 062711. [Google Scholar] [CrossRef]
  30. Bartschat, K.; Madison, D.H. Non-statistical branching ratios for excitation of the (np5(n+1)s)1,3P0,1,2 states in noble gases. J. Phys. B At. Mol. Opt. Phys. 1992, 25, 4619. [Google Scholar] [CrossRef]
  31. Zatsarinny, O.; Bartschat, K. Electron impact excitation of neon at intermediate energies. Phys. Rev. A 2012, 86, 022717. [Google Scholar] [CrossRef]
  32. Zatsarinny, O.; Bartschat, K. Benchmark calculations for near-threshold electron-impact excitation of krypton and xenon atoms. J. Phys. B At. Mol. Opt. Phys. 2010, 43, 074031. [Google Scholar] [CrossRef]
  33. Jönsson, P.; He, X.; Froese Fischer, C.; Grant, I.P. The grasp2K relativistic atomic structure package. Comp. Phys. Commun. 2007, 177, 597–622. [Google Scholar] [CrossRef]
  34. Khakoo, M.A.; Beckmann, C.E.; Trajmar, S.; Csanak, G. Electron-impact excitation of the ns[3/2]o and n[1/2]o levels of Ne, Ar, Kr and Xe. J. Phys. B At. Mol. Opt. Phys. 1994, 27, 3159. [Google Scholar] [CrossRef]
  35. Brunt, J.H.; King, G.C.; Read, F.H. Resonance structure in elastic electron scattering from helium, neon and argon. J. Phys. B At. Mol. Phys. 1977, 10, 1289. [Google Scholar] [CrossRef]
  36. Gopalan, A.; Bömmels, J.; Götte, S.; Landwehr, A.; Franz, K.; Ruf, M.W.; Hotop, H.; Bartschat, K. A novel electron scattering apparatus combining a laser photoelectron source and a triply differentially pumped supersonic beam target: Characterization and results for the He(1s2s2) resonance. Eur. Phys. J. D 2002, 22, 17. [Google Scholar] [CrossRef]
  37. Khakoo, M.A.; Keane, K.; Campbell, C.; Guzman, N.; Hazlett, K. Low energy elastic electron scattering from ethylene. J. Phys. B At. Mol. Opt. Phys. 2007, 40, 3601. [Google Scholar] [CrossRef]
  38. Hargreaves, L.R.; Campbell, C.; Khakoo, M.A.; Zatsarinny, O.; Bartschat, K. Unusual angular momentum transfer in electron-impact excitation of neon. Phys. Rev. A 2012, 85, 050701(R). [Google Scholar] [CrossRef] [Green Version]
Figure 1. Electron energy loss spectrum of Kr taken at E 0 = 13.5 eV and θ = 130 . The data are the red dots, and the line is a linear least squares fit using a multi-Gaussian instrumental line profile centered at the empirical E L values given by Moore [10] and also listed in Table 2. The features m = 1–4 are discussed in the text following Equations (1)–(4).
Figure 1. Electron energy loss spectrum of Kr taken at E 0 = 13.5 eV and θ = 130 . The data are the red dots, and the line is a linear least squares fit using a multi-Gaussian instrumental line profile centered at the empirical E L values given by Moore [10] and also listed in Table 2. The features m = 1–4 are discussed in the text following Equations (1)–(4).
Atoms 09 00061 g001
Figure 2. Ratio r for various E 0 values as a function of θ ; ( Atoms 09 00061 i001) present experimental work; ( Atoms 09 00061 i002) past experimental work of Guo et al. [18]; ( Atoms 09 00061 i003), DSBR-31 calculation.
Figure 2. Ratio r for various E 0 values as a function of θ ; ( Atoms 09 00061 i001) present experimental work; ( Atoms 09 00061 i002) past experimental work of Guo et al. [18]; ( Atoms 09 00061 i003), DSBR-31 calculation.
Atoms 09 00061 g002
Figure 3. Same as Figure 2 but for the r ratio.
Figure 3. Same as Figure 2 but for the r ratio.
Atoms 09 00061 g003
Figure 4. Same as Figure 2 but for the r ratio.
Figure 4. Same as Figure 2 but for the r ratio.
Atoms 09 00061 g004
Figure 5. Same as Figure 2 but for the r ratio.
Figure 5. Same as Figure 2 but for the r ratio.
Atoms 09 00061 g005
Table 1. Intermediate-coupling coefficients for the 4p 6 and 4p 5 5s states of Kr in a multi-configuration expansion, produced with the Cowan code [8] by Fontes [9]. E L values from the Cowan code and from Moore [10] are also listed. The notation [ K ] J (see also below) indicates the total electronic angular momentum of the ( 4 p 5 ) 2 P K ionic core. After coupling to the outer orbital (here 5s), this results in the total electronic angular momentum J of the state. Refer to the text for discussions.
Table 1. Intermediate-coupling coefficients for the 4p 6 and 4p 5 5s states of Kr in a multi-configuration expansion, produced with the Cowan code [8] by Fontes [9]. E L values from the Cowan code and from Moore [10] are also listed. The notation [ K ] J (see also below) indicates the total electronic angular momentum of the ( 4 p 5 ) 2 P K ionic core. After coupling to the outer orbital (here 5s), this results in the total electronic angular momentum J of the state. Refer to the text for discussions.
Level #ConfigurationIntermediate Coupling E L (eV)
CowanMoore
1 4 p 6 0.998 ( 4 p 6 ) 1 S 0 0.000.00
2 4 p 5 5 s [ 3 / 2 ] 2 0.998 ( 4 p 5 5 s ) 3 P 2 9.889.92
3 4 p 5 5 s [ 3 / 2 ] 1 0.67 ( 4 p 5 5 s ) 3 P 1 + 0.74 ( 4 p 5 5 s ) 1 P 1 9.9710.03
4 4 p 5 5 s [ 1 / 2 ] 0 0.996 ( 4 p 5 5 s ) 3 P 0 10.5010.56
5 4 p 5 5 s [ 1 / 2 ] 1 0.74 ( 4 p 5 5 s ) 3 P 1 0.67 ( 4 p 5 5 s ) 1 P 1 10.5610.64
Table 2. Transmission parameters that affect the ratios r, r , r , and  r . The E L values, listed in decreasing order for the 5 s [ 3 / 2 ] 2 , 5 s [ 3 / 2 ] 1 , 5 s [ 1 / 2 ] 0 , and  5 s [ 1 / 2 ] 1 states, are taken from the recommended values of Moore [10]. See text for a detailed discussion.
Table 2. Transmission parameters that affect the ratios r, r , r , and  r . The E L values, listed in decreasing order for the 5 s [ 3 / 2 ] 2 , 5 s [ 3 / 2 ] 1 , 5 s [ 1 / 2 ] 0 , and  5 s [ 1 / 2 ] 1 states, are taken from the recommended values of Moore [10]. See text for a detailed discussion.
E 0 (eV) →11.512.013.515.020.0
E L (eV) ↓ E R (eV) ↓
9.8761.6242.1243.6245.12410.124
9.9691.5312.0313.5315.03110.031
10.4951.0051.5053.0054.5059.505
10.5660.9341.4342.9344.4349.434
Δ E R / E R ¯
r0.4709010.3411410.1867550.128570.06307
r 0.4843810.3445890.1846870.1261490.061341
r 0.0147390.0111910.0064990.0045790.002307
r 0.0183080.0120790.0059770.0039710.001874
Table 3. Ratio r for Kr taken at E 0 = 11.5  eV, with one standard deviation uncertainty. See text for discussion.
Table 3. Ratio r for Kr taken at E 0 = 11.5  eV, with one standard deviation uncertainty. See text for discussion.
θ ( ) r Err
200.2380.018
250.2420.019
300.1950.015
350.1740.014
400.2550.020
450.2900.023
500.3590.028
600.6790.053
701.250.10
802.370.19
902.660.22
1001.700.10
1100.9700.055
1200.5200.030
1300.3700.019
Table 4. Ratios r, r , r , and  r for Kr taken at E 0 = 12.0  eV with one standard deviation error. See text for discussion.
Table 4. Ratios r, r , r , and  r for Kr taken at E 0 = 12.0  eV with one standard deviation error. See text for discussion.
θ   ( ) rErr r Err r Err r Err
2010.31.32.220.220.3830.0390.0820.009
2510.81.52.180.220.4920.0500.0990.012
3011.31.42.120.210.5070.0510.0960.010
3510.91.42.000.200.4850.0490.0880.010
4010.81.41.670.170.5920.0600.0910.010
459.561.201.430.140.6450.0650.0960.010
507.590.951.360.140.6650.0670.1190.013
605.890.741.220.120.7820.0800.1620.017
703.620.451.040.110.9730.1000.2810.030
802.460.310.9940.1031.130.120.4580.048
902.730.341.060.111.080.110.4200.044
1003.460.431.170.120.930.090.3140.033
1104.470.561.310.130.760.080.2240.023
1206.750.851.330.130.650.070.1290.014
13012.71.71.580.160.540.060.0680.008
Table 5. Same as Table 4 but for E 0 = 13.5  eV. See text for discussion.
Table 5. Same as Table 4 but for E 0 = 13.5  eV. See text for discussion.
θ   ( ) rErr r Err r Err r Err
155.490.991.690.220.1110.0120.0340.004
204.950.801.410.180.1720.0170.0490.005
255.650.931.350.180.2530.0260.0610.006
305.560.911.260.170.3390.0340.0780.008
356.141.011.220.160.4180.0430.0830.009
406.141.001.260.170.4780.0490.0990.010
457.801.301.320.170.5130.0520.0870.009
507.501.241.410.190.5140.0520.0970.011
608.281.381.480.200.5910.0600.1070.012
708.151.431.6490.2230.660.070.1340.016
807.541.291.670.220.700.070.1560.018
906.041.021.540.210.690.070.1760.020
1004.980.831.420.190.570.060.1640.018
1104.840.801.360.180.530.050.1490.016
1205.870.971.380.180.540.060.1270.014
1306.81.11.440.190.530.050.1110.012
Table 6. Same as Table 4 but for E 0 = 15.0  eV. See text for discussion.
Table 6. Same as Table 4 but for E 0 = 15.0  eV. See text for discussion.
θ   ( ) rErr r Err r Err r Err
155.400.781.810.090.0580.0060.0190.002
205.920.871.760.090.0970.0100.0290.004
254.980.711.690.090.1600.0170.0540.007
356.350.951.510.090.3470.0370.0830.011
405.500.781.480.090.4040.0430.1090.014
456.010.871.440.080.4040.0430.0970.012
506.420.921.500.090.4460.0470.1050.013
605.840.821.450.090.5100.0540.1300.016
705.600.781.500.090.5770.0610.1540.019
805.680.821.450.090.6590.0710.1690.022
905.240.791.260.090.7030.0780.1700.023
1004.930.741.290.080.5810.0640.1550.020
1104.400.621.230.080.5530.0600.1550.019
1204.880.671.230.070.5390.0570.1360.016
1305.220.701.210.070.6200.0650.1420.016
Table 7. Same as Table 4 but for E 0 = 20.0  eV. See text for discussion.
Table 7. Same as Table 4 but for E 0 = 20.0  eV. See text for discussion.
θ   ( ) rErr r Err r Err r Err
105.510.731.490.080.0860.0090.0230.003
155.690.751.470.080.1460.0150.0380.004
205.100.671.420.070.2250.0230.0620.007
255.980.761.320.070.2780.0280.0610.007
305.610.731.270.070.2900.0300.0660.007
356.120.821.270.070.2640.0270.0550.006
404.970.751.300.070.2370.0260.0620.008
455.290.751.240.070.1940.0200.0460.006
504.870.801.200.070.1960.0220.0500.007
604.350.701.230.070.2010.0230.0570.008
704.680.691.190.070.2590.0280.0660.009
805.030.751.040.060.3870.0420.0800.011
905.280.811.130.070.4520.0500.0970.013
1005.250.851.110.070.4710.0540.1020.015
1105.580.941.090.080.5130.0590.1000.015
1204.850.691.010.060.5810.0620.1210.015
1304.740.650.9640.0560.6240.0660.1280.015
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sakaamini, A.; Faure, J.-B.; Khakoo, M.A.; Zatsarinny, O.I.; Bartschat, K. Benchmark Angle-Differential Cross-Section Ratios for Excitation of the 4p5s Configuration in Krypton. Atoms 2021, 9, 61. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms9030061

AMA Style

Sakaamini A, Faure J-B, Khakoo MA, Zatsarinny OI, Bartschat K. Benchmark Angle-Differential Cross-Section Ratios for Excitation of the 4p5s Configuration in Krypton. Atoms. 2021; 9(3):61. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms9030061

Chicago/Turabian Style

Sakaamini, Ahmad, Jean-Baptiste Faure, Murtadha A. Khakoo, Oleg I. Zatsarinny, and Klaus Bartschat. 2021. "Benchmark Angle-Differential Cross-Section Ratios for Excitation of the 4p5s Configuration in Krypton" Atoms 9, no. 3: 61. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms9030061

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