Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomistic simulation”

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 541 records · Page 30

Multi-fidelity machine-learning with uncertainty quantification and Bayesian optimization for materials design: Application to ternary random alloys

We present a scale-bridging approach based on a multi-fidelity (MF) machine-learning (ML) framework leveraging Gaussian processes (GP) to fuse atomistic computational model predictions across multiple levels of fidelity. Through the posterior variance of the MFGP, our framework naturally enables uncertainty quantification, providing estimates of confidence in the predictions. We used density functional theory as high-fidelity prediction, while a ML interatomic potential is used as low-fidelity prediction. Practical materials’ design efficiency is demonstrated by reproducing the ternary composition dependence of a quantity of interest (bulk modulus) across the full aluminum–niobium–titanium ternary random alloy composition space. The MFGP is then coupled to a Bayesian optimization procedure, and the computational efficiency of this approach is demonstrated by performing an on-the-fly search for the global optimum of bulk modulus in the ternary composition space. The framework presented in this manuscript is the first application of MFGP to atomistic materials simulations fusing predictions between density functional theory and classical interatomic potential calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A density-functional theory study of the Al/AlOx/Al tunnel junction

The aluminum oxide tunnel junction is a key component of the majority of superconducting quantum devices. For high-quality, reproducible, and scalably manufacturable qubits, the ability to fabricate Josephson junctions (JJs) with a targeted critical current and high uniformity is essential. In this study, we use first-principles modeling to assess fundamental aspects of the atomic structure of both amorphous and crystalline aluminum oxide tunnel junctions and relate the structure to predicted performance metrics. We use modified ab initio molecular dynamics to develop realistic models of the tunnel junction, from which interface roughness and local thickness fluctuations are analyzed in an unbiased manner by training a neural network to identify the boundary between metal and oxide. We show that the effective thickness of the insulating part of the junction can be different from the apparent physical thickness. We calculate the rate of Cooper pair tunneling for the atomically resolved electrostatic potential using direct numerical solution in 3D, which shows a channeling effect that impacts the junction critical current. The predicted critical current is a useful JJ design parameter that can be accessed from the ab initio calculations without fitting parameters. To assess the limits of uniformity and fabrication choices (e.g., oxidation vs epitaxy), we compare the amorphous junctions to crystalline models, which show order of magnitude more efficient tunneling compared to the amorphous case, underlining the connection between atomistic structure and Cooper pair tunneling efficiency. Further, this work provides a foundation for ab initio materials design and evaluation to help accelerate future development of improved tunnel junctions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Geometry optimization speedup through a geodesic approach to internal coordinates

We present a new geodesic-based method for geometry optimization in a basis set of redundant internal coordinates. Overall, our method updates the molecular geometry by following the geodesic generated by a displacement vector on the internal coordinate manifold, which dramatically reduces the number of steps required to converge to a minimum. Our method can be implemented in any existing optimization code, requiring only implementation of derivatives of the Wilson B-matrix and the ability to numerically solve an ordinary differential equation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Violation of the Cauchy–Born rule in multi-principal element alloys

Multi-principal element alloys are a novel class of materials that are formed by combining multiple elements in high concentrations and show exceptional properties compared to conventional alloys. These alloys have high configurational entropy due to inherent atomic disorder. The Cauchy–Born rule is a popular homogenization method for linking atomistic to continuum length scales. In this Letter, we use ab initio density functional theory calculations to report that the Cauchy–Born rule, which holds in conventional alloys in the absence of defects, is not valid in multi-principal element alloys. Further, the violation of the Cauchy–Born rule in these alloys is attributed to the presence of atomic disorder due to which the deformations are inhomogeneous. Our results also show the deviations in deformation are related to the magnitude of the stretch and shear by power laws.

36 MATERIALS SCIENCE↗

RANGE: A robust adaptive nature-inspired global explorer of potential energy surfaces

With the growing demand for realistic representations of chemical structures and the advent of exascale computing, the intelligent sampling of potential energy surfaces and efficient identification of global minima have become more essential but also more feasible. Building on prior studies demonstrating the efficiency of the Artificial Bee Colony (ABC) swarm intelligence algorithm, we report a hybrid metaheuristic framework that integrates the adaptive exploration capabilities of ABC coupled with the exploitation strengths of genetic algorithms (GA) in a scalable, Python-based implementation. The resulting tool, RANGE (Robust Adaptive Nature-inspired Global Explorer), provides seamless interfaces to multiple potential energy evaluators, either directly or via widely used Python libraries, and is designed for high-performance computing environments. We describe the implementation details of RANGE and evaluate its performance, relative to ABC- or GA-alone based algorithms, on a variety of chemical systems, including molecular clusters and heterogeneous surfaces. In conclusion, our results demonstrate RANGE’s efficiency, robustness, and broad applicability in addressing challenging global optimization problems in computational chemistry and materials science.

Algorithms and data structure↗

Magnus force induced magnetic diode effect in skyrmion systems

We show that skyrmions can exhibit a “magnetic diode effect,” where there is a nonreciprocal response in the transport when the magnetic field is reversed. This effect can be achieved for skyrmions moving in channels with a sawtooth potential on one side and a reversed sawtooth potential on the other side. We consider the cases of both spin-transfer torque (STT) and spin–orbit torque (SOT). When the magnetic field is held fixed, the velocity response of the skyrmion is the same for current applied in either direction for both torques, so there is no current diode effect. When the magnetic field is reversed, under STT driving the velocity of the skyrmion reverses and its absolute value changes. Under SOT driving, the velocity remains in the same direction but drops to a much lower value, resulting in negative differential conductivity. For a fixed current, we find a nonreciprocal skyrmion velocity as a function of the applied field’s sign, in analogy to the velocity–current curves observed in the usual diode effect. The nonreciprocity is generated by the Magnus force, which causes skyrmions to interact preferentially with one side of the channel. Since the channel sides have opposite asymmetry, a positive magnetic field can cause the skyrmion to interact with the “hard” asymmetry side of the channel, while a negative magnetic field causes the skyrmion to interact with the “easy” asymmetry side. This geometry could be used to create new kinds of magnetic-field-induced diode effects that can be harnessed in new types of skyrmion-based devices.

36 MATERIALS SCIENCE↗

Maximizing efficiency of dataset compression for machine learning potentials with information theory

Machine learning interatomic potentials (MLIPs) balance high accuracy and lower costs compared to density functional theory calculations, but their performance often depends on the size and diversity of training datasets. Large datasets improve model accuracy and generalization but are computationally expensive to produce and train on, while smaller datasets risk discarding rare but important atomic environments and compromising MLIP accuracy/reliability. Here, we develop an information-theoretical framework to quantify the efficiency of dataset compression methods and propose an algorithm that maximizes this efficiency. By framing atomistic dataset compression as an instance of the minimum set cover (MSC) problem over atom-centered environments, our method identifies the smallest subset of structures that contains as much information as possible from the original dataset while pruning redundant information. The approach is extensively demonstrated on the GAP-20 and TM23 datasets and validated on 64 varied datasets from the ColabFit repository. Across all cases, MSC consistently retains outliers, preserves dataset diversity, and reproduces the long-tail distributions of forces even at high compression rates, outperforming other subsampling methods. Furthermore, MLIPs trained on MSC-compressed datasets exhibit reduced error for out-of-distribution data even in low-data regimes. We explain these results using an outlier analysis and show that such quantitative conclusions could not be achieved with conventional dimensionality reduction methods. The algorithm is implemented in the open-source QUESTS package and can be used for several tasks in atomistic modeling, from data subsampling, outlier detection, and training improved MLIPs at a lower cost.

36 MATERIALS SCIENCE↗

Data efficiency and extrapolation trends in neural network interatomic potentials

Abstract Recently, key architectural advances have been proposed for neural network interatomic potentials (NNIPs), such as incorporating message-passing networks, equivariance, or many-body expansion terms. Although modern NNIP models exhibit small differences in test accuracy, this metric is still considered the main target when developing new NNIP architectures. In this work, we show how architectural and optimization choices influence the generalization of NNIPs, revealing trends in molecular dynamics (MD) stability, data efficiency, and loss landscapes. Using the 3BPA dataset, we uncover trends in NNIP errors and robustness to noise, showing these metrics are insufficient to predict MD stability in the high-accuracy regime. With a large-scale study on NequIP, MACE, and their optimizers, we show that our metric of loss entropy predicts out-of-distribution error and data efficiency despite being computed only on the training set. This work provides a deep learning justification for probing extrapolation and can inform the development of next-generation NNIPs.

36 MATERIALS SCIENCE↗

Topological energy barrier for skyrmion lattice formation in MnSi

We report the direct measurement of the topological skyrmion energy barrier through a hysteresis of the skyrmion lattice in the chiral magnet MnSi. Measurements were made using small-angle neutron scattering with a custom-built resistive coil to allow for high-precision minor hysteresis loops. The experimental data were analyzed using an adapted Preisach model to quantify the energy barrier for skyrmion formation and corroborated by the minimum-energy path analysis based on atomistic spin simulations. Here, we reveal that the skyrmion lattice in MnSi forms from the conical phase progressively in small domains, each of which consisting of hundreds of skyrmions, and with an activation barrier of several eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Triple junction solute segregation in Al-based polycrystals

Solute segregation is a crucial means of stabilizing nanostructured alloys, and at very small grain sizes, this requires consideration of triple junctions (TJs), which attain a meaningful volume fraction and thus become relevant for bulk material behavior. Here, in this study, the solute segregation spectra for grain boundary (GB) and TJ sites are calculated for a large number of dilute Al-based binary systems with available interatomic potentials. A defect-identifying algorithm is applied to quantify the average GB thickness and classify the intergranular site spectra into GB and junction subspectra. The algorithm is also applied to a hybrid electronic-level database for GB segregation of various solutes in Al, yielding polycrystalline TJ solute segregation spectra from first principles. The results suggest that TJ segregation is alloy or interatomic potential dependent and can exhibit either boundary or junction preference. With these spectra as inputs, the spectral GB segregation model gives quantitative predictions of segregation as a function of grain size, temperature, and total solute concentration, suitable for alloy screening and design.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

OpenFerro v0.1.0

OpenFerro is a Python package for on-lattice atomistic dynamics simulation of ferroic materials. OpenFerro is based on JAX, a high-performance linear algebra package supporting auto-differentiation and GPU acceleration. OpenFerro is designed to minimize the effort required to build on-lattice Hamiltonian models, and to perform molecular dynamics (MD) and Landau-Lifshitz-Gilbert simulations. Unlike existing codes, OpenFerro provides a unified interface to model different types of local order parameters.

Xie, Pinchen [Lawrence Berkeley National Laborator↗

Thermodynamics of grain-boundary phases

The past decade has seen a significant increase in research efforts aimed at understanding the thermodynamics of low-dimensional phases existing in many materials systems, ranging from two-dimensional materials to core regions of extended defects in crystalline solids. We review the current status of theoretical, computational, and experimental research on the “defect phases,” focusing on grain boundaries (GBs) in elemental and multicomponent polycrystalline materials. After reviewing the generalized concept of a phase of any dimensionality, we discuss recent progress in atomistic computer simulations of GB phase transformations and phase coexistences, including the observation of one-dimensional defects separating GB phases (defects in defects). Computational predictions compare well with experimental observations of multiple GB phases and segregation-induced phase transformations. An intriguing open question of GB thermodynamics is whether the GB free energy can be driven to a zero value by increasing solute segregation. We review recent efforts to understand this ultimate thermodynamic stabilization of GB phases and the possible polycrystalline microstructures that may arise. An outlook for future research in the field is discussed.

Materials science↗

Multi-Scale Modeling Framework for Mercury Biogeochemistry

Multi-Scale modeling of mercury (Hg) geochemical speciation and reactions has been performed by integrating atomistic quantum chemical calculations with continuum scale speciation models. Major progress has been made in the improvement of quantum chemical models to calculate critical thermodynamic data for Hg complexes in aquatic environments. Rapid and reliable quantum chemical approaches have been developed for calculating acid dissociation constants (pK a ) and stability constants (log K), with calculated mean unsigned errors of 0.5 and 1.5 log units, respectively for ligand molecules and Hg complexes. At the continuum scale, systematic analysis of uncertainty propagation in mercury (Hg) speciation modeling has been conducted and was used to identify environmental conditions under which thermodynamic constant uncertainties are significant and recommended to be accounted for. The integrated framework for multi-scale modeling of mercury geochemistry is open to the research community through the web-based multiscale modeling aqueous speciation resource, AQUA-MER. The improved quantum chemical approaches for thermodynamic constant calculations are accessible through AQUA-MER and can be used to provide the missing constants in the continuum scale speciation calculations. In addition to low molecular mass Hg complex speciation, modeling natural aquatic environments also involve the transport of high molecular weight dissolved organic matter (DOM) in reactive flows simultaneously with equilibrium and kinetic reactions. To this end, atomistic MD simulations were performed to capture the details of aggregation, mechanisms and distribution of functional groups in DOM at the molecular level. The elemental composition and calculated bulk properties of the DOM models are in close agreement with experimental measurements. A travel-time based reactive transport model in the hyporheic zone of stream corridors was established for the multicomponent Hg-DOM-S system and implemented through PFLOTRAN.

54 ENVIRONMENTAL SCIENCES↗

Nanoscale Bio-engineering Solutions for Space Exploration: The Nanopore Sequencer

Characterization of biological systems at the molecular level and extraction of essential information for nano-engineering design to guide the nano-fabrication of solid-state sensors and molecular identification devices is a computational challenge. The alpha hemolysin protein ion channel is used as a model system for structural analysis of nucleic acids like DNA. Applied voltage draws a DNA strand and surrounding ionic solution through the biological nanopore. The subunits in the DNA strand block ion flow by differing amounts. Atomistic scale simulations are employed using NASA supercomputers to study DNA translocation, with the aim to enhance single DNA subunit identification. Compared to protein channels, solid-state nanopores offer a better temporal control of the translocation of DNA and the possibility to easily tune its chemistry to increase the signal resolution. Potential applications for NASA missions, besides real-time genome sequencing include astronaut health, life detection and decoding of various genomes.

Stolc, Viktor↗

Nanoscale Bioengineering Solutions for Space Exploration the Nanopore Sequencer

Characterization of biological systems at the molecular level and extraction of essential information for nano-engineering design to guide the nano-fabrication of solid-state sensors and molecular identification devices is a computational challenge. The alpha hemolysin protein ion channel is used as a model system for structural analysis of nucleic acids like DNA. Applied voltage draws a DNA strand and surrounding ionic solution through the biological nanopore. The subunits in the DNA strand block ion flow by differing amounts. Atomistic scale simulations are employed using NASA supercomputers to study DNA translocation. with the aim to enhance single DNA subunit identification. Compared to protein channels, solid-state nanopores offer a better temporal control of the translocation of DNA and the possibility to easily tune its chemistry to increase the signal resolution. Potential applications for NASA missions, besides real-time genome sequencing include astronaut health, life detection and decoding of various genomes. http://phenomrph.arc.nasa.gov/index.php

Ioana, Cozmuta↗

Interfacial Atomistic Mechanisms of Lithium Metal Stripping and Plating in Solid–State Batteries

All-solid-state batteries based on a Li metal anode represent a promising next-generation energy storage system, but are currently limited by low current density and short cycle life. Further research to improve the Li metal anode is impeded by the lack of understanding in its failure mechanisms at lithium–solid interfaces, in particular, the fundamental atomistic processes responsible for interface failure. Here, using large-scale molecular dynamics simulations, the first atomistic modeling study of lithium stripping and plating on a solid electrolyte is performed by explicitly considering key fundamental atomistic processes and interface atomistic structures. In the simulations, the interface failure initiated with the formation of nano-sized pores, and how interface structures, lithium diffusion, adhesion energy, and applied pressure affect interface failure during Li cycling are observed. By systematically varying the parameters of solid-state lithium cells in the simulations, the parameter space of applied pressures and interfacial adhesion energies that inhibit interface failure during cycling are mapped to guide selection of solid-state cells. Furthermore, this study establishes the atomistic modeling for Li stripping and plating, and predicts optimal solid interfaces and new strategies for the future research and development of solid-state Li-metal batteries.

25 ENERGY STORAGE↗