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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Measured space environmental effects to LDEF during retrieval

The Interim Operational Contamination Monitor (IOCM) is an attached shuttle payload to the Long Duration Exposure Facility (LDEF) that was used on 2 earlier flights to quantify the contamination deposited during the course of the missions. The IOCM can characterize by direct measurement the deposition of molecular and particulate contamination during any phase of flight. Also, the IOCM actively measures the optical property changes of thermal control surfaces by calorimetry, the flux of the ambient atomic oxygen environment, the incident solar flux, and the absolute ambient pressure in the payload bay. The IOCM also provides a structure and sample holders for the exposure of passive material samples to the space environment, e.g., thermal cycling, atomic oxygen, and micrometeoroids and/or debris. Some of the more salient results from the flight suggests that the payload bay was slightly contaminated during both the prelaunch phase of the mission and after the deployment of the SYNCOM IV payload. Measurements during the postflight phases, i.e., ferry flight and deintegration processing, show negligible to very low mass deposition. The results of the data are discussed in depth.

Maag, Carl R.↗
DOE OSTI2024

Insight into molecular basis and dynamics of full-length CRaf kinase in cellular signaling mechanisms

Raf kinases play key roles in signal transduction in cells for regulating proliferation, differentiation, and survival. Despite decades of research into functions and dynamics of Raf kinases with respect to other cytosolic proteins, understanding Raf kinases is limited by the lack of their full-length structures at the atomic resolution. Here, we present the first model of the full-length CRaf kinase obtained from artificial intelligence/machine learning algorithms with a converging ensemble of structures simulated by large-scale temperature replica exchange simulations. Our model is validated by comparing simulated structures with the latest cryo-EM structure detailing close contacts among three key domains and regions of the CRaf. Our simulations identify potentially new epitopes of intramolecule interactions within the CRaf and reveal a dynamical nature of CRaf kinases, in which the three domains can move back and forth relative to each other for regulatory dynamics. The dynamic conformations are then used in a docking algorithm to shed insight into the paradoxical effect caused by vemurafenib in comparison with a paradox breaker PLX7904. In this study, we propose a model of Raf-heterodimer/KRas-dimer as a signalosome based on the dynamics of the full-length CRaf.

59 BASIC BIOLOGICAL SCIENCES↗
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 OSTI2022

Impact of Remdesivir Incorporation along the Primer Strand on SARS-CoV-2 RNA-Dependent RNA Polymerase

Remdesivir was the first antiviral drug that received emergency use authorization from the United States Food and Drug Administration and is now formally approved to treat COVID-19. Remdesivir is a nucleotide analogue that targets the RNA-dependent RNA polymerase (RdRp) of coronaviruses, including SARS-CoV-2. The solution of multiple RdRp structures has been one of the main axes of research in the race against the SARS-CoV-2 virus. Several hypotheses of the mechanism of inhibition of RdRp by remdesivir have been proposed, although open questions remain. This work uses molecular dynamics simulations to explore the impact of remdesivir and two analogues as incoming nucleotides and of up to four incorporations of remdesivir along the primer strand on RdRp. The simulation results suggest that the overall structure and the dynamical behavior of RdRp are destabilized by remdesivir and the two analogues in the incoming position. The incorporation of remdesivir along the primer strand impacts specific non-bonded interactions between the nascent RNA and the polymerase subunit, as well as the overall dynamical networks on RdRp. The strongest impact on the structure and dynamics are observed after three incorporations, when remdesivir is located at position -A3, in agreement with previously reported experimental and computational results. Our results provide atomic-level details of the role played by remdesivir on the disruption of RNA synthesis by RdRp and the main drivers of these disruptions.

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

Minimum Free-Energy Shapes of Ag Nanocrystals: Vacuum vs Solution

Here, we use two variants of replica-exchange molecular dynamics (MD) simulations, parallel tempering MD and partial replica exchange MD, to probe the minimum free-energy shapes of Ag nanocrystals containing 100–200 atoms in a vacuum, ethylene glycol (EG) solvent, and EG solvent with a PVP polymer containing 100 repeat units. Our simulations reveal a shape intermediate between a Dh and an Ih, a Dh-Ih, that has distinct structural signatures and magic sizes. We find several prominent features associated with entropy: pure FCC nanocrystals are less common than FCC crystals containing stacking faults, and crystals with the minimum potential energy are not always preferred over the range of relevant temperatures. The shapes of the nanocrystals in solution are influenced by the chemical identities of the solution-phase molecules. Comparing Ag nanocrystal shapes in EG to those in an EG+PVP solution, we find more icosahedra in EG and more decahedra in EG+PVP across all of the nanocrystal sizes probed in this study. At certain critical sizes, nanocrystal shapes can change dramatically with the addition and removal of a single atom or with a change in temperature at a fixed size. The information in our study could be useful in efforts to devise processing routes to achieve selective nanocrystal shapes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular Fe-N 4 Moieties Coupled with Atomic Co-N 4 Sites Toward Improved Oxygen Reduction Performance

Research on high-efficiency and cost-efficient catalysts for oxygen reduction reaction (ORR) is still a vital but challenging issue for commercializing metal–air batteries. Herein, a single-molecule/atom hybrid catalyst is developed to boost the ORR, in which iron phthalocyanine molecules containing molecular Fe-N 4 moieties couple with atomic Co-N 4 sites on the surface of polyhedral carbon. Density functional theory calculations reveal that face-to-face laminated construction of Fe-N 4 and Co-N 4 in the hybrid catalyst can effectively modulate the electronic structure of active iron atoms and reduce the energy barrier of the rate-determining step for ORR. As a result, this hybrid catalyst demonstrates excellent ORR performance, featuring a half-wave potential of 0.904 V, a peak power density of 238.3 mW cm -2 for zinc–air battery, and outstanding electrocatalytic stability. Here, this work offers a distinctive and robust molecular/atomic engineering approach to creating efficient electrocatalysts, advancing the fields of metal–air batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Deuteration Effects on the Physical and Optoelectronic Properties of Donor–Acceptor Conjugated Polymers

The significant differences in scattering cross sections between deuterium and protium are unique to neutron scattering techniques and have been a long-standing area of interest within the neutron scattering community. Researchers have explored selective deuteration to manipulate scattering contrast in soft matter systems, leading to the widespread use of deuterium labeling in materials development. As deuteration changes the atomic mass, it alters physical properties such as molecular volume, polarizability, and polarity, which in turn may affect noncovalent interactions and crystal ordering. Despite previous studies, there remains a limited understanding of how deuteration impacts donor–acceptor (DA) conjugated polymers. To address this, we synthesized deuterated DPP polymers and systematically investigated the effects of side-chain deuteration on their thermal stability, crystal packing, morphology, and optoelectronic properties. We found that deuteration increased the melting and crystallization temperatures of DPP polymers, although it did not significantly alter their morphology, molecular packing, or charge mobility. These properties were assessed by using atomic force microscopy (AFM), X-ray scattering, and thin-film transistor device measurements, respectively, for DPP polymers. Our work shows that deuterium labeling could be a powerful method for controlling scattering length density, enabling neutrons to study the structure and dynamics of conjugated polymers without impacting their electronic performance.

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

Controllable Formation of Threefold-Coordinated Oxygen in Graphene by Low-Energy Ion Implantation

The atomically precise engineering of impurities in graphene and the understanding of their structural and carrier-dependent electronic properties at the nanoscale are crucial for advancing graphene-based nanoelectronics, catalysis, and energy technologies. Here, we demonstrate controllable incorporation of the elusive 3-fold-coordinated O substitutions into graphene using low-energy O + ion implantation under ultrahigh-vacuum conditions. By combining high-resolution scanning tunneling microscopy and spectroscopy (STM/S), bond-resolved noncontact atomic force microscopy techniques, and density functional theory (DFT) calculations, we resolve both the structural and electronic properties of the O-related defects. The STM/S measurements, corroborated by DFT calculations, uncover a characteristic impurity state that is energetically pinned to the Dirac point across different charge-carrier doping regimes. Molecular dynamics simulations further reveal the distribution of implantation-induced configurations and identify the formation of 3-fold-coordinated O dopants. Furthermore, this work provides a viable route to incorporate 3-fold-coordinated O dopants and opens new opportunities for controlled defect engineering in graphene.

3-fold-coordinated oxygen↗
DOE OSTI2024

Elucidating and contrasting the mechanisms for Mg and Ca sulfate ion-pair formation with multi-level embedded quantum mechanics/molecular dynamics simulations

Here, solutions and minerals containing sulfate (SO 4 2- ), and Ca 2+ and Mg 2+ cations, are ubiquitous throughout the lithosphere and are significant components of seawater, thus presenting a prototypical system for the study of strong electrolytes and crystal nucleation mechanisms. However, despite their relative abundance, key questions remain unanswered about the most fundamental atomic-level steps of their mineralization pathways and aqueous dynamics. Here, we carry out enhanced sampling multi-level molecular dynamics (MD) embedded correlated wavefunction theory simulations to elucidate ion-pairing mechanisms for Mg–SO 4 and Ca–SO 4 in concentrated aqueous solution, accurately capturing effects arising from both structural dynamics and electron exchange–correlation. We predict contact-ion-pair formation to be barrierless and highly exoergic for Ca–SO 4 , in agreement with its minimal solubility, whereas for Mg–SO 4 , solvent-shared and contact ion pairs have similar free energies, qualitatively consistent with its higher solubility. Finally, we demonstrate that brief high-temperature pre-equilibration may be utilized to accelerate convergence of free energies in blue-moon-ensemble enhanced-sampling MD.

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

Theory of strong-field sequential double ionization of polyatomic molecules

Understanding of strong-field sequential double ionization (SDI) of molecules by a highly intense infrared laser pulse could be the key to observing electron motion in molecules in the attosecond to femtosecond timescale. Based on a novel density matrix approach for SDI (DM-SDI), a recent theoretical study has shown that SDI can be used as a probe to monitor the changes in vibronic coherence in homonuclear diatomic molecules. In this article, we extend the DM-SDI model to general molecules that could possess permanent dipole moments and arbitrary symmetry. We apply the model to SDI of a water molecule and identify the formation pathway of individual dication states. We further deduce that kinetic energy release spectra from two- and three-body fragments could carry the signature of vibronic coherence between the lowest two states of the water cation. Furthermore, our results suggest that observables from the SDI probe can be interpreted intuitively with only the knowledge of electronic structures of populated ionic states, making the SDI probe to be a highly desirable probing scheme for vibronic coherence in generic molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Quantum Many-Body Theory from a Solution of the N -Representability Problem

Here, in this study, we present a many-body theory based on a solution of the N-representability problem in which the ground-state two-particle reduced density matrix (2-RDM) is determined directly without the many-particle wave function. We derive an equation that re-expresses physical constraints on higher-order RDMs to generate direct constraints on the 2-RDM, which are required for its derivation from an N-particle density matrix, known as N-representability conditions. The approach produces a complete hierarchy of 2-RDM constraints that do not depend explicitly upon the higher RDMs or the wave function. By using the two-particle part of a unitary decomposition of higher order constraint matrices, we can solve the energy minimization by semidefinite programming in a form where the low-rank structure of these matrices can be potentially exploited. We illustrate by computing the ground-state electronic energy and properties of the H 8 ring.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2011

Discovery of Interstellar Anions in Cepheus and Auriga

We report the detection of microwave emission lines from the hydrocarbon anion C6H(-) and its parent neutral C6H in the star-forming region LI251 A (in Cepheus), and the pre-stellar core LI512 (in Auriga). The carbon chain-bearing species C4H, HC3N, HC5N, HC7N, and C3S are also detected in large abundances. The observations of L1251A constitute the first detections of anions and long-chain polyynes and cyanopolyynes (with more than five carbon atoms) in the Cepheus Flare star-forming region, and the first detection of anions in the vicinity of a protostar outside of the Taurus molecular cloud complex, indicating a possible wider importance for anions in the chemistry of star formation. Rotational excitation temperatures have been derived from the HC3N hyperfine structure lines and are found to be 6.2 K for L1251A and 8.7 K for LI5l2. The anion-to-neutral ratios are 3.6% and 4.1%, respectively, which are within the range of values previously observed in the interstellar medium, and suggest a relative uniformity in the processes governing anion abundances in different dense interstellar clouds. This research contributes toward the growing body of evidence that carbon chain anions are relatively abundant in interstellar clouds throughout the Galaxy, but especially in the regions of relatively high density and high depletion surrounding pre-stellar cores and young, embedded protostars.

Cordiner, M. A.↗
DOE OSTI2019

Molecular-Level Design of Heterogeneous Chiral Catalysts

A recent major focus of catalyst research has been on improving the selectivity or atom efficiency of catalytic processes. Perhaps the most subtle and difficult form of selectivity to control is enantioselectivity, a catalyst’s ability to preferentially yield one enantiomer of a chiral product. Since two enantiomers are identical in all respects except for their mirror reflection symmetry, both are formed in non-chiral environments at exactly equal rates; only by introducing some type of chiral modification to the system can the rate of formation of one enantiomer be favored over the other. This has been successfully achieved for homogeneous-phase, organometallic catalysts, but the empirical development of heterogeneous-phase analogs has proven elusive. The design of such heterogeneous-phase catalysts would obviate the need for inefficient purification steps to remove the heavy-metal catalysts. The goal of the proposed research by the Chiral Catalysis Group is to understand the interactions of prochiral reactants on chirally modified surfaces that lead to the preferential formation of one enantiomer over the other. There are two general concepts for designing heterogeneous-phase enantioselective catalysts. The first adopts the strategy used in homogeneous phase, by adding a chiral modifier to the surface. This approach, while conceptually simple has challenges associated with modifying an extended surface. The second concept, which is unique to extended surfaces, is the idea of creating surfaces that are themselves naturally chiral. A fundamental understanding of the interactions that control enantioselectivity will enable the Chiral Catalysis Group to develop the concepts necessary to design novel, heterogenous catalytic materials. Enantioselectivity is also central to catalytic processes that underpin the pharmaceutical and agrochemical industries. With funding from the Department of Energy, the Chiral Catalysis Group has played a leading role in developing a fundamental understanding of chiral catalysis and of the adsorbate structures, interactions and kinetic effects that control heterogenous-phase enantioselectivity. The proposed research program is organized conceptually around three types of enantiospecific interactions: interactions with naturally chiral surfaces; those with modified surfaces categorized as chirally templated surfaces; and one-to-one interactions of prochiral molecules with chiral surface modifiers. Our specific goals are to: (1) determine the structures of chiral surfaces, prochiral reactants, chiral products and chiral modifiers; (2) understand the enantiospecific interactions that occur between prochiral reactants and chiral products with chirally modified surfaces; (3) develop methods for detection of enantiospecific surface phenomena; and (4) elucidate enantioselective reaction pathways and the interactions that control them. The range of expertise within the Chiral Catalysis Group includes molecular-scale theory, surface science experiments, state-of-the-art spectroscopic, kinetic and atomic-scale imaging capabilities carried out in vacuum and under more-realistic liquid- and gas-phase environments complemented by catalytic measurements. By fully integrating such a breadth of approaches and extensive expertise, the Chiral Catalysis Group will, for the first time, be able to understand, rationally design and then control active sites for optimal enantioselectivity. In addition to having a significant impact on the field of chiral catalysis, the proposed research program will further our general understanding of the molecular-scale origins of catalytic selectivity.

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

Department of Energy Project DE-SC0012198: Probing Electron-Molecule Scattering with HHG: Theory and Experiment (Final Report)

The study of the motion of atoms in molecules is important to understanding many areas of physical and life sciences. Such motion occurs on many different times scales, with electronic motion occurring on a sub-femtosecond time scale, and simple vibrational motion in tens to hundreds of femtoseconds. One way to follow such processes in real time is by the use of shortpulsed lasers, and in particular by studying time-resolved photoionization and the related process of high-harmonic generation (HHG). Thus there has been much effort to develop the tools necessary to probe molecular systems using short pulse lasers and understanding the sensitivity of the different possible probes to the time dependent geometric structure as well as the electronic structure of molecules. Our research has particularly focused on the connection between high-field processes and the more traditional weak field photoionization processes. Strong field and weak field processes can be connected through models that involve the same matrix elements. We have demonstrated in our study of HHG from SF 6 that the spectrum is sensitive to the interplay between the angular dependence of the ionization step and recombination step. In our study of rescattering spectroscopy, we have shown that with a combination of experiment and theory, we can use this high-field spectroscopy to determine molecular structure in molecules such as C 2 H 4 . We have also developed new computational tools based on overset grids to enable studies on larger molecular systems which use much more robust numerical approaches so that the resulting code can be a tool that non-specialists can use to study related systems.

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

Crystallographic and temperature effects in low-energy collisions for plasma–material interactions

The interaction of plasma with materials is critical for the fundamental understanding of non-equilibrium processes and their wide application. Recent achievements in fusion energy research emphasize the importance of this problem. Because modelling and predicting plasma–material interactions (PMI) require considering a tremendous number of single plasma particle–surface interaction events, this type of modelling was, until recently, possible only within a binary collision approximation (BCA). The BCA approach considers materials as uniform temperature-insensitive media. The research presented here addresses the PMI problem within an atomistic-based approach for the first time at a statistically significant level. Approximately 10 5 molecular dynamics trajectories were generated for 100 eV deuterium ions interacting with a tungsten surface. The research demonstrated the critical importance of incorporating the discrete lattice structure of matter into the model. Deuterium penetration depth and fraction of backscattered deuterium ions strongly depend on the surface orientation, the impact incident ion directions, and the ions’ initial positions. On average, the BCA-based calculations underestimate the penetration depth by a factor of two or more, and the fraction of backscattered atoms is overestimated by a factor of two or more. As a result, increasing the material temperature from 500 to 3000 K reduces penetration depth by about 30% and increases the fraction of backscattered atoms by 250%.

36 MATERIALS SCIENCE↗
NASA NTRS2015

Effects of Atomic-Scale Structure on the Fracture Properties of Amorphous Carbon - Carbon Nanotube Composites

The fracture of carbon materials is a complex process, the understanding of which is critical to the development of next generation high performance materials. While quantum mechanical (QM) calculations are the most accurate way to model fracture, the fracture behavior of many carbon-based composite engineering materials, such as carbon nanotube (CNT) composites, is a multi-scale process that occurs on time and length scales beyond the practical limitations of QM methods. The Reax Force Field (ReaxFF) is capable of predicting mechanical properties involving strong deformation, bond breaking and bond formation in the classical molecular dynamics framework. This has been achieved by adding to the potential energy function a bond-order term that varies continuously with distance. The use of an empirical bond order potential, such as ReaxFF, enables the simulation of failure in molecular systems that are several orders of magnitude larger than would be possible in QM techniques. In this work, the fracture behavior of an amorphous carbon (AC) matrix reinforced with CNTs was modeled using molecular dynamics with the ReaxFF reactive forcefield. Care was taken to select the appropriate simulation parameters, which can be different from those required when using traditional fixed-bond force fields. The effect of CNT arrangement was investigated with three systems: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. For each arrangement, covalent bonds are added between the CNTs and AC, with crosslink fractions ranging from 0-25% of the interfacial CNT atoms. The SWNT and MWNT array systems represent ideal cases with evenly spaced CNTs; the SWNT bundle system represents a more realistic case because, in practice, van der Waals interactions lead to the agglomeration of CNTs into bundles. The simulation results will serve as guidance in setting experimental processing conditions to optimize the mechanical properties of CNT composites.

Jensen, Benjamin D.↗
DOE OSTI2026

Environmental Contributions to Proton Sharing in Protein Low-Barrier Hydrogen Bonds

Hydrogen bonds (H-bonds) are central to biomolecular structure and dynamics. Although H-bonds are typically characterized by well-defined proton positions, proton delocalization has been proposed to play a role in facilitating enzyme catalysis and allostery in some systems. Experimentally locating protons is difficult, hampering the study of proton mobility in H-bonds. We used neutron crystallography, atomic resolution X-ray bond length analysis, and large quantum mechanics/molecular mechanics-Born–Oppenheimer molecular dynamics (QM/MM-BOMD) simulations to comprehensively characterize the shared proton/deuteron in a Glu–Asp low-barrier hydrogen bond (LBHB) in the bacterial protein YajL that is a conventional H-bond in the homologous disease-associated human protein DJ-1. X-ray bond length analysis of protiated and perdeuterated DJ-1 and YajL shows no significant effect of deuteron substitution on these carboxylic acid-carboxylate H-bonds but does reveal an effect at the active site glutamic acid near a cysteine thiolate. Residues in an H-bonded network that might favor LBHB formation in YajL were interrogated by the mutation of homologous residues in DJ-1. A distal DJ-1 substitution increases proton delocalization in the Glu–Asp H-bond, demonstrating that mutations within extended H-bond networks can modulate proton transfer barriers in carboxylic acid-carboxylate H-bonds. In addition, proton mobility in the H-bond is correlated with dimer-spanning motions in the QM/MM-BOMD simulations of YajL and DJ-1. Our results show that proton delocalization can be tuned using combined bioinformatic, structural, and computational information, opening the possibility of using engineered proton delocalization as a probe of H-bonding environments and as a tool to test hypotheses about LBHB function.

Lin, Jiusheng [University of Nebraska, Lincoln, NE↗
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↗