ThunderBoltz DSMC Simulation of Pulsed Discharges with Disorder-Induced Heating
Explore the source record for details and available documents.
Park, Ryan↗
Engineering topics
Zammit, Mark Christian
Publications and source records attributed to Zammit, Mark Christian.
Ground and Excited Electronic Structure Analysis of CH Diatomic Molecule and Its Astrophysical Importance [Slides]
Ground and excited PECs and DMCs of CH are reported. The aug-cc-pV6Z basis set is a good accuracy-efficiency compromise to study C-H system. CASSCF results are in reasonable agreement with the experiment. High-level MRCI analysis are vital for reaching for the experimental accuracy.
79 ASTRONOMY AND ASTROPHYSICS↗
Ultraviolet H 2 luminescence in molecular clouds induced by cosmic rays
Context. Galactic cosmic rays (CRs) play a crucial role in ionisation, dissociation, and excitation processes within dense cloud regions where UV radiation is absorbed by dust grains and gas species. CRs regulate the abundance of ions and radicals, leading to the formation of more and more complex molecular species, and determine the charge distribution on dust grains. A quantitative analysis of these effects is essential for understanding the dynamical and chemical evolution of star-forming regions. Aims. The CR-induced photon flux has a significant impact on the evolution of the dense molecular medium in its gas and dust components. This study evaluates the flux of UV photons generated by CRs to calculate the photon-induced dissociation and ionisation rates of a vast number of atomic and molecular species, as well as the integrated UV photon flux. Methods. To achieve these goals, we took advantage of recent developments in the determination of the spectra of secondary electrons, in the calculation of state-resolved excitation cross sections of H 2 by electron impact, and of photodissociation and photoionisation cross sections. Results. We calculated the H 2 level population of each rovibrational level of the X, B, C, B', D, B", D', and a states. We then computed the UV photon spectrum of H 2 in its line and continuum components between 72 and 700 nm, with unprecedented accuracy, as a function of the CR spectrum incident on a molecular cloud, the H 2 column density, the isomeric H 2 composition, and the dust properties. The resulting photodissociation and photoionisation rates are, on average, lower than previous determinations by a factor of about 2, with deviations of up to a factor of 5 for the photodissociation of species such as AlH, C 2 H 2 , C 2 H 3 , C 3 H 3 , LiH, N 2 , NaCl, NaH, O$^+_2$, S 2 , SiH, l-C 4 , and l-C 5 H. A special focus is given to the photoionisation rates of H 2 , HF, and N 2 , as well as to the photodissociation of H 2 , which we find to be orders of magnitude higher than previous estimates. We give parameterisations for both the photorates and the integrated UV photon flux as a function of the CR ionisation rate, which implicitly depends on the H 2 column density, as well as the dust properties.
79 ASTRONOMY AND ASTROPHYSICS↗
Photoionization from the Ground and Excited Vibrational States of H$^+_2$ and Its Deuterated Isotopologues
Photoionization cross sections and rate coefficients have been calculated for all bound vibrational levels of the 1sσ g state of H$^+_2$, HD+, and D$^+_2$. The Born–Oppenheimer approximation is employed in our calculation of vibrationally resolved photoionization cross sections. Vibrationally resolved and local thermal equilibrium photoionization rate coefficients are presented for photon temperatures less than 50,000 K and are found to be several orders of magnitude larger than previous results in the literature. Analytic fits for the vibrationally resolved and local thermal equilibrium photoionization rate coefficients are provided. Near-threshold oscillations in the vibrationally resolved photoionization are observed. A benchmark set of photoionization cross sections are presented. Fixed-nuclei photoionization cross sections are calculated using two-center true continuum wave functions and are verified by comparison with previous calculations and are found to be in excellent agreement in all cases. Data files for our set of benchmark cross sections, rate coefficients, and fitting parameters for H$^+_2$, HD + , and D$^+_2$ are available on Zenodo under an open-source Creative Commons Attribution license at doi:10.5281/zenodo.8304061.
74 ATOMIC AND MOLECULAR PHYSICS↗
Elastic scattering and rotational excitation of H 2 by electron impact: Convergent close-coupling calculations
We apply the adiabatic-nuclei molecular convergent close-coupling (MCCC) method to the study of elastic scattering and rotational excitation of H 2 by 0.01–20-eV electrons. Integral cross sections are presented for all rotational transitions with |ΔN| ≤ 2 and N i = 0–31 within the v = 0 vibrational level and differential cross sections for a selection of transitions. Agreement with the available measurements and previous calculations is mostly excellent, depending on the transition and incident energy. Finally, we suggest possible reasons, and argue for the accuracy of the MCCC data.
74 ATOMIC AND MOLECULAR PHYSICS↗
Complete collision data set for electrons scattering on molecular hydrogen and its isotopologues: IV. Vibrationally-resolved ionization of the ground and excited electronic states
Here, we present a comprehensive set of vibrationally-resolved cross sections for electron-impact ionization of molecular hydrogen and its isotopologues (H 2 , D 2 , T 2 , HD, HT, and DT) in both the ground and excited electronic states. We apply the adiabatic-nuclei molecular convergent close-coupling (MCCC) method to calculate cross sections from threshold to 1000 eV for ionization of the ground and excited vibrational levels of the X 1 Σ$^{+}_{g}$, B 1 Σ$^{+}_{u}$, C 1 π u , EF 1 Σ$^{+}_{g}$, a 3 Σ$^{+}_{g}$, and c 3 π u electronic states, representing all states with united-atoms-limit principle quantum number n=1–2. The cross sections are presented in graphical form and provided as both numerical values and analytic fit functions in supplementary data files. The data can also be downloaded from the MCCC database at mccc-db.org.
74 ATOMIC AND MOLECULAR PHYSICS↗
Convergent close-coupling calculations of electron scattering on HeH +
Here, we use the molecular convergent close-coupling (MCCC) method to perform calculations of 10–1000 eV electron scattering on the ground state of HeH + . Cross sections are presented for excitation of the n = 2 – 3 singlet and triplet states (where n is the united-atoms-limit principle quantum number), as well as ionization. We also present cross sections for He + and H + ion production following dissociative excitation and ionization. The He + production cross section is compared with the measurements of Lecointre, Jureta, Urbain, and Defrance [J. Lecointre, J. J. Jureta, X. Urbain, and P. Defrance, J. Phys. B At. Mol. Opt. Phys. 47, 015203 (2014)]. We find that the MCCC results are up to 30% higher than experiment.
74 ATOMIC AND MOLECULAR PHYSICS↗
Complete collision data set for electrons scattering on molecular hydrogen and its isotopologues: III. Vibrational excitation via electronic excitation and radiative decay
We present state-resolved cross sections for vibrational excitation via electronic excitation followed by radiative decay (ERD), for electrons scattering on all bound vibrational levels of the ground electronic state ($X\hspace{1.5 mm} ^1Σ_g^+$) of molecular hydrogen and its isotopologues (H 2 , HD, HT, D 2 , DT and T 2 ). We consider excitation of the singlet $\textit{n}$ = 2–3 states (where n refers to the united-atoms-limit principle quantum number) and account for all possible decay pathways back to the bound vibrational levels of the ground electronic state ($X\hspace{1.5 mm} ^1Σ_g^+$) to produce an estimate for the ERD cross sections. A selection of the results are presented in graphical form and a full data set is provided as both numerical values and analytic fit functions in supplementary data files. The uncertainty in the cross sections is estimated to be 11%, except for scattering on the highest bound vibrational level of HD, HT, D 2 , and DT, where the uncertainty is 21%. The data can be downloaded from the MCCC database at mccc-db.org.
74 ATOMIC AND MOLECULAR PHYSICS↗
Anisotropic angular scattering models of elastic electron-neutral collisions for Monte Carlo plasma simulations
Many laboratory and industrial plasma applications require accurate modeling techniques to understand the interplay between microscopic and macroscopic processes. A prime example of this interplay is how particle and Monte Carlo (MC) simulation codes describe angular scattering of electrons following elastic scattering events. The forward peaked nature of high energy electron elastic scattering is relatively trivial to accurately describe in plasma simulations. However, for lower energy collisions, which produce near isotropic or backward peaked differential cross sections, there is not a strong consensus among the plasma modeling community on how to best describe these angular scattering trends. Here, in this study, we propose a systematic method to approximate the aforementioned non-trivial angular scattering behavior with a formula that can be readily implemented in particle-in-cell (PIC) and/or MC plasma simulation codes. The present approach is specifically applied to fusion relevant atomic hydrogen and helium, as well as for molecular hydrogen, and results are also applicable to the atomic isotopes and homonuclear molecular isotopologues of these species. Comparisons between the present angular distribution function and benchmark scattering data were used to validate the proposed models. In addition, two-term Boltzmann calculations and PIC direct simulation MC simulations revealed that the proposed angular distribution function is accurate, agreeing very well with benchmark convergent close-coupling scattering calculations, and electron transport measurements. These studies confirmed that the present angular distribution function model can be utilized without the need of renormalization to the momentum transfer cross section (as opposed to using the elastic scattering integrated cross section), which has been suggested by several studies in order to correct for deficient angular scattering models, and to agree with transport measurements. Hence, the present anisotropic angular scattering model can be utilized to accurately model the momentum transfer as well as the electron trajectories of elastic collisions.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Introduction to spectral line shape theory
Spectral line-shape models are an important part of understanding high-energy-density (HED) plasmas. Models are needed for calculating opacity of materials and can serve as diagnostics for astrophysical and laboratory plasmas. However, much of the literature on line shapes is directed toward specialists. This perspective makes it difficult for non-specialists to enter the field. We have two broad goals with this topical review. First, we aim to give information so that others in HED physics may better understand the current field. Additionally, this first goal may help guide future experiments to test different aspects of the theory. Second, we provide an introduction for those who might be interested in line-shape theory, and enough materials to be able to navigate the field and the literature. We give a high-level overview of line broadening process, as well as dive into the formalism, available methods, and approximations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Cosmic rays in molecular clouds probed by H 2 rovibrational lines
Context. Low-energy cosmic rays (<1 TeV) play a fundamental role in the chemical and dynamical evolution of molecular clouds, as they control the ionisation, dissociation, and excitation of H 2 . Their characterisation is therefore important both for the interpretation of observations and for the development of theoretical models. However, the methods used so far for estimating the cosmic-ray ionisation rate in molecular clouds have several limitations due to uncertainties in the adopted chemical networks. Aims. Here, we refine and extend a previously proposed method to estimate the cosmic-ray ionisation rate in molecular clouds by observing rovibrational transitions of H 2 at near-infrared wavelengths, which are mainly excited by secondary cosmic-ray electrons. Methods. Combining models of interstellar cosmic-ray propagation and attenuation in molecular clouds with the rigorous calculation of the expected secondary electron spectrum and updated electron-H 2 excitation cross sections, we derive the intensity of the four H 2 rovibrational transitions observable in cold dense gas: (1—0)O(2), (1—0)Q(2), (1—0)S(0), and (1—0)O(4). Results. The proposed method allows the estimation of the cosmic-ray ionisation rate for a given observed line intensity and H 2 column density. We are also able to deduce the shape of the low-energy cosmic-ray proton spectrum impinging upon the molecular cloud. In addition, we present a look-up plot and a web-based application that can be used to constrain the low-energy spectral slope of the interstellar cosmic-ray proton spectrum. We finally comment on the capability of the James Webb Space Telescope to detect these near-infrared H 2 lines, which will make it possible to derive, for the first time, spatial variation in the cosmic-ray ionisation rate in dense gas. Besides the implications for the interpretation of the chemical-dynamic evolution of a molecular cloud, it will finally be possible to test competing models of cosmic-ray propagation and attenuation in the interstellar medium, as well as compare cosmic-ray spectra in different Galactic regions
79 ASTRONOMY AND ASTROPHYSICS↗
Understanding how minority relativistic electron populations may dominate charge state balance and radiative cooling of a post-thermal quench tokamak plasma
Minority relativistic electron populations can occur in a range of complex plasmas. Of specific interest is when runaway electrons form among the presence of high-atomic-number ion species in a tokamak plasma discharge. It has been recently demonstrated that ion charge state distributions and radiation losses at low bulk electron temperatures can be dominated by relativistic electrons, even though their density is orders of magnitude lower. Furthermore, this was attributed to the relativistic enhancement of electron impact inelastic cross sections. In this work, we provide a closer inspection of the atomic physics underpinning this effect. We also demonstrate the consequences of runaway enhanced scattering on post-disruption tokamak fusion discharges with neon and argon impurities present. Effects on charge state distributions, radiation and spectral characteristics, and reduced-order modeling considerations are discussed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
All-Order Full-Coulomb Quantum Spectral Line-Shape Calculations
Understanding how atoms interact with hot dense matter is essential for astrophysical and laboratory plasmas. Interactions in high-density plasmas broaden spectral lines, providing a rare window into interactions that govern, for example, radiation transport in stars. However, up to now, spectral line-shape theories employed at least one of three common approximations: second-order Taylor treatment of broadening operator, dipole-only interactions between atom and plasma, and classical treatment of perturbing electrons. In this Letter, we remove all three approximations simultaneously for the first time and test the importance for two applications: neutral hydrogen and highly ionized magnesium and oxygen. Finally, we found 15%–50% change in the spectral line widths, which are sufficient to impact applications including white-dwarf mass determination, stellar-opacity research, and laboratory plasma diagnostics.
79 ASTRONOMY AND ASTROPHYSICS↗
Linear polarization fractions of Fulcher- α fluorescence in electron collisions with H 2
We apply the molecular convergent close-coupling method to the calculation of linear polarization fractions for Fulcher-band fluorescence following electron-impact excitation of the H 2 d 3 Π u state. Furthermore, the results exhibit the opposite threshold behavior compared to the only previous calculations, but are in agreement with the most recent measurements for the Q(1), R(1), and Q(3) transitions.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Isotopic and vibrational-level dependence of H 2 dissociation by electron impact
The low-energy electron-impact dissociation of molecular hydrogen has been a source of disagreement between various calculations and measurements for decades. Excitation of the ground state of H 2 into the dissociative b 3 $Σ$$$^{+}_{u}$ state is now well understood, with the most recent measurements being in excellent agreement with the molecular convergent close-coupling (MCCC) calculations of both integral and differential cross sections [Zawadzki et al., Phys. Rev. A 98, 062704 (2018)]. However, in the absence of similar measurements for vibrationally excited or isotopically substituted H 2 , cross sections for dissociation of these species must be determined by theory alone. We have identified large discrepancies between MCCC calculations and the recommended R-matrix cross sections for dissociation of vibrationally excited H 2 , D 2 , T 2 , HD, HT, and DT [Trevisan et al., Plasma Phys. Contr. Fusion 44 1263 (2002); Plasma Phys. Contr. Fusion 44, 2217 (2002)], with disagreement in both the isotope effect and dependence on initial vibrational level. Furthermore, we investigate the source of the discrepancies, and discuss the consequences for plasma models, which have incorporated the previously recommended data.
74 ATOMIC AND MOLECULAR PHYSICS↗
Complete collision data set for electrons scattering on molecular hydrogen and its isotopologues: II. Fully vibrationally-resolved electronic excitation of the isotopologues of H 2 (X 1 $Σ^{+}_{g}$)
Here, we present a comprehensive set of vibrationally-resolved cross sections for electron-impact electronic excitation of the isotopologues of molecular hydrogen (D 2 , T 2 , HD, HT, and DT) initially in the ground electronic state. We apply the adiabatic-nuclei molecular convergent close-coupling (MCCC) method to calculate cross sections from threshold to 500 eV for excitation of all bound vibrational levels and dissociative excitation of the B 1 $Σ^{+}_{g}$, C 1 Π u , EF 1 $Σ^{+}_{g}$, B' 1 $Σ^{+}_{g}$, GK 1 $Σ^{+}_{g}$, I 1 Π g , J 1 Δ g , D 1 Π u , H 1 $Σ^{+}_{g}$, b 1 $Σ^{+}_{g}$, c 3 Π u , a 3 $Σ^{+}_{g}$, e 3 $Σ^{+}_{u}$, d 3 Π u , h 3 $Σ^{+}_{g}$, g 3 $Σ^{+}_{g}$, i 3 Π g , and j 3 Δ g electronic states from all bound vibrational levels of the ground electronic (X 1 $Σ^{+}_{g}$) state. Including the previously-published MCCC e-H2 cross sections the data set contains cross sections for over 60,000 electronic and vibrational transitions. The cross sections are presented in graphical form and provided as both numerical values and analytic fit functions in supplementary data files. The data can also be downloaded from the MCCC database at mccc-db.org.
74 ATOMIC AND MOLECULAR PHYSICS↗
Application of molecular convergent close-coupling cross sections in a collisional radiative model for the triplet system of molecular hydrogen
Collisional radiative (CR) models for molecular hydrogen are of high relevance for performing qualitative and quantitative analysis of excited-state population densities measured in plasmas or predicting the dependence of plasma emission on parameter variations. Although the development of such models for H 2 started decades ago, major uncertainties still exist regarding the most important set of input parameters, namely the cross sections for electron-impact excitation. The deviations between cross sections from different datasets are particularly pronounced in the energy region close to the threshold energy, strongly increasing the uncertainty of CR models applied to low-temperature plasmas. This paper presents experimental validation of a set of newly calculated non ro-vibrationally resolved electron-impact cross sections calculated for the triplet system of H 2 using the molecular convergent close-coupling method in the adiabatic-nuclei formulation. These cross sections are implemented into a CR model based on the flexible solver Yacora. A first comparison of CR calculations with the different datasets to experimentally-determined population densities is performed at a planar ICP discharge for varying pressure (between 1 and 10 Pa) and RF power (between 700 and 1100 W). For the experimentally-accessible electron temperature and density range (2.5–10 eV and 1.8–3.3 × 10 16 m –3 , respectively), very good agreement between the model and experiment is obtained using the new data set, in contrast to previously used cross sections.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Complete collision data set for electrons scattering on molecular hydrogen and its isotopologues: I. Fully vibrationally-resolved electronic excitation of H 2 $X^1Σ^+_g)$
Here, we present a comprehensive set of vibrationally-resolved cross sections for electron-impact electronic excitation of molecular hydrogen suitable for implementation in collisional-radiative models. The adiabatic-nuclei molecular convergent close-coupling method is used to calculate cross sections for excitation of all bound vibrational levels and dissociative excitation of the B 1 Σ u + , C 1 Π u , E F 1 Σ g + , B ′ 1 Σ u + , G K 1 Σ g + , I 1 Π g , J 1 Δ g , D 1 Π u , H 1 Σ g + , b 3 Σ u + , c 3 Π u , a 3 Σ g + , e 3 Σ u + , d 3 Π u , h 3 Σ g + , g 3 Σ g + , i 3 Π g , and j 3 Δ g electronic states from all –14 bound vibrational levels of the ground electronic ($\text{X}$ $^1Σ^+_g$) state. The data set consists of cross sections from threshold to 500 eV for over 5000 transitions, representing all possible electronic and vibrational transitions between the state and the –3 singlet and triplet states (where refers to the united-atoms-limit principle quantum number). The cross sections are presented in graphical form and provided as both numerical values and analytic fit functions in supplementary data files.
74 ATOMIC AND MOLECULAR PHYSICS↗