Search NASA⌕ Search

SEARCH · Search NASA

Results for “Eam”

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 19 records

Graph-EAM: An Interpretable and Efficient Graph Neural Network Potential Framework

The development of deep learning interatomic potentials has enabled efficient and accurate computations in quantum chemistry and materials science, circumventing computationally expensive ab initio calculations. However, the huge number of learnable parameters in deep learning models and their complex architectures hinder physical interpretability and affect the robustness of the derived potential. In this work, we propose graph-EAM, a lightweight graph neural network (GNN) inspired by the empirical embedded atom method to model the interatomic potential of single-element structures. Four material systems: platinum, niobium, silicon, and amorphous-carbon, for which quantum simulation data sets are publicly available, are examined to demonstrate that graph-EAM can achieve high energy and force prediction accuracy-comparable or better than existing state-of-the-art machine learning models-with much fewer parameters. It is also shown that the explicit inclusion of the angular information via three-body atomic density increases the prediction accuracy. In conclusion, the accuracy and efficiency of potentials obtained from graph-EAM can help accelerate the molecular dynamics simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Disordering and H-embrittlement of Pb(110) surface using embedded atom method and molecular dynamics

The dynamics and structure of the crystal melt interface are of great importance in crystal growth, melting, and other properties of the solid-liquid interface. The Embedded Atom Method (EAM), a semiempirical method developed by Daw and Baskes for determining the energetics of atoms in a bulk environment, is a useful method in the computer simulation. The EAM functions of PB in conjunction with the MD code are employed to study melting and disordering of the Pb(100) surface. The EAM functions are fitted to the bulk experimental data at zero temperature, and it is not clear if it will behave correctly at higher temperatures. In fact, it is shown that use of a similar type of potential at higher temperatures might result in a negative coefficient of thermal expansion. The primary purpose here is to perform MD simulation combined with the EAM functions of PB to study melting of the Pb(100) surface. This provides an accurate test of the EAM functions at higher temperatures. In particular, we calculate the two dimensional structure factor, the two dimensional distribution function, and density profile as a function of temperature up to the melting point of lead.

Karimi, Majid↗

Energetics and structural properties of twist grain boundaries in Cu

Structural and energetics properties of atoms near a grain boundary are of great importance from theoretical and experimental standpoints. From various experimental work it is concluded that diffusion at low temperatures at polycrystalline materials take place near grain boundary. Experimental and theoretical results also indicate changes of up to 70 percent in physical properties near a grain boundary. The Embedded Atom Method (EAM) calculations on structural properties of Au twist grain boundaries are in quite good agreement with their experimental counterparts. The EAM is believed to predict reliable values for the single vacancy formation energy as well as migration energy. However, it is not clear whether the EAM functions which are fitted to the bulk properties of a perfect crystalline solid can produce reliable results on grain boundaries. One of the objectives of this work is to construct the EAM functions for Cu and use them in conjunction with the molecular static simulation to study structures and energetics of atoms near twist grain boundaries in Cu. This provides tests of the EAM functions near a grain boundary. In particular, we determine structure, single vacancy formation energy, migration energy, single vacancy activation energy, and interlayer spacing as a function of distance from grain boundary. Our results are compared with the available experimental and theoretical results from grain boundaries and bulk.

Karimi, Majid↗

The Evolution of Tropospheric Temperature Field and its Relationship With The Onset of Asian Summer Monsoon

The mean state and year-to-year variations of the tropospheric temperature fields and their relationship with the establishment of the summertime East Asian monsoon (EAM) and the Indian monsoon (INM) are studied using the NCEP reanalysis data of 15 years (1982-1996). The results show that the seasonal shift of the South Asian High in the upper troposphere and the establishment of the EAM and the INM are closely related to the seasonal warming which causes a reversal of the meridional gradient of upper tropospheric mean temperature over the monsoon regions. On the average of 15 years, the reversal time of the temperature gradient in the EAM region (INM region) is concurrent with (one pentad earlier than) the onset time of the summer monsoon. In most years of the 15-year period, the reversal of temperature gradient coincides or precedes the onset time of the summer monsoon in both the EAM region and the INM region. The results suggest an important role of thermal processes on the establishment of the Asian monsoon. The contributors to the upper tropospheric warming over the EAM region are the strong horizontal warm advection and the diabetic heating against the adiabatic cooling due to upward motion. In the INM region, strong adiabatic heating by subsidence and the diabetic heating are major warming processes against the strong horizontal cold advection related to the persistent northwestlies to the southwestern periphery of the Tibetan Plateau. It appears that the early or late establishment of the Asian summer monsoon is not directly related to the differential warming near the surface.

He, H.↗

Atomistic Modeling of the Interphase Structure between Alumina and Al12Fe2Cr Quasicrystal Approximant

Molecular dynamics simulations are performed to study the stability of an fcc metallic interlayer phase between an iron aluminide quasicrystal approximant (QCA) phase (84.2 Al, 5.7 Fe, 5.4 Ni, and 4.7 at.% Cr) and aluminum oxides (α-Al 2 O 3 and γ-Al 2 O 3 ). The presence of the interlayer phase was experimentally observed in some regions of the inner diameter surface of a TPBAR cladding tube. The simulations employed existing many-body (MEAM and EAM) interatomic potentials to describe atomic interactions in the interlayer and QCA phases, and Buckingham+Coulomb pair interactions for the alumina and between the alumina and the QCA (or fcc). Interface systems of QCA/alumina, fec/alumina, and QCA/fcc are modeled and the fracture energy of the interfaces is calculated. Based on the fracture energies, for the γ-Al 2 O 3 , both MEAM and EAM suggest that the fcc interlayer phase would not form. For the α-Al 2 O 3 , the MEAM predicts that the fcc phase would form, however, the EAM predicts that it forms only if the α-Al 2 O 3 is Al-terminated at the interface. For the O-terminated α-Al 2 O 3 , the EAM predicts the fcc phase would not form, consistent with the results for γ-Al 2 O 3 .

36 MATERIALS SCIENCE↗

Diffusion on Cu surfaces

Understanding surface diffusion is essential in understanding surface phenomena, such as crystal growth, thin film growth, corrosion, physisorption, and chemisorption. Because of its importance, various experimental and theoretical efforts have been directed to understand this phenomena. The Field Ion Microscope (FIM) has been the major experimental tool for studying surface diffusion. FIM have been employed by various research groups to study surface diffusion of adatoms. Because of limitations of the FIM, such studies are only limited to a few surfaces: nickel, platinum, aluminum, iridium, tungsten, and rhodium. From the theoretical standpoint, various atomistic simulations are performed to study surface diffusion. In most of these calculations the Embedded Atom Method (EAM) along with the molecular static (MS) simulation are utilized. The EAM is a semi-empirical approach for modeling the interatomic interactions. The MS simulation is a technique for minimizing the total energy of a system of particles with respect to the positions of its particles. One of the objectives of this work is to develop the EAM functions for Cu and use them in conjunction with the molecular static (MS) simulation to study diffusion of a Cu atom on a perfect as well as stepped Cu(100) surfaces. This will provide a test of the validity of the EAM functions on Cu(100) surface and near the stepped environments. In particular, we construct a terrace-ledge-kink (TLK) model and calculate the migration energies of an atom on a terrace, near a ledge site, near a kink site, and going over a descending step. We have also calculated formation energies of an atom on the bare surface, a vacancy in the surface, a stepped surface, and a stepped-kink surface. Our results are compared with the available experimental and theoretical results.

Karimi, Majid↗

Sensory-Feedback Exoskeletal Arm Controller

An electromechanical exoskeletal arm apparatus has been designed for use in controlling a remote robotic manipulator arm. The apparatus, called a force-feedback exoskeleton arm master (F-EAM) is comfortable to wear and easy to don and doff. It provides control signals from the wearer s arm to a robot arm or a computer simulator (e.g., a virtual-reality system); it also provides force and torque feedback from sensors on the robot arm or from the computer simulator to the wearer s arm. The F-EAM enables the wearer to make the robot arm gently touch objects and finely manipulate them without exerting excessive forces. The F-EAM features a lightweight design in which the motors and gear heads that generate force and torque feedback are made smaller than they ordinarily would be: this is achieved by driving the motors to power levels greater than would ordinarily be used in order to obtain higher torques, and by providing active liquid cooling of the motors to prevent overheating at the high drive levels. The F-EAM (see figure) includes an assembly that resembles a backpack and is worn like a backpack, plus an exoskeletal arm mechanism. The FEAM has five degrees of freedom (DOFs) that correspond to those of the human arm: 1. The first DOF is that of the side-to-side rotation of the upper arm about the shoulder (rotation about axis 1). The reflected torque for this DOF is provided by motor 1 via drum 1 and a planar four-bar linkage. 2. The second DOF is that of the up-and-down rotation of the arm about the shoulder. The reflected torque for this DOF is provided by motor 2 via drum 2. 3. The third DOF is that of twisting of the upper arm about its longitudinal axis. This DOF is implemented in a cable remote-center mechanism (CRCM). The reflected torque for this DOF is provided by motor 3, which drives the upper-arm cuff and the mechanism below it. A bladder inflatable by gas or liquid is placed between the cuff and the wearer s upper arm to compensate for misalignment between the exoskeletal mechanism and the shoulder. 4. The fourth DOF is that of flexion and extension of the elbow. The reflected torque for this DOF is provided by motor 4 and drum 4, which are mounted on a bracket that can slide longitudinally by a pin-and-slot engagement with the upper-arm cuff to compensate for slight variations in the position of the kinematic center of the elbow. Attached to drum 4 is an adapter plate to which is attached a CRCM for the lower arm. 5. The lower-arm CRCM implements the fifth DOF, which is the twist of the forearm about its longitudinal axis. Motor 5 provides the reflected torque for this DOF by driving the lower-arm cuff. A rod transmits twist and torsion between the lower-arm cuff and the hand cuff. With this system, the motion of the wearer s joints and the reflected torques applied to these joints can be measured and controlled in a relatively simple manner. This is because the anthropomorphic design of the mechanism imitates the kinematics of the human arm, eliminating the need for kinematic conversion of joint-torque and joint-angle data.

An, Bin↗

Region and cloud regime dependence of parametric sensitivity in E3SM atmosphere model

Abstract The Department of Energy (DOE)’s Energy Exascale Earth System Model (E3SM), including its atmosphere model (EAM), has many relatively new features. In a previous study we conducted a systematic parametric sensitivity analysis for EAM based on short, perturbed parameter ensemble (PPE) simulations, mainly focusing on global mean climate features and metrics. While parameter values in global climate models are generally invariant in space and time, model response to parameters perturbation may vary by regions and climate regimes, which motivates the need to better understand the EAM model behaviors and physics at regional scale and process level. In this study, using the same set of PPE simulations and a similar sensitivity analysis framework, we identify parameters that cause largest sensitivities over different regions and compare model responses in fast atmospheric processes to the parameters across different cloud regimes for several important cloud-related fidelity metrics. We find that cloud forcing has opposite response to some parameters over mid-latitude vs. tropical land. We also analyze how the parametric sensitivity varies as stratocumulus transitions to shallow convection and to deep convection over ocean. Low cloud forcing and shortwave cloud forcing in the subtropical eastern Pacific are most sensitive to the parameters controlling the width of the probability density function (PDF) of the subgrid vertical velocity ( w’ ) ( gamma ) and the damping of the w’ skewness ( c8 ) near the coast but become more sensitive to the parameter affecting the damping of the w’ variance ( c1 ) further offshore. Detailed interpretation of the spatial dependence of parametric sensitivity is provided. We also investigate how the parametric sensitivity evolves with prediction duration. This study improves our process-level understanding of cloud physics and parameterization and provides insights for developing more advanced regime-aware parameterization schemes in global climate model.

54 ENVIRONMENTAL SCIENCES↗

Atmospheric excitation of nonseasonal polar motion

Analysis of nonseasonal polar motion excitation and atmospheric mass equatorial angular momentum (EAM) over land for the period 1980-1989 reveals a clear pattern of high power and correlation during the northern hemisphere (NH) winter followed by low power and correlation during the NH summer. A special case of this pattern occurs for longer than 14 months (from January 1987 to March 1988) when the correlation throughout the NH summer remains statistically significant. During this epoch an average of 72% of the nonseasonal polar motion excitation power at frequencies between -30 and +12 cycles/yr linearly related to atmospheric EAM over land. During the southern hemisphere winter there is significant correlation between the atmospheric EAM over midlatitude southern oceans and polar motion excitation indicating the existence of a dynamic atmosphere-ocean excitation. The atmospheric excitation power is too small to explain the large correlation during the NH winter. The effects of winds probably account for the deficit in power. The implication of these results is that there are two main excitation sources each dominant at different seasons. Atmospheric mass redistribution over land forces polar motion during the NH winter, and a dynamic atmpshere-ocean response is important during the SH winter.

Kuehne, John↗

Unlocking naphtha from polyolefins using Ni-based hydrocracking catalysts

Naphtha (C 5 -C 12 alkanes) is a platform feedstock in the petroleum industry and an ideal target for reutilizing plastic waste at scale. Hydrocracking has emerged as a promising technology for deconstructing polyolefins to naphtha. However, contemporary catalysts rely on Pt to achieve high rates, while earth-abundant metals (EAM) perform poorly. Herein, we develop high-performance Ni/BEA catalysts for polyolefin hydrocracking, achieving complete low-density polyethylene (LDPE) deconstruction with 80 % maximum naphtha yield within 12 h at 250 °C, 60 bar H 2 , and a catalyst-to-polymer ratio of 1:100. These catalysts are versatile, accommodating virgin resins and commercial polymer products, and are directly reusable and regenerable. Comparative analysis highlights that 5Ni/BEA(25) exhibits the highest naphtha productivity of 11.4 $\huge{(}$$\frac{g_{Naphtha}}{g_{cat}{•h}}$$\huge{)}$, surpassing previously reported Pt- and EAM-based catalysts by 2.5–6.2x. Technoeconomic and life-cycle analyses reveal naphtha from LDPE hydrocracking is economically and environmentally competitive to fossil-fuel-derived naphtha. In conclusion, these advancements pave the way for the industrial implementation of earth-abundant Ni-based catalysts in polyolefin hydrocracking, followed by retrofitting steam crackers to address plastics waste at scale.

Circular Economy↗

A liquid drop model for embedded atom method cluster energies

Minimum energy configurations for homonuclear clusters containing from two to twenty-two atoms of six metals, Ag, Au, Cu, Ni, Pd, and Pt have been calculated using the Embedded Atom Method (EAM). The average energy per atom as a function of cluster size has been fit to a liquid drop model, giving estimates of the surface and curvature energies. The liquid drop model gives a good representation of the relationship between average energy and cluster size. As a test the resulting surface energies are compared to EAM surface energy calculations for various low-index crystal faces with reasonable agreement.

Finley, C. W.↗

Characterizing diurnal and seasonal cycles in monsoon systems from TRMM and CEOP observations

The CEOP Inter-Monsoon Study (CIMS) is one of the two main science drivers of CEOP that aims to (a) provide better understanding of fundamental physical processes in monsoon regions around the world, and (b) demonstrate the synergy and utility of CEOP data in providing a pathway for model physics evaluation and improvement. As the data collection phase for EOP-3 and EOP-4 is being completed, two full annual cycles (2003-2004) of research-quality data sets from satellites, reference sites, and model output location time series (MOLTS) have been processed and made available for data analyses and model validation studies. This article presents preliminary results of a CIMS study aimed at the characterization and intercomparison of all major monsoon systems. The CEOP reference site data proved its value in such exercises by being a powerful tool to cross-validate the TRMM data, and to intercompare with multi-model results in ongoing work. We use 6 years (1998-2003) of pentad CEOP/TRMM data with 2deg x 2.5deg latitude-longitude grid, over the domain of interests to define the monsoon climatological diurnal and annual cycles for the East Asian Monsoon (EAM), the South Asian Monsoon (SAM), the West Africa Monsoon (WAM), the North America/Mexican Monsoon (NAM), the South American Summer Monsoon (SASM) and the Australian Monsoon (AUM). As noted, the TRMM data used in the study were cross-validated using CEOP reference site data, where applicable. Results show that the observed diurnal cycle of rain peaked around late afternoon over monsoon land, and early morning over the oceans. The diurnal cycles in models tend to peak 2-3 hours earlier than observed. The seasonal cycles of the EAM and SAM show the strongest continentality, i.e, strong control by continental processes away from the ITCZ. The WAM, and the AUM shows the less continentality, i.e, strong control by the oceanic ITCZ.

Lau, William K. M.↗

Characterizing Diurnal and Seasonal Cycles in Monsoon Systems from TRMM and CEOP Observations

The CEOP Inter-Monsoon Study (CIMS) is one of the two main science drivers of CEOP that aims to (a) provide better understanding of fundamental physical processes in monsoon regions around the world, and (b) demonstrate the synergy and utility of CEOP data in providing a pathway for model physics evaluation and improvement. As the data collection phase for EOP-3 and EOP-4 is being completed, two full annual cycles (2003-2004) of research-quality data sets from satellites, reference sites, and model output location time series (MOLTS) have been processed and made available for data analyses and model validation studies. This article presents preliminary results of a CIMS study aimed at the characterization and intercomparison of all major monsoon systems. The CEOP reference site data proved its value in such exercises by being a powerful tool to cross-validate the TRMM data, and to intercompare with multi-model results in ongoing work. We use 6 years (1998-2003) of pentad CEOP/TRMM data with 2 deg x 2.5 deg. latitude-longitude grid, over the domain of interests to define the monsoon climatological diurnal and annual cycles for the East Asian Monsoon (EAM), the South Asian Monsoon (SAM), the West Africa Monsoon (WAM), the North America/Mexican Monsoon (NAM), the South American Summer Monsoon (SASM) and the Australian Monsoon (AUM). As noted, the TRMM data used in the study were cross-validated using CEOP reference site data, where applicable. Results show that the observed diurnal cycle of rain peaked around late afternoon over monsoon land, and early morning over the oceans. The diurnal cycles in models tend to peak 2-3 hours earlier than observed. The seasonal cycles of the EAM and SAM show the strongest continentality, i.e, strong control by continental processes away from the ITCZ. The WAM, and the AUM shows the less continentality, i.e, strong control by the oceanic ITCZ.

Lau, William K. M.↗

Studying performance portability of LAMMPS across diverse GPU-based platforms

The molecular dynamics simulation software, LAMMPS, utilizes the Kokkos acceleration library to port computation to a diverse set of architectures including those based on GPU accelerators. In addition to Kokkos, LAMMPS contains a vast code base that leverages the CUDA application programming interface using library functions such as cuFFT, CUDA's fast-fourier transform (FFT) library, and, more recently, also support for AMD's Heterogeneous Interface for Portability (HIP) that is rapidly growing. While preparing LAMMPS tests for the AMD GPU-based test system precursors to Frontier, we investigated several strategies for accelerating LAMMPS on AMD GPUs, using the AMD Instinct MI100 and MI250X. In this work, we integrated the HIP FFT library, hipFFT, into the particle-particle particle-mesh (PPPM) long-range solver, which allowed the porting of PPPM calculations to the GPUs. Kokkos behavior on the MI100 and MI250X was also investigated through the package kokkos command of LAMMPS, targeting communication, memory usage, and particle grid decomposition. The Tersoff, Reax, Lennard-Jones (LJ), EAM, Granular, and PPPM potentials were investigated in this effort, and results from these experiments are provided. In conclusion, the selected potentials were run on Spock (AMD Instinct MI100), Crusher (AMD Instinct MI250X), AFW HPC11 (NVIDIA A100) and Summit (NVIDIA V100), for comparison. Operational roofline models were constructed and analyzed for the Tersoff, Reax, and Lennard–Jones potentials on Crusher and Summit.

97 MATHEMATICS AND COMPUTING↗

Evolution of dislocation loops in irradiated α-Uranium: An atomistically-informed cluster dynamics investigation

An atomistically-informed mean field cluster dynamics model has been developed to investigate the nucleation and growth of prismatic loops in irradiated a-Uranium. TEM analysis of neutron irradiated a-Uranium shows the evolution of self-interstitial atom and vacancy loops on (010) and (100) crystallographic planes, respectively, resulting in an anisotropic lattice swelling of its face-centered orthorhombic crystal. To provide model parameters, the crystallography of loops and the binding energy of point defects to these loops were studied using an angular dependent EAM potential and classical molecular dynamics (MD) simulations. Furthermore, using the bond-boost hyperdynamics method, the anisotropic diffusion of interstitials and vacancies in a-Uranium was investigated. Here, the mechanisms of point defect diffusion and the associated migration energies were reported and compared with previous DFT studies. The energetics and kinetic quantities mentioned above were used in the cluster dynamics model to predict the cluster density as a function of dose rate, dose and temperature and the results were compared to the reported neutron irradiation experiments. The model predictions reveal an accumulation of small sized vacancy loops along with a population of large and growing self-interstitial loops, which closely corresponds to the TEM observations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Probing multi-dimensional composition spaces in search of strong metallic alloys

Refractory complex concentrated alloys (RCCA) offer exceptionally high-temperature strength compared to pure metals and dilute alloys, but predictive theory for RCCA design is lacking. We present large-scale molecular Dynamics (MD) simulations of crystal plasticity to explore alloy compositions for maximum mechanical strength, focusing on Fe-Ta-W and Nb-Ta-Mo-W alloy families modeled with Embedded Atom Model (EAM) and Spectral Neighbor Analysis Potentials (SNAP). To efficiently guide the search for strong alloy compositions, we employ iterative optimization using Gaussian process regression. Many simulated RCCA compositions exhibit pronounced cocktail strengthening, with strengths surpassing their strongest constituent metal, tungsten. Contrary to expectations, the highest strength is found on binary edges of the RCCA composition space. Detailed analyses of atomistic simulations reveal that, similar to pure BCC metals, plastic response in RCCA is primarily governed by screw dislocations. However, at large strains, dislocation multiplication and interactions (Taylor hardening) become the dominant mechanisms contributing to RCCA strength.

Materials science↗

A review of displacement cascade simulations using molecular dynamics emphasizing interatomic potentials for TPBAR components

This review explores molecular dynamics simulations for studying radiation damage in Tritium Producing Burnable Absorber Rod (TPBAR) materials, emphasizing the role of interatomic potentials in displacement cascades. Recent machine learning potentials (MLPs), trained on quantum data, enhance prediction accuracy over traditional models like EAM. We highlight temperature, PKA energy, and composition effects on damage evolution in TPBAR components, recommending suitable potentials and discussing advancements for materials in extreme radiation environments.

36 MATERIALS SCIENCE↗

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗