Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Rotational excitation and de-excitation of magnesium mono-sulphide (MgS) by He collisions

ABSTRACT Magnesium mono-sulphide (MgS) plays a crucial role in astrochemical processes, particularly in the interstellar medium (ISM), where metal-sulphur chemistry influences molecular formation and evolution. This work presents a detailed study of the rotational excitation and de-excitation of MgS through collisions with helium (He) atoms, the second most abundant species in the ISM after hydrogen, which includes both atomic (H) and molecular forms (H2). The focus on MgS–He collisions arises from He's high abundance, chemical inertness, and simpler electronic structure, which make it well suited for quantum scattering calculations. These characteristics establish He as an ideal candidate for initial studies, providing fundamental data for future investigations involving H2. The study uses quantum scattering methods to calculate the collisional rate coefficients over a broad temperature range. These rates are critical for interpreting observational data on MgS and predicting its abundance in space. The interaction potential between MgS and He is calculated using the rigid rotor approximation and the Jacobi coordinate system, employing the CCSD(T)-F12a/aug-cc-pVTZ method for accurate two-dimensional potential energy surface. The study explores the anisotropic nature of the MgS–He interaction, which favours odd Δj rotational transitions at low collision energies. The inelastic cross-sections for rotational transitions involving up to 16 rotational levels of MgS were computed up to 1000 cm−1, enabling the calculation of rate coefficients up to 150 K for Δj = ±1, ±2, and ± 3 rotational transitions. The results show that Δj = 1 transitions dominate at low temperatures, while Δj = 2 transitions become more significant at higher temperatures. This study provides valuable data for interpreting future astrophysical observations of MgS. The findings also propose new rotational transitions for MgS detection in space, enhancing our ability to track and study this molecule in various cosmic environments.

Hendaoui, Hamza (ORCID:0000000218641872)↗
DOE OSTI2025

Parton distribution functions from scalar light-front parton gas model

Here, we propose an application of a microcanonical ensemble with light-front kinematics to model the phase-space distribution of relativistic constituents of a bound state. These constituents denoted by partons are treated as classical spin-zero particles confined inside the bound state with inter-parton collisions as their only interaction. The microcanonical molecular dynamics ensemble is applied to obtain the phase-space distribution of such a thermodynamic system. We sample this phase-space distribution using Monte Carlo algorithms to obtain the parton distribution functions (PDFs) in scenarios with 3, 4, and 5 identical partons. In addition PDFs when a selected number of massless partons are mixed with 3 massive partons are also presented.

Microcanonical ensemble↗
DOE OSTI2024

Direct observation of ultrafast cluster dynamics in supercritical carbon dioxide using X-ray Photon Correlation Spectroscopy

Supercritical fluids exhibit distinct thermodynamic and transport properties, making them of particular interest for a wide range of scientific and engineering applications. These anomalous properties emerge from structural heterogeneities due to the formation of molecular clusters at conditions above the critical point. While the static behavior of these clusters and their effects on the thermodynamic response functions have been recognized, the relation between the ultrafast cluster dynamics and transport properties remains elusive. By measuring the intermediate scattering function in carbon dioxide at conditions near the critical point with X-ray photon correlation spectroscopy, we directly capture the cross-over dynamics between 4 and 13 picoseconds, revealing the transition between ballistic and diffusive motion. Complementary analysis using large-scale molecular dynamics simulations reveals that this behavior arises from collisions between unbound molecules and clusters. This study provides direct evidence of the ultrafast momentum exchange between clusters, which has significant impact on transport properties, solvation processes, and reaction kinetics in supercritical fluids.

carbon capture and storage↗