Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic properties”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Examination of Short- and Long-Range Atomic Order Nanocrystalline SiC and Diamond by Powder Diffraction Methods

The real atomic structure of nanocrystals determines unique, key properties of the materials. Determination of the structure presents a challenge due to inherent limitations of standard powder diffraction techniques when applied to nanocrystals. Alternate methodology of the structural analysis of nanocrystals (several nanometers in size) based on Bragg-like scattering and called the "apparent lattice parameter" (alp) is proposed. Application of the alp methodology to examination of the core-shell model of nanocrystals will be presented. The results of application of the alp method to structural analysis of several nanopowders were complemented by those obtained by determination of the Atomic Pair Distribution Function, PDF. Based on synchrotron and neutron diffraction data measured in a large diffraction vector of up to Q = 25 Angstroms(exp -1), the surface stresses in nanocrystalline diamond and SiC were evaluated.

Palosz, B.↗
DOE OSTI2025

Quantum annealing-assisted lattice optimization

High Entropy Alloys (HEAs) have drawn great interest due to their exceptional properties compared to conventional materials. The configuration of HEA system is considered a key to their superior properties, but exhausting all possible configurations of atom coordinates and species to find the ground energy state is extremely challenging. In this work, we proposed a quantum annealing-assisted lattice optimization (QALO) algorithm, which is an active learning framework that integrates the Field-aware Factorization Machine (FFM) as the surrogate model for lattice energy prediction, Quantum Annealing (QA) as an optimizer and Machine Learning Potential (MLP) for ground truth energy calculation. By applying our algorithm to the NbMoTaW alloy, we reproduced the Nb depletion and W enrichment observed in bulk HEA. We found our optimized HEAs to have superior mechanical properties compared to the randomly generated alloy configurations. Our algorithm highlights the potential of quantum computing in materials design and discovery, laying a foundation for further exploring and optimizing structure-property relationships.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Large Single Crystals of V 1– x Mo x O 2 from a Two-Step Chemical Vapor Transport Synthesis

This article details the synthesis of single crystals with the chemical formulas V 1-x Mo x O 2 ( 0 ≤ x ≤ 0.6) using the chemical vapor transport method at high temperatures (850 to 950 °C). These compounds are of interest due to the special relationship between their structure and electronic properties. Understanding this relationship will help guide future efforts in designing energy-saving devices. While some of these compounds have been synthesized before, this is the first time their synthesis has been reported for such a wide range of compounds in the form of bulk single crystals. The single crystal form is required for detailed characterization of the compounds. The materials were characterized chemically and structurally in order to validate the precise chemical compositions that were reported. Characterization methods included x-ray diffraction, inductively-coupled plasma atomic emission spectroscopy, and electron dispersive spectroscopy. The synthesis method is tunable, and the introduction of multiple steps both overcome the previous chemical challenges precluding synthesis of high molybdenum content crystals and allow for the creation of larger single crystals that were otherwise possible. The crystals synthesized this way are suitable for specialized measurements for determining the precise physical properties and atomic structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Learning from metastable symmetric-tilt grain boundaries using physics-based descriptors

Grain boundaries (GBs) govern critical properties of polycrystalline materials. Although significant advancements have been made in characterizing minimum energy and ordered GBs, real GBs are seldom found in such well-defined states. This diversity of atomic arrangements in metastable states makes it challenging to establish structure-property relationships with physical insights. Here, to address this challenge, we use data-driven methods to explore these relationships and examine the underlying physics. In this study, we utilize a large atomistic database (~5000) of minimum energy and metastable states of symmetric-tilt copper GBs, combined with physically motivated local atomic environment (LAE) descriptors [strain functional descriptors (SFDs)], to predict GB properties and gain physical insights. Our regression models exhibit robust predictive capabilities using only 19 descriptors, generalizing to atomic environments in nanocrystals. A significant highlight of our work is the integration of an unsupervised method with SFDs to elucidate LAEs at GBs and their role in determining properties. The model, trained on these minimum energy and metastable GBs using SFDs, predicts the properties of unseen nanocrystals with good accuracy. Our research underscores the role of a physics-based representation of LAEs and the efficacy of data-driven methods in establishing GB structure-property relationships.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Probing individual single atom electrocatalyst sites by advanced analytical scanning transmission electron microscopy

Single atom electrocatalysts (SAEs) are promising next-generation materials for promoting a variety of important reactions, such as the oxygen reduction, nitrogen reduction, and CO 2 reduction reactions. While bulk characterization techniques such as X-ray absorption spectroscopy and Mössbauer spectroscopy have significantly enhanced our understanding of these catalysts, direct probing of individual single metal atom sites at the atomic scale is necessary to understand local variations in the properties of these sites and accelerate design and synthesis of improved SAEs. Aberration-corrected scanning transmission electron microscopy (STEM) has become a powerful tool for providing this type of atomic-scale information about SAE metal sites. These sites are typically unstable under the electron beam, however, which, in combination with conventional acquisition methods and detectors, has limited the type and quantity of information obtainable by spectroscopic STEM techniques. Here, we map multiple individual SAE metal sites in a nitrogen-doped carbon containing atomically dispersed Fe and Re (FeReNC) at the atomic scale by direct electron detection electron energy-loss spectroscopy (EELS). Direct electron detection provides an improved signal-to-noise ratio over conventional scintillator-based detectors and enables detection and real space localization of weak signals. In addition, we demonstrate an automated method for identification of metal atom positions, placement of the probe on these sites, and simultaneous EELS and energy dispersive X-ray spectroscopic (EDS) signal acquisition. This simultaneous acquisition of EELS and EDS provides access to the composition and bonding of a wide range of SAE metal sites. In this study, focusing the probe directly on the metal sites also increases the relevant data acquisition rate by more than an order of magnitude over two-dimensional mapping, enabling improved statistical measurements of site properties. The versatility, sensitivity, and speed that these techniques provide enhances our ability to probe the local elemental and chemical environment of a large number of individual SAE metal site structures at the atomic scale, enabling an improved understanding of the variations in the local properties of these electrocatalysts to be gained. As a result, significantly increased information about individual metal sites will be available to future electrochemical studies through these techniques, accelerating the development of advanced SAEs.

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

Synergistic Combination of Living Ring-Opening Metathesis Polymerization and Atom Transfer Radical Polymerization to Synthesize Structurally Tailored and Engineered Macromolecular Networks

Structurally tailored and engineered macromolecular (STEM) networks are attractive materials for soft robotics, stretchable electronics, tissue engineering, and 3D printing due to their tunable properties. To date, STEM networks have been synthesized by atom transfer radical polymerization (ATRP) or the combination of reversible addition–fragmentation chain-transfer (RAFT) polymerization and ATRP. RAFT polymerization could have limited selectivity with ATRP inimer sites that can participate in radical-transfer processes. On the other hand, living ring-opening metathesis polymerization (ROMP) can produce a polymeric network with latent ATRP initiator sites in high selectivity. Herein, for the first time, we report the syntheses of STEM zero-generation (STEM-0) networks using a monomer, a cross-linker, and an ATRP/ROMP inimer via living ROMP, followed by their modification using a second monomer via ATRP to synthesize STEM first-generation (STEM-1) networks. The mechanical property and swelling capacity analyses of these networks were carried out. A change in mechanical properties and swelling capacity of these networks was observed due to their structural modification.

Absorption↗
DOE OSTI2020

Characterization of organic aerosol across the global remote troposphere: a comparison of ATom measurements and global chemistry models

The spatial distribution and properties of submicron organic aerosol (OA) are among the key sources of uncertainty in our understanding of aerosol effects on climate. Uncertainties are particularly large over remote regions of the free troposphere and Southern Ocean, where very few data have been available and where OA predictions from AeroCom Phase II global models span 2 to 3 orders of magnitude, greatly exceeding the model spread over source regions. The (nearly) pole-to-pole vertical distribution of non-refractory aerosols was measured with an aerosol mass spectrometer onboard the NASA DC-8 aircraft as part of the Atmospheric Tomography (ATom) mission during the Northern Hemisphere summer (August 2016) and winter (February 2017). This study presents the first extensive characterization of OA mass concentrations and their level of oxidation in the remote atmosphere. OA and sulfate are the major contributors by mass to submicron aerosols in the remote troposphere, together with sea salt in the marine boundary layer. Sulfate was dominant in the lower stratosphere. OA concentrations have a strong seasonal and zonal variability, with the highest levels measured in the lower troposphere in the summer and over the regions influenced by biomass burning from Africa (up to 10 µg sm -3 ). Lower concentrations (~0.1–0.3 µg sm -3 ) are observed in the northern middle and high latitudes and very low concentrations (<0.1 µg sm -3 ) in the southern middle and high latitudes. The ATom dataset is used to evaluate predictions of eight current global chemistry models that implement a variety of commonly used representations of OA sources and chemistry, as well as of the AeroCom-II ensemble. The current model ensemble captures the average vertical and spatial distribution of measured OA concentrations, and the spread of the individual models remains within a factor of 5. These results are significantly improved over the AeroCom-II model ensemble, which shows large overestimations over these regions. However, some of the improved agreement with observations occurs for the wrong reasons, as models have the tendency to greatly overestimate the primary OA fraction and underestimate the secondary fraction. Measured OA in the remote free troposphere is highly oxygenated, with organic aerosol to organic carbon (OA/OC) ratios of ~2.2–2.8, and is 30 %–60 % more oxygenated than in current models, which can lead to significant errors in OA concentrations. The model–measurement comparisons presented here support the concept of a more dynamic OA system as proposed by Hodzic et al. (2016), with enhanced removal of primary OA and a stronger production of secondary OA in global models needed to provide better agreement with observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

First-principles and cluster expansion study of the effect of magnetism on short-range order in Fe–Ni–Cr austenitic stainless steels

Short-range order (SRO), the regular and predictable arrangement of atoms over short distances, alters the mechanical properties of technologically relevant structural materials such as medium/high entropy alloys and austenitic stainless steels. In this study, we present a generalized spin cluster expansion (CE) model and show that magnetism is a primary factor influencing the level of SRO present in austenitic Fe-Ni-Cr alloys. The spin CE consists of a chemical cluster expansion combined with an Ising model for Fe-Ni-Cr austenitic alloys. It explicitly accounts for local magnetic exchange interactions, thereby capturing the effects of finite temperature magnetism on SRO. Model parameters are obtained by fitting to a first-principles data set comprising both chemically and magnetically diverse FCC configurations. The magnitude of the magnetic exchange interactions are found to be comparable to the chemical interactions. Compared to a conventional implicit magnetism CE built from only magnetic ground state configurations, the spin CE shows improved performance on several experimental benchmarks over a broad spectrum of compositions, particularly at higher temperatures due to the explicit treatment of magnetic disorder. We find that SRO is strongly influenced by alloy Cr content, since Cr atoms prefer to align antiferromagnetically with nearest neighbors but become magnetically frustrated with increasing Cr concentration. Using the spin CE, we predict that increasing the Cr concentration in typical austenitic stainless steels promotes the formation of SRO and increases order-disorder transition temperatures. Furthermore, this study underscores the significance of considering magnetic interactions explicitly when exploring the thermodynamic properties of complex transition metal alloys. It also highlights guidelines for customizing SRO through adjustments of alloy composition.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Computational and Experimental Determination of the Properties, Structure, and Stability of Peptoid Nanosheets and Nanotubes

Peptoids (N-substituted glycines) are a group of highly controllable peptidomimetic polymers. Amphiphilic diblock peptoids have been engineered to assemble crystalline nanospheres, nanofibrils, nanosheets, and nanotubes with biochemical, biomedical, and bioengineering applications. The mechanical properties of peptoid nanoaggregates and their relationship to the emergent self-assembled morphologies have been relatively unexplored and are critical for the rational design of peptoid nanomaterials. In this work, we consider a family of amphiphilic diblock peptoids consisting of a prototypical tube-former (Nbrpm6Nc6, a NH 2 -capped hydrophobic block of six N-((4-bromophenyl)methyl)glycine residues conjugated to a polar NH 3 (CH 2 ) 5 CO tail), a prototypical sheet-former (Nbrpe6Nc6, where the hydrophobic block comprises six N-((4-bromophenyl)ethyl)glycine residues), and an intermediate sequence that forms mixed structures ((NbrpeNbrpm)3Nc6). We combine all-atom molecular dynamics simulations and atomic force microscopy to determine the mechanical properties of the self-assembled 2D crystalline nanosheets and relate these properties to the observed self-assembled morphologies. We find good agreement between our computational predictions and experimental measurements of Young’s modulus of crystalline nanosheets. A computational analysis of the bending modulus along the two axes of the planar crystalline nanosheets reveals bending to be more favorable along the axis in which the peptoids stack by interdigitation of the side chains compared to that in which they form columnar crystals with π-stacked side chains. Here, we construct molecular models of nanotubes of the Nbrpm6Nc6 tube-forming peptoid and predict a stability optimum in good agreement with experimental measurements. A theoretical model of nanotube stability suggests that this optimum is a free energy minimum corresponding to a “Goldilocks” tube radius at which capillary wave fluctuations in the tube wall are minimized.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomic-Scale Surface Studies of Bulk Metallic Glasses. Final Report

Bulk metallic glasses (BMGs) are of both scientific and technological interest because the absence of periodic atomic arrangements provides them with unique physical, chemical, and mechanical properties. Their high strength, superior elasticity, and an ability to be easily formed into virtually unlimited shapes with feature sizes from centimeters to Angstroms by blow molding and thermoplastic forming makes them an attractive choice for more and more practical applications and products. Due to their complex internal structure, however, experiments that yield insight into their exact atomic arrangements have been scarce. As a result, glass physics is one of the last remaining unexplored fields of materials science despite the scientific and technological importance of glasses in our daily lives, and the question how to characterize and control matter away from equilibrium, as glasses are, was listed as one of five Grand Challenges in a recent DoE report. The main reason for the slow progress in glass physics is the lack of experimental tools that enable access to atomic-scale structural and behavioral information for disordered materials. Such atomic-scale knowledge is mandatory to establish structure-property relationships that ultimately could allow to custom-design alloys featuring specific desired characteristics. With no such relationships available, theory development in glass remains basic and the few that exist are often untested. The aim of this research was to enable progress in our understanding of BMGs by developing a new approach that will allow a meaningful application of local surface science methods to specially prepared BMG samples to obtain a wealth of quantitative information on their atomic arrangements. Key was the availability of specially prepared samples whose surfaces feature large atomically flat terraces despite being entirely amorphous, which we have produced from a Pt 57.5 Cu 14.7 Ni 5.3 P 22.5 alloy (‘Pt-BMG’) both under ambient conditions as well as in ultrahigh vacuum using a unique setup that has been specially developed within this grant. Our approach starts with the in-situ preparation of oxide crystals that are terminated by large terraces, from which exact mirror images out of BMG will be produced using thermoplastic forming (TPF); for the research within this grant, we have successfully used (001)-oriented SrTiO 3 single crystal surfaces as well as (100)-, (110)-, and (111)-oriented single crystals made from LaAlO 3 . Since the resulting BMG replicas display all features of the original crystal with sub-Angstrom fidelity, thereby mimicking the original crystal’s termination by atomically flat terraces without being crystalline themselves, they are ideally suited for further investigation. The following atomic-scale local studies were then carried within this proposal: (i) high-resolution surface imaging and local spectroscopy using scanning probe microscopy, which showed disordered atom-like features and revealed changes in the gradient of the local surface potential on a 1-2 nm length scale; (ii) characterization of atomic-scale plastic flow with affected volumes as low as 1000 atoms, which showed local hardness near or above the theoretically predicted maximum and, once plastic deformation was initiated, homogeneous flow of the atoms involved; (iii) characterization of surface relaxation processes and the onset of crystallization induced by annealing, which showed that upon heating over the material’s glass transition temperature, the surface rearranges and relaxes towards a more stable, denser packed glass, which increases on-terrace surface roughness, while surface tension smoothens step edges; and (iv) studies that investigate the dependence of the material’s mechanical properties and structure on processing parameters, revealing that relaxed glasses get denser, harder, and more elastic. In combination, this information allows to combine structural models with mechanical properties and preparation history, thereby facilitating the development of preparation-structure-property relationships for metallic glasses. With the availability of such information, bulk metallic glasses can be further optimized to be used in more and more applications in industry.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electron Tomography and Machine Learning for Understanding the Highly Ordered Structure of Leafhopper Brochosomes

Insects known as leafhoppers (Hemiptera: Cicadellidae) produce hierarchically structured nanoparticles known as brochosomes that are exuded and applied to the insect cuticle, thereby providing camouflage and anti-wetting properties to aid insect survival. Although the physical properties of brochosomes are thought to depend on the leafhopper species, the structure–function relationships governing brochosome behavior are not fully understood. Brochosomes have complex hierarchical structures and morphological heterogeneity across species, due to which a multimodal characterization approach is required to effectively elucidate their nanoscale structure and properties. In this work, we study the structural and mechanical properties of brochosomes using a combination of atomic force microscopy (AFM), electron microscopy (EM), electron tomography, and machine learning (ML)-based quantification of large and complex scanning electron microscopy (SEM) image data sets. This suite of techniques allows for the characterization of internal and external brochosome structures, and ML-based image analysis methods of large data sets reveal correlations in the structure across several leafhopper species. Our results show that brochosomes are relatively rigid hollow spheres with characteristic dimensions and morphologies that depend on leafhopper species. Nanomechanical mapping AFM is used to determine a characteristic compression modulus for brochosomes on the order of 1–3 GPa, which is consistent with crystalline proteins. Altogether, this work provides an improved understanding of the structural and mechanical properties of leafhopper brochosomes using a new set of ML-based image classification tools that can be broadly applied to nanostructured biological materials.

Chemical structure↗
DOE OSTI2024

Spectroscopy-guided discovery of three-dimensional structures of disordered materials with diffusion models

Spectroscopy techniques such as x-ray absorption near edge structure (XANES) provide valuable insights into the atomic structures of materials, yet the inverse prediction of precise structures from spectroscopic data remains a formidable challenge. In this study, we introduce a framework that combines generative artificial intelligence models with XANES spectroscopy to predict three-dimensional atomic structures of disordered systems, using amorphous carbon (a-C) as a model system. In this work, we introduce a new framework based on the diffusion model, a recent generative machine learning method, to predict 3D structures of disordered materials from a target property. For demonstration, we apply the model to identify the atomic structures of a-C as a representative material system from the target XANES spectra. We show that conditional generation guided by XANES spectra reproduces key features of the target structures. Furthermore, we show that our model can steer the generative process to tailor atomic arrangements for a specific XANES spectrum. Finally, our generative model exhibits a remarkable scale-agnostic property, thereby enabling generation of realistic, large-scale structures through learning from a small-scale dataset (i.e. with small unit cells). Our work represents a significant stride in bridging the gap between materials characterization and atomic structure determination; in addition, it can be leveraged for materials discovery in exploring various material properties as targeted.

36 MATERIALS SCIENCE↗
NASA NTRS2008

Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Data Processing, Sky Maps, and Basic Results

The Wilkinson Microwave Anisotropy Probe (WMAP) is a Medium-Class Explorer (MIDEX) satellite aimed at elucidating cosmology through full-sky observations of the cosmic microwave background (CMB). The WMAP full-sky maps of the temperature and polarization anisotropy in five frequency bands provide our most accurate view to date of conditions in the early universe. The multi-frequency data facilitate the separation of the CMB signal from foreground emission arising both from our Galaxy and from extragalactic sources. The CMB angular power spectrum derived from these maps exhibits a highly coherent acoustic peak structure which makes it possible to extract a wealth of information about the composition and history of the universe. as well as the processes that seeded the fluctuations. WMAP data have played a key role in establishing ACDM as the new standard model of cosmology (Bennett et al. 2003: Spergel et al. 2003; Hinshaw et al. 2007: Spergel et al. 2007): a flat universe dominated by dark energy, supplemented by dark matter and atoms with density fluctuations seeded by a Gaussian, adiabatic, nearly scale invariant process. The basic properties of this universe are determined by five numbers: the density of matter, the density of atoms. the age of the universe (or equivalently, the Hubble constant today), the amplitude of the initial fluctuations, and their scale dependence. By accurately measuring the first few peaks in the angular power spectrum, WMAP data have enabled the following accomplishments: Showing the dark matter must be non-baryonic and interact only weakly with atoms and radiation. The WMAP measurement of the dark matter density puts important constraints on supersymmetric dark matter models and on the properties of other dark matter candidates. With five years of data and a better determination of our beam response, this measurement has been significantly improved. Precise determination of the density of atoms in the universe. The agreement between the atomic density derived from WMAP and the density inferred from the deuterium abundance is an important test of the standard big bang model. Determination of the acoustic scale at redshift z = 1090. Similarly, the recent measurement of baryon acoustic oscillations (BAO) in the galaxy power spectrum (Eisenstein et al. 2005) has determined the acoustic scale at redshift z approx. 0.35. When combined, these standard rulers accurately measure the geometry of the universe and the properties of the dark energy. These data require a nearly flat universe dominated by dark energy consistent with a cosmological constant. Precise determination of the Hubble Constant, in conjunction with BAO observations. Even when allowing curvature (Omega(sub 0) does not equal 1) and a free dark energy equation of state (w does not equal -1), the acoustic data determine the Hubble constant to within 3%. The measured value is in excellent agreement with independent results from the Hubble Key Project (Freedman et al. 2001), providing yet another important consistency test for the standard model. Significant constraint of the basic properties of the primordial fluctuations. The anti-correlation seen in the temperature/polarization (TE) correlation spectrum on 4deg scales implies that the fluctuations are primarily adiabatic and rule out defect models and isocurvature models as the primary source of fluctuations (Peiris et al. 2003).

Weiland, J.L.↗
DOE OSTI2025

Accelerating the discovery of low-energy structure configurations: A computational approach that integrates first-principles calculations, Monte Carlo sampling, and Machine Learning

Finding Minimum Energy Configurations (MECs) is essential in fields such as physics, chemistry, and materials science, as they represent the most stable states of the systems. In particular, identifying such MECs in multi-component alloys considered candidate PFMs is key because it determines the most stable arrangement of atoms within the alloy, directly influencing its phase stability, structural integrity, and thermo-mechanical properties. However, since the search space grows exponentially with the number of atoms considered, obtaining such MECs using computationally expensive first-principles DFT calculations often results in a cumbersome task. To escape the above compromise between physical fidelity and computational efficiency, we have developed a novel physics-based data-driven approach that combines Monte Carlo sampling, first-principles DFT calculations, and Machine Learning to accelerate the discovery of MECs in multi-component alloys. More specifically, we have leveraged well-established Cluster Expansion (CE) techniques with Local Outlier Factor models to establish strategies that enhance the reliability of the CE method. In this work, we demonstrated the capabilities of the proposed approach for the particular case of a tungsten-based quaternary high-entropy alloy. However, the method is applicable to other types of alloys and enables a wide range of applications.

36 MATERIALS SCIENCE↗
NASA NTRS2019

High-Performance Computing Optimization for Aladyn – Adaptive Neural Network Molecular Dynamics Mini-Application

This report provides a description and performance evaluation of the optimization techniques for high performance computing (HPC) implementation of the open source Computational Materials mini-application Aladyn (https://github.com/nasa/aladyn). Aladyn is a basic molecular dynamics code written in FORTRAN 2003, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. While achieving orders of magnitude faster computational performance than DFT, the ANN-based approach was still very computationally demanding compared to the conventional approach of using empirically fitted energy functions. After its initial development, Aladyn was evaluated and optimized by experts at the NASA Advanced Supercomputing (NAS) division to exploit modern supercomputer architectures. The code has been optimized for execution on multicore central processing units (CPUs), including Intel® Skylake microarchitecture, and on graphic accelerators, such as Nvidia® V100 graphic processing units (GPUs), using Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) programming interfaces. The optimization achieved a speedup of 4.7 times the baseline version on CPU performance and an additional 2.4 times on CPU+GPU performance. Atomistic computer simulations are a fundamental tool in materials research to model material properties form physics-based first principles. Atomic interaction, governed by Quantum Mechanics (QM) require sophisticated and highly computationally demanding mathematical models to calculate [1]. Classical methods use approximate functional forms, empirically fitted through a set of variable parameters to emulate atomic energies as direct functions of atomic coordinates [2]. While empirical potentials are computationally much simpler, allowing simulations of large-scale systems of up to a trillion (1012) atoms [3], they are substantially less accurate compared to quantum calculations and applicable only to very specific atomic configurations or predefined crystallographic phases. A recently suggested approach is to use heuristic machine learning methods [4], such as those based on Adaptive Neural Networks (ANNs) to predict atomic energies, after being trained on a sufficiently large database of QM-calculated structures [5,6]. This approach reduces significantly the computational complexity, allowing for simulations of orders of magnitude larger systems compared to QM-based methods without compromising accuracy. Still, compared to classical methods using empirical energy functions, ANN methods remain two- to three orders of magnitude more computationally demanding. Hence, the computational cost of simulations, together with the need for extensive training of ANNs, still makes the practical implementation of ANN-based methods quite challenging. The purpose of the Aladyn mini-application software [7], available as open source at https://github.com/nasa/aladyn, is to be a testbed for exploring possible optimization strategies to develop highly scalable parallel algorithms for ANN-based atomistic simulations. Aladyn is aimed at utilizing the architecture of the high-end modern highperformance computing (HPC) hardware based on multicore central processing units (CPUs) equipped with graphic processing unit (GPU) accelerators. Specifically, the goal is to optimize the performance on a single HPC compute node, before implementing scaling to multi-node parallelization using message passing interface (MPI). At the same time, the open source code of Aladyn can serve as a training model for students and professors in academia.

Yamakov, Vesselin I.↗
DOE OSTI2023

Lattice dynamics, sound velocities, and atomic environments of szomolnokite at high pressure

Complex mixtures of sulfates, silicates, and ice have been observed in a variety of planetary environments on Earth, Mars and the icy satellites of the solar system. Characterizing the properties of the corresponding compositional endmembers is important for understanding the interiors of a range of planetary bodies in which these phases are observed. To measure the electronic and vibrational properties of the pure ferrous iron endmember of the kieserite group, szomolnokite, (FeSO 4 ∙H 2 O), we have performed synchrotron 57 Fe nuclear resonant inelastic and forward scattering experiments in the diamond-anvil cell up to 14.5 GPa. This pressure range covers depths within Earth’s interior relevant to sulfur cycling in subduction zones and the range of pressures expected within icy satellite interiors. We find evidence of crystal lattice softening, changes in elastic properties, and changes in the electric field gradients of iron atoms associated with two structural transitions occurring within the experimental pressure range. Here we apply these findings to icy satellite interiors, including discussion of elastic properties, modeling of ice-sulfate aggregates, and implications for tidal observations.

58 GEOSCIENCES↗