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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

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

Energy-Energy Correlators for Precision QCD

In this contribution to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021) we review recent progress in the evaluation and application of the Energy-Energy Correlator (EEC) event shape observable in $e^+e^⁻$ annihilation, hadronic collisions, and deep inelastic scattering. The importance of EEC as a precision probe of the perturbative and non perturbative aspects of QCD dynamics is emphasized. It can be used to extract the strong coupling constant and to constrain TMD distribution functions. Closely related energy-correlation shape variables have also been used to tag boosted objects produced in high energy collisions. The opportunities to study EEC at the future Electron-Ion Collider are also highlighted.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Multinucleon transfer in the interaction of 977 MeV and 1143 MeV Hg 204 with Pb 208

A previous study of symmetric collisions of massive nuclei has shown that current models of multinucleon transfer (MNT) reactions do not adequately describe the transfer product yields. To gain further insight into this problem, we have measured the yields of MNT products in the interaction of 977 (E/A = 4.79 MeV) and 1143 MeV (E/A = 5.60 MeV) 204 Hg with 208 Pb. We find that the yield of multinucleon transfer products are similar in these two reactions and are substantially lower than those observed in the reaction of 1257 MeV (E/A = 6.16 MeV) 204 Hg+ 198 Pt. We compare our measurements with the predictions of the GRAZING-F, dinuclear systems (DNS), and improved quantum molecular dynamics (ImQMD) models. We report that for the observed isotopes of the elements Au, Hg, Tl, Pb, and Bi, the measured values of the MNT cross sections are orders of magnitude larger than the predicted values. Furthermore, the various models predict the formation of nuclides near the N = 126 shell, which are not observed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Initial position optimization in molecular dynamics simulations for a Coulomb system

A new algorithm for molecular dynamics (MD) simulations is developed to optimize plasma particle distributions at given initial temperatures. By combining velocity scaling and reassignment, the method effectively eliminates the initial rise and oscillation in temperatures observed with randomly distributed positions. These rises and oscillations are undesired numerical artifacts observed in conventional plasma MD simulations, arising from unoptimized particle positions. The algorithm demonstrates temperature relaxation without initial rises or oscillations, as well as precise flow velocity relaxation, enabling accurate measurement of relaxation times. The code is accelerated using graphics processing units for parallel processing, enhancing the study of plasma dynamics. The proposed method for distributing physically valid particles in MD simulations enables accurate studies of intrinsic collision processes in plasmas, including the dynamics of strongly coupled plasmas, plasma–wave interactions, and transport phenomena in magnetized plasmas. The paper concludes with a discussion of potential applications and future enhancements to the algorithm.

Jo, Jawon (ORCID:0009000924193285)↗
DOE OSTI2020

The charge reduction rate for multiply charged polymer ions via ion–ion recombination at atmospheric pressure

The charge reduction of multiply charged macromolecular ions via recombination with small ions in the gas phase is commonly employed to modulate the charge on macromolecules prior to mass spectrometric and mobility analyses. We employ a recently developed continuum-Molecular Dynamics (MD) calculation approach to determine the recombination rate coefficient of multiply charged (1 to 7 excess positive charged) polyethylene glycol ions (mass of 4600 Da) with smaller singly charged anions, modeled as NO 2 – ions. The continuum-MD approach accounts explicitly for the influence of the background gas on the recombination process, accounts explicitly for ion translational, vibrational, and rotational motion, and enables recombination rate coefficient calculation in nitrogen near atmospheric pressure, wherein neither low pressure nor high pressure recombination theories are strictly applicable. Continuum-MD simulations yield recombination rate coefficients near 3.9 × 10 –14 m 3 s –1 for singly charged ions, increasing to 3.0 × 10 –11 m 3 s –1 for the +7 ion. Pre-existing collision rate coefficient expressions for rigid ions are found to be within a factor of 2–5 of calculations for all charge states, but their use requires knowledge of an appropriate collision distance, which is not well-defined for flexible polymer ions. Continuum-MD-inferred rate coefficients are incorporated into a model of charge reduction, and the charge state distribution versus anion concentration determined with it is compared to charge reduction measurements made via atmospheric pressure differential mobility analysis. Furthermore, good agreement is observed between simulations and experiments; although results are extremely sensitive to the recombination rate coefficients, experimental results are bound by models utilizing rates within a factor of 2 (0.5–2.0×) of the continuum-MD rates.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Study of the interaction between $Ξ$ baryons and light mesons via femtoscopy at the LHC

Meson-baryon systems with strangeness content provide a unique laboratory for investigating the strong interaction and testing theoretical models of hadron structure and dynamics. In this work, the measured correlation functions for oppositely charged $Ξ$-K and $Ξ$ − 𝜋 pairs obtained in high-multiplicity pp collisions at $\sqrt{𝑠}$ = 13 TeV at the LHC are presented. For the first time, high-precision data on the $Ξ$-K interaction are delivered at small relative momenta. The scattering lengths, extracted via the Lednický–Lyuboshits expression of the pair wavefunction, indicate a repulsive and a shallow attractive strong interaction for the $Ξ$-K and $Ξ$ − 𝜋 systems, respectively. The $Ξ$(1620) and $Ξ$(1690) states are observed in the $Ξ$ − 𝜋 correlation function and their properties, mass and width, are determined. These measurements are in agreement with other available results. Such high-precision data can help refine the understanding of these resonant states, provide stronger constraints for chirally motivated potentials, and address the key challenge of describing the coupled-channel dynamics that may give rise to molecular configurations .

Femtoscopy↗
DOE OSTI2022

Impact of fragment formation on shear viscosity in the nuclear liquid-gas phase transition region

Within the improved quantum molecular dynamics (ImQMD) model we follow the evolution of nuclear matter for planar Couette flow in a periodic box. We focus on the region of liquid-gas phase transition and extract the shear viscosity coefficient from the local stress tensor, directly following the viscosity definition. Here, by switching on and off the mean field and thus inducing the phase transition, we are able to observe the impact of clumping in the phase transition region on the viscosity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Full radiator-perturber interaction in computer simulations of hydrogenic spectral line broadening by plasmas

Spectral line broadening by plasmas can be computed by solving the equation of motion for the dipole of the radiating system perturbed by a fluctuating potential obtained from computer simulations. Such calculations have relied on the multipole expansion of the radiator-plasma interaction often keeping only the dipole term. With increasing density, however, higher multipoles as well as plasma perturbers overlapping the bound electron wave functions are expected to become important. For hydrogenic systems, the atomic matrix elements of the full Coulomb and screened Coulomb interactions are given by analytical formulas. Using these results, a computer simulation approach that accounts for the full radiator-plasma interaction is developed. One benefit is the removal of inherent strong collision divergences in the multipole expansion approximation. Furthermore, it yields the plasma polarization shift produced by perturbers penetrating the wave function of the radiator bound electrons. Here, the model was applied to hydrogenic argon Ly-α , Ly- β , and Ly-$\gamma$ spectral lines in a dense argon plasma at free electron densities of 10 24 or 10 25 cm ⁻3 and temperature of $800$ $eV$ relevant to plasma diagnostic techniques for inertial confinement fusion implosions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Searching for solitons in heavy-ion reactions near the fermi energy

Solitons are special shape-conserving hydrodynamical solutions that appear in many areas of physics. Here, we explore the existence of such solutions in microscopic descriptions of the heavy ion reaction 12 C + 28 Si → 12 C*+ 7α in the range E/A = 10 – 65 MeV/u. After recognizing the centrality of the collision and time-reversibility as fundamental requirements for the presence of solitonic 12 C, we utilize the Hybrid α-Cluster model for our analysis with a novel methodology. Our results suggest soliton production for E/A = 25 – 45 MeV/u in the forward direction, with a total cross section being at least in the order of a few μb. This, apart from being encouraging for possible experimental studies, maybe connected to possible toroid states of 28 Si at low angular momenta.

Heavy ion reactions↗
DOE OSTI2021

Ab initio electron dynamics in high electric fields: Accurate prediction of velocity-field curves

Electron dynamics in external electric fields governs the behavior of solid-state electronic devices. First-principles calculations enable precise predictions of charge transport in low electric fields. However, studies of high-field electron dynamics remain elusive due to a lack of accurate and broadly applicable methods. In this work, we develop an efficient approach to solve the real-time Boltzmann transport equation with both the electric field term and ab initio electron-phonon collisions. These simulations provide field-dependent electronic distributions in the time domain, allowing us to investigate both transient and steady-state transport in electric fields ranging from low to high (> 10 kV/cm). The broad capabilities of our approach are shown by computing nonequilibrium electron occupations and velocity-field curves in Si, GaAs, and graphene, obtaining results in quantitative agreement with experiment. Our approach sheds light on microscopic details of transport in high electric fields, including the dominant scattering mechanisms and valley occupation dynamics. Our results demonstrate quantitatively accurate calculations of electron dynamics in low to high electric fields, with broad application to power and micro-electronics, optoelectronics, and sensing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Derivation of a Multi-Species Mix Model

Multi-species flow takes place when two or more gases, plasmas or fluids are flowing and interacting with one another through particle collisions, chemical reactions, turbulence, and thermonuclear processes; this list is not exhaustive. It has a wide range of applications in geological, biological and engineering systems. In what follows species refers to different species of atoms or molecules that may or may not be partially or fully ionized. The interactions between these species takes place at an atomic or molecular level. The species are electrically neutral. These are very important points, because the coupling between species that are interacting on a mesoscale level or are immiscible would be very different than what is being considered in this work. In addition, in section 4, we present a heuristic argument to enhance the collisional mean free path to incorporate the effects of turbulent eddies by use of Prandtl’s assumption that the integral scale length of a turbulent boundary layer is proportional to the width of the turbulent layer. This leads to coupling terms that are non-Newtonian, because they are a function of the dynamics of the flow. It is also assumed that the turbulent stress terms are much smaller than the species partial pressures. This may be a good assumption if the species are hot, dense plasmas that exist in ICF targets.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Shapes of rotating normal fluid He-3 versus superfluid He-4 droplets in molecular beams

Previous single-pulse extreme ultraviolet and X-ray coherent diffraction studies revealed that superfluid 4He droplets obtained in free jet expansion acquire sizable angular momentum, resulting in significant centrifugal distortion. Similar experiments with normal fluid 3He droplets may help elucidating the origin of the of the large degree of rotational excitation and highlight similarities and differences of dynamics in normal and superfluid droplets. Here, we present the first comparison of the shapes of isolated 3He and 4He droplets following expansion of the corresponding fluids in vacuum at temperatures as low as ~ 2 K. Large 3He and 4He droplets with average radii of ~160 nm and ~350 nm, respectively, were produced. We find that the majority of the 3He droplets in the beam correspond to rotating oblate spheroids with reduced average angular momentum ($\Lambda$) and reduced angular velocities ($\Omega$) similar to that of 4He droplets. Given the different physical nature of 3He and 4He, this similarity in $\Lambda$ and $\Omega$ may be surprising and suggest that similar mechanisms induce rotation regardless of the isotope. We hypothesized that the observed distribution of droplet sizes and angular momenta stem from processes in the dense region close to the nozzle. In this region, the significant velocity spread and collisions between the droplets induce excessive rotation followed by droplet fission. The process may repeat itself several times before the droplets enter the collision-fee high vacuum region further downstream.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Nonadiabatically Driven Quantum Interference Effects in the Ultracold K + KRb → Rb + K 2 Chemical Reaction

The K + KRb → Rb + K 2 chemical reaction is the first ultracold atom–diatom chemical reaction for which experimental results have been reported for temperatures below 1 μK more than a decade ago. The reaction occurs through coupling with an excited electronic state that is accessible even in the ultracold limit. A previous quantum dynamics study, excluding nonadiabatic effects, has reported a rate coefficient that is about 35% below the experimental value. Here, we report the first nonadiabatic quantum dynamics study of this reaction and obtain rate coefficients in better agreement with experiments. Our results show that short-range dynamics mediated by coupling with the excited electronic state introduces quantum interference effects that influence both the state-to-state rate coefficients and the overall reaction rates.

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

SDSS-IV MaNGA: global properties of kinematically misaligned galaxies

ABSTRACT We select 456 gas–star kinematically misaligned galaxies from the internal Product Launch-10 of MaNGA survey, including 74 star-forming (SF), 136 green-valley (GV), and 206 quiescent (QS) galaxies. We find that the distributions of difference between gas and star position angles for galaxies have three local peaks at ∼0°, 90°, and 180°. The fraction of misaligned galaxies peaks at log (M*/M⊙) ∼ 10.5 and declines to both low- and high-mass end. This fraction decreases monotonically with increasing star formation rate and specific star formation rate. We compare the global parameters including gas kinematic asymmetry Vasym, H i detection rate and mass fraction of molecular gas, effective radius Re, Sérsic index n as well as spin parameter $\lambda _{R_e}$ between misaligned galaxies and their control samples. We find that the misaligned galaxies have lower H i detection rate and molecular gas mass fraction, smaller size, higher Sérsic index, and lower spin parameters than their control samples. The SF and GV misaligned galaxies are more asymmetric in gas velocity fields than their controls. These pieces of observational evidence point to the gas accretion scenario followed by angular momentum redistribution from gas–gas collision, leading to gas inflow and central star formation for the SF and GV misaligned galaxies. We propose three possible origins of the misaligned QS galaxies: (1) external gas accretion, (2) merger, and (3) GV misaligned galaxies evolve into QS galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Reaction-induced departures from continuum Navier–Stokes turbulence

Reactive hydrodynamic turbulence is an inherently multiscale phenomenon, characterized by the separation between energy-containing, viscous, and molecular length and time scales. The separation between the viscous scale (the Kolmogorov scale) and the molecular mean free path ostensibly justifies a macroscopic description of reactive turbulence via the Navier–Stokes (NS) equations. However, here we use molecular-level simulations to demonstrate that exothermic bimolecular reactions can cause the NS description of turbulence to break down in the near-continuum regime. Sufficiently energetic heat-releasing reactive collisions strongly distort the Maxwell–Boltzmann velocity distribution function, modifying not only the macroscopic chemical rate law but the kinetic-energy-transfer processes as well. This translational nonequilibrium ultimately introduces significant departures from the NS description in the kinetic energy spectra at scales orders of magnitude larger than both the molecular mean free path and the Kolmogorov length scale. These departures prove substantial enough to meaningfully alter integrated quantities, including the overall turbulence kinetic energy itself.

molecular gas dynamics↗
DOE OSTI2024

Multi-level Monte Carlo methods in chemical applications with Lennard-Jones potentials and other landscapes with isolated singularities

We describe and compare outcomes of various Multi-Level Monte Carlo (MLMC) method variants, motivated by the potential of improved computational efficiency over rejection based Monte Carlo, which scales poorly with problem dimension. With an eye toward its application to computational chemical physics, we test MLMC's ability to sample trajectories on two problems — a familiar double-well potential, with known stationary distributions, and a Lennard-Jones solid potential (a Galton Board). By sampling Brownian motion trajectories, we are able to compute expectations of observable averages. These multi-basin potential energy problems capture the essence of the challenges with using MLMC, namely, maintaining correspondence of sample paths as time-resolution is varied. Addressing this challenge properly can lead to MLMC significantly outperforming standard Monte Carlo path sampling. We describe the essence of this problem and suggest strategies that circumvent diverging multilevel sample paths for an important class of problems. In the tests we also compare the computational cost of several, “adaptive,” variants of MLMC. Our results demonstrate that MLMC overcomes the collision, time scale limitation of the more familiar Brownian path MC samplers, and our implementation provides tunable error thresholds, making MLMC a promising candidate for application to larger and more complex molecular systems.

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

Detecting low-energy interactions and the effects of energy accumulation in materials

Elusive as dark matter particles are, they are not the only entities that can produce small energy releases deep inside the most sensitive detectors. Cosmogenic and residual radioactivity as well as other factors can produce and slowly accumulate long-living excitations and defect configurations in materials. Unsteady and avalanche-like releases of accumulated energy can limit the dynamic range of detectors and the sensitivity of experiments. This type of mechanism, while not widely discussed in the dark matter and neutrino detectors community, got a lot of attention in discussions about systems with energy flow in the context of nonequilibrium thermodynamic. Here this paper explores the avalanche hypothesis using published data on low-energy background in detectors and available condensed-matter physics information on excitations, defects, and energetic molecules which can be produced in detector materials and examine ways to use these models to understand the origin of low-energy background in dark matter and other detectors and for suppression of these backgrounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Isotopic Substitution Reveals the Importance of Aluminate Diffusion Dynamics in Gibbsite (Al(OH) 3 ) Crystallization from Alkaline Aqueous Solution

Understanding molecular-scale factors governing the precipitation of aluminum hydroxides, such as gibbsite, under alkaline conditions is important for the formation of laterite deposits, as well as aluminum processing. However, mechanisms enabling tetrahedral aluminate ions to assemble into octahedral sites of the gibbsite lattice remain unclear. Formation of oligomeric complexes has been hypothesized as a critical intermediate step. Here, we report a study of gibbsite solubility in highly alkaline solutions using deuterium substitution to probe equilibrium and kinetic factors that could affect oligomeric intermediate formation, including the reactivity and the diffusivity of aluminate ions. When substituting sodium hydroxide with sodium deuteroxide, solution analysis shows a nearly 40% and 50% increase in gibbsite solubility in 2.4 and 3.3 mol·kg –1 total sodium solutions, respectively. Raman spectroscopy indicated that monomeric aluminate ions are the predominant species in solution irrespective of deuteration. However, both 27 Al and 23 Na nuclear magnetic resonance (NMR) spectroscopy revealed significant differences in chemical shifts in deuterated solutions, and analysis of 1 H, 23 Na, and 27 Al diffusion coefficients with pulsed-field gradient, stimulated echo NMR spectroscopy shows a decrease in ion diffusivity with increasing total deuterium in solution, consistent with the notion of increasing strength in the hydrogen/deuterium bonding network. In conclusion, because the relative change in 1 H diffusion coefficients are commensurate with the difference in apparent solubility constants, there is evidence for an oligomeric intermediate whose steady-state concentration is maintained by the collision frequency of aluminate ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗