Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular collisions”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simulations of classical three-body thermalization in one dimension

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of pointlike particles, typically the dominant process for reaching thermal equilibrium in higher-dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential V ij ∝ 1/|x i –x j | n or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature T toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low-density limit. For the inverse power-law potential, the rate of thermalization is proportional to ρ 2 ⁢T$\frac{1}{2}$ – $\frac{1}{n}$, where ρ is the number density. Furthermore, the corresponding proportionality constant decreases with increasing n.

1-dimensional systems↗
DOE OSTI2021

Simulations of ion heating due to ion-acoustic instabilities in presheaths

Particle-in-cell, direct simulation Monte Carlo simulations reveal that ion-acoustic instabilities excited in presheaths can cause significant ion heating. Ion-acoustic instabilities are excited by the ion flow toward a sheath when the neutral gas pressure is small enough and the electron temperature is large enough. A series of 1D simulations were conducted in which neutral plasma (electrons and ions) was uniformly sourced with an ion temperature of 0.026 eV and different electron temperatures (0.1 eV–50 eV). Ion heating was observed when the electron-to-ion temperature ratio exceeded the minimum value predicted by linear response theory to excite ion-acoustic instabilities at the sheath edge (T e /T i ≈ 28). When this threshold was exceeded, the temperature equilibration rate between ions and electrons rapidly increased near the sheath so that the local temperature ratio did not significantly exceed the threshold for instability. This resulted in significant ion heating near the sheath edge, which also extended back into the bulk plasma; presumably due to wave reflection from the sheath. Furthermore, this ion-acoustic wave heating mechanism was found to decrease for higher neutral pressures, where ion-neutral collisions damp the ion-acoustic waves and ion heating is instead dominated by inelastic collisions in the presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Vibrationally resolved cross sections for proton collisions with hydrogen molecules

The semi-classical molecular convergent close-coupling approach has been applied to proton collisions with molecular hydrogen in the ground electronic and vibrational state. Cross sections for electron loss and electron capture agree well with the experimental data in the intermediate energy range where coupling between electronic reaction channels is strong and previously available calculations produce inconsistent results. Results for elastic scattering and total excitation are compared with effective one-electron coupled-channel calculations, showing that accurate target structure makes an important difference to these cross sections at intermediate energies. The excitation cross sections are calculated using both the fixed-nuclei and adiabatic-nuclei approximations for transitions from the $X^1Σ^+_{\textrm{g}}$ state to the $B^1Σ^+_{\textrm{u}}$, $B'$$^1Σ^+_{\textrm{u}}$, $B''$$^1Σ^+_{\textrm{u}}$, $EF$$^1Σ^+_{\textrm{g}}$, $GK$$^1Σ^+_{\textrm{g}}$, $H$$^1Σ^+_{\textrm{g}}$, $C$$^1Π_{\textrm{u}}$, $D$$^1Π_{\textrm{u}}$, $D'$$^1Π_{\textrm{u}}$, $I$$^1Π_{\textrm{u}}$, and $J$$^1Δ_{\textrm{g}}$ states for proton energies from to keV. We find significant differences between the present ab initio calculations and the equivelocity scaled electron-scattering data that is currently used for collisional-radiative modelling. Furthermore, the adiabatic-nuclei calculations enable us to resolve the final vibrational level after excitation, producing a complete set of cross sections for excitation of ground-state molecular hydrogen for all electronic states up to $n = 3$, where $n$ is the united-atoms-limit principle quantum number.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

A comprehensive conductivity model for drift and micro-tearing modes

The parallel electrical conductivity is a crucial parameter in the study of the linear stability of drift-modes like the micro-tearing mode (MTM). The conductivity enters by closing the electromagnetic tearing layer equations. Recent progress in the understanding of the pedestal suggests that the MTM could play an important role in its structure and evolution. For this reason, we revisit and improve previous model conductivities. This parameter is generally derived from the linearized drift kinetic equation. In the past literature, it has been computed using either simplified collision operators or neglecting the spatial dependence away from the rational surface. Here, a fully consistent expression for the conductivity that would accurately model the pedestal has not been available. By applying a novel variational procedure and with the full Fokker Plank collision operator, including electron–electron collisions, we compute a closed expression for the parallel electrical conductivity in the form of a rational function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Binding the Power of Cycloaddition and Cross-Coupling in a Single Mechanism: An Unexpected Bending Journey to Radical Chemistry of Butadiynyl with Conjugated Dienes

What if an experiment could combine the power of cycloaddition and cross-coupling with the in situ formation of an aromatic molecule in a single collision? Crossed molecular beam experiments augmented with electronic structure and statistical calculations provided compelling evidence on a novel radical route involving 1,3-butadiynyl (HCCCC; X 2 Σ + ) radicals synthesizing (substituted) arylacetylenes in the gas phase upon reactions with 1,3-butadiene (CH 2 CHCHCH 2 ; X 1 Ag) and 2-methyl-1,3-butadiene (isoprene; CH 2 C(CH 3 )CHCH 2 ; X 1 A'). Furthermore, this elegant mechanism de facto merges two previously disconnected concepts of cross-coupling and cycloaddition–aromatization in a single collision event via the formation of two new C(sp 2 )–C(sp 2 ) bonds and bending the 180° moiety of the linear 1,3-butadiynyl radical out of the ordinary by 60° to 120°. In addition to its importance to fundamental organic chemistry, this unconventional mechanism links two previously separated routes of gas-phase molecular mass growth processes of polyacetylenes and polycyclic aromatic hydrocarbons (PAHs), respectively, in low-temperature environments such as in cold molecular clouds like the Taurus Molecular Cloud (TMC-1) and in hydrocarbon-rich atmospheres of planets and their moons such as Titan, what revises the established understanding of low-temperature molecular mass growth processes in the Universe.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Impact of penetrating collisions of plasma ions on spectral line shapes

Spectral-line-broadening models have been moving towards including full Coulomb interactions between the atom and plasma, replacing the commonly used dipole approximation. The effects of the full Coulomb interaction have been thoroughly explored for plasma electrons, resulting in redshifts of spectral lines in high-energy-density plasmas. We explore the impact of a full Coulomb treatment on ion broadening. Penetrating collisions due to ions do not significantly impact the linewidths. The most significant aspect is the appearance of quasimolecular resonances in the far line wings, such as those previously observed in white-dwarf spectra. We identify several problems with the existing models. The most direct implementation causes the quasimolecular features to appear at the wrong photon energies. It is possible to include the detailed molecular structure, but this currently requires the Anderson-Talman approximation, which ignores the N-body properties of the plasma. We find that N-body effects can substantially broaden these quasimolecular features and even shift them, depending on the plasma conditions. We conclude that penetration of ions into the spatial extent of the radiator wave function does not strongly affect the usual diagnostics and may have a moderate to weak effect on opacity. Furthermore, the Rosseland mean opacity is weighted towards low-opacity regions of a spectrum, where these quasimolecular features are found.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Supervised machine learning-based multivariate regression of parallel closures for a high-collisionality deuterium-carbon plasma

Many plasmas of interest in laboratory experiments and space consist of multiple ion species. In tokamak edge plasmas, for instance, ionized impurities expelled from the vessel wall influence plasma transport. When describing multi-species plasmas using fluid equations, we need accurate closure relations to close the set of fluid equations. In this study, we introduce the development of fitting formulas for parallel closures using supervised machine learning, in conjunction with the recent closure theory, considering multi-ion collisions and arbitrary ion temperatures. We apply this approach to a high-collisionality deuterium-carbon plasma and demonstrate its effectiveness. As a result, the machine learning-based method for developing practical and accurate closures can be extended to a wider range of plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Multi-fluid and kinetic models of partially ionized magnetic reconnection

Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomenon in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet–Parker scaling, where neutrals are unimportant, or a fully coupled Sweet–Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfvén speed, $v^*_A$. In this work, these differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $ν_{ni}L/v^*_A$ ≳1, where $ν_{ni}$ is the neutral–ion collision frequency and L is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $ν_{ni}L/v^*_A$ <1, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Diffusion regime of electron–electron collisions in weakly ionized plasmas

We consider weakly ionized plasma where an elastic scattering of electrons on neutrals changes the individual acts and the rate of electron-electron collisions significantly. In this case, the kinetics of electron thermalization is very different from that in fully ionized plasma. Here, the colliding electrons do not move freely. They, instead, diffuse because of fast scattering on neutrals. We demonstrate how a proper account of this diffusion enables one to estimate the characteristic time of electron thermalization. We also present a rigorous derivation of the kinetic equation for electrons by using Bogolyubov’s method based on Liouville equations for multi-particle distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The plasma–sheath transition and Bohm criterion in a high recycling divertor

The high recycling regime of a divertor is characterized by high plasma particle fluxes and low temperature at the target, where a strong hydrogen recirculation loop exists. Atomic processes in the high recycling regime, such as ion–neutral friction and radiation, can affect the plasma momentum and energy transport in the sheath transition region. Here, in this work, the plasma–sheath transition near a high recycling wall is investigated. The Bohm speed, which constraints the ion exit flow speed, is evaluated from a transport model that accounts for the effect of the anisotropic transport and atomic collisions in the transition layer. A first principles kinetic code vector particle-in-cell with the atomic collision package is used to investigate a 1D self-consistent slab plasma with a high recycling boundary for the tungsten and carbon divertors. The results demonstrate the accuracy of the Bohm speed model in predicting the ion exit flow speed in the transition region, as well as the reduction of the Bohm speed due to the ion–neutral friction. The effect of different wall materials, tungsten, and carbon, on the Bohm speed and near-wall plasma profile is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Enforcing detailed balance in the Borgnakke–Larsen redistribution method with temperature dependent relaxation models

For decades, it has been observed that the commonly used Borgnakke–Larsen method for energy redistribution in Direct Simulation Monte Carlo codes fails to satisfy the principle of detailed balance when coupled to a wide variety of temperature dependent relaxation models, while seemingly satisfying detailed balance when coupled to others. Many attempts have been made to remedy the issue, yet much ambiguity remains, and no consensus appears in the literature regarding the root cause of the intermittent compatibility of the Borgnakke–Larsen method with temperature dependent relaxation models. This paper alleviates that ambiguity by presenting a rigorous theoretical derivation of the Borgnakke–Larsen method's requirement for satisfying detailed balance. Specifically, it is shown that the Borgnakke–Larsen method maintains detailed balance if and only if the probability of internal-energy exchange during a collision depends only on collision invariants (e.g., total energy). Here, the consequences of this result are explored in the context of several published definitions of relaxation temperature, including translational, total, and cell-averaged temperatures. Of particular note, it is shown that cell-averaged temperatures, which have been widely discussed in the literature as a way to ensure equilibrium is reached, also fail in a similar, although less dramatic, fashion when the aforementioned relationship is not enforced. The developed theory can be used when implementing existing or new relaxation models and will ensure that detailed balance is satisfied.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Multi-scale dynamics of atmospheric-pressure discharges ignited over liquid electrodes

A multiphase computational model of atmospheric-pressure gas discharges ignited over liquid electrodes is developed. The model takes into account both the liquid electrode motion under applied electric fields and plasma generation in the gas phase. The influence of the applied voltage and liquid properties (density, viscosity, and surface tension) on the liquid surface deformation is analyzed. It is shown that the shape of dynamic cones formed on the liquid surface differs from the static Taylor's cones. The influence of the liquid surface protrusions on gas breakdown dynamics is demonstrated. It is shown that the breakdown develops in two stages: first, a fast ionization wave propagates from the cathode to the anode, and then the cathode sheath collapses due to secondary electron emission from the liquid surface. It is shown that dynamics of the multiphase system containing gas, plasma, and liquid states is characterized by three disparate time scales: the fast electron time scale, the intermediate ion time scale, and the slow liquid dynamics time scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Impact of a minority relativistic electron tail interacting with a thermal plasma containing high-atomic-number impurities

A minority relativistic electron component can arise in both laboratory and naturally-occurring plasmas. In the presence of high-atomic-number ion species, the ion charge state distribution at low bulk electron temperature can be dominated by relativistic electrons, even though their density is orders of magnitude lower. This is due to the relativistic enhancement of the collisional excitation and ionization cross sections. The resulting charge state effect can dramatically impact the radiative power loss rate and the related Bethe stopping power of relativistic electrons in a dilute plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

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↗
DOE OSTI2022

Microwave breakdown of sub-atmospheric argon gas in the vicinity of a microgap

The microwave electric field at breakdown for argon gas in a microgap (generally, gap size <1 mm) over a wide range of pressures and gap sizes was found to exhibit two behavioral pressure regimes. A high microwave electric field was established in the gap of a re-entrant mode resonator. By slowly raising the microwave electric field level, a breakdown level is reached where impact ionization induces low pressure argon gas to breakdown preferentially outside of the gap. At high pressure, breakdown is found to occur preferentially inside the gap, but only after the pressure is high enough for impact ionization to produce breakdown there. Impact ionization has a lower threshold electric field than the boundary-controlled processes, such as secondary electron emission, that would eventually produce breakdown inside the microgap if impact ionization did not occur first. Thus, the lower pressure regime is characterized by breakdown outside of the gap, while in the high-pressure regime, breakdown occurs inside the gap, in both cases by impact ionization. The transition pressure P T between these two pressure regimes depends on the gap size d gap , roughly as P T ∝ d gap –2 . Consistent with the model developed here, this corresponds to the full-amplitude of electron oscillation in the gap field being commensurate with the gap size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Evaluation of a collisional radiative model for electron temperature determination in hydrogen plasma

A collisional-radiative (CR) model that extracts the electron temperature, T e , of hydrogen plasmas from Balmer-line-ratio measurements is examined for the plasma electron density, n e , and T e ranges of 10 10 –10 15 cm –3 and 5–500 eV, respectively. The CR code, developed and implemented in Python, has a forward component that computes the densities of excited states up to n = 15 as functions of T e , n e , and the molecular-to-atomic neutral ratio r(H 2 /H). The backward component provides ne and r(H 2 /H) as functions of the Balmer ratios to predict the T e . The model assumes Maxwellian electrons. Furthermore, the density profiles of the electrons and of the molecular and atomic hydrogen neutrals are shown to be of great importance, as is the accuracy of the line-ratio measurement method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nonlocal suppression of Biermann battery magnetic-field generation for arbitrary atomic numbers and magnetization

The Biermann battery term of magnetohydrodynamics (MHD) generates a magnetic field where electron density gradients and electron temperature gradients are perpendicular to one another. Kinetic simulations and experiments have shown that the rate of magnetic-field generation is lower than Biermann when the electron mean free path becomes comparable to or greater than the temperature gradient scale length, known as the nonlocal regime. We investigate the nonlocal suppression of the Biermann term using simplified Fokker–Planck simulations covering a wide range of parameters. We provide the first fit for nonlocal Biermann suppression that has physically accurate behavior for small and large values of a suitable nonlocality parameter, valid for an arbitrary atomic number, and that includes the effect of magnetization on nonlocality. Lastly, the fit is intended to provide an approximate method to account for reduced magnetic-field generation in MHD codes and theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗