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 109 records · Page 6

Recent progress on in-medium properties of heavy mesons from finite-temperature EFTs

Mesons with heavy flavor content are an exceptional probe of the hot QCD medium produced in heavy-ion collisions. In the past few years, significant progress has been made toward describing the modification of the properties of heavy mesons in the hadronic phase at finite temperature. Ground-state and excited-state thermal spectral properties can be computed within a self-consistent many-body approach that employs appropriate hadron-hadron effective interactions, providing a unique opportunity to confront hadronic Effective Field Theory predictions with recent and forthcoming lattice QCD simulations and experimental data. In this article, we revisit the application of the imaginary-time formalism to extend the calculation of unitarized scattering amplitudes from the vacuum to finite temperature. These methods allow us to obtain the ground-state thermal spectral functions. The thermal properties of the excited states that are dynamically generated within the molecular picture are also directly accessible. We present here the results of this approach for the open-charm and open-bottom sectors. We also analyze how the heavy-flavor transport properties, which are strongly correlated to experimental observables in heavy-ion collisions, are modified in hot matter. In particular, transport coefficients can be computed using an off-shell kinetic theory that is fully consistent with the effective theory describing the scattering processes. The results of this procedure for both charm and bottom transport coefficients are briefly discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Radiation induced athermal diffusivity in uranium mononitride

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuel considered for light water reactors and advanced reactor designs. Properties like self- and fission gas diffusivity need to be better understood, given that they influence key fuel performance phenomena such as fission gas swelling and release. In particular, the radiation induced athermal (D 3 ) diffusivity remains challenging to accurately predict and has only been sparsely characterized in UN, despite its importance as it likely governs diffusion at the low temperatures this high-thermal-conductivity fuel form may operate. Molecular Dynamics simulations are used to estimate the mean square displacement induced by a primary knock-on atom (PKA) with a given kinetic energy. These results are combined with the PKA energy distributions obtained from binary collision approximation calculations to obtain the displacement due to a particular fission fragment. Finally, this is combined with experimental fission fragment yields to determine the displacement due to an average fission event and, thus, express the athermal diffusivity as a function of the fission rate density. These results are in excellent agreement with available experimental data. In conclusion, a particular importance is given to the understanding and the quantification of the variability of these results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Viscosity of the magnetized strongly coupled one-component plasma

The viscosity tensor of the magnetized one-component plasma, consisting of five independent shear viscosity coefficients, a bulk viscosity coefficient, and a cross coefficient, is computed using equilibrium molecular dynamics simulations and the Green-Kubo relations. A broad range of Coulomb coupling and magnetization strength conditions are studied. Magnetization is found to strongly influence the shear viscosity coefficients when the gyrofrequency exceeds the Coulomb collision frequency. Three regimes are identified as the Coulomb coupling strength and magnetization strength are varied. Here, the Green-Kubo relations are used to separate kinetic and potential energy contributions to each viscosity coefficient, showing how each contribution depends upon the magnetization strength. The shear viscosity coefficient associated with the component of the pressure tensor parallel to the magnetic field, and the two coefficients associated with the component perpendicular to the magnetic field, are all found to merge to a common value at strong Coulomb coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pair momentum dependence of a tilted source in heavy-ion collisions

In noncentral heavy-ion collisions, the particle-emitting source can be tilted away from the beam direction, an effect that becomes particularly significant at collision energies of a few GeV and lower. This phenomenon, manifest itself in many observables such as directed flow, polarization, and vorticity, is therefore important to investigate. Here, in this paper, we study the consistency between the tilt extracted directly from the freeze-out distribution of pions and the tilt parameter obtained using the azimuthally sensitive femtoscopy (asHBT) method. Using the Ultrarelativistic Quantum Molecular Dynamics model, we demonstrate a strong dependence of the tilt parameter extracted with asHBT on the momentum of the particle pair. Considering the experimental challenges in accessing low particle momenta—where the tilt parameter extracted with asHBT closely matches the tilt of the freeze-out distribution of pions—we propose an exponential extrapolation method to obtain the tilt of the entire freeze-out distribution. This approach aims to enhance the accuracy of experimental measurements of tilt in noncentral heavy-ion collisions.

particle correlations & fluctuations↗

X-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula

Here, we study ab initio approaches for calculating x-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula that expresses the inelastic contribution in terms of the dielectric function. We study the electronic dynamic structure factor computed from the Mermin dielectric function using an ab initio electron-ion collision frequency in comparison to computations using a linear-response time-dependent density functional theory (LR-TDDFT) framework for hydrogen and beryllium and investigate the dispersion of free-free and bound-free contributions to the scattering signal. A separate treatment of these contributions, where only the free-free part follows the Mermin dispersion, shows good agreement with LR-TDDFT results for ambient-density beryllium, but breaks down for highly compressed matter where the bound states become pressure ionized. LR-TDDFT is used to reanalyze x-ray Thomson scattering experiments on beryllium demonstrating strong deviations from the plasma conditions inferred with traditional analytic models at small scattering angles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Photoinitiated Reactions of Molecules and Radicals in Molecular Beams

The UV photochemistry of organic molecules is a fundamental process that governs reactions in the atmosphere, synthetic chemistry, processing of organic aerosols, and biological damage in living tissues. In many organic molecules, the ensuing evolution involves pathways that are in competition, giving rise to different products, isomerization, coupling to other electronic states, and secondary reactive collisions. Because photodissociation is usually fast (picoseconds to microseconds), these processes are far from equilibrium and are controlled by kinetic competition and dynamical forces. The work accomplished was focused primarily on the photochemistry of alpha-keto carboxylic acids, a group of acids produced from natural sources, which are implicated in aerosol formation and biological processes. In the atmosphere, they are destroyed mainly by solar radiation. Studying their photochemistry has been surprisingly difficult because of the complexity of their excited electronic states, and the effect of collisions and secondary reactions. The second project investigated the lifetime of excited electronic states of pyrazine and picoline and had both experimental and theoretical aspects. The work focused on a model molecule, pyruvic acid (PA), because it was hypothesized that it should be a good source of the unstable carbene, methylhydroxycarbene (MHC). PA has an internal hydrogen bond that controls the evolution of its decomposition. The decomposition was studied following excitation to two of its lowest excited states reached by laser irradiation at 351 and 193 nm. The photodissociation dynamics in the absence and presence of collisions was monitored and compared. Understanding the UV photochemistry requires the use of complementary experimental approaches. Two methods were enlisted, which together generated a comprehensive and detailed set of results: (i) The time-sliced velocity map imaging (SVMI) instrument at USC was exploited to determine kinetic energy release of fragments, fast dissociation timescales, and internal state distributions of fragments for which Resonance Enhanced Multiphoton Ionization (REMPI) schemes exist; and (ii) The multiplexed photoionization mass spectrometer (MPIMS) setup developed at the Sandia Combustion Research Facility was used for product discovery, achieved by exploiting tunable narrowband VUV radiation at the Advanced Light Source (ALS), and to follow in real time their subsequent unimolecular and bimolecular reactions. The experimental conditions ranged from collisionless molecular beams to study nascent products to flow reactors of variable pressures to study subsequent bimolecular reactions of the products. The goals stated above were successfully accomplished. By exciting PA to its lowest excited state, MHC was identified as the only primary product and its isomerization to vinylalcohol and acetaldehyde was directly observed in real time. Moreover, we reported the first bimolecular reaction of MHC, which was with acetaldehyde, and identified its reaction product. This paves the way for the study of other reactions of this important carbene intermediate. When excitation was carried out at 193 nm, which imparted much higher energy to PA, many more products were directly observed in real time, some deriving from three-body dissociation. Quantitative branching ratios and secondary reactions of radical products were also determined and analyzed. The studies of pyrazine and picoline focused on the second excited state of these molecules and showed (via ionization studies) that their excited states were very short lived (<100 fs) because of efficient couplings to lower electronic states via vibronic coupling. Theoretical work on modeling the absorption spectrum and decay mechanisms of the excited states are in progress in collaboration with theoreticians.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Molecular dynamics simulations of radiation cascade evolution near cellular dislocation structures in additively manufactured stainless steels

The dislocation dense cellular structures observed in additively manufactured (AM) stainless steels can allow these materials to have enhanced strength and ductility compared to conventionally manufactured materials. Effective design of radiation resistant materials often necessitates creating microstructures containing of a high density of dislocations that act as sinks for radiation induced defects. This work uses molecular dynamics simulations to study the impact of pre-existing defects on radiation damage in stainless steel 316L fabricated by the laser powder bed fusion process. The evolution of synthetically generated dislocations in a crystal system in response to multiple radiation collision cascades suggest that dislocation dense regions reduce the probability of surviving point defects forming clusters. This finding seems to result from the dislocation cores’ ability to absorb point defects, thereby limiting the formation of clusters during the recombination phase. However, after successive cascades, the ability of the dislocation entanglement to prevent defect formation is diminished and defect cluster formation trends to an equilibrium. Accordingly, it is hypothesized that in an experimental setting the dislocation cells will initially act as neutral sinks for point defects, potentially delaying the onset of radiation damage effects; but, these structures will be degraded due to radiation-enhanced diffusion and lose their effectiveness at doses above 1.5–2 displacements per atom. Nevertheless, the ability to produce complex and feature-specific microstructures suggests promise for the design of radiation tolerant materials using AM methods.

36 MATERIALS SCIENCE↗

Calculation of self-diffusion coefficients in supercritical carbon dioxide using mean force kinetic theory

This paper presents an application of mean force kinetic theory (MFT) to the calculation of the self-diffusivity of CO 2 in the supercritical fluid regime. Two modifications to the typical application of MFT are employed to allow its application to a system of molecular species. Therefore, the first is the assumption that the inter-particle potential of mean force can be obtained from the molecule center-of-mass pair correlation function, which in the case of CO 2 is the C–C pair correlation function. The second is a new definition of the Enskog factor that describes the effect of correlations at the surface of the collision volume. The new definition retains the physical picture that this quantity represents a local density increase, resulting from particle correlations, relative to that in the zero density homogeneous fluid limit. These calculations are facilitated by the calculation of pair correlation functions from molecular dynamics (MD) simulations using the FEPM2 molecular CO2 model. The self-diffusivity calculated from theory is in good agreement with that from MD simulations up to and slightly beyond the density at the location of the Frenkel line. The calculation is compared with and is found to perform similarly well to other commonly used models but has a greater potential for application to systems of mixed species and to systems of particles with long range interatomic potentials due to electrostatic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The CO Emission in the Taffy Galaxies (UGC 12914/15) at 60 pc Resolution. I. The Battle for Star Formation in the Turbulent Taffy Bridge

We present Atacama Large Millimeter/submillimeter Array observations at a spatial resolution of $0^{"}_{.}2$ (60 pc) of CO emission from the Taffy galaxies (UGC 12914/5). The observations are compared with narrowband Paα, mid-IR, radio continuum and X-ray imaging, plus optical spectroscopy. The galaxies have undergone a recent head-on collision, creating a massive gaseous bridge that is known to be highly turbulent. The bridge contains a complex web of narrow molecular filaments and clumps. The majority of the filaments are devoid of star formation, and fall significantly below the Kennicutt–Schmidt relationship for normal galaxies, especially for the numerous regions undetected in Paα emission. Within the loosely connected filaments and clumps of gas we find regions of high velocity dispersion that appear gravitationally unbound for a wide range of likely values of X CO . Like the "Firecracker" region in the Antennae system, they would require extremely high external dynamical or thermal pressure to stop them dissipating rapidly on short crossing timescales of 2–5 Myr. We suggest that the clouds may be transient structures within a highly turbulent multiphase medium that is strongly suppressing star formation. Despite the overall turbulence in the system, stars seem to have formed in compact hotspots within a kiloparsec-sized extragalactic H ii region, where the molecular gas has a lower velocity dispersion than elsewhere, and shows evidence for a collision with an ionized gas cloud. Like the shocked gas in the Stephan's Quintet group, the conditions in the Taffy bridge shows how difficult it is to form stars within a turbulent, multiphase, gas.

79 ASTRONOMY AND ASTROPHYSICS↗

Alloying effects on low-energy recoil events in concentrated solid-solution alloys

Alloying elements at low-concentrations into pure metals is usually adopted to improve their mechanical properties and irradiation performance. Recently developed single-phase concentrated solid-solution alloys (CSAs) comprised of two or more elements, all at high concentrations, have demonstrated good radiation resistance. CSAs are characterized by their extreme disordered states that arise from the random arrangement of different elements and accompanied random local lattice distortions. In this work, we investigate low-energy recoil events in Ni 0.5 X 0.5 (X = Fe and Co) CSAs and Ni 0.8 X 0.2 (X = Fe, Co, Cr, and Pd) CSAs using ab initio molecular dynamics simulations to understand the effects of different disorder on defect production in CSAs. The threshold displacement energies are determined along three high-symmetric directions by randomly choosing independent primary knock-on atoms in each direction. As expected, the threshold energies in Ni and its CSAs are anisotropic, with the highest values found in the [111] direction. The calculated threshold energies of Fe in NiFe are smaller than those in NiCo and Ni, especially along the [111] direction. An inspection of the atomic trajectories inside the collision cascade reveals that the effect of chemical disorder outweighs the site-to-site lattice distortions in determining the threshold energies. Especially, different interaction properties between elements due to their different electronic structures are responsible for the observed different threshold energies. Finally, the local environment dependence of threshold energies suggests that local elemental arrangement can be used to understand and predict threshold energies in disordered alloys.

36 MATERIALS SCIENCE↗

Molecular dynamics simulations of reflection and sputtering behavior of boron under deuterium ion irradiation

Boronization is a commonly used method of wall conditioning in fusion reactors. The application of boron films to the plasma-facing materials results in enhanced plasma performance due to the reduction of intrinsic impurities. This is primarily driven by a reduction in oxygen content that is chemically trapped in the boron film. The reactive nature of these boron films also raises questions concerning interactions with hydrogen isotopes. In this work, boron-deuterium interactions were studied using molecular dynamics (MD). Reactive force field potentials were used to model the chemical interactions between B and D. An amorphous boron substrate was irradiated by D atoms at varying incident energies, 10 eV < E i < 150 eV and angles, 0° < α < 85°. The reflection probability was calculated and compared to results from the commonly-used binary collision approximation (BCA) method. This comparison found that the BCA underestimated the reflection probability at E i < 35 eV and α > 45 ∘ . The source of this discrepancy was found to be the surface binding energy model. The BCA calculation with an isotropic surface binding energy model was more closely aligned to the MD result. This, in combination with a correction function based on the MD results allows for corrections to the reflection probability of deuterium impinging on boron surfaces. The sputtering of the substrate material was also studied. While this study did not contain sufficient events to quantitatively describe the sputtering behavior, some qualitative results emerged: namely, chemical sputtering of B and D-containing molecules (BD, BD 2 , BD 3 ) at low ( < 20 eV) incident deuterium energies. This result suggests that chemical sputtering could be a significant factor in limiting boron coating lifetime when exposed to lower ion energies, such as those in detached plasmas. The results show that chemical interactions should be taken into account when modeling the interactions between D ions and B surfaces.

Boron↗

Thermodynamics, dynamics, and structure of supercritical water at extreme conditions

Molecular dynamics (MD) simulations to understand the thermodynamic, dynamic, and structural changes in supercritical water across the Frenkel line and the melting line have been performed. The two-phase thermodynamic model [J. Phys. Chem. B(2010), 114(24), 8191-8198] and the velocity autocorrelation functions are used to locate the Frenkel line and to calculate the thermodynamic and dynamic properties. The Frenkel lines obtained from the two-phase thermodynamic model and the velocity autocorrelation criterion do not agree with each other. Structural characteristics and the translational diffusion dynamics of water suggest that this inconsistency could arise from the two oscillatory modes in water, which are associated with the bending of hydrogen bonds and intermolecular collisions inside the first coordination shell. Finally, the overall results lead us to conclude that the universality of the Frenkel line as a dynamic crossover line from rigid to nonrigid fluids is preserved in water.

74 ATOMIC AND MOLECULAR PHYSICS↗

Observation of Universal Quench Dynamics and Townes Soliton Formation from Modulational Instability in Two-Dimensional Bose Gases

We experimentally study universal nonequilibrium dynamics of two-dimensional atomic Bose gases quenched from repulsive to attractive interactions. Here, we observe the manifestation of modulational instability that, instead of causing collapse, fragments a large two-dimensional superfluid into multiple wave packets universally around a threshold atom number necessary for the formation of Townes solitons. We confirm that the density distributions of quench-induced solitary waves are in excellent agreement with the stationary Townes profiles. Furthermore, our density measurements in the space and time domain reveal detailed information about this dynamical process, from the hyperbolic growth of density waves, the formation of solitons, to the subsequent collision and collapse dynamics, demonstrating multiple universal behaviors in an attractive many-body system in association with the formation of a quasistationary state.

74 ATOMIC AND MOLECULAR PHYSICS↗

Dynamics of resonant low-energy electron attachment to ethanol-producing hydroxide anions

Here, the dynamics of dissociative electron attachment to ethanol is experimentally investigated at the Feshbach resonance formed with incident electron energies near 9.5 eV. Highly differential laboratory-frame momentum distributions of OH – fragments are measured for a series of energies spanning the resonance width, using the velocity-map-imaging technique. The OH – kinetic-energy distribution indicates that the C-O breaking dissociation process could either be a three-body dissociation or a two-body dissociation with significant rovibrational excited fragments. The small, but significant, anisotropy in the OH – angular distribution provides signatures of the molecular symmetry of the associated resonant state under the axial recoil approximation, which assumes the dissociation is much faster than any rotation of the dissociation axis. Within these assumptions, the 9.5-eV Feshbach resonance can be assigned to the electronic transition from the ($10a'$) orbital with its ground-state C s symmetry to the empty ($4a''$) level, involving the simultaneous electron attachment. This dynamics could be a model for C-O dissociation in larger alcohols and ethers.

74 ATOMIC AND MOLECULAR PHYSICS↗

Exploring the impacts of conformer selection methods on ion mobility collision cross section predictions

The prediction of structure dependent molecular properties, such as collision cross sections as measured using ion mobility spectrometry, are crucially dependent on the selection of the correct population of molecular conformers. Here, we report an in-depth evaluation of multiple conformation selection techniques, including simple averaging, Boltzmann weighting, lowest energy selection, low energy threshold reductions, and similarity reduction. Generating 50,000 conformers each for 18 molecules, we used the In Silico Chemical Library Engine (ISiCLE) to calculate the collision cross sections for the entire dataset. First, we employed Monte Carlo simulations to understand the variability between conformer structures as generated using simulated annealing. Then we employed Monte Carlo simulations to the aforementioned conformer selection techniques applied on the simulated molecular property—the ion mobility collision cross section. Based on our analyses, we found Boltzmann weighting to be a good tradeoff between precision and theoretical accuracy. Combining multiple techniques revealed that energy thresholds and root-mean-squared deviation-based similarity reductions can save considerable computational expense while maintaining property prediction accuracy. Molecular dynamic conformer generation tools like AMBER can continue to generate new lowest energy conformers even after tens of thousands of generations, decreasing precision between runs. This reduced precision can be ameliorated and theoretical accuracy increased by running density functional theory geometry optimization on carefully selected conformers.

Nielson, Felicity F.↗

Pressure Effects on the Relaxation of an Excited Ethane Molecule in High-Pressure Bath Gases

Here, we use molecular dynamics to calculate the rotational and vibrational energy relaxation of C 2 H 6 in Ar, Kr, and Xe bath gases over a pressure range of 10 to 400 atm and at temperatures of 300 K and 800 K. The C 2 H 6 is instantaneously excited by 80 kcal/mol randomly distributed into both vibrational and rotational modes. The computed relaxation rates show little sensitivity to the identity of the noble gas in the bath. Vibrational relaxation rates show a non-linear pressure dependence at 300 K. At 800 K the reduced range of bath gas densities covered by the range of pressures do not yet show any non-linearity in the pressure dependence. Rotational relaxation is characterized with two relaxation rates. The slower rate is comparable to the vibrational relaxation rate. The faster rate has a linear pressure dependence at 300 K but an irregular, nonlinear pressure dependence at 800 K. To understand this, a model was developed based on approximating the periodic box used in the molecular dynamics simulations by an equal-volume collection of cubes where each cube is sized to allow only single occupancy by the noble gas or the molecule. Combinatorial statistics then leads to a pressure and temperature dependent analytic distribution of the bath gas species the molecule encounters in a collision. This distribution, the dissociation energy of molecule/bath gas complexes and bath gas clusters, and the computed energy release per collision combine to show that only at 300 K is the energy release sufficient to dissociate likely complexes and clusters. This suggests that persistent and pressure-dependent clusters and complexes at 800 K may be responsible for the non-linear pressure dependence of rotational relaxation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗