Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomic & molecular structure”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Implementing a neural network interatomic model with performance portability for emerging exascale architectures

The two main thrusts of computational science are increasingly accurate predictions and faster calculations; to this end, the zeitgeist in molecular dynamics (MD) simulations is pursuing machine learned and data driven interatomic models, e.g. neural network potentials, and novel hardware architectures, e.g. GPUs. Current implementations of neural network potentials are orders of magnitude slower than traditional interatomic models and while looming exascale computing offers the ability to run large, accurate simulations with these models, achieving portable performance for MD with new and varied exascale hardware requires rethinking traditional algorithms, using novel data structures, and library solutions. We re-implement a neural network interatomic model in CabanaMD, an MD proxy application, built on libraries developed for performance portability. Our implementation shows significantly improved thread scaling in this complex kernel as compared to a current LAMMPS implementation, across both strong and weak scaling. Our single-source solution enables simulations up to 20 million atoms on a single CPU node and 4 million atoms with improved performance on a single GPU. Furthermore, we also explore parallelism and data layout choices (using flexible data structures called AoSoAs) and their effect on performance, seeing up to ~50% and ~5% improvements in performance on a GPU by choosing the right level of parallelism and data layout respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

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

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

Refining T c Prediction in Hydrides via Symbolic‐Regression‐Enhanced Electron‐Localization‐Function‐Based Descriptors

Hydrogen‐based materials are able to possess extremely high superconducting critical temperatures, T c s , due to hydrogen's low atomic mass and strong electron–phonon interaction. Recently, a descriptor based on the Electron Localization Function (ELF) has enabled the rapid estimation of the T c of hydrogen‐containing compounds from electronic networking properties, but its applicability has been limited by the small size and homogeneity of the training dataset used. Herein, the model is re‐examined, compiling a publicly available combined dataset of 244 binary and ternary hydride superconductors. The analysis shows that though ELF‐based networking remains a valuable descriptor, its predictive power declines with increasing compositional complexity. However, by introducing the molecularity index, defined as the highest value of the ELF at which two hydrogen atoms connect, and applying symbolic regression, the accuracy of the predictions can be substantially enhanced. These results establish a more robust framework for assessing superconductivity in hydride materials, facilitating accelerated screening of novel candidates through integration with crystal structure prediction methods or high‐throughput searches.

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

Graphene nanoribbons for quantum electronics

Graphene nanoribbons (GNRs) are a family of one-dimensional (1D) materials with a graphitic lattice structure. GNRs possess high mobility and current-carrying capability, sizeable bandgap and versatile electronic properties, which make them promising candidates for quantum electronic applications. In the past 5 years, progress has been made towards atomically precise bottom-up synthesis of GNRs and heterojunctions that provide an ideal platform for functional molecular devices, as well as successful production of semiconducting GNR arrays on insulating substrates potentially useful for large-scale digital circuits. With further development, GNRs can be envisioned as a competitive candidate material in future quantum information sciences. In this Perspective, we discuss recent progress in GNR research and identify key challenges and new directions likely to develop in the near future.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Preferential interactions of primary amine-terminated quantum dots with membrane domain boundaries and lipid rafts revealed with nanometer resolution

The initial interactions of engineered nanoparticles (NPs) with living cells are governed by physicochemical properties of the NP and the molecular composition and structure of the cell membrane. Eukaryotic cell membranes contain lipid rafts – liquid-ordered nanodomains involved in membrane trafficking and molecular signaling. However, the impact of these membrane structures on cellular interactions of NPs remains unclear. Here we investigate the role of membrane domains in the interactions of primary amine-terminated quantum dots (Qdots) with liquid-ordered domains or lipid rafts in model membranes and intact cells, respectively. Using correlative atomic force and fluorescence microscopy, we found that the Qdots preferentially localized to boundaries between liquid-ordered and liquid-disordered phases in supported bilayers. The Qdots also induced holes at these phase boundaries. Using super resolution fluorescence microscopy (STORM), we found that the Qdots preferentially co-localized with lipid rafts in the membrane of intact trout gill epithelial cells – a model cell type for environmental exposures. Our observations uncovered preferential interactions of amine-terminated Qdots with liquid-ordered domains and their boundaries, possibly due to membrane curvature at phase boundaries creating energetically favorable sites for NP interactions. The preferential interaction of the Qdots with lipid rafts supports their potential internalization via lipid raft-mediated endocytosis and interactions with raft-resident signaling molecules.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Many-body reactive force field development for carbon condensation in C/O systems under extreme conditions

In this paper, we describe the development of a reactive force field for C/O systems under extreme temperatures and pressures, based on the many-body Chebyshev Interaction Model for Efficient Simulation (ChIMES). The resulting model, which targets carbon condensation under thermodynamic conditions of 6500 K and 2.5 g cm –3 , affords a balance between model accuracy, complexity, and training set generation expense. We show that the model recovers much of the accuracy of density functional theory for the prediction of structure, dynamics, and chemistry when applied to dissociative condensed phase systems at 1:1 and 1:2 C:O ratios, as well as molten carbon. Our C/O modeling approach exhibits a 104 increase in efficiency for the same system size (i.e., 128 atoms) and a linear system size scalability over standard quantum molecular dynamics methods, allowing the simulation of significantly larger systems than previously possible. We find that the model captures the condensed-phase reaction-coupled formation of carbon clusters implied by recent experiments, and that this process is susceptible to strong finite size effects. Overall, we find the present ChIMES model to be well suited for studying chemical processes and cluster formation at pressures and temperatures typical of shock waves. We expect that the present C/O modeling paradigm can serve as a template for the development of a broader high pressure–high temperature force-field for condensed phase chemistry in organic materials.

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

Epitaxial PbGeSe thin films and their photoluminescence in the mid-wave infrared

PbSe is a narrow bandgap IV–VI compound semiconductor with application in mid-wave infrared optoelectronics, thermoelectrics, and quantum devices. Alkaline-earth or rare-earth elements such as Sr and Eu can substitute Pb to widen the bandgap of PbSe in heterostructure devices, but they come with challenges such as deteriorating optical and electronic properties, even in dilute concentrations due to their dissimilar atomic nature. We substitute Pb instead with column IV Ge and assess the potential of rocksalt phase PbGeSe as a wider bandgap semiconductor in thin films grown by molecular beam epitaxy on GaAs substrates. Low sticking of GeSe adatoms requires synthesis temperatures below 260 °C to incorporate Ge, but this yields poor structural and compositional uniformity as determined by x-ray diffraction. Consequently, as-grown films in the range Pb0.94Ge0.06Se–Pb0.83Ge0.17Se (6%–17% Ge) show much less bandgap widening in photoluminescence than prior work on bulk crystals using absorption. We observe that post-growth rapid thermal annealing at temperatures of 375–450 °C improves the crystal quality and recovers bandgap widening. Rapid interdiffusion of Ge during annealing, however, remains a challenge in harnessing such PbGeSe materials for compositionally sharp heterostructures. Annealed 17% Ge films emit light at 3–3.1 μm with a minimal shift in wavelength vs temperature. These samples are wider in bandgap than PbSe films by 55 meV at room temperature, and the widening increases to 160 meV at 80 K, thanks to sharply different dependence of bandgap on temperature in PbSe vs PbGeSe.

Xiao, Kelly (ORCID:0000000321000838)↗
DOE OSTI2020

Dimensionality-Mediated Semimetal-Semiconductor Transition in Ultrathin PtTe 2 Films

Platinum ditelluride (PtTe2), a type-II Dirac semimetal, remains semimetallic in ultrathin films down to just two tri-atomic layers (TLs) with a negative gap of -0.36 eV. Further reduction of the film thickness to a single TL induces a Lifshitz electronic transition to a semiconductor with a large positive gap of +0.79 eV. This transition is evidenced by experimental band structure mapping of films prepared by layer-resolved molecular beam epitaxy, and by comparing the data to first-principles calculations using a hybrid functional. The results demonstrate a novel electronic transition at the two-dimensional limit through film thickness control.

36 MATERIALS SCIENCE↗
NASA NTRS1998

X-Ray Diffraction Analysis of NLO Crystals: Traditional Applications and More New Opportunities

Single crystal X-ray diffraction analysis is one of the more important methods for the molecular and crystal structure determination of matter and therefore it has a great importance in material science including design and engineering of different compounds with non-linear optical (NLO) properties. It was shown in our previous publications that this method provides unique information about molecular structure of NLO compounds, their crystal symmetry and crystal packing arrays, molecular conformation and geometries and many other structural and electronic characteristics that are important for understanding the nature of NLO properties of solids. A very new application of the X-ray diffraction method is related to analysis of the electron density distribution p(r) in crystals and some of its characteristics (atomic and group charges, dipole and higher multipole moments, etc.), that may be obtained directly form the diffraction measurements. In the present work, we will discuss our preliminary low temperature high-resolution X-ray data for the m-nitroaniline (mNA) single crystal (VI). This is one of the "classical" organic NLO materials and electron density distribution analysis in this simple compound has a great scientific interest.

Antipin, Mikhail Yu.↗
DOE OSTI2021

AFM Slip Length Measurements for Water at Selected Phyllosilicate Surfaces

Most reported slip length measurements have been made at the surfaces of synthetic materials and modified synthetic materials. In contrast, few slip length measurements at the surface of unmodified natural mineral surfaces have been reported. In this regard, flow at the silica face surfaces of the phyllosilicate minerals, talc and mica, was considered. A slip boundary condition was expected at the nonpolar hydrophobic silica surface of talc leading to enhanced flow, and a no-slip boundary condition was expected at the hydrophilic silica surface of mica. Atomic force microscopy (AFM) slip length measurements were made at the talc and mica surfaces. The slip length results for the hydrophobic silica surface of talc were contrasted to the results for the hydrophilic silica surface of mica (no-slip flow). The results are discussed based on molecular dynamics simulations (MDS), as reported in the literature, and AFM images of surface nanobubbles. For nonpolar hydrophobic surfaces (such as talc), it is doubtful that the MDS interfacial water structure and the water exclusion zone (3.2 Å) account for the AFM slip flow with slip lengths as great as 95 nm. Rather, a better explanation for the AFM slip flow condition is based on reduced interfacial viscosity due to the presence of dissolved gas and the accommodation of pancake nanobubbles at the talc surface having a height dimension of magnitude similar to the slip length.

42 ENGINEERING↗
DOE OSTI2024

Evidence of short chains in liquid sulfur

High energy x-ray pair distribution function measurements show the average coordination number of the first shell in liquid sulfur is 1.86 ± 0.04 across the λ-transition, not precisely 2.0 as widely accepted. This indicates that upon melting, liquid sulfur does not comprise solely of S 8 rings but also possesses a significant number of short chains. Intensities of the pre-peak and first diffraction peak of the x-ray structure factor and third peak height of the pair distribution function all show deviations at the λ-transition temperature T λ , associated with the break-up of S 8 rings and the start of oligomer polymerization. A significant number of non-bonded or loosely bonded “interstitial atoms,” with an average coordination number of 0.20 ± 0.005, are also observed in the so-called “forbidden zone” between the first and second shells upon melting. The number of interstitial atoms is found to decrease to a minimum at the λ-transition, but the majority persist into the high temperature polymerized liquid. Furthermore, the existence of short chains and nearby interstitial atoms represent the two main factors required to initiate the S 8 -ring to chain transition, as proposed by recent molecular dynamics simulations.

Chemical bonding↗
DOE OSTI2023

Exploring NaCl-PuCl 3 molten salts with machine learning interatomic potentials and graph theory

Actinide molten salts are the basis of the liquid fuels used in molten salt reactors. Due to the inherent difficulties associated with high temperature and hazardous conditions, experimental investigations of fundamental properties of these materials are usually challenging. In this work, we describe the structure and transport of NaCl-PuCl 3 mixtures using computational techniques. Three compositions were considered (16, 25, and 36 mol% PuCl 3 ) over a temperature range (730 – 1257K) using ab initio molecular dynamics, which provided the necessary data sets for training machine learned interatomic potentials. Further, molecular dynamics simulations based on these potentials were then used to determine structure and transport properties. A substantial change was noted in the structure factor when increasing the PuCl 3 content from 25 to 36 mol%. This change is linked to the aggregation of larger Pu 3+ clusters. In addition, the similarity of the atomic environments of metal cations in molten salt systems to their solid states counterparts was investigated using an unsupervised learning technique. Finally, graph theory was employed to explore the structure and size of actinide networks. Consistent with the structure factor, a dense Pu 3+ intermolecular structure is observed within the 36 mol% PuCl 3 mixture. The structure of cation-cation inter-junctions is also discussed. In all cases, the diffusion of Pu 3+ is significantly lower than that of Na + and Cl - .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Trigonal Bipyramidal V 3+ Complex as an Optically Addressable Molecular Qubit Candidate

Synthetic chemistry enables a bottom-up approach to quantum information science, where atoms can be deterministically positioned in a quantum bit or qubit. Two key requirements to realize quantum technologies are qubit initialization and read-out. By imbuing molecular spins with optical initialization and readout mechanisms, analogous to solid-state defects, molecules could be integrated into existing quantum infrastructure. To mimic the electronic structure of optically addressable defect sites, we designed the spin-triplet, V 3+ complex, (C 6 F 5 ) 3 trenVCN t Bu (1).We measured the static spin properties as well as the spin coherence time of1demonstrating coherent control of this spin qubit with a 240 GHz electron paramagnetic resonance spectrometer powered by a free electron laser. We found that1exhibited narrow, near-infrared photoluminescence (PL) from a spin-singlet excited state. Using variable magnetic field PL spectroscopy, we resolved emission into each of the ground-state spin sublevels, a crucial component for spin-selective optical initialization and readout. This work demonstrates that trigonally symmetric, heteroleptic V 3+ complexes are candidates for optical spin addressability.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

UV-Induced Reaction Pathways in Bromoform Probed with Ultrafast Electron Diffraction

For many chemical reactions, it remains notoriously difficult to predict and experimentally determine the rates and branching ratios between different reaction channels. This is particularly the case for reactions involving short-lived intermediates, whose observation requires ultrafast methods. The UV photochemistry of bromoform (CHBr 3 ) is among the most intensely studied photoreactions. Yet, a detailed understanding of the chemical pathways leading to the production of atomic Br and molecular Br 2 fragments has proven challenging. In particular, the role of isomerization and/or roaming and their competition with direct C–Br bond scission has been a matter of continued debate. Here, in this work, gas-phase ultrafast megaelectronvolt electron diffraction (MeV-UED) is used to directly study structural dynamics in bromoform after single 267 nm photon excitation with femtosecond temporal resolution. The results show unambiguously that isomerization contributes significantly to the early stages of the UV photochemistry of bromoform. In addition to direct C–Br bond breaking within <200 fs, formation of iso-CHBr 3 (Br-CH-Br-Br) is observed on the same time scale and with an isomer lifetime of >1.1 ps. The branching ratio between direct dissociation and isomerization is determined to be 0.4 ± 0.2:0.6 ± 0.2, i.e., approximately 60% of molecules undergo isomerization within the first few hundred femtoseconds after UV excitation. The structure and time of formation of iso-CHBr 3 compare favorably with the results of an ab initio molecular dynamics simulation. The lifetime and interatomic distances of the isomer are consistent with the involvement of a roaming reaction mechanism.

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

Real-space atomic dynamics in metallic liquids investigated by inelastic neutron scattering

Understanding the dynamics of liquids at the atomic level remains a major challenge. Even though viscosity is one of the most fundamental properties of liquids, its atomistic origin is not fully elucidated. Through inelastic neutron scattering experiment on levitated metallic liquid droplets, the time-dependent pair correlation function, the Van Hove function, was determined for Z⁢r 50 ⁢C⁢u 50 and Z⁢r 80 ⁢P⁢t 20 liquids at various temperatures. The time for change in local atomic connectivity, τ LC , which is the timescale of atomic bond cutting and forming, is estimated based on the exponential decay of the nearest neighbor peak of the Van Hove function. At high temperatures above the crossover temperature T A , τ LC is equal to the Maxwell relaxation time, τ M =η/G ∞ , where η is the macroscopic shear viscosity and G ∞ is the high-frequency shear modulus. Below T A the ratio of τ M /τ LC increases with decreasing temperature, indicating increased atomic cooperativity as predicted by molecular dynamics simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Normal mode description of phases of matter: Application to heat capacity

Understanding thermodynamics in liquids at the atomic level is challenging because of strong atomic interactions and lack of spatial symmetry. Recent prior theoretical works have focused on describing heat capacity of liquids in terms of phonon-like excitations but often rely on fitting factors and assumptions. In this work, we propose characterizing various phases in terms of instantaneous normal modes (INMs) of structural snapshots from molecular dynamics simulations of single-element systems over wide ranges of temperature and pressure. We use the INMs to build a mode-level microscopic description of heat capacity and demonstrate that heat capacity of liquids can be described by a combination of both solidlike and gaslike degrees of freedom, leading to a more unified framework to fundamentally describe heat capacity of all three phases of matter: solid, liquid, and gas. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Hydration-driven stabilization and volume collapse of grain boundaries in Mg2SiO4 forsterite predicted by first-principles simulations

Grain boundaries in mantle minerals are of critical importance to geophysical and geochemical processes of the Earth’s interior. One of the fundamental issues is to understand how the water (H 2 O) component influences the properties of grain boundaries in silicate materials. Here, we report the results of the structure and stability of several tilt grain boundaries in Mg 2 SiO 4 forsterite over the pressure range 0 to 15 GPa using density functional theory-based first-principles simulations. Here, the results suggest greater energetic stability and hydration-driven volume collapse (negative excess volume) at zero pressure for the majority of hydrous grain boundaries relative to the anhydrous (dry) ones. All the hydrous grain boundaries become increasingly favorable at elevated pressures as the calculated hydration enthalpy systematically decreases with increasing pressure. The hydrous components at the interfacial regions are predominantly in the hydroxyl form and, to a lesser extent, in the molecular H 2 O form. Their calculated ratio ranges from 1.6 to 8.7 among the different grain boundary configurations. Our structural analysis also reveals that the hydroxyls are bound to either both Mg and Si or to Mg only. In comparison, the molecular species are bound only to Mg sites. Besides direct oxygen-hydrogen bonding, intermolecular hydrogen bonding becomes important with compression. On the basis of our results, we suggest that local atomic rearrangements caused by dissociative adsorption of water facilitate efficient compaction of the boundary interfaces, which, in turn, results in greater relative stability of hydrous grain boundaries. This means that water prefers to be incorporated within the grain boundaries over the bulk of silicate materials.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗