Debye-screening effect on electron-impact excitation of helium-like Al11+and Fe24+ions

Yu-Long Ma(馬玉龍) Ling Liu(劉玲) Lu-You Xie(頡錄有) Yong Wu(吳勇) Deng-Hong Zhang(張登紅)Chen-Zhong Dong(董晨鐘) Yi-Zhi Qu(屈一至) and Jian-Guo Wang(王建國(guó))

1College of Physics and Electronic Engineering,Northwest Normal University,Lanzhou 730070,China

2Key Laboratory of Computational Physics,Institute of Applied Physics and Computational Mathematics,Beijing 100088,China

3Center for Applied Physics and Technology,Peking University,Beijing 100084,China

4School of Optoelectronics,University of Chinese Academy of Sciences,Beijing 100049,China

Keywords: Debye plasmas, electron-impact excitation, Debye-H¨uckel model potential, He-like Al11+ and Fe24+ions

1. Introduction

Electron-impact excitation (EIE) processes in various plasmas have been subject to continued interest due to their potential application. For example, EIE cross sections are required for calculating the spectral line intensities and especially level populations of nonlocal-thermodynamicequilibrium (NLTE) plasmas.[1,2]Line emissions from the excited states produced by the EIE process provide valuable information on plasma parameters such as temperature and density.[3,4]Theoretical studies on electron-ion collision processes are needed to consider screening effects due to plasma environments,covering both weakly and strongly coupled plasmas. The electron-ion interaction in the strongly coupled plasmas is often described by the ion-sphere model potential.[5-8]For weakly coupled (or Debye) plasmas, the Debye-H¨uckel (DH) potential[9-12]of an ion of positive chargeZwhich interacts with an electron is given by

When a continuum electron interacts with ions that are immersed in Debye plasmas,two types of DH potentials[12,13]are employed to describe the interaction between the continuum electron and bound electron of the target ion, referred to as the continuum-bound(CB)interaction for simplicity, in an electron-ion collision. The first one is a non-spherical DH(NSDH)potential[12]

whereD=(kBTe/4πe2ne)1/2is the Debye screening length,Teandneare the electron temperature and density, respectively,andkBis the Boltzmann constant. It is noted that when the Coulomb coupling parameterΓ=e2/ˉakBTe?1 the Debyescreening potential(1)is valid,where ˉa=(3/4πne)1/3is the average distance between particles.

Due to its spherical symmetry,the SDH potential is well suited to atomic calculations. Therefore, a significant number of theoretical studies have also been devoted to electron- and proton-impact excitation processes based on the SDH potential. The first attempt at using the SDH model was made by Hattonet al.[16]for the study of inelastic scattering of H-like ions in a dense plasma. Later, the Debye-screened impact cross sections were calculated for H-like ions by using the Born approximation and close-coupling methods,[7]semiclassical impact-parameter formalism[17]and relativistic distorted-wave (RDW) method,[18]as well as for He-like ions by using the RDW method.[19-21]Moreover, the Debyescreening effect on the electron-impact ionization of H-like ions[22]and multielectron Ne4+and Au47+ions has been studied with the SDH potential.[23]

Yoon and Jung[13]have performed a comparison of the Debye-screened behavior between the SDH and NSDH potentials for the EIE process of H-like ions. They suggested that the NSDH model is physically more reliable, but only calculated the transition probabilities. To the best of our knowledge, such comparisons have not yet been performed for the EIE cross sections of multielectron ions. Therefore, the motivation of the present study is to present a detailed analysis and comparison of cross sections between NSDH and SDH model potentials that are used to describe the Debye-screened CB interaction for the EIE of multielectron ions.

Apart from the plasma screening on the CB interaction,plasma screening imposed on the target ion will also affect the EIE cross sections. Zhanget al.[24-26]studied the Debyescreening effect on low-energy electron-H scattering using the R-matrix method with pseudostates. Zammitet al.[27-29]also calculated the excitation and ionization cross sections of electron-H/He-atom collisions in Debye plasmas by the convergent close-coupling (CCC) method. They have naturally included the Debye screening on the target ion when the DH potentials for the electron-electron and electron-nucleus interactions are adopted in the(N+1)-electron Hamiltonian,but the influence of the Debye screening imposed on the target ion on the cross sections has not been discussed. Therefore, the influence of the Debye screening imposed on the target ion on the EIE cross section are investigated in the present work.

We shall restrict our discussions to the simplest multielectron system of He-like ions such as Al11+and Fe24+,but it can be extended to other complex ions straightforwardly. Specifically, the dipole-allowed transition 1s21S→1s2p1P is discussed,as this line of atomic systems embedded in plasmas is often investigated to test various models that have been developed to describe the plasma screening effects.[30,31]To obtain the EIE cross sections, we performed fully relativistic multiconfiguration Dirac-Fock(MCDF)and RDW calculations.

2. Theoretical methods

We briefly present the RDW method for calculating the EIE cross sections(please refer to Refs.[32-34]for more details). The EIE cross sectionσiffrom initial state ‘i’ to final state‘f’is expressed in terms of the collision strengthΩifas[33]

whereκiandκfare the initial and final relativistic quantum numbers for the continuum electron, respectively.Ψrepresents the antisymmetric wavefunctions for the total (N+1)-electron system with total angular momentumJT,which consists of the incident(or scattered)continuum electron and the bound target wavefunction.

In the practical RDW calculations,the Coulomb potential between electrons is expanded in the partial-wave form,which is given by

In Eqs. (4) and (5), the bound target wavefunctions are obtained by means of the MCDF method,in which the atomic state wavefunction is represented as a linear combination of configuration state wavefunctions.[36]In order to include the plasma-screening effect on the target ion (changes of target wavefunctions due to Debye screening), the MCDF calculations were performed by the GRASP2K program[37]with some modifications,[35]in which the DH potential is incorporated into the relativistic Hamiltonian. The MCDF method including the Debye-screening effect has been described in detail by Xieet al.[35]as well as Saha and Fritzsche.[38]The distortedwave wavefunctions of the continuum electron are obtained by solving the coupled Dirac equations numerically.[36]The numerical procedures are accomplished by using the component COWF of RATIP,[39]in which the two-body Coulomb potential (6) has been modified by replacement with the screened ones (7) and (8) for the SDH and NSDH potentials, respectively. Moreover,the Debye-screened electron-nucleus potential is considered asVnuc(r)=-(Z/r)exp(-r/D). In addition,the asymptotic continuum wavefunctions are determined numerically for the boundary conditions, based on the WKB method.[40]

3. Results and discussion

The Yukawa-type screened potential (1) has been often utilized for the theoretical studies of atomic structures and collision dynamics in weakly coupled plasmas, but Chang and Fang[41]pointed out that the Debye-screening length should be longer than the average size of a free atom, as the mean radius of the 2p orbital in an N5+ion was estimated to be～0.4 a.u. (atomic units). Moreover, the DH potential supports only a finite number of bound states, which results in a migration of the bound state into the continuum whenDis less than the critical screening length.[42]According to the work of Xieet al.,[35]the critical screening length was determined to be 0.93 a.u. for the state 1s2p1P in a He-like C4+ion. Therefore, the range of the screening length was chosen as from 1 a.u.to ∞so that the 1s2p1P in Al11+and Fe24+ions is still the bound state. Due to weak electron correlations in the present cases of highly charged He-like ions, the singleconfiguration approximation is employed.

保證數(shù)據(jù)的可控性，數(shù)億量級(jí)的數(shù)據(jù)需要實(shí)時(shí)處理，對(duì)于現(xiàn)在的計(jì)算機(jī)硬件來(lái)說(shuō)是一個(gè)極大的挑戰(zhàn)，如何搭建一個(gè)高性能，高并發(fā)的實(shí)時(shí)數(shù)據(jù)分析系統(tǒng)，并在突發(fā)情況，自然災(zāi)害發(fā)生的情況下保證數(shù)據(jù)不會(huì)丟失，是智慧城市發(fā)展的第二個(gè)需要面臨的問(wèn)題。

3.1. Test of numerical methods

In order to verify the numerical methods and computing codes, the present EIE cross sections for the transition 1s→2p in an Ne9+ion are compared with the distortedwave calculations[23]forD=3 a.u. with the SDH potential in Fig. 1, and a good agreement is obtained. As a check of our numerical methods, which includes the Debye screening on both the target ions and the CB interaction with NSDH potential, our RDW results are compared with the 110-state CCC calculations[28]for the transition 1s21S→1s2p1P in He in Fig.2. It is noted that our purpose is not to discuss the accuracy of the RDW method for calculating the EIE cross section of a neutral He atom, but to characterize the plasmascreening effect. It is found that at the high incident electron energy of 800 eV,both of the present unscreened and screened(D=10 a.u.) calculations nearly agree with the CCC calculations,which supports the reliability of our theoretical method and codes used to include Debye-screening effects on EIE processes.

3.2. Transition energy and oscillator strength

Table 1 shows the Debye-screening effect on the transition energy and oscillator strength for the dipole-allowed transition 1s21S→1s2p1P in Al13+and Fe24+ions. Generally,when the screening lengthDis sufficiently greater than the radius of the ions, the influence of the Debye screening can be neglected. For the case ofD →∞, the present transition energy agrees well with the NIST value[43]with differences of 0.09% and 0.11% for Al13+and Fe24+ions, respectively.Moreover,the relative differences between our unscreened oscillator strengths and the NIST data[43]are determined to be 5.37%and 0.57%for Al13+and Fe24+ions,respectively. This indicates the reliability of the single-configuration approximation employed in the present calculations. As expected,when the screening lengthDis comparable to the radius of ions,the Debye-screening effect strongly affects the atomic properties.Thus,in Table 1 the screening effect makes significant reductions in the transition energy and oscillator strength whenDis smaller than about 5 a.u.,while the screening effect is very small for the case ofD ＞10 a.u.

Fig. 1. Comparison of the present unscreened and screened EIE cross sections for the transition 1s →2p of Ne9+ ion with those of the distorted-wave(DW)calculations by Pindzola et al.[23] The screened results for D=3 a.u.consider the bound-continuum interaction by the SDH potential.

Fig. 2. Comparison of the present unscreened and screened EIE cross sections for the transition 1s21S →1s2p 1P in He with those of the 110-state CCC calculations by Zammit et al.[28] The screened results for D=10 a.u.consider the plasma screening on both the target ion and interaction between the continuum electron and bound electron with the NSDH potential.

Table 1. Transition energy(eV)and oscillator strength for the transition 1s21S →1s2p 1P in Al13+ and Fe24+ ions for screening length D(a.u.).

3.3. Cross section

Next, we investigate the Debye-screening effect on the EIE cross sections for the transition 1s21S→1s2p1P in He-like Al11+and Fe24+ions. In Fig. 3, the variations of the EIE cross section at incident energies of 3000 eV and 15000 eV are illustrated for Al11+and Fe24+ions, respectively, with SDH and NSDH potentials. The variation Δσis defined as the difference between the cross section with the pure Coulomb interaction and that with the Debye-screened case,i.e., Δσ=σ(D →∞)-σ(D ＞0). As shown in Fig.3,since Δσis always positive, the Debye-screening effect decreases the EIE cross sections. Moreover, the Δσbecomes more obvious with decreasingD. It is noted that variations Δσof the SDH potential are dramatically greater than those of the NSDH potential,except that the Debye screening is weak in the large-Dcases. Thus, the SDH model can significantly overestimate the Debye-screening effects.

When the Debye-screening effect on the target ion is considered, the EIE cross sections decrease further, which is demonstrated in Fig. 3. However, for aboutD ＞5 a.u., the variations of the EIE cross section are almost the same as in the case in which the Debye screening on the target ion is not included. This indicates that the screening imposed on the target ion has negligible contributions to the EIE cross section when the plasma screening is too weak to affect the target ions for cases ofD ＞5 a.u.

We define the ratio

whereσion(D)andσee(D)represent the contributions from the Debye screening on the target ion and CB interaction,respectively, for reducing the EIE cross section. Here, the screened CB interaction is described by the NSDH model. In Fig. 4,since the ratioRis always less than one for all the screening length, the screening on the CB interaction is dominant for affecting (reducing) the EIE cross sections for the transition 1s21S→1s2p1P in Al11+and Fe24+ions at the incident energies of 3000 eV and 15000 eV, respectively. Moreover,as the screening length decreases, the influence of the Debye screening on the target ion increases(the ratioRincreases).As shown in Fig.4,in the case ofD=1 a.u.,R～0.24 and 0.18 for Al11+and Fe24+ions,respectively;hence,it indicates that the influences of screening on the target ion should be included for strong plasma screening. On the other hand,for the screening lengthD=10 a.u., the ratio is estimated to beR～0.05 and 0.03 for the Al11+and Fe24+cases, respectively, which implies that the reductions of the EIE cross sections are almost determined by the screening on the CB interaction. This also results in the fact that the influences of screening on the target ion are negligible for sufficiently weak screening strength,which will be demonstrated further in Fig.5.

In Fig.5,the screened(D=10 a.u.) and unscreened EIE cross sections are shown for the transition 1s21S→1s2p1P in He-like Al11+ions at different incident energies. In the process of the electron-ion collision,the cross sections at the threshold are finite since the long-range Coulomb interaction has contributions to the large separation between the continuum electron and target ion. Thus, the EIE cross sections at the threshold should be more significantly reduced due to the short-range character of the DH potential with the exponential screening factor [see Eq. (1)]. However, in Fig. 5 the reductions of the EIE cross section due to the SDH potential are basically the same for incident energy ranging from the threshold to high energy. For the NSDH potential, it clearly demonstrates that the reductions of the cross section at the threshold are indeed larger than those at high energies. Hence,the NSDH potential could reasonably describe the Debyescreening effects on the EIE process in Debye plasmas, as expected since this potential results in a reasonable changing behavior of the cross section with incident energy.

Fig.3. Variation of the EIE cross section,Δσ =σ(D →∞)-σ(D ＞0),for the transition 1s21S →1s2p 1P in Al11+(a)and Fe24+(b)ions,respectively.SHD and NSDH[see Eqs.(2)-(3)and text]indicate that the CB interaction is described by the SDH and NSDH potentials, respectively. NSDH with screened target ion means that Debye screening on the target ion is included as well.

Fig.4. Ratio R[see Eq.(9)and text]for the transition 1s21S →1s2p 1P in Al11+ and Fe24+ ions,respectively.

Moreover, as mentioned above, in Fig. 3 the ratioRat the incident electron energy of 3000 eV is shown for a Helike Al11+ion, which is estimated to be～0.05 for the case ofD= 10 a.u., indicating that the screening on the target ion results in negligible reductions of the cross sections compared to the screening on the CB interaction. As illustrated in Fig.5,this is also true for the energy range from threshold to 10000 eV, since there are no significant differences between the cross sections with and without inclusions of the screening on the target ion forD=10 a.u.

In Fig. 6, the EIE cross sections are displayed for the dipole-allowed transition 1s21S→1s2p1P in He-like Al11+and Fe24+ions withD=∞,5,2,and 1 a.u.For the unscreened EIE cross sections, our results agree with the available theoretical calculations by Kiefferet al.,[44]Hakelet al.[45]and Zhanget al.[33]The EIE cross sections are decreased due to Debye-screening effects,which become distinct with decreasingD.For example,at the incident energy of 3400 eV,the EIE cross section is reduced by about 13.4%, 30.2%, and 50.9%forD=5, 2, and 1 a.u., respectively, for He-like Al11+ions.Moreover, the screening effect decreases the cross section by about 3.1%,7.8%,and 15.6%forD=5,2,and 1 a.u.,respectively,near the threshold of 7500 eV for He-like Fe24+ions.

Fig. 5. Unscreened (black line) and screened (D=10 a.u.) EIE cross sections for the transition 1s21S →1s2p 1P in Al11+ ions. SHD and NSDH[see Eqs. (2)-(3) and text] indicate that the CB interaction is described by the SDH and NSDH potentials,respectively. NSDH with screened target ion means that Debye screening on the target ion is included as well.

Fig. 6. EIE cross sections for the transition 1s21S →1s2p 1P in Al11+ (a)and Fe24+ (b)ions with the screening length D=∞,5,2,and 1 a.u. Present unscreened cross sections are compared with the available RDW results by Kieffer et al.,[44] Zhang et al.,[33] and Hakel et al.[45]

4. Conclusion

We have studied the EIE processes for the dipole-allowed transition 1s21S→1s2p1P in He-like Al11+and Fe24+ions embedded in weakly coupled(Debye)plasmas using the MCDF and RDW methods with the DH potential. In order to account for the Debye-screening effect on the CB interaction in the EIE process,two model potentials(SDH and NSDH)are employed,which result in reductions of the EIE cross section.However,the reduced magnitude of the EIE cross sections for the SDH potential is far greater than that for the NSDH potential. Therefore,the well-suited SDH potential obviously overestimates the Debye-screening effects on the EIE processes for multielectron He-like ions. Moreover, the Debye-screening effect on the target ion also decreases the EIE cross sections.Our results reveal that the influence of the screening on the CB interaction is dominant for the reduction of the EIE cross section,while the screening on the target ion could contribute to reducing the cross sections further as the screening strength increases. However, for cases of sufficiently weak screening strength, the contributions of the screening on the target ion are negligible.

Finally, the EIE cross sections are presented with the inclusion of the Debye-screening effect on both the target ion and the CB interaction. It is found that the Debyescreening effects obviously reduce the EIE cross sections,which clearly indicates that atomic processes are perturbed in a non-negligible way in plasmas. Therefore,it is highly desirable to investigate environmental effects on plasma kinetics,as well as on absorption and emission spectra.

Acknowledgments

Project supported by the Science Challenge Project(Grant No. TZ2016001), the National Key Research and Development Program of China(Grants Nos.2017YFA0403200 and 2017YFA0402300),the Funds for Innovative Fundamental Research Group Project of Gansu Province, China (Grant No.20JR5RA541),and the National Natural Science Foundation of China(Grants Nos.11774037 and 12064041).