Towards a spectroscopically accurate set of potentials for heavy hydride laser cooling candidates: effective core potential calculations of BaH

BaH (and its isotopomers) is an attractive molecular candidate for laser cooling to ultracold temperatures and a potential precursor for the production of ultracold gases of hydrogen and deuterium. The theoretical challenge is to simulate the laser cooling cycle as reliably as possible and this paper addresses the generation of a highly accurate ab initio $^{2}\Sigma^+$ potential for such studies. The performance of various basis sets within the multi-reference configuration-interaction (MRCI) approximation with the Davidson correction (MRCI+Q)is tested and taken to the Complete Basis Set (CBS) limit. It is shown that the calculated molecular constants using a 46 electron Effective Core-Potential (ECP) and even-tempered augmented polarized core-valence basis sets (aug-pCV$n$Z-PP, n= 4 and 5) but only including three active electrons in the MRCI calculation are in excellent agreement with the available experimental values. The predicted dissociation energy De for the X$^2\Sigma^+$ state (extrapolated to the CBS limit) is 16895.12 cm$^{-1}$ (2.094 eV), which agrees within 0.1$\%$ of a revised experimental value of <16910.6 cm$^{-1}$, while the calculated re is within 0.03 pm of the experimental result.

Electron-impact excitation of W 44+ and W 45+

Plans to employ tungsten in the divertor region of the International Thermonuclear Experimental Reactor require radiative and collisional data for modelling x-ray emissions of highly ionized stages of tungsten. In an earlier paper, we reported on the results of fully relativistic R -matrix calculations for W 46+ that included the effects of radiation damping on the resonance contributions. In this paper, we present the results of similar fully relativistic, radiatively damped R -matrix calculations for W 44+ and W 45+ . Radiation damping is found to be small for W 45+ , but is appreciable for many of the excitations from the ground and metastable levels of W 44+ . Rates from the present calculations will be combined with those from the calculations for W 46+ and employed for collisional-radiative modelling for these ions.

Relativistic radiatively damped R-matrix calculation of the electron-impact excitation of W46+

The current design plans for the International Thermonuclear\nExperimental Reactor ( ITER) call for tungsten to be employed for\ncertain plasma facing components in the divertor region. Thus, accurate\natomic collision data are needed for emission modelling of tungsten.\nElectron-impact excitation and radiative rates are of particular\nimportance for Ni-like W, since this ion emits some of the most intense\nspectral lines of all ionization stages. We report on a fully\nrelativistic 115-level R-matrix calculations of W46+, which includes the\neffects of radiation damping. Although radiation damping is very\nimportant in most highly ionized species, its effects are reduced in\nthis case because of the closed-shell Ni-like ground state. The rates\nfrom these relativistic atomic calculations will be employed for\ncollisional-radiative modelling of this ion.

Electron-impact excitation of Fe 4+ : an intermediate-coupling R -matrix calculation

For a number of years, there has been a major effort to calculate electron-impact excitation data for every ion stage of iron embodied by the ongoing efforts of the IRON project by Hummer et al (1993 Astron. Astrophys. 279 298). Due to the complexity of the targets, calculations for the lower stages of ionization have been limited to either intermediate-coupling calculations within the ground configurations or LS -coupling calculations of the ground and excited configurations. However, accurate excitation data between individual levels within both the ground and excited configurations of the low charge-state ions are urgently required for applications to both astrophysical and laboratory plasmas. Here we report on the results of the first intermediate-coupling R -matrix calculation of electron-impact excitation for Fe 4+ for which the close-coupling (CC) expansion includes not only those levels of the 3d 4 ground configuration, but also the levels of the 3d 3 4s, 3d 3 4p, 3d 3 4d and 3d 2 4s 2 excited configurations. With 359 levels in the CC expansion and over 2400 scattering channels for many of the J Π partial waves, this represents the largest electron–ion scattering calculation to date and it was performed on massively parallel computers using a recently developed set of relativistic parallel R -matrix programs.

Rydberg structure associated with core-excited autoionizing states of the LiH radical

We have investigated inner-shell excitation of the LiH + molecular ion by electron impact within several different collision models to delineate Rydberg autoionizing resonance structure associated with the LiH + (1σ2σ 2 2 Σ + ) core-excited threshold. The minimal representation requires only the retention of the 1σ and 2σ molecular orbitals, in which the core-excited state involves the promotion of a single electron into the 2σ orbital. This model is extended to include two further representations, in which both the 3σ and 4σ orbitals obtained from a self-consistent field calculation improve target representation, correlation and support additional autoionization channels. This affects the autoionization widths and to a lesser degree the positions of the LiH (1σ2σ 2 n s, n p 1,3 Σ + ) resonance series. Comparing our work with calculations on the counterpart atomic Be system assists in the assignment of the core-excited molecular resonance states. The results from our investigation provide helpful insights into the study of inner-shell transitions produced by electron or photon impact in more complex diatomic molecules.

Electron-impact ionization of diatomic molecules using a configuration-average distorted-wave method

Electron-impact ionization cross sections for diatomic molecules are calculated in a configuration-average distorted-wave method. Core bound orbitals for the molecular ion are calculated using a single-configuration self-consistent-field method based on a linear combination of Slater-type orbitals. The core bound orbitals are then transformed onto a two-dimensional (r,θ) numerical lattice from which a Hartree potential with local exchange is constructed. The single-particle Schrödinger equation is then solved for the valence bound orbital and continuum distorted-wave orbitals with S-matrix boundary conditions. Total cross section results for H2 and N2 are compared with those from semiempirical calculations and experimental measurements.

Electron-Impact Uncertainty Analysis and its Impact on Certain Temperature Diagnostics

In this study we calculate the electron-impact uncertainties in atomic data for direct ionization and recombination and investigate the role of these uncertainties on spectral diagnostics. We outline a systematic approach to assigning meaningful uncertainties that vary with electron temperature. Once these uncertainty parameters have been evaluated, we can then calculate the uncertainties on key diagnostics through a Monte Carlo routine, using the Astrophysical Emission Code (APEC) [Smith et al. 2001]. We incorporate these uncertainties into well known temperature diagnostics, such as the Lyman alpha versus resonance line ratio and the G ratio. We compare these calculations to a study performed by [Testa et al. 2004], where significant discrepancies in the two diagnostic ratios were observed. We conclude that while the atomic physics uncertainties play a noticeable role in the discrepancies observed by Testa, they do not explain all of them. This indicates that there is another physical process occurring in the system that is not being taken into account. This work is supported in part by the National Science Foundation REU and Department of Defense ASSURE programs under NSF Grant no. 1262851 and by the Smithsonian Institution.

Non-equilibrium modeling of the FE XVII 3C/3D Line Ratio in an Intense X-ray Free-electron Laser Excited Plasma

Recent measurements using an X-ray Free Electron Laser (XFEL) and an Electron Beam Ion Trap at the Linac Coherent Light Source facility highlighted large discrepancies between the observed and theoretical values for the Fe XVII 3C/3D line intensity ratio. This result raised the question of whether the theoretical oscillator strengths may be significantly in error, due to insufficiencies in the atomic structure calculations. We present time-dependent spectral modeling of this experiment and show that non-equilibrium effects can dramatically reduce the predicted 3C/3D line intensity ratio, compared with that obtained by simply taking the ratio of oscillator strengths. Once these non-equilibrium effects are accounted for, the measured line intensity ratio can be used to determine a revised value for the 3C/3D oscillator strength ratio, giving a range from 3.0 to 3.5. We also provide a framework to narrow this range further, if more precise information about the pulse parameters can be determined. We discuss the implications of the new results for the use of Fe XVII spectral features as astrophysical diagnostics and investigate the importance of time-dependent effects in interpreting XFEL-excited plasmas.

Electron-impact excitation of Li to high principal quantum numbers

The time-dependent close-coupling method is used to calculate electron-impact excitation cross sections for the Li(2s)--{\textgreater}Li(nl) and Li(2p)--{\textgreater}Li(nl) transitions at incident energies just above the ionization threshold. The implementation of the time-dependent close-coupling method on a nonuniform lattice allows the study of continuum-coupling effects in excitations to high principal quantum number, i.e., n{\textless}=10. Good agreement is found with R-matrix with pseudostates calculations, which also include continuum-coupling effects, for excitations to low principal quantum number, i.e., n{\textless}=4. Poor agreement is found with standard distorted-wave calculations for excitations to all principal quantum numbers, with differences still at the 50% level for n=10. We are able to give guidance as to the accuracy expected in the n3 extrapolation of nonperturbative close-coupling calculations of low n cross sections and rate coefficients.

Electron-impact ionization of Al2+

We report on a nonperturbative R-matrix with PseudoStates (RMPS) calculation for the electron-impact ionization cross section of the ground state of Al2+. We include both the direct ionization of the 3s and 2p subshells and the indirect ionization from the 2p subshell. This calculation, thus, includes extra decay channels for the indirect-ionization process not included in previous RMPS calculations. This lowers the total-ionization cross section, resulting in closer agreement with the most recent experimental measurements. This calculation also shows better agreement with the position and height of the resonant-excitation double autoionization features seen in the experiment.

Electron-impact excitation of Ar+: an improved determination of Ar impurity influx in tokamaks

Electron-impact scattering data for argon and its ions continue to be of interest in studies of magnetically confined plasmas. In an earlier paper, Griffin et al (1997 J. Phys. B: At. Mol. Opt. Phys. 30 3543) employed the results of 28-term and 40-term R-matrix calculations of electron-impact excitation in Ar+ to carry out a collisional-radiative modelling study of the impurity influx of argon in tokamaks. We have now completed a 452-term R-matrix with pseudo-states (RMPS) calculation of electron-impact excitation for Ar+ in order to provide more accurate excitation data; using these improved data, we have repeated the modelling studies presented in the earlier paper. We compare our excitation data, as well as the results of the collisional radiative calculations, with those arising from the 40-term R-matrix calculation and find significant differences.

Vibrational excitation of the N2+ first negative (0,0), (1,0) and (2,0) bands by electron impact: a theoretical study using the R-matrix approach

Ab initio cross section calculations for vibronic excitation using the R -matrix approach have been performed on the N 2 + molecular ion complex. A three-state close-coupling expansion is used where the electronic target states; X 2 g + , A 2 u and B 2 u + of the molecular cation are represented by a valence configuration-interaction approximation. A non-adiabatic approximation is invoked to study vibronic excitation for the first three negative bands, (0,0), (1,0) and (2,0) of the X-B transition (B 2 u + v ´ X 2 g + v ´´ ) of N 2 + . Fixed-nuclei and non-adiabatic cross section results are compared with the available experimental data for the (0,0) band and the breakdown of the adiabatic fixed-nuclei approximation is clearly evident for the vibronic excitation of the (1,0) and (2,0) bands in this molecular ion complex.

Electron-impact single ionization of Mg and Al+

Theory and experiment are compared for the electron-impact single ionization of Mg and Al+. Nonpertur- bative R matrix with pseudostates RMPS and time-dependent close-coupling TDCC calculations have been carried out that exhibit large reductions from perturbative distorted-wave results of 38% for Mg and 20% for Al+. Experimental single-ionization data available for Mg and Al+ are in reasonable accord with distorted-wave data and lie substantially above the new theoretical results. Rate coefficients, necessary for the collisional- radiative modeling of Mg and Al plasmas were generated from the RMPS ionization cross sections. In the collisional-ionization region near the ionization threshold, the resulting rates were found to be up to two times lower for Mg and three times lower for Al+ than the rates generated from experimental data.

Electron-impact excitation and ionization of W 3+ for the determination of tungsten influx in a fusion plasma

Tungsten will be employed as a plasma facing material in the ITER fusion reactor under construction in Cadarache, France; therefore, there is a significant need for accurate electron-impact excitation and ionization data for the ions of tungsten. We report on the results of extensive calculations of ionization and excitation for W 3+ that are intended to provide the atomic data needed for the determination of impurity influx diagnostics of tungsten in several existing tokamak reactors. The electron-impact excitation rate coefficients for this study were determined using the relativistic R -matrix method. The contribution to direct electron-impact ionization was determined using the distorted-wave approximation, the accuracy of which was verified by an R -matrix with pseudo states calculation. Contributions to total ionization from excitation autoionization were also generated from the relativistic R -matrix method. These results were then employed to calculate values of ionization per emitted photon, or SXB ratios, for four carefully selected spectral lines; these data will allow the determination of impurity influx from tungsten facing surfaces. For the range of densities of importance in the edge region of a tokamak reactor, these SXB ratios are found to be nearly independent of electron density but vary significantly with electron temperature.

Electron-impact excitation of beryllium and its ions

Inelastic electron scattering from light atomic species is of fundamental importance and has significant applications in fusion-plasma modeling. Therefore, it is of interest to apply advanced nonperturbative, close-coupling methods to the determination of electron-impact excitation for these atoms. Here we present the results of R matrix with pseudostate (RMPS) calculations of electron-impact excitation cross sections through the n=4 terms in Be, Be+, Be2+, and Be3+. In order to determine the effects of coupling of the bound states to the target continuum in these species, we compare the RMPS results with those from standard R-matrix calculations. In addition, we have performed time-dependent close-coupling calculations for excitation from the ground and the metastable terms of Be+ and the metastable term of Be3+. In general, these results are found to agree with those from our RMPS calculations. The full set of data resulting from this work is now available on the Oak Ridge National Laboratory Controlled Fusion Atomic Data Center web site, and will be employed for collisional-radiative modeling of Be in magnetically confined plasmas.