Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamics of 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Transport model comparison studies of intermediate-energy heavy-ion collisions

Transport models are the main method to obtain physics information on the nuclear equation of state and in-medium properties of particles from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. To this end, calculations under controlled conditions of physical input and set-up were performed with various participating codes. These included both calculations of nuclear matter in a box with periodic boundary conditions, which test separately selected ingredients of a transport code, and more realistic calculations of heavy-ion collisions. Over the years, six studies have been performed within this project. In this intermediate review, we summarize and discuss the present status of the project. We also provide condensed descriptions of the 26 participating codes, which contributed to some part of the project. These include the major codes in use today. After a compact description of the underlying transport approaches, we review the main results of the studies completed so far. They show, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known under the names of Boltzmann–Uehling–Uhlenbeck (BUU) and Quantum Molecular Dynamics (QMD) type codes. However, when the codes were compared in full heavy-ion collisions using different physical models, as recently for pion production, they still yielded substantially different results. This calls for further comparisons of heavy-ion collisions with controlled models and of box comparisons of important ingredients, like momentum-dependent fields, which are currently underway. Furthermore, our evaluation studies often indicate improved strategies in performing transport simulations and thus can provide guidance to code developers. Results of transport simulations of heavy-ion collisions from a given code will have more significance if the code can be validated against benchmark calculations such as the ones summarized in this review.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1979

Semiclassical theory of electronically nonadiabatic transitions in molecular collision processes

An introductory account of the semiclassical theory of the S-matrix for molecular collision processes is presented, with special emphasis on electronically nonadiabatic transitions. This theory is based on the incorporation of classical mechanics with quantum superposition, and in practice makes use of the analytic continuation of classical mechanics into the complex space of time domain. The relevant concepts of molecular scattering theory and related dynamical models are described and the formalism is developed and illustrated with simple examples - collinear collision of the A+BC type. The theory is then extended to include the effects of laser-induced nonadiabatic transitions. Two bound continuum processes collisional ionization and collision-induced emission also amenable to the same general semiclassical treatment are discussed.

Lam, K. S.↗
DOE OSTI2022

Functionalization of Electrodes with Tunable [EMIM] x [Cl] x +1 – Ionic Liquid Clusters for Electrochemical Separations

Functionalization of electrodes with clusters of hydrophobic molecules may improve the energy efficiency and selectivity of electrochemical separations by modulating the desolvation process occurring at the interface. Ionic liquids (IL), which have a broad range of compositions and properties, are potential candidates for controlling the transport, desolvation, and adsorption of target ions at electrochemical interfaces. We report a joint experimental and theoretical investigation of the structure, stability, and selective adsorption properties of the IL clusters 1-ethyl-3-methylimidazolium chloride [EMIM] x [Cl] x+1 - (x = 1 – 10) and demonstrate their ability to adsorb and separate ions from solution. The structure and stability of the IL clusters are determined experimentally using high-mass-resolution electrospray ionization mass spectrometry, collision-induced dissociation, and negative ion photoelectron spectroscopy. Global optimization theory and ab initio molecular dynamics simulations provide molecular-level insight into the bonding and structural fluxionality of these species. Ion soft landing is used to selectively functionalize the surface of highly oriented pyrolytic graphite (HOPG) working electrodes with [EMIM] 1 [Cl] 2 - , [EMIM] 3 [Cl] 4 - , and [EMIM] 5 [Cl] 6 - clusters. Kelvin probe microscopy provides insight into the relative stability of the clusters on HOPG and their effect on the work function of IL-functionalized electrodes. Cyclic voltammetry measurements reveal irreversible adsorption of Fe(CN) 6 4-/3- anions during redox cycling, while electrochemical impedance spectroscopy indicates a substantial decrease in the electron transfer resistance of the IL-functionalized electrodes due to adsorption of Fe(CN) 6 4-/3- . Overall, our findings demonstrate that IL clusters with different size and stoichiometry may be used to increase the efficiency of electrochemical separations, opening new horizons in selective electrode functionalization.

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

Formation of a single quasicrystal upon collision of multiple grains

Abstract Quasicrystals exhibit long-range order but lack translational symmetry. When grown as single crystals, they possess distinctive and unusual properties owing to the absence of grain boundaries. Unfortunately, conventional methods such as bulk crystal growth or thin film deposition only allow us to synthesize either polycrystalline quasicrystals or quasicrystals that are at most a few centimeters in size. Here, we reveal through real-time and 3D imaging the formation of a single decagonal quasicrystal arising from a hard collision between multiple growing quasicrystals in an Al-Co-Ni liquid. Through corresponding molecular dynamics simulations, we examine the underlying kinetics of quasicrystal coalescence and investigate the effects of initial misorientation between the growing quasicrystalline grains on the formation of grain boundaries. At small misorientation, coalescence occurs following rigid rotation that is facilitated by phasons. Our joint experimental-computational discovery paves the way toward fabrication of single, large-scale quasicrystals for novel applications.

36 MATERIALS SCIENCE↗
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 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↗
NASA NTRS1988

Rotational and translational effects in collisions of electronically excited diatomic hydrides

Collisional quenching and vibrational energy proceed competitively with rotational energy transfer for several excited states of the diatomic radicals OH, NH, and CH. This occurs for a wide variety of molecular collision partners. This phenomenon permits the examination of the influence of rotational motion on the collision dynamics of these theoretically tractable species. Measurements can also be made as a function of temperature, i.e., collision velocity. In OH (sup 2 sigma +), both vibrational transfer and quenching are found to decrease with an increase in rotational level, while quenching decreases with increasing temperature. This behavior indicates that for OH, anisotropic attractive forces govern the entrance channel dynamics for these collisions. The quenching of NH (sup 3 pi sub i) by many (although not all) collision partners also decreases with increasing rotational and translational energy, and NH (sup 1 pi) behaves much like OH (sup 2 sigma +). However, the quenching of CH (sup 2 delta) appears to decrease with increasing rotation but increases with increasing temperature, suggesting in this case anisotropic forces involving a barrier or repulsive wall. Such similarities and differences should furnish useful comparisons with both simple and detailed theoretical pictures of the appropriate collision dynamics.

Crosley, David R.↗
DOE OSTI2021

Non-adiabatic quantum interference in the ultracold Li + LiNa → Li 2 + Na reaction

Electronically non-adiabatic effects play an important role in many chemical reactions. However, how these effects manifest in cold and ultracold chemistry remains largely unexplored. Here for the first time we present from first principles the non-adiabatic quantum dynamics of the reactive scattering between ultracold alkali-metal LiNa molecules and Li atoms. We show that non-adiabatic dynamics induces quantum interference effects that dramatically alter the ultracold rotationally resolved reaction rate coefficients. The interference effect arises from the conical intersection between the ground and an excited electronic state that is energetically accessible even for ultracold collisions. These unique interference effects might be exploited for quantum control applications such as a quantum molecular switch. The non-adiabatic dynamics are based on full-dimensional ab initio potential energy surfaces for the two electronic states that includes the non-adiabatic couplings and an accurate treatment of the long-range interactions. A statistical analysis of rotational populations of the Li2 product reveals a Poisson distribution implying the underlying classical dynamics are chaotic. The Poisson distribution is robust and amenable to experimental verification and appears to be a universal property of ultracold reactions involving alkali metal dimers.

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

Crystallographic and temperature effects in low-energy collisions for plasma–material interactions

The interaction of plasma with materials is critical for the fundamental understanding of non-equilibrium processes and their wide application. Recent achievements in fusion energy research emphasize the importance of this problem. Because modelling and predicting plasma–material interactions (PMI) require considering a tremendous number of single plasma particle–surface interaction events, this type of modelling was, until recently, possible only within a binary collision approximation (BCA). The BCA approach considers materials as uniform temperature-insensitive media. The research presented here addresses the PMI problem within an atomistic-based approach for the first time at a statistically significant level. Approximately 10 5 molecular dynamics trajectories were generated for 100 eV deuterium ions interacting with a tungsten surface. The research demonstrated the critical importance of incorporating the discrete lattice structure of matter into the model. Deuterium penetration depth and fraction of backscattered deuterium ions strongly depend on the surface orientation, the impact incident ion directions, and the ions’ initial positions. On average, the BCA-based calculations underestimate the penetration depth by a factor of two or more, and the fraction of backscattered atoms is overestimated by a factor of two or more. As a result, increasing the material temperature from 500 to 3000 K reduces penetration depth by about 30% and increases the fraction of backscattered atoms by 250%.

36 MATERIALS SCIENCE↗
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

Experimental investigation of abnormal transverse flow enhancement of α particles in heavy-ion collisions

The mass dependence of the transverse flow for Z = 1–5 fragments from the collisions of 40 Ar + 27 Al, 40 Ar + 48 Ti, and 40 Ar + 58 Ni at 47 MeV/nucleon is investigated experimentally in this article. The transverse flow values are determined using the in-plane components of the fragment transverse momenta, where three conventional methods, i.e., the kinetic flow tensor method, the transverse momentum analysis method, and the azimuthal correlation method, are applied to reconstruct the reaction plane in an event-by-event basis. It is demonstrated from the comparison of the present experimental mass dependent flow measurements and the model simulations using an improved antisymmetrized molecular dynamics model that the experimentally observed abnormal α transverse flow enhancement is closely related to the reaction plane reconstruction procedure in the flow extraction. We further investigate the physical existence of the abnormal α flow behavior using a two-particle azimuthal correlation method, which allows us to provide the relative flow magnitude information with an identification of fragment charge number without the knowledge of the reaction plane differing from the three conventional methods. It is found that the relative flow magnitudes deduced from the two-particle azimuthal correlation functions with an identification of Z, with the correction for the recoil effect imposed by the momentum conservation, show a monotonically increasing trend as a function of fragment charge number, with no exception of the α flow enhancement. These results, in addition to those from the improved antisymmetrized molecular dynamics model simulations, definitely provide experimental evidences for the inexistence of the abnormal α flow behavior in the heavy-ion collisions at the present incident energy region in nature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Constraining nucleon effective masses with flow and stopping observables from the S π RIT experiment

Properties of the nuclear equation of state (EoS) can be probed by measuring the dynamical properties of nucleus-nucleus collisions. In this study, we present the directed flow (υ 1 ), elliptic flow (υ 2 ) and stopping (VarXZ) measured in fixed target Sn + Sn collisions at 270 AMeV with the SπRIT Time Projection Chamber. We perform Bayesian analyses in which EoS parameters are varied simultaneously within the Improved Quantum Molecular Dynamics-Skyrme (ImQMD-Sky) transport code to obtain a multivariate correlated constraint. The varied parameters include symmetry energy, S 0 , and slope of the symmetry energy, L, at saturation density, isoscalar effective mass, $m^{⁎}_{S}$/$m_{N}$ , isovector effective mass, $m^{⁎}_{υ}$/$m_{N}$ and the in-medium cross-section enhancement factor η. We find that the flow and VarXZ observables are sensitive to the splitting of proton and neutron effective masses and the in-medium cross-section. Comparisons of ImQMD-Sky predictions to the SπRIT data suggest a narrow range of preferred values for $m^{⁎}_{S}$/$m_{N}$, $m^{⁎}_{υ}$/$m_{N}$ and η.

Astronomy & Astrophysics↗
DOE OSTI2026

Non-adiabatic quantum interference and complex formation in ultracold collisions of Rb with KRb

Ultracold elastic collisions of 87 Rb with 40 K 87 Rb in its ground vibrational and rotational state are investigated using a first principles based theoretical methodology. Full-dimensional ab initio computed potential energy surfaces are reported that include the two lowest-lying electronic states, their conical intersection, non-adiabatic couplings and an accurate long-range behavior. A numerically exact time-independent quantum dynamics method in hyperspherical coordinates is used to compute the elastic scattering cross sections, rate coefficients and collision lifetime spectrum. The quantum scattering calculations include all degrees of freedom and treat both electronic states and their non-adiabatic couplings using a two-state diabatic representation. The theoretically computed elastic rate coefficient is in good agreement with the recently reported experimental value. Significant non-adiabatic quantum interference effects are shown to originate from the unique properties of ultracold collisions and the geometric phase associated with the conical intersection. A high-resolution collision energy grid is used to investigate the origin of the experimentally reported long-lived 3-body collision complexes.

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

State-to-state rate coefficients for HCS+ in rotationally inelastic collisions with H2 at low temperatures

ABSTRACT HCS+ ions have been detected in several regions of the interstellar medium (ISM), but an accurate determination of the chemical-physical conditions in the molecular clouds where this molecule is observed requires detailed knowledge of the collisional rate coefficients with the most common colliders in those environments. In this work, we study the dynamics of rotationally inelastic collisions of HCS+ + H2 at low temperature, and report, for the first time, a set of rate coefficients for this system. We used a recently developed potential energy surface for the HCS+–H2 van der Waals complex and computed state-to-state rotational rate coefficients for the lower rotational states of HCS+ in collision with both para- and ortho-H2, analysing the influence of the computed rate coefficients on the determination of critical densities. Additionally, the computed rate coefficients are compared with those obtained by scaling the ones from HCS+ in collision with He (an approximation that is sometimes used when data is lacking), and large differences are found. Furthermore, the approximation of using the rates for the HCO+ + H2 collision as a rough approximation for those of the HCS+ + H2 system is also evaluated. Finally, the complete set of de-excitation rate coefficients for the lowest 30 rotational states of HCS+ by collision with H2 is reported from 5 to 100 K.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Directed gas phase preparation of ethynylallene (H 2 CCCHCCH; X 1 A') via the crossed molecular beam reaction of the methylidyne radical (CH; X 2 Π) with vinylacetylene (H 2 CCHCCH; X 1 A')

The gas-phase bimolecular reaction of the methylidyne (CH; X 2 Π) radical with vinylacetylene (H 2 CCHCCH; X 1 A') was conducted at a collision energy of 20.3 kJ mol –1 under single collision conditions exploiting the crossed molecular beam experimental results merged with ab initio electronic structure calculations and ab initio molecular dynamics (AIMD) simulations. The laboratory data reveal that the bimolecular reaction proceeds barrierlessly via indirect scattering dynamics through long-lived C 5 H 5 reaction intermediate(s) ultimately dissociating to C 5 H 4 isomers along with atomic hydrogen with the latter predominantly originating from the vinylacetylene reactant as confirmed by the isotopic substitution experiments in the D1-methylidyne–vinylacetylene reaction. Combined with ab initio calculations of the potential energy surface (PES) and statistical Rice–Ramsperger–Kassel–Marcus (RRKM) calculations, the experimental determined reaction energy of –146 ± 26 kJ mol –1 along with the distribution minimum of T(θ) at 90° and isotopic substitution experiments suggest ethynylallene (p1; ΔrG = –230 ± 4 kJ mol –1 ) as the dominant product. The ethynylallene (p1) may be formed with extensive rovibrational excitation, which would result in a lower maximum translational energy. Further, AIMD simulations reveal that the reaction dynamics leads to p1 (ethynylallene, 75%) plus atomic hydrogen with the dominant initial complex being i1 formed by methylidyne radical addition to the double C[double bond, length as m-dash]C bond in vinylacetylene. Altogether, combining the crossed molecular beam experimental results with ab initio electronic structure calculations and ab initio molecular dynamics (AIMD) simulations, ethynylallene (p1) is expected to represent the dominant product in the reaction of the methylidyne (CH; X 2 Π) radical with vinylacetylene (H 2 CCHCCH; X 1 A').

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

Sputtering Erosion in the Ion Thruster

During the first phase of this research, the sputtering yields of molybdenum by low energy (100 eV and higher) xenon ions were measured by using the methods of secondary neutral mass spectrometry (SNMS) and Rutherford backscattering spectrometry (RBS). However, the measured sputtering yields were found to be far too low to explain the sputtering erosions observed in the long-duration tests of ion thrusters. The only difference between the sputtering yield measurement experiments and the ion thruster tests was that the later are conducted at high ion fluences. Hence, a study was initiated to investigate if any linkage exists between high ion fluence and an enhanced sputtering yield. The objective of this research is to gain an understanding of the causes of the discrepancies between the sputtering rates of molybdenum grids in an ion thruster and those measured from our experiments. We are developing a molecular dynamics simulation technique for studying low-energy xenon ion interactions with molybdenum. It is difficult to determine collision sequences analytically for primary ions below the 200 eV energy range where the ion energy is too low to be able to employ a random cascade model with confidence and it is too high to have to consider only single collision at or near the surface. At these low energies, the range of primary ions is about 1 to 2 nm from the surface and it takes less than 4 collisions on the average to get an ion to degrade to such an energy that it can no longer migrate. The fine details of atomic motion during the sputtering process are revealed through computer simulation schemes. By using an appropriate interatomic potential, the positions and velocities of the incident ion together with a sufficient number of target atoms are determined in small time steps. Hence, it allows one to study the evolution of damages in the target and its effect on the sputtering yield. We are at the preliminary stages of setting up the simulation program.

Ray, Pradosh K.↗
NASA NTRS1997

Data Base for CFD Validation

The flow behind the shock wave formed around objects which fly at hypervelocity behaves differently from that of a perfect gas. Molecules become vibrationally excited, dissociated, and ionized. The hot gas may emit or absorb radiation. When the atoms produced by dissociation reach the wall surface, chemical reactions, including recombination, may occur. The thermochemical phenomena of vibration, dissociation, ionization, surface chemical reaction, and radiation are referred to commonly as high-temperature real-gas phenomena. The phenomena cause changes in the dynamic behavior of the flow and the surface pressure and heat transfer distribution around the object. The character of a real gas is described by the internal degrees of freedom and state of constituent molecules; nitrogen and oxygen for air. The internal energy states, rotation, vibration and electronic, of the molecules are excited and, in the limit, the molecular bonds are exceeded and the gas dissociated into atomic and, possibly, ionic constituents. The process of energy transfer causing excitation, dissociation and recombination is a rate process controlled by particle collisions. Binary, two-body, collisions are sufficient to cause internal excitation, dissociation and ionization while three-body collisions are required to recombine the particles into molecular constituents. If the rates of energy transfer are fast with respect to the local fluid dynamic time scale the gas is in, or nearly in, equilibrium. If the energy transfer rates are very slow the gas can be described as frozen. In all other instances, wherein any of the energy exchange rates are comparable to the local fluid time scale, the gas will be thermally or chemically reacting and out of equilibrium. Real gas thermochemical nonequilibrium processes are important in the determination of aerodynamic heating; both convective (including wall catalytic effects) and radiative heating. To illustrate this we consider the hypervelocity flow over a bluff body typical of an atmospheric entry vehicle or an aerospace transfer vehicle (ASTV). The qualitative aspects of a hypersonic flow field over a bluff body are discussed in two parts, forebody and afterbody, with attention to which particular physical effects must be included in an analysis. This will indicate what type of numerical modeling will be adequate in each region of the flow. A bluff forebody flow field is dominated by the presence of the strong bow shock wave and the consequent heating, and chemical reaction of the gas. At high altitude hypersonic flight conditions the thermal excitation and chemical reaction of the gas occur slowly enough that a significant portion of the flow field is in a state of thermochemical nonequilibrium. A second important effect is the presence of the thick boundary layer along the forebody surface. In this region there are large thermal and chemical species gradients due to the interaction of the gas with the wall. Also at high altitudes the shock wave and the boundary layer may become so thick that they merge; in this case the entire shock layer is dominated by viscous effects.

Deiwert, George S.↗