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 OSTI2023

Hydrogen irradiation-driven computational surface chemistry of lithium oxide and hydroxide

We have investigated, using molecular dynamics, the surface chemistry of hydrogen incident on the amorphous and crystalline lithium oxide and lithium hydroxide surfaces upon being slowed down by a collision cascade and retained in the amorphous surface of either Li 2 O or LiOH. We looked for the bonding of H to the resident structures in the surface to understand a possible chain of chemical reactions that can lead to surface transformation upon H atom impact. Our findings, using Density-Functional Theory (DFT) trained ReaxFF force field/electronegativity equalization method potentials, stress the importance of inclusion of polarization in the dynamics of a Li–O–H system, which is also illustrated by DFT energy minimization and quantum–classical molecular dynamics using tight binding DFT. The resulting polar-covalent chemistry of the studied systems is complex and very sensitive to the instantaneous positions of all atoms as well as the ratio of concentrations of various resident atoms in the surface.

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

Diffusion coefficients in the envelopes of white dwarfs

The diffusion of elements is a key process in understanding the unusual surface composition of white dwarfs and their spectral evolution. The diffusion coefficients of Paquette et al. have been widely used to model diffusion in white dwarfs. In this work, we perform new calculations of the coefficients of interdiffusion and ionic thermal diffusion with (1) a more advanced model that uses a recent modification of the calculation of the collision integrals that is more suitable for the partially ionized, partially degenerate, and moderately coupled plasma and (2) classical molecular dynamics. The coefficients are evaluated for silicon and calcium in white dwarf envelopes of hydrogen and helium. A comparison of our results with Paquette et al. shows that the latter systematically underestimates the coefficient of interdiffusion yet provides reliable estimates for the relatively weakly coupled plasmas found in nearly all types of stars, as well as in white dwarfs with hydrogen envelopes. In white dwarfs with cool helium envelopes (${T}_{\mathrm{eff}}$ < 15,000 K), the difference grows to more than a factor of two. We also explored the effect of the ionization model used to determine the charges of the ions and found that it can be a substantial source of discrepancy between different calculations. Finally, we consider the relative diffusion timescales of Si and Ca in the context of the pollution of white dwarf photospheres by accreted planetesimals and find factor of ≳3 differences between calculations based on Paquette et al. and our model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Dislocation dynamics in polycrystalline materials

An extension of the dislocation dynamics method to polycrystalline materials is presented with emphasis on handling interactions between dislocations and grain boundaries. The key new features of the method include: (1) a grain boundary detection algorithm, (2) new dislocation emission criteria based on power dissipation and area growth, (3) mobility of dislocations at grain boundaries, and (4) extension of topological operations e.g. mesh adaption, collision, and dissociation to handle dislocations interacting with grain boundaries. Furthermore, verification and validation examples are executed to compare our calculations with existing analytic methods, molecular dynamics simulations, and experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multichannel dynamics in the OH+ n -butane reaction revealed by crossed-beam slice imaging and quasiclassical trajectory calculations

Multidimensional reactions present various channels that can exhibit very different dynamics and give products of varying subsequent reactivity. Here, we present a combination of experiment and theory to reveal the dynamics of hydrogen abstraction by OH radical at primary and secondary sites in n-butane at a collision energy of 8 kcal/mol. Crossed molecular beam slice imaging experiments unequivocally probe the secondary abstraction channel showing backward angular distributions with mild energy release to product translation, which are accurately captured by trajectory calculations using a specific-reaction-parameter Hamiltonian. Additionally, experiments containing both reaction channels indicate a less marked backward character in the angular distribution, whose origin is shown by trajectory calculations to appear as an evolution toward more sideways scattering from the secondary to primary channel. While the two channels have markedly different angular distributions, their energy release is largely comparable, showing ample energy release into the water product. The synergistic combination of crossed-beam imaging and trajectories opens the door to detailed reaction-dynamics studies of chemical reactions with ever-increasing complexity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2022

Silicon nanocluster anion-argon cation recombination via hybrid continuum-molecular dynamics calculations

Prediction of the charge distribution on particles in an aerosol is critical not only in electrical mobility based characterization methods, but also in understanding the role that charging plays in particle growth in gas phase synthesis reactors. The latter is particularly important in non-thermal plasma synthesis reactors, wherein nanoclusters form and grow from vapor phase precursors in a high electron and high ion density environment. In plasmas, because free electrons are much less massive and much more mobile than positive gas ions, nanocluster charge distributions are biased negative. Furthermore, neutral or even positively charged nanoclusters may exist, depending on the rate of nanocluster-gas ion recombination, and nanocluster-nanocluster collisions may greatly contribute to nanocluster growth if not Coulombically suppressed.

42 ENGINEERING↗
DOE OSTI2026

Dynamics of iodine geminate recombination in supercritical xenon solvent: Caging effect

Understanding the dynamics of chemical reactions in solutions is vital, as their rates and kinetics are significantly affected by the solvent environment. Supercritical solvents offer extensive applications in chemical reactions by enabling the manipulation of the solution environment. Here, in this study, we investigate the geminate recombination of iodine in a supercritical xenon solvent by using ReaxFF-based molecular dynamics simulations. Our findings reveal that the highest iodine recombination rate occurs near supercritical conditions, while lower-pressure conditions lead to reduced collision rates and unstable recombination, and higher-pressure conditions hinder iodine diffusion, resulting in a lower recombination rate. Our analysis shows that the xenon local density at the time of recombination is at least 2.5 times higher than the global density, confirming the presence of xenon clusters surrounding the Iodine atoms. This observation is further supported by coordination number analysis, which confirms an elevated xenon local density during recombination. In addition, the correlation between the total energy of xenon atoms within a cluster and recombined iodine atoms underscores the kinetic energy transfer process, validating the occurrence of geminate recombination. The excess kinetic energy from the recombining iodine atoms is transferred to the surrounding xenon atoms. Our examination of geminate recombination demonstrates that iodine atoms confined within xenon clusters—whether through manual insertion of atoms or the fast dissociation of an iodine molecule within xenon clusters—are more likely to recombine as primary geminate recombination. However, extending the iodine molecule dissociation time allows iodine atoms to diffuse out of the cluster, and the recombination to shift toward secondary geminate recombination.

Cage effect↗
DOE OSTI2024

Multinucleon transfer in 48 Ti + 48 Ti collisions at 11.5 MeV/nucleon

In this work, experimental data of projectile fragment distributions from a preliminary experimental study of the reaction of 11.5 MeV/nucleon 48 Ti on 48 Ti analyzed with the MARS recoil separator at the Cyclotron Institute of Texas A&M University are presented. Production cross sections, momentum distributions and excitation-energy distributions are extracted and compared with calculations with the Deep Inelastic Transfer (DIT) model and the Constrained Molecular Dynamics model (CoMD). Despite the limited extent of the present dataset, an overall agreement of the models with the data is found pointing at the prevailing nucleon-exchange character of the collisions at this energy.

Souliotis, Georgios↗
DOE OSTI2024

Spinodal enhancement of fluctuations in nucleus-nucleus collisions

Subensemble Acceptance Method (SAM) [1, 2] is an essential link between measured event-by-event fluctuations and their grand canonical theoretical predictions such as lattice QCD. The method allows quantifying the global conservation law effects in fluctuations. In its basic formulation, SAM requires a sufficiently large system such as created in central nucleus-nucleus collisions and sufficient space-momentum correlations. Directly in the spinodal region of the First Order Phase Transition (FOPT) different approximations should be used that account for finite size effects. Thus, we present the generalization of SAM applicable in both the pure phases, metastable and unstable regions of the phase diagram [3]. Obtained analytic formulas indicate the enhancement of fluctuations due to crossing the spinodal region of FOPT and are tested using molecular dynamics simulations. A rather good agreement is observed. Using transport model calculations with interaction potential we show that the spinodal enhancement of fluctuations survives till the later stages of collision via the memory effect [4]. However, at low collision energies the space-momentum correlation is not strong enough for this signal to be transferred to second and third order cumulants measured in momentum subspace. This result agrees well with recent HADES data on proton number fluctuations at $\sqrt{S_{NN}}$ = 2.4 GeV which are found to be consistent with the binomial momentum space acceptance [5].

Poberezhnyuk, Roman↗
DOE OSTI2025

From the laboratory to space: unveiling isomeric diversity of C 5 H 2 in the reaction of tricarbon (C 3 , X 1 Σ g + ) with the vinyl radical (C 2 H 3 , X 2 A′)

By connecting laboratory dynamics with cosmic observables, this work highlights the critical role of reactions between highly reactive species in shaping the molecular inventory of the interstellar medium and opens new windows into the spectroscopically elusive corners of astrochemical complexity. The gas phase formation of distinct C 5 H 2 isomers is explored through the bimolecular reaction of tricarbon (C 3 , X 1 Σ + g ) with the vinyl radical (C 2 H 3 , X 2 A′) at a collision energy of 44 ± 1 kJ mol −1 employing the crossed molecular beam technique augmented by electronic structure and Rice–Ramsperger–Kassel–Marcus (RRKM) calculations. This barrierless and exoergic reaction follows indirect dynamics and is initiated by the addition of tricarbon to the radical center of the vinyl radical forming a Cs symmetric doublet collisional complex (CCCCHCH 2 ). Subsequent low-barrier isomerization steps culminate in the resonantly stabilized 2,4-pentadiynyl-1 radical (CHCCCCH 2 ), which decomposes via atomic hydrogen loss. Statistical calculations identify linear, triplet pentadiynylidene (p2, X 3 Σ − g ) as the dominant product, while singlet carbenes ethynylcyclopropenylidene (p1, X 1 A′), pentatetraenylidene (p3, X 1 A 1 ), and ethynylpropadienylidene (p4, X 1 A′) are formed with lower branching ratios. The least stable isomer, 2-cyclopropen-1-ylidenethenylidene (‘eiffelene’; p5, X 1 A 1 ), remains thermodynamically feasible, but exhibits negligible branching ratios. Two isomers detected in TMC-1 to date (p1 and p3) possess significant dipole moments making them amenable to radio telescopic observations, whereas linear pentadiynylidene (p2; D ∞h ) is only traceable via infrared spectroscopy or through its cyanopentadiynylidene derivative (HCCCCCCN). This study highlights the isomer diversity accessed in the low temperature hydrocarbon chemistry of barrierless and exoergic bimolecular reactions involving two unstable, reactants in cold molecular clouds.

Medvedkov, Iakov A. [University of Hawai'i at Mano↗
DOE OSTI2020

Stratification of polymer mixtures in drying droplets: Hydrodynamics and diffusion

In this work, we study the evaporation-induced stratification of a mixture of short and long polymer chains in a drying droplet using molecular simulations. We systematically investigate the effects of hydrodynamic interactions (HI) on this process by comparing hybrid simulations accounting for HI between polymers through the multiparticle collision dynamics technique with free-draining Langevin dynamics simulations neglecting the same. We find that the dried supraparticle morphologies are homogeneous when HI are included but are stratified in core–shell structures (with the short polymers forming the shell) when HI are neglected. The simulation methodology unambiguously attributes this difference to the treatment of the solvent in the two models. We rationalize the presence (or absence) of stratification by measuring phenomenological multicomponent diffusion coefficients for the polymer mixtures. The diffusion coefficients show the importance of not only solvent backflow but also HI between polymers in controlling the dried supraparticle morphology.

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

A crossed molecular beams and computational study on the unusual reactivity of banana bonds of cyclopropane (c-C 3 H 6 ; X 1 A$^{'}_{1}$ ) through insertion by ground state carbon atoms (C( 3 P j ))

The mechanism and chemical dynamics of the reaction of ground electronic state atomic carbon C( 3 P j ) with cyclopropane c-C 3 H 6 (X 1 A$^{'}_{1}$) have been explored by combining crossed molecular beams experiments with electronic structure calculations of the pertinent triplet C 4 H 6 potential energy surface and statistical computations of product branching ratios under single-collision conditions. The experimental findings suggest that the reaction proceeds via indirect scattering dynamics through triplet C 4 H 6 reaction intermediate(s) leading to C 4 H 5 product(s) plus atomic hydrogen via a tight exit transition state, with the overall reaction exoergicity evaluated as 231 ± 52 kJ mol –1 . The calculations indicate that C( 3 P j ) can easily insert into one of the three equivalent C–C ‘banana’ bonds of cyclopropane overcoming a low barrier of only 2 kJ mol –1 following the formation of a van der Waals reactant complex stabilized by 15 kJ mol –1 . The carbon atom insertion into one of the six C–H bonds is also feasible via a slightly higher barrier of 5 kJ mol –1 . Here, these results highlight an unusual reactivity of cyclopropane's banana C–C bonds, which behave more like unsaturated C–C bonds with a π-character than saturated σ C–C bonds, which are known to be generally unreactive toward the ground electronic state atomic carbon such as in ethane (C 2 H 6 ). The statistical theory predicts the overall product branching ratios at the experimental collision energy as 50% for 1-butyn-4-yl, 33% for 1,3-butadien-2-yl, i-C 4 H 5 , and 11% for 1,3-butadien-1-yl, n-C 4 H 5 , with i-C 4 H 5 (230 kJ mol –1 below the reactants) favored by the C–C insertion providing the best match with the experimentally observed reaction exoergicity. The C( 3 P j ) + c-C 3 H 6 reaction is predicted to be a source of C 4 H 5 radicals under the conditions where its low entrance barriers can be overcome, such as in planetary atmospheres or in circumstellar envelopes but not in cold molecular clouds. Both i- and n-C 4 H 5 can further react with acetylene eventually producing the first aromatic ring and hence, the reaction of the atomic carbon with c-C 3 H 6 can be considered as an initial step toward the formation of benzene.

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

Exploring the chemical dynamics of phenanthrene (C 14 H 10 ) formation via the bimolecular gas-phase reaction of the phenylethynyl radical (C 6 H 5 CC) with benzene (C 6 H 6 )

The exploration of the fundamental formation mechanisms of polycyclic aromatic hydrocarbons (PAHs) is crucial for the understanding of molecular mass growth processes leading to two- and three-dimensional carbonaceous nanostructures (nanosheets, graphenes, nanotubes, buckyballs) in extraterrestrial environments (circumstellar envelopes, planetary nebulae, molecular clouds) and combustion systems. While key studies have been conducted exploiting traditional, high-temperature mechanisms such as the hydrogen abstraction–acetylene addition (HACA) and phenyl addition–dehydrocyclization (PAC) pathways, the complexity of extreme environments highlights the necessity of investigating chemically diverse mass growth reaction mechanisms leading to PAHs. Employing the crossed molecular beams technique coupled with electronic structure calculations, we report on the gas-phase synthesis of phenanthrene (C 14 H 10 )—a three-ring, 14π benzenoid PAH—via a phenylethynyl addition–cyclization–aromatization mechanism, featuring bimolecular reactions of the phenylethynyl radical (C 6 H 5 CC, X 2 A 1 ) with benzene (C 6 H 6 ) under single collision conditions. The dynamics involve a phenylethynyl radical addition to benzene without entrance barrier leading eventually to phenanthrene via indirect scattering dynamics through C 14 H 11 intermediates. The barrierless nature of reaction allows rapid access to phenanthrene in low-temperature environments such as cold molecular clouds which can reach temperatures as low as 10 K. Furthermore, this mechanism constitutes a unique, low-temperature framework for the formation of PAHs as building blocks in molecular mass growth processes to carbonaceous nanostructures in extraterrestrial environments thus affording critical insight into the low-temperature hydrocarbon chemistry in our universe.

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

Beam energy dependence of net- Λ fluctuations measured by the STAR experiment at the BNL Relativistic Heavy Ion Collider

The measurements of particle multiplicity distributions have generated considerable interest in understanding the fluctuations of conserved quantum numbers in the quantum chromodynamics (QCD) hadronization regime, in particular near a possible critical point and near the chemical freeze-out. Net-protons and net-kaons have been used as proxies for the net-baryon number and net-strangeness, respectively. We report the measurement of efficiency- and centrality-bin width-corrected cumulant ratios (C 2 /C 1 , C 3 /C 2 ) of net- Λ distributions, in the context of both strangeness and baryon number conservation, as a function of collision energy, centrality, and rapidity. The results are for Au + Au collisions at five beam energies ( √s NN = 19.6 , 27, 39, 62.4, and 200 GeV) recorded with the Solenoidal Tracker at RHIC (STAR). We compare our results to the Poisson and negative binomial (NBD) expectations, as well as to ultrarelativistic quantum molecular dynamics (UrQMD) and hadron resonance gas (HRG) model predictions. Both NBD and Poisson baselines agree with data within the statistical and systematic uncertainties. UrQMD describes the measured net- Λ C 1 and C 3 at 200 GeV reasonably well but deviates from C 2 , and the deviation increases as a function of collision energy. The ratios of the measured cumulants show no features of critical fluctuations. The chemical freeze-out temperatures extracted from a recent HRG calculation, which was successfully used to describe the net-proton, net-kaon, and net-charge data, indicate Λ freeze-out conditions similar to those of kaons. However, large deviations are found when comparing with temperatures obtained from net-proton fluctuations. The net- Λ cumulants show a weak but finite dependence on the rapidity coverage in the acceptance of the detector, which can be attributed to quantum number conservation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Elucidating the chemical dynamics of the elementary reactions of the 1-propynyl radical (CH 3 CC; X 2 A 1 ) with 2-methylpropene ((CH 3 ) 2 CCH 2 ; X 1 A 1 )

Exploiting the crossed molecular beam technique, we studied the reaction of the 1-propynyl radical (CH 3 CC; X 2 A 1 ) with 2-methylpropene (isobutylene; (CH 3 ) 2 CCH 2 ; X 1 A 1 ) at a collision energy of 38 ± 3 kJ mol –1 . The experimental results along with ab initio and statistical calculations revealed that the reaction has no entrance barrier and proceeds via indirect scattering dynamics involving C 7 H 11 intermediates with lifetimes longer than their rotation period(s). The reaction is initiated by the addition of the 1-propynyl radical with its radical center to the π-electron density at the C1 and/or C2 position in 2-methylpropene. Further, the C 7 H 11 intermediate formed from the C1 addition either emits atomic hydrogen or undergoes isomerization via [1,2-H] shift from the CH 3 or CH 2 group prior to atomic hydrogen loss preferentially leading to 1,2,4-trimethylvinylacetylene (2-methylhex-2-en-4-yne) as the dominant product. The molecular structures of the collisional complexes promote hydrogen atom loss channels. RRKM results show that hydrogen elimination channels dominate in this reaction, with a branching ratio exceeding 70%. Since the reaction of the 1-propynyl radical with 2-methylpropene has no entrance barrier, is exoergic, and all transition states involved are located below the energy of the separated reactants, bimolecular collisions are feasible to form trimethylsubstituted 1,3-enyne (p1) via a single collision event even at temperatures as low as 10 K prevailing in cold molecular clouds such as G+0.693. As a result, the formation of trimethylsubstituted vinylacetylene could serve as the starting point of fundamental molecular mass growth processes leading to di- and trimethylsubstituted naphthalenes via the HAVA mechanism.

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

Collective behavior of “flexicles”

In recent years the functionality of synthetic active microparticles has edged even closer to that of their biological counterparts. However, we still lack the understanding needed to recreate at the microscale key features of autonomous behavior exhibited by microorganisms or swarms of macroscopic robots. In this study, we propose a model for a three-dimensional deformable cellular composite particle consisting of self-propelled rod-shaped colloids confined within a flexible vesicle—representing a superstructure we call a “flexicle” that couples particle deformation to the internal dynamics of the internal active components. Using molecular dynamics simulations, we investigate the collective behavior of dense systems composed of many flexicles. We show that individual flexicles exhibit shape changes upon collisions with other flexicles that lead to rearrangements of the internal active rods, which slows flexicle motion. This shape deformability gives rise to a diverse set of motility-induced phase separation phenomena and the spontaneous flow of flexicles reminiscent of the migration of cells in dense tissues. Our findings establish a foundation for designing responsive, cell-like active particles and developing strategies for controlling swarm migration and other autonomous swarm behaviors at cellular and colloidal scales.

Science & Technology - Other Topics↗
DOE OSTI2023

Semirelativistic antisymmetrized molecular dynamics for the production of energetic neutrons in intermediate-energy heavy-ion reactions

Relativistic corrections have been made in the nonrelativistic antisymmetrized molecular dynamics (amd) simulations to apply to the high-energy neutron production in the 12 C + 12 C and 16 O + 12 C collisions at incident energies of 290 and 400 MeV/nucleon. The corrections are made in kinematics alone and neither nucleon-nucleon inelastic scatterings nor meson productions are taken into account, and amd with the relativistic corrections is called semirelativistic amd. The three-nucleon collision (3NC) and Fermi boost in the collision processes are taken into account in the nonrelativistic AMD. Since the relativistic corrections tend to compensate in each other, the difference between the semirelativistic and nonrelativistic results become small. High-energy tails of the available experimental neutron double-differential cross sections, especially at larger angles, are well reproduced by amd with the 3NC term both with nonrelativistic and semirelativistic simulations. Furthermore, these results indicate that the high-energy neutrons are dominantly produced by the 3NC process in this incident energy range.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Probability Density Function for the spatial and intensity distribution of neutron-induced defects in Silicon

The ability to model semiconductor device degradation under neutron irradiation depends upon having a robust modeling capability for the neutron-induced collision cascades as well as a means to analytically fit the resulting probability distributions of defect production and ionizing energy deposition for purposes of extrapolation to low-probability, high-consequence scenarios. In this paper, the widely-utilized binary collision approximation codes MARLOWE and SRIM are deployed in conjunction with a critical examination of their parameterizations as benchmarked against higher-fidelity molecular dynamics simulations. A simple 3-parameter form described by the Generalized Logistic Distribution is shown to be a good fit to Frenkel pair and ionization intensity distributions in bulk silicon. The BCA codes are then applied to simulate cascades in 5 nm layers of a representative gate-all-around nanosheet transistor, where joint probability distributions of threshold levels of damage to multiple layers are evaluated.

36 MATERIALS SCIENCE↗