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 127 records · Page 7

One Collision—Two Substituents: Gas–Phase Preparation of Xylenes under Single–Collision Conditions

The fundamental reaction pathways to the simplest dialkylsubstituted aromatics—xylenes (C 6 H 4 (CH 3 ) 2 )—in high-temperature combustion flames and in low-temperature extraterrestrial environments are still unknown, but critical to understand the chemistry and molecular mass growth processes in these extreme environments. Exploiting crossed molecular beam experiments augmented by state-of-the-art electronic structure and statistical calculations, this study uncovers a previously elusive, facile gas-phase synthesis of xylenes through an isomer-selective reaction of 1-propynyl (methylethynyl, CH 3 CC) with 2-methyl-1,3-butadiene (isoprene, C 5 H 8 ). The reaction dynamics are driven by a barrierless addition of the radical to the diene moiety of 2-methyl-1,3-butadiene followed by extensive isomerization (hydrogen shifts, cyclization) prior to unimolecular decomposition accompanied by aromatization via atomic hydrogen loss. This overall exoergic reaction affords a preparation of xylenes not only in high-temperature environments such as in combustion flames and around circumstellar envelopes of carbon-rich Asymptotic Giant Branch (AGB) stars, but also in low-temperature cold molecular clouds (10 K) and in hydrocarbon-rich atmospheres of planets and their moons such as Triton and Titan. Furthermore, our study established a hitherto unknown gas-phase route to xylenes and potentially more complex, disubstituted benzenes via a single collision event highlighting the significance of an alkyl-substituted ethynyl-mediated preparation of aromatic molecules in our Universe.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamical control in a prethermalized molecular ultracold plasma: Local dissipation drives global relaxation

Prethermalization occurs as an important phase in the dynamics of isolated many-body systems when coupled degrees of freedom relax in a subspace separated from a state of complete thermodynamic equilibrium by a gap in energy or other conserved quantity. Slow equilibration from a prethermal phase can localize the dynamics of natural and model systems despite high dimensionality and limited disorder. Here, we report the signature of an enduring prethermal regime of arrested relaxation in the molecular ultracold plasma that forms following the avalanche of a state-selected Rydberg gas of nitric oxide. For a wide range of initial conditions, this system enters a critical phase in which a density of NO + and electrons balances a population of Rydberg molecules. Electron collisions mix orbital angular momentum, scattering Rydberg electrons to states of very high ℓ . The rapid predissociation of molecules in states of low- ℓ purifies this high- ℓ ensemble. The angular momentum barrier that separates NO + ions and electrons creates an extraordinary gap between the plasma states of n ≈ ℓ , with measured n > 200 and penetrating states of ℓ = 0 , 1 , and 2. Evolution to a statistically equilibrated state of N and O atoms cannot occur without Rydberg electron penetration, and this gap blocks relaxation for a millisecond or more. A weak radiofrequency (RF) field drives ℓ -mixing electron collisions that erase this gap, causing wholesale dissipation. Remarkably, a local quantum-state transition promoting an exceedingly small fraction of the molecules in the prethermalized ensemble to a predissociative state also acts with a global effect, driving the entire system to a dissipative equilibrium. Using the Lindblad master equation, we illustrate qualitatively similar dynamics for a toy model of an open quantum system that consists of a localized set of spins on which dissipation acts at a single site as a gateway to equilibrium. Published by the American Physical Society 2025

Wang, Ruoxi↗

Gas–Phase Synthesis of 3–Vinylcyclopropene via the Crossed Beam Reaction of the Methylidyne Radical (CH; X 2 Π) with 1,3–Butadiene (CH 2 CHCHCH 2 ; X 1 Ag)

The crossed molecular beam reactions of the methylidyne radical (CH; X 2 Π) with 1,3-butadiene (CH 2 CHCHCH 2 ; X 1 Ag) along with their (partially) deuterated counterparts were performed at collision energies of 20.8 kJ mol -1 under single collision conditions. Combining our laboratory data with ab initio calculations, we reveal that the methylidyne radical may add barrierlessly to the terminal carbon atom and/or carbon-carbon double bond of 1,3-butadiene, leading to doublet C 5 H 7 intermediates with life times longer than the rotation periods. Furthermore, these collision complexes undergo non-statistical unimolecular decomposition through hydrogen atom emission yielding the cyclic cis- and trans-3-vinyl-cyclopropene products with reaction exoergicities of 119±42 kJ mol -1 . Since this reaction is barrierless, exoergic, and all transition states are located below the energy of the separated reactants, these cyclic C 5 H 6 products are predicted to be accessed even in low-temperature environments, such as in hydrocarbon-rich atmospheres of planets and cold molecular clouds such as TMC-1.

3-vinylcyclopropene↗

Selectivity in gas–liquid interactions: Molecular beam scattering of CD 4 and ND 3 from an aqueous flat liquid jet

The dynamics of polar and nonpolar molecules colliding with an aqueous surface are characterized by scattering molecular beams of deuterated methane and ammonia, CD 4 and ND 3 (E i = 28.9 and 30.3 kJ mol −1 , respectively), from a flat liquid jet of cold salty water (8 m LiBr, 230 K). Translational energy distributions of scattered species collected as a function of collision geometry probe both impulsive scattering (IS) and thermal desorption (TD) mechanisms. Here, we find that CD 4 scattering is dominated by IS and exhibits a super-specular angular distribution. The fraction of TD scattering events is notably smaller for cold salty water than for dodecane, consistent with a higher free energy of solvation for CD 4 in the water jet. In contrast, no scattering signal is seen for ND 3 from the water jet, a result attributed to the high solubility and efficient protonation of ND 3 in liquid water. The IS channel for CD 4 was analyzed using a soft-sphere model, yielding a higher internal energy (E int ) and lower effective surface mass (m eff ) than was seen for Ne/water; the higher value of E int is attributed to rotational excitation of the scattered CD 4 . These findings demonstrate that the outcomes of a gas–liquid collision—scattering trajectory, surface adherence, and energy transfer—are directed at the molecular level by both the gaseous scatterer and liquid surface.

Foreman, Madison M. [University of California, Ber↗

Quasiclassical sampling and Wigner sampling of initial vibrational coordinates and momenta for polyatomic molecules in Monte Carlo molecular dynamics simulations

In a quasiclassical trajectory simulation, the vibrational modes are initialised with quantised vibrational energies, but vibrational phases are sampled by Monte Carlo. This requires an algorithm to assign coordinates and momenta to the various atoms. In this work, we present two methods for implementing this for nonrotating polyatomic molecules, namely, fixed-energy vibrational-state-selected initial conditions and thermal initial conditions. We also present a method for initiating classical trajectories with a ground-state Wigner distribution. These vibrational treatments are sufficient to initialise trajectories for unimolecular processes, and we also show how they can be applied to simulate bimolecular collision processes. The treatments of unimolecular and bimolecular collision processes are available in two Python codes called wigner_state_selected.py and bimolecular_collision.py, respectively, which will generate initial condition files that are recognisable by the SHARC and SHARC-MN computer programs for dynamics calculations. Both codes are available as standalone programs, as well as being included in SHARC-MN, and they will be included in future versions of SHARC. Here, the methods implemented in these codes are mostly also available in the ANT computer program, and those that are not available in ANT will be incorporated in future versions of ANT.

Wigner distribution↗

Collisional-energy-cascade model for nonthermal velocity distributions of neutral atoms in plasmas

Nonthermal velocity distributions with much greater tails than a Maxwellian have been observed for radical atoms in plasmas for a long time. Historically, such velocity distributions have been modeled by a two-temperature Maxwell distribution. In this paper, I propose a model based on collisional energy cascade, which has been studied in the field of granular materials. In the collisional energy cascade, a particle ensemble undergoes energy input at the high-energy region, entropy production by elastic collisions among particles, and energy dissipation. For radical atoms, energy input may be caused by the Franck-Condon energy of molecular dissociation or charge-exchange collision with hot ions, and the input energy is eventually dissipated by collisions with the walls. Here, I show that the steady-state velocity distribution in the collisional energy cascade is approximated by the generalized Mittag-Leffler distribution, which is a one-parameter extension of the Maxwell distribution. This parameter indicates the degree of the nonthermality and is related to the relative importance of energy dissipation over entropy production. This model is compared with a direct molecular dynamics simulation for simplified gaseous systems with energy input, as well as some experimentally observed velocity distributions of light radicals in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High-order cumulants and correlation functions near the critical point from molecular dynamics

We present a systematic investigation of particle-number fluctuations in the crossover region near the critical end point of a first-order phase transition using molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our prior studies to third- and fourth-order cumulants in both coordinate- and momentum-space acceptances and integrated correlation functions (factorial cumulants). We find that, even near the critical point, non-Gaussian cumulants equilibrate on timescales comparable to those of the second-order cumulants, but show stronger finite-size effects. The presence of interactions and of the critical point leads to strong deviations of the cumulants from the ideal-gas baseline in coordinate space; these deviations are expected to persist in momentum space in the presence of collective expansion. In particular, the kurtosis becomes strongly negative, κσ 2 ≪ − 1 , on the crossover side of the critical point. However, this signal is significantly diluted once an efficiency cut used to distinguish protons from baryons is applied, leading to |κσ 2 | ≲ 1 even in the presence of the critical point. We discuss our results in the context of ongoing measurements of proton-number cumulants in heavy-ion collisions in RHIC-BES-II.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Sculpting Non-Machinable Porous Materials

The goal of this project was to establish the scientific basis for subtractive manufacturing (sculpting) of porous materials using ion beams. Using the catcher foil technique, complemented with molecular dynamics simulations, we studied the dependence of the sputtering yield of Ag and Au nanofoams as a function of foam monolith density, ion mass, and the beam incidence angle, and compared these dependencies to those measured for full-density solids. We found that the foam sputtering yield was considerably lower than for the full-density solid. The sputtering yield increased with beam incidence angle, and the open structure of the foam suppressed the development of collision cascades. We also found that the sputtering yield decreased with increasing foam density. The simulation revealed a significant redeposition of sputtered atoms in the foam.

36 MATERIALS SCIENCE↗

Enhanced dilepton emission from a phase transition in dense matter

It is demonstrated that the presence of a phase transition in heavy-ion collisions, at beam energies that probe dense quantum chromodynamics (QCD) matter, leads to a significant enhancement of the dilepton yield at low invariant mass and requires the embedding of low beam energies per produced pion due to the extended emission time. In addition, the temperature of low-mass dileptons shows a modest decrease due to the mixed phase. The emission of dileptons in the SIS18–SIS100 beam energy range is studied by augmenting the ultra-relativistic quantum molecular dynamics (UrQMD) transport model with a realistic density-dependent equation of state, as well as two different phase transitions. This is achieved by extending the molecular dynamics interaction part of the UrQMD model to a density-dependent interaction potential with a high-density minimum, leading to a phase transition and metastable coexisting high-density states. Together with a high-precision measurement, these simulations will be able to constrain the existence of a phase transition in QCD up to densities of several times nuclear saturation density.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Radiation-Induced Defect Formation Kinetics in Inconel–Cu Multimetallic Layered Composites

This study investigates the stability of Inconel–Cu Multimetallic Layered Composites (MMLCs) in nuclear reactor applications using Molecular Dynamics simulations. The focus is on understanding the underlying mechanisms governing the properties of MMLCs for advanced nuclear reactors, specifically, the mechanochemistry of the interface between Inconel and copper alloys. The selection of Inconel–Cu MMLCs is primarily due to copper’s superior thermal conductivity, enhancing heat management within reactors by preventing hotspots and ensuring uniform temperature distribution. This research examines Incoloy 800H and two Inconel variants (718 and 625), assessing their stability at 1000 K after exposure to 10 keV collision cascades up to 0.12 dpa. Notable findings include defect clustering on the {1 2 0} family of planes of Inconel and Cu, with Stacking Faults and Lomer–Cottrell locks on the Inconel side.

Ramesh, Rajesh↗

Statistical analysis of correlations in the x-ray induced Coulomb explosion of iodopyridine

Coulomb explosion imaging is a promising experimental tool to study individual molecules. Here, in this work simulation data for the x-ray induced Coulomb explosion of 2-iodopyridine is analyzed and the involved fragmentation dynamics are described. It is found that particular final ion momenta show correlations that reflect a collision of two atoms during the explosion. Variations of the proton kinetic energies can be associated with variations in the charge build-up time. Covariances of the forces during the explosion can be utilized to simplify the description of the dynamics in reduced dimensionality using only four collective coordinates.

74 ATOMIC AND MOLECULAR PHYSICS↗

New experimental approach to understanding the chemical reactivity of oxide surfaces

Metal oxides have been an attractive option for a range of applications, including hydrogen sensors, microelectronics, and catalysis, due to their reactivity and tunability. The properties of metal oxides can vary greatly on their precise surface structure; however, few surface science techniques can achieve atomistic-level determinations of surface structure, and fewer yet can do so for insulator surfaces. Low energy ion beam analysis offers a potential insulator-compatible solution to characterizing the surface structure of metal oxides. As a feasibility study, we apply low energy ion beam analysis to investigate the surface structure of a magnetite single crystal, Fe 3 O 4 (100). We obtain multi-angle maps using both forward-scattering low energy ion scattering (LEIS) and backscattering impact-collision ion scattering spectroscopy (ICISS). Both sets of experimental maps have intensity patterns that reflect the symmetries of the Fe 3 O 4 (100) surface structure. However, analytical interpretation of these intensity patterns to extract details of the surface structure is significantly more complex than previous LEIS and ICISS structural studies of one-component metal crystals, which had far more symmetries to exploit. To gain further insight into the surface structure, we model our experimental measurements with ion-trajectory tracing simulations using molecular dynamics. Our simulations provide a qualitative indication that our experimental measurements agree better with a subsurface cation vacancy model than a distorted bulk model.

36 MATERIALS SCIENCE↗

Strain-Induced Growth of Twisted Bilayers during the Coalescence of Monolayer MoS 2 Crystals

Tailoring the grain boundaries (GBs) and twist angles between two-dimensional (2D) crystals are two crucial synthetic challenges to deterministically enable envisioned applications such as moiré excitons, emerging magnetism, or single-photon emission. Here, we reveal how twisted 2D bilayers can be synthesized from the collision and coalescence of two growing monolayer MoS 2 crystals during chemical vapor deposition. The twisted bilayer (TB) moiré angles are found to preserve the misorientation angle (θ) of the colliding crystals. The shapes of the TB regions are rationalized by a kink propagation model that predicts the GB formed by the coalescing crystals. Optical spectroscopy measurements reveal a θ-dependent long-range strain in crystals with stitched grain boundaries and a sharp (θ > 20°) threshold for the appearance of TBs, which relieves this strain. Reactive molecular dynamics simulations explain this strain from the continued growth of the crystals during coalescence due to the insertion of atoms at unsaturated defects along the GB, a process that self-terminates when the defects become saturated. The simulations also reproduce atomic-resolution electron microscopy observations of faceting along the GB, which is shown to arise from the growth-induced long-range strain. These facets align with the axes of the colliding crystals to provide favorable nucleation sites for second-layer growth of a TB with twist angles that preserve the misorientation angle θ. This interplay between strain generation and aligned nucleation provides a synthetic pathway for the growth of TBs with deterministic angles.

2D materials↗

Attosecond delays in X-ray molecular ionization

The photoelectric effect is not truly instantaneous, but exhibits attosecond delays that can reveal complex molecular dynamics. Sub-femtosecond duration light pulses provide the requisite tools to resolve the dynamics of photoionization. Accordingly, the past decade has produced a large volume of work on photoionization delays following single photon absorption of an extreme ultraviolet (XUV) photon. However, the measurement of time-resolved core-level photoionization remained out of reach. The required x-ray photon energies needed for core-level photoionization were not available with attosecond tabletop sources. We have now measured the x-ray photoemission delay of core-level electrons, and here report unexpectedly large delays, ranging up to 700 attoseconds in NO near the oxygen K-shell threshold. These measurements exploit attosecond soft x-ray pulses from a free-electron laser (XFEL) to scan across the entire region near the K-shell threshold. Furthermore, we find the delay spectrum is richly modulated, suggesting several contributions including transient trapping of the photoelectron due to shape resonances, collisions with the Auger-Meitner electron that is emitted in the rapid non-radiative relaxation of the molecule, and multi-electron scattering effects. Finally, the results demonstrate how x-ray attosecond experiments, supported by comprehensive theoretical modelling, can unravel the complex correlated dynamics of core-level photoionization.

74 ATOMIC AND MOLECULAR PHYSICS↗

Pulsed laser synthesis of mesoporous metal chalcogenide thin films

Mesoporous films of the metal chalcogenide B-FeSe were grown on MgO substrates by KrF pulsed laser deposition (PLD) in an argon background. At 100 mTorr, gated intensified charge-coupled device imaging and ion probe measurements showed that the plasma plume responsible for crystal growth initially comprised three components, with distinct expansion velocities. Plume interactions with the substrate heater and ablation target gave rise to complex dynamics, including collisions between the charged leading edge—rebounding between the substrate and the target—and slower-moving species in the plume interior. Film growth was dominated by species with kinetic energies ≤0.5 eV/atom. X-ray reflectivity revealed that films grown in this environment—with a substrate temperature of 350 ° C, a laser fluence of 1.0 J cm −2 , and a 7.5 mm 2 spot area—formed a porous framework with 15% porosity. Atomic force microscopy showed surface features that suggest pore sizes below 100 nm. X-ray diffraction indicated that the porous films were epitaxial with respect to the substrate and likely grew by oriented-attachment of gas-phase molecular clusters or very small nanoparticles, in contrast to the conventional epitaxy of vacuum films from atomic constituents. The in-plane orientation of the mesoporous films was B-FeSe [100]||[110] MgO, attributed to the soft landing of pre-formed crystallites on the MgO substrates, where protruding Se rows of B-FeSe aligned with corrugations of the MgO surface. In conclusion, this work implies that growth of candidate electrocatalyst materials by PLD in inert gas background may allow mesoporous frameworks with a single crystallographic orientation that expose specific crystal facets for electrochemical reactions and active site engineering.

cluster epitaxy↗

Application of machine learning in the determination of impact parameter in the 132 Sn+ 124 Sn system

Here, 132 Sn + 124 Sn collisions at a beam energy of 270 MeV/nucleon were performed at the Radioactive Isotope Beam Factory (RIBF) in RIKEN to investigate the nuclear equation of state. Reconstructing the impact parameter is one of the important tasks in the experiment as it relates to many observable. In this work, we employ three commonly used algorithms in machine learning, the artificial neural network (ANN), the convolutional neural network (CNN), and the light gradient boosting machine (LightGBM), to determine the impact parameter by analyzing either the charged particle spectra or several features simulated with events from the ultrarelativistic quantum molecular dynamics (UrQMD) model. To closely imitate experimental data and investigate the generalizability of the trained machine learning algorithms, incompressibility of nuclear equation of state and the in-medium nucleon-nucleon cross sections are varied in the UrQMD model to generate the training data. The mean absolute error Δb between the true and the predicted impact parameter is smaller than 0.45 fm if training and testing sets are sampled from the UrQMD model with the same parameter set. However, if training and testing sets are sampled with different parameter sets, Δb would increase to 0.8 fm. The generalizability of the trained machine learning algorithms suggests that these machine learning algorithms can be used reliably to reconstruct the impact parameter in experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamics of Product Branching from Radical Intermediates in Elementary Combustion Reactions (Final Technical Report)

While the total rate constant for many elementary reactions is well-characterized, understanding the product branching in complex reactions presents a formidable challenge. To gain an incisive probe of such reactions, our experiments investigate the dynamics of the product channels that arise from transient radical intermediates along the bimolecular reaction coordinates. Our recent published work uses the methodology developed in my group in the last fifteen years, using both imaging and scattering apparatuses. The experiments generate a particular isomeric form of an unstable radical intermediate along a bimolecular reaction coordinate and study the branching between the ensuing product channels of the energized radical as a function of its internal rotational and vibrational energy under collision-less conditions. The experiments use a combination of: 1) measurement of product velocity and angular distributions in a crossed laser-molecular beam apparatus, with electron bombardment detection in my lab in Chicago or 2) with tunable vacuum ultraviolet photoionization detection at Taiwan's National Synchrotron Radiation Research Center (NSRRC), and 3) velocity map imaging using state-selective REMPI and single photon VUV ionization of radical intermediates and reaction products. We also employ tunable VUV photoionization detection in our imaging apparatus, using difference frequency four-wave mixing to produce photoionization light tunable from 8 to 10.8 eV. This year of support funded our work to photolytically generate the radical intermediate important in the OH + propene reaction when OH adds to an end carbon. Our other major result this past year was on the photoionization detection of vinoxy radicals. Our experiments at the NSRRC definitively show that vinoxy dissociatively ionizes to both CH 3 + and HCO + , and establishes the absolute partial photoionization cross section of vinoxy to CH 3 + at 10.5 and 11.44 eV. This should allow experimentalists detecting vinoxy radicals produced in bimolecular reactions such as OH + propene to get good branching fractions to that product channel. Finally we initiated studies of the photodissociation of 1-bromo-2-chloroethane and 1,1-bromochloroethane at 193 nm at the NSRRC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimized collision-specific parameters for binary mixtures of nitrogen, oxygen, argon, and helium

Recently proposed collision-specific parameters for direct simulation Monte Carlo simulations are tested for binary mixtures of nitrogen, oxygen, and argon. Near ambient conditions, the traditional collision-averaged parameters are highly accurate, whereas the collision-specific parameters are not. The simulated transport using the collision-averaged parameters for mixtures with helium, however, is found to be inaccurate. Therefore, we propose a novel method to determine molecular parameters by combining the Chapman–Enskog theory with empirical mixing rules and experimental data. The optimized parameters are highly accurate for the binary mixtures of nitrogen, oxygen, and argon and greatly improve the simulated transport for the helium mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗