skip to main content

This content will become publicly available on March 1, 2023

Title: Generalizations of the R-Matrix Method to the Treatment of the Interaction of Short-Pulse Electromagnetic Radiation with Atoms
Since its initial development in the 1970s by Phil Burke and his collaborators, the R-matrix theory and associated computer codes have become the method of choice for the calculation of accurate data for general electron–atom/ion/molecule collision and photoionization processes. The use of a non-orthogonal set of orbitals based on B-splines, now called the B-spline R-matrix (BSR) approach, was pioneered by Zatsarinny. It has considerably extended the flexibility of the approach and improved particularly the treatment of complex many-electron atomic and ionic targets, for which accurate data are needed in many modelling applications for processes involving low-temperature plasmas. Both the original R-matrix approach and the BSR method have been extended to the interaction of short, intense electromagnetic (EM) radiation with atoms and molecules. Here, we provide an overview of the theoretical tools that were required to facilitate the extension of the theory to the time domain. As an example of a practical application, we show results for two-photon ionization of argon by intense short-pulse extreme ultraviolet radiation.
; ;
Award ID(s):
1834740 2110023 1803844
Publication Date:
Journal Name:
Sponsoring Org:
National Science Foundation
More Like this
  1. Analytic fits to the recommended electron-impact excitation and ionization cross sections for Be I are presented. The lowest 19 terms of configurations 2snl (n≤4) and 2p2 terms below the first ionization limit are considered. The fits are based on the accurate calculations with the convergent close coupling (CCC) method as well as the B-spline R-matrix (BSR) approach. The fitted cross sections provide rate coefficients that are believed to approximate the original data within 10% with very few exceptions. The oscillator strengths for the dipole-allowed transitions between all the considered states are calculated with the relativistic multi-configuration Dirac–Hartree–Fock (MCDHF) approach andmore »compared with the CCC and BSR results. This comparison shows a very good agreement except for a handful of cases with likely strong cancellations.« less
  2. Context. Older models of Galactic chemical evolution (GCE) predict [K/Fe] ratios as much as 1 dex lower than those inferred from stellar observations. Abundances of potassium are mainly based on analyses of the 7698 Å resonance line, and the discrepancy between GCE models and observations is in part caused by the assumption of local thermodynamic equilibrium (LTE) in spectroscopic analyses. Aims. We study the statistical equilibrium of K  I , focusing on the non-LTE effects on the 7698 Å line. We aim to determine how non-LTE abundances of potassium can improve the analysis of its chemical evolution, and help to constrain themore »yields of GCE models. Methods. We construct a new model K  I atom that employs the most up-to-date atomic data. In particular, we calculate and present inelastic e+K collisional excitation cross-sections from the convergent close-coupling (CCC) and the B -Spline R -matrix (BSR) methods, and H+K collisions from the two-electron model (LCAO). We constructed a fine, extended grid of non-LTE abundance corrections based on 1D MARCS models that span 4000 < T eff ∕K < 8000, 0.50 < log g < 5.00, − 5.00 < [Fe/H] < + 0.50, and applied the corrections to potassium abundances extracted from the literature. Results. In concordance with previous studies, we find severe non-LTE effects in the 7698 Å line. The line is stronger in non-LTE and the abundance corrections can reach approximately − 0.7 dex for solar-metallicity stars such as Procyon. We determine potassium abundances in six benchmark stars, and obtain consistent results from different optical lines. We explore the effects of atmospheric inhomogeneity by computing for the first time a full 3D non-LTE stellar spectrum of K  I lines for a test star. We find that 3D modeling is necessary to predict a correct shape of the resonance 7698 Å line, but the line strength is similar to that found in 1D non-LTE. Conclusions. Our non-LTE abundance corrections reduce the scatter and change the cosmic trends of literature potassium abundances. In the regime [Fe/H] ≲−1.0 the non-LTE abundances show a good agreement with the GCE model with yields from rotating massive stars. The reduced scatter of the non-LTE corrected abundances of a sample of solar twins shows that line-by-line differential analysis techniques cannot fully compensate for systematic LTE modelling errors; the scatter introduced by such errors introduces a spurious dispersion to K evolution.« less
  3. Abstract This work reports large-scale calculations of electron excitation effective collision strengths and transition rates for a wide range of Sc ii spectral lines for astrophysical analysis and modeling. The present results are important for reliable abundance determinations in various astrophysical objects, including metal-poor stars, H ii regions, and gaseous nebulae. Accurate descriptions of the target wave functions and adequate accounts of the various interactions between the target levels are of primary importance for calculations of collision and radiative parameters. The target wave functions have been determined by a combination of the multiconfiguration Hartree–Fock and B-spline box-based close-coupling methods, togethermore »with the nonorthogonal orbitals technique. The calculations of the collision strengths have been performed using the B-spline Breit–Pauli R-matrix method. The close-coupling expansion includes 145 fine-structure levels of Sc ii belonging to the terms of the 3 p 6 3 d 2 , 3 p 6 3 d 4 l ( l = 0–3), 3 p 6 3 d 5 l ( l = 0–3), 3 p 6 3 d 6 s , 3 p 6 4 s 2 , 3 p 6 4 s 4 l ( l = 0–3), 3 p 6 4 s 5 l ( l = 0–1), and 3 p 6 4 p 2 configurations. The effective collision strengths are reported as a function of electron temperature in the range from 10 3 to 10 5 K. The collision and radiative rates are reported for all of the possible transitions between the 145 fine-structure levels. Striking discrepancies exist with the previous R-matrix calculations of the effective collision strengths for the majority of the transitions, indicating possible systematic errors in these calculations. Thus, there is a need for accurate calculations to reduce the uncertainties in the atomic data. The likely uncertainties in our effective collision strengths and radiative parameters have been assessed by means of comparisons with other collision calculations and available experimental radiative parameters.« less
  4. Partial and total photoionization cross sections of iron-peak elements are important for the determination of abundances in late-type stars and nebular objects. We have investigated photoionization of neutral chromium from the ground and excited states in the low energy region from the first ionization threshold at 6.77 eV to 30 eV. Accurate descriptions of the initial bound states of Cr I and the final residual Cr II ionic states have been obtained in the multiconfiguration Hartree-Fock method together with adjustable configuration expansions and term-dependent non-orthogonal orbitals. The B-spline R-matrix method has been used for the calculation of photoionization cross sections.more »The 194 LS final ionic states of Cr II 3d44s, 3d34s2, 3d5, 3d44p, and 3d34s4p principal configurations have been included in the close-coupling expansion. The inclusion of all terms of these configurations has significant impact on the near-threshold resonance structures as well as on the nonresonant background cross sections. Total photoionization cross sections from the ground 3d54sa7S and excited 3d54sa5S, 3d44s2a5D, 3d54pz5P, and 3d44s4py5P states of Cr I have been compared with other available R-matrix calculation to estimate the likely uncertainties in photoionization cross sections. We analyzed the partial photoionization cross sections for leaving the residual ion in various states to identify the important scattering channels, and noted that 3d electron ionization channel becomes dominant at higher energies.« less
  5. We have applied the full-relativistic Dirac B-Spline R-matrix method to obtain cross sections for electron scattering from ytterbium atoms. The results are compared with those obtained from a semi-relativistic (Breit-Pauli) model-potential approach and the few available experimental data.