We present an efficient numerical method to solve the time-dependent Schrödinger equation in the single-active electron picture for atoms interacting with intense optical laser fields. Our approach is based on a non-uniform radial grid with smoothly increasing steps for the electron distance from the residual ion. We study the accuracy and efficiency of the method, as well as its applicability to investigate strong-field ionization phenomena, the process of high-order harmonic generation, and the dynamics of highly excited Rydberg states.
Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
-
Free, publicly-accessible full text available July 1, 2025
-
Abstract The active-particle number density is a key parameter for plasma material processing, space propulsion, and plasma-assisted combustion. The traditional actinometry method focuses on measuring the density of the atoms in the ground state, but there is a lack of an effective optical emission spectroscopy method to measure intra-shell excited-state densities. The latter atoms have chemical selectivity and higher energy, and they can easily change the material morphology as well as the ionization and combustion paths. In this work, we present a novel state-resolved actinometry (SRA) method, supported by a krypton line-ratio method for the electron temperature and density, to measure the number densities of nitrogen atoms in the ground and intra-shell excited states. The SRA method is based on a collisional-radiative model, considering the kinetics of atomic nitrogen and krypton including their excited states. The densities measured by our method are compared with those obtained from a dissociative model in a miniature electron cyclotron resonance (ECR) plasma source. Furthermore, the saturation effect, in which the electron density remains constant due to the microwave propagation in an ECR plasma once the power reaches a certain value, is used to verify the electron density measured by the line-ratio method. An ionization balance model is also presented to examine the measured electron temperature. All the values obtained with the different methods are in good agreement with each other, and hence a set of verified rate coefficient data used in our method can be provided. A novel concept, the ‘excited-state system’, is presented to quickly build an optical diagnostic method based on the analysis of quantum number propensity and selection rules.
Free, publicly-accessible full text available May 1, 2025 -
Abstract The ionization fraction is a key figure of merit for optimizing the performance of plasma device. This work presents an optical emission spectroscopy (OES) method to determine the ionization fraction in low-temperature xenon plasma. The emission line-ratio of xenon ionic and atomic 6
p –6s transitions is used in this method. A comprehensive collisional-radiative model developed in our previous work is employed to describe the relationship between the line-ratios and the plasma parameters. It is found that some special line-ratios have a sensitive relationship to the ionization fraction, e.g. the ratio of the 460.30 nm line and 828.01 nm lines. These line-ratios are selected for the diagnostic method. The method is demonstrated in a magnetized discharge chamber. The axially-resolved emission spectra of the ionization chamber are measured, and from those the ionization fraction along the chamber axis is determined via the OES method. The axially-resolved ionization fraction is found to be dependent on the magnetic field and agrees well with those obtained from a Langmuir probe. In the experiment, the probe is overheated under some conditions, possibly due to the bombardment by energetic particles. In this case, no results can be obtained from the probe, while the OES method can still obtain reasonable results. Combined with optical tomography and spectral imaging technology, the OES method can also provide the spatial distribution of the ionization fraction, which is needed for revealing the discharge mechanisms of plasma devices. -
Within the framework of the improved quantitative rescattering (QRS) model, we simulate the correlated two-electron momentum distributions (CMDs) for nonsequential double ionization (NSDI) of Ar by elliptically polarized laser pulses with a wavelength of 788 nm at an intensity of 0.7 × 1014W/cm2for the ellipticities ranging from 0 to 0.3. Only the CMDs for recollision excitation with subsequent ionization (RESI) are calculated and the contribution from recollision direct ionization is neglected. According to the QRS model, the CMD for RESI can be factorized as a product of the parallel momentum distribution (PMD) for the first released electron after recollision and the PMD for the second electron ionized from an excited state of the parent ion. The PMD for the first electron is obtained from the laser-free differential cross sections for electron impact excitation of Ar+calculated using state-of-the-art many-electron
R -matrix theory while that for the second electron is evaluated by solving the time-dependent Schrödinger equation. The results show that the CMDs for all the ellipticities considered here exhibit distinct anticorrelated back-to-back emission of the electrons along the major polarization direction, and the anticorrelation is more pronounced with increasing ellipticity. It is found that anticorrelation is attributed to the pattern of the PMD for the second electron ionized from the excited state that, in turn, is caused by the delayed recollision time with respect to the instant of the external field crossing. Our work shows that both the ionization potential of the excited parent ion and the laser intensity play important roles in the process. -
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.more » « less
-
We reinvestigate a key process in electron-atom collision physics, the elastic scattering of electrons from helium atoms. Specifically, results from a special-purpose relativistic polarized-orbital method, which is designed to treat elastic scattering only, are compared with those from a very extensive, fully ab initio, general-purpose B-spline R-matrix (close-coupling) code.more » « less