Search NASA⌕ Search

SEARCH · Search NASA

Results for “MOLECULAR COLLISION”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Quantum scattering of HC 5 N and para -H 2 on a new potential energy surface

In the interstellar medium (ISM), non-local thermodynamic equilibrium situations are common due to low density, and one needs to consider the effect of molecular collisions in order to interpret the observations. Among the species detected in the ISM, cyanopolyynes, with the general molecular formula HC 2n+1 N (n = 1, 2, …), are characterized by large dipole moments and small rotational constants and constitute an indispensable class of candidates for the sensitive tracers of local density and temperature. We present a study of the collisional (de-) excitation of HC 5 N by para -H 2 (p-H 2 ) in its ground rotational state, namely HC 5 N ( j 1 ) + H 2 ( j 2 = 0) → HC 5 N (j$_1^′$) + H2 (j$_2^′$ = 0), where j 1 (or j$_1^′$) and j 2 (or j$_2^′$) denote the initial (or final) rotational quantum numbers of HC 5 N and H 2 , respectively. We performed the quantum scattering calculations at low collision energy using a new four-dimensional ab initio potential energy surface. In the regime where p-H 2 remains in its rotational ground state, converged cross sections did not require including excited rotational states of p-H 2 in the rotational basis. State-to-state cross sections were computed by means of the quantum-mechanical close-coupling (CC) method and the coupled states (CS) approximation, and rate coefficients for the first 61 levels of HC 5 N were computed for the first time up to 20 K with the CC approach and up to 50 K with the CS method. CC and CS results were found to agree well at temperatures up to 20 K. Finally, these data should allow a more accurate derivation of the HC 5 N abundance in molecular clouds.

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

Ab initio investigation of a hypersonic double cone experiment

This article presents a direct molecular simulation (DMS) of a reactive Mach 8.2 oxygen flow over a double cone geometry. The free stream conditions and article configuration generate a flow with thermal and chemical nonequilibrium, which are common attributes of hypersonic flight. This scenario was first studied experimentally at Calspan University of Buffalo Research Center’s test facility. DMS is a particle method that uses quantum mechanically derived interaction potentials to simulate molecular collisions within a flow field. Since these interaction potentials are the only modeling inputs used in the simulation, all flow features can solely be attributed to the ab initio potential energy surfaces. Hence, providing a comparison of a hypersonic ground test and numerical data anchored to quantum mechanics. The fundamental nature of DMS is leveraged to investigate molecular level mechanisms prevalent in the flow, and comparisons with lower fidelity simulations are presented to highlight the role of these first principles calculations as benchmark solutions.

Science & Technology - Other Topics↗
DOE OSTI2026

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗
DOE OSTI2024

Utilizing Quantum Cascade Lasers for Ultranarrow Velocity Resolution and Quantum-State Selectivity in Molecular Beam Scattering and Spectroscopy

Here, we demonstrate the capability of a narrow linewidth quantum cascade laser (QCL) to selectively excite a very narrow velocity range of nitric oxide (σ ≤ 7(3) m/s) with a pure ro-vibrational quantum state. By implementing a counter-propagating geometry, the molecules are selectively excited according to the Doppler shift of the ro-vibrational transition frequency such that the velocity width associated with the excited molecules depends only on the QCL linewidth. We demonstrate a velocity distribution limited by the effective linewidth of our free-running QCL (Γ = 3.2 MHz). Our development provides a cost-effective, flexible approach to resolve quantum-state selective chemical dynamics with excellent velocity resolution in a wide variety of molecules with infrared-active transitions. This technique has been formulated to provide ultrahigh collisional energy resolution in molecular beams to delineate final quantum-state product pairs in studies of molecular collisions.

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

Experimental validation of a collision-radiation dataset for molecular hydrogen in plasmas

Quantitative spectroscopy of molecular hydrogen has generated substantial demand, leading to the accumulation of diverse elementary process data encompassing radiative transitions, electron-impact transitions, predissociations, and quenching. However, their rates currently available are still sparse, and there are inconsistencies among those proposed by different authors. In this study, we demonstrate an experimental validation of such a molecular dataset by composing a collisional-radiative model (CRM) for molecular hydrogen and comparing experimentally obtained vibronic populations across multiple levels. From the population kinetics of molecular hydrogen, the importance of each elementary process in various parameter space is studied. In low-density plasmas (electron density ne≲1017 m−3) the excitation rates from the ground states and radiative decay rates, both of which have been reported previously, determine the excited state population. The inconsistency in the excitation rates affects the population distribution the most significantly in this parameter space. However, in higher density plasmas (ne≳1018 m−3), the excitation rates from excited states become important, which have never been reported in the literature, and may need to be approximated in some way. In order to validate these molecular datasets and approximated rates, we carried out experimental observations for two different hydrogen plasmas; a low-density radio frequency heated plasma (ne≈1016 m−3) and the Large Helical Device (LHD) divertor plasma (ne≳1018 m−3). The visible emission lines from EF1Σg+, HH¯1Σg+, D1Πu±, GK1Σg+, I1Πg±, J1Δg±, h3Σg+, e3Σu+, d3Πu±,g3Σg+, i3Πg±, and j3Δg± states were observed simultaneously and their population distributions were obtained from their intensities. We compared the observed population distributions with the CRM prediction, in particular the CRM with the rates compiled by Janev et al., Miles et al., and those calculated with the molecular convergent close-coupling (MCCC) method. The MCCC prediction gives the best agreement with the experiment, particularly for the emission from the low-density plasma. However, the population distribution in the LHD divertor shows a worse agreement with the CRM than those from low-density plasma, indicating the necessity of the precise excitation rates from excited states. We also found that the rates for the electron attachment is inconsistent with experimental results. This requires further investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 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 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 OSTI2025

Enhanced shear stabilization of turbulence in NSTX

In studying a particular non-stationary NSTX L-mode plasma, we observed unexpectedly high levels of flux—first with quasilinear (TGLF) modeling and subsequently with nonlinear gyrokinetic simulations. Upon more detailed analysis, a novel confinement regime was discovered in which a modest increase in E x B shear (beyond baseline experimental estimates) rapidly reduced turbulent transport to levels consistent with power balance. This modest increase is plausible given the errors inherent to the estimation of shearing rates, and the added complexity of the non-stationary (time-dependent) power balance. Remarkably, an additional small increase in shear yields the familiar ion-neoclassical transport level with what appears to be the onset of high-k electron transport only. Although analyses using the TGLF-SAT2 model successfully capture numerous parametric dependencies of this plasma, TGLF does not reproduce the rapid E x B stabilization seen in CGYRO. We believe the results presented should help to better characterize the nonlinear physics of spherical tokamak confinement regimes, provide useful ST datasets for reduced model development, and motivate more accurate experimental diagnosis of E x B shearing rates.

Atomic and molecular collisions↗
DOE OSTI2025

Anomalies in the Electronic Stopping of Slow Antiprotons in LiF

We present first-principles theoretical calculations for the electronic stopping power (SP) of both protons and antiprotons in LiF. Our results show the presence of the Barkas effect: a higher stopping for positively charged particles than their negatively charged antiparticles. In contrast, a previous study has predicted an anti-Barkas effect (higher stopping for negative charges) at low velocity [Qi, Bruneval and Maliyov, Phys. Rev. Lett. 128, 043401 (2022)]. Here, we explain this discrepancy by showing that this anti-Barkas effect appears for highly symmetric trajectories and disappears when considering trajectories that better reproduce the experimental setup. Our low-velocity results show that the SP of both protons and antiproton vanish for velocities under 0.1 a.u.

3-dimensional systems↗
DOE OSTI2026

Close-coupling approach to electron scattering with multielectron targets

Momentum-space close coupling calculations of electron scattering require removal of spurious and unphysical solutions. Here, we demonstrate here that removal of these solutions involve regulator operators that enforce Pauli exclusion selection rules in addition to removing spurious solutions. The form of a regular operator for e-H scattering has already been established, but a general extension to the multielectron case has been elusive. Here we present a general method for scattering on multielectron targets, atoms, or molecules, ensuring that the scattering solutions obey Pauli-exclusion selection rules. The regulator operator is obtained by finding the null space vectors of the 𝑁+1 electrons of the projectile and target atom scattering system. We demonstrate that this general procedure reduces to the e-H result and provide examples for He- and Li-like targets as well as guidance for implementation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Mapping Rydberg States of H 2 with the Halfium R-Matrix Method

In this article, we use the Halfium R-matrix method to investigate the Rydberg states of the H 2 molecule up to n = 20, filling the gap above the low-lying bound states already calculated with configuration interaction packages. Moreover, we show that the use of Quantum Defect Theory scaling laws, allows for a comprehensive analysis of the regular patterns resulting from the coupling between Rydberg series and doubly excited states. The results should open the door for more efficient quasi-diabatization of the potential energy curves which is required for calculating cross sections and rate coefficients of the (e + H 2 + ) collisional processes, involved in the plasma modeling for fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Resistive Coatings for High-performance, Low-background MCPs Operating Across Broad Temperature Ranges and at Cryogenic Temperatures

Microchannel plate with improved thermo-electric properties are a high risk, high payoff development undertaken by a consortium of effort that links the Argonne National Lab, the Space-Science Lab at UC Berkeley and the small businesses (Incom Inc.) that will commercialize the advanced technology in open MCPs and LAPPDs. This development will satisfy the needs for new instrumentation for homeland security (non-proliferation) sensors to screen vehicles and cargo for Special Nuclear Materials (SNMs) and scientific detectors for astrophysics, electron microscopy, time-of-flight mass spectrometry, molecular and atomic collision studies, and fluorescence imaging applications in biotechnology and medical imaging products including positron emission tomography (PET scanning).

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Competing ionization and dissociation: Extension of the energy-dependent frame transformation to the gerade symmetry of H 2

This article solves two major tasks that frequently arise in the theory of electron collisions with a target molecular cation. First, it extends the energy-dependent frame transformation (EDFT) treatment, which is needed to map fixed-nuclei electron-molecule scattering matrices into an energy-dependent laboratory-frame scattering matrix with vibrational channel indices. The EDFT mapping can now be carried out even when the target molecule possesses multiple low-energy potential curves, significantly transcending previous applications. Second, it implements a method to extract the rest of the full laboratory-frame scattering matrix, i.e., the columns and rows describing input and/or output dissociation channels. The treatment is benchmarked in this article against the essentially exact solution of a refined two-dimensional model of the singlet gerade Σ symmetry of H 2 . Our tests demonstrate that the theory accurately maps fixed-nuclei scattering information, of the type provided by existing electron-molecule computer codes, into a laboratory-frame scattering matrix that includes both ionization and dissociation. Furthermore, this treatment can provide a general framework applicable to a broad class of electron collision processes involving diatomic target ions, suitable for an accurate description of challenging processes such as dissociative recombination.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Formation of 1 H -Phenalene (C 13 H 10 ) in the Taurus Molecular Cloud via Methylidyne Addition-Cyclization-Aromatization (MACA)

The formation of 1H-phenalene (C 13 H 10 ) in cold molecular clouds, such as the Taurus Molecular Cloud-1 (TMC-1), presents a significant challenge to traditional astrochemical models, which predominantly suggest high-temperature pathways for polycyclic aromatic hydrocarbon (PAH) formation. In this study, we explore computationally the Methylidyne Addition-Cyclization-Aromatization (MACA) mechanism as a viable, barrierless pathway for phenalene synthesis under low-temperature conditions. Through electronic structure calculations and Rice–Ramsperger–Kassel–Marcus (RRKM) statistical methods, we demonstrate that the reaction of 1-vinylnaphthalene (C 10 H 7 C 2 H 3 ) with the methylidyne radical (CH) leads to the formation of 1H-phenalene via a bimolecular reaction, a process that is exoergic and without entrance barrier. The MACA mechanism facilitates the growth of the aromatic carbon backbone via a [5 + 1] ring annulation, providing a new insight into PAH formation in cold molecular clouds. Notably, the MACA mechanism has previously been shown to form indene (C9H8), which was detected in TMC-1 as well, via a [4 + 1] annulation, demonstrating its potential to produce a variety of complex PAHs by addition of a five- and six-membered ring to a benzene moiety via [4 + 1] and [5 + 1] annulation, respectively. As a result, this work highlights the importance of barrierless, exoergic reactions involving MACA in the synthesis of complex aromatic molecules in space, expanding our physicochemical understanding of carbon-rich chemistry in cold molecular clouds.

Aromatic compounds↗
DOE OSTI2025

Gas-phase synthesis of naphthalene through an unconventional thermal alkyne–alkene [2 + 2] cycloaddition mechanism

Exotic cycloaddition entrance channels were discovered for the bimolecular gas-phase reactions of the phenylethynyl radical (C 6 H 5 CC, X 2 A 1 ) with ethylene-d 4 (C 2 D 4 ) and propylene (C 3 H 6 ) as explored under single-collision conditions utilizing the crossed molecular beams technique combined with electronic structure and statistical calculations. Connecting the concepts of barrierless entrance channels, excited states, and facile non-photochemically activated cycloadditions, the reaction pathway features an unconventional thermal [2 + 2] cycloaddition forming a four-membered ring collision complex followed by multiple isomerizations prior to unimolecular decomposition via atomic hydrogen loss to (un)substituted naphthalenes—naphthalene-d 4 (C 10 H 4 D 4 ) and 1-/2-methylnaphthalene (C 11 H 10 ). The small energy gap between the singly-occupied a 1 highest occupied molecular orbital (HOMO) with a σ-character and the underlying doubly-occupied b 1 molecular orbital with a π-character allows a facile promotion of an electron. This in turn enables a versatile low-temperature reactivity of phenylethynyl, where the end-on and side-on barrierless approaches of ethylene are due to its interaction with the σ and π orbitals, respectively, thus suggesting this mechanism as a possible method for tuning substituents in polycyclic aromatic hydrocarbon (PAH) formation and highlighting its versatility as a probe of fundamental carbon chemistry via counterintuitive cycloaddition reactions under single-collision conditions.

Goettl, Shane J. [University of Hawaii at Manoa, H↗
DOE OSTI2025

Bayesian inference of nuclear incompressibility from collective flow in mid-central Au+Au collisions at 400–1500 MeV/nucleon

The incompressibility K of symmetric nuclear matter (SNM) is determined through a Bayesian analysis of collective flow data from Au + Au collisions at beam energies $E = 400 -1500$ MeV/nucleon. This analysis utilizes a Gaussian process (GP) emulator applied to the isospin-dependent quantum molecular dynamics (IQMD) model for heavy-ion collisions, both with and without incorporating the momentum dependence of the single-nucleon potentials. Specifically, at the 68% confidence level, using rapidity and transverse velocity dependence of proton elliptic flow data with and without consideration of the momentum dependence, the inferred incompressibility values are $K=188.9^{+2.9}_{-4.5}$ MeV and $256.1^{+8.2}_{-8.7}$ MeV at $E = 400$ MeV/nucleon, respectively. When the transverse momentum dependence of proton-like directed flow data is included, the inferred incompressibility values become $K=222.3^{+9.0}_{-9.9}$ MeV and $K=285.5^{+6.7}_{-7.3}$ MeV, respectively. Furthermore, we found that the value of K derived from observables of proton elliptic flow increases with beam energy. Finally, this indicates that the equation of state (EoS) of nuclear matter hardens at higher densities and temperatures in reactions with higher beam energies.

Bayesian inference↗