Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “force field”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Comparing Classical and Machine Learning Force Fields for Modeling Deformation of Metal–Organic Frameworks Relevant for Direct Air Capture

Deformation of metal–organic frameworks (MOFs) induced by adsorbate molecules can affect adsorption properties such as capacity and selectivity, but most computational studies of MOFs assume framework rigidity to simplify calculations. Although flexible force fields (FFs) for MOFs have been parametrized for specific materials, the generality of FFs for reliably modeling adsorbate-induced deformation to accuracy nearing that of density functional theory (DFT) has not been established. This work confirms using DFT calculations that adsorbate-induced deformation can affect CO 2 and H 2 O adsorption energies in a considerable fraction of MOFs promising for direct air capture (DAC). We then benchmark the efficacy of several general-purpose FFs in describing adsorbate-induced deformation for DAC against DFT. Our results show that current classical FFs are insufficient for describing MOF deformation, especially in cases of interest for DAC where strong interactions exist between adsorbed molecules and MOF frameworks. Some emerging machine learning force fields (MLFFs) we tested, particularly CHGNet, MACE-MP-0, and Equiformer V2, appear to be more promising than the classical FF for emulating the deformation behavior described by DFT. The best performing FF (CHGNet), however, fails to achieve the accuracy required for practical predictions with a mean absolute adsorption energy error of 0.124 eV.

adsorption↗
DOE OSTI2026

Unsupervised atomic data mining via multi-kernel graph autoencoders for machine learning force fields

Constructing a chemically diverse dataset while avoiding sampling bias is critical to training efficient and generalizable force fields. However, in computational chemistry and materials science, many common dataset generation techniques are prone to oversampling regions of the potential energy surface. Furthermore, these regions can be difficult to identify and isolate from each other or may not align well with human intuition, making it challenging to systematically remove bias in the dataset. While traditional clustering and pruning (down-sampling) approaches can be useful for this, they can often lead to information loss or a failure to properly identify distinct regions of the potential energy surface due to difficulties associated with the high dimensionality of atomic descriptors. In this work, we introduce the Multi-kernel Edge Attention-based Graph Autoencoder (MEAGraph) model, an unsupervised approach for analyzing atomic datasets. MEAGraph combines multiple linear kernel transformations with attention-based message passing to capture geometric sensitivity and enable effective dataset pruning without relying on labels or extensive training. Demonstrated applications on niobium, tantalum, and iron datasets show that MEAGraph efficiently groups similar atomic environments, allowing for the use of basic pruning techniques for removing sampling bias. This approach provides an effective method for representation learning and clustering that can be used for data analysis, outlier detection, and dataset optimization.

Materials science↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A Transferable Force Field for Simulating Adsorption in Metal–Organic Frameworks with Open Metal Sites Based on the 12–6–4 Lennard-Jones Potential

Metal−organic frameworks (MOFs) that contain coordinatively unsaturated open metal sites (OMSs) provide strong host− guest interactions, making them promising sorbents for low-concentration gas adsorption applications such as direct air capture and atmospheric water harvesting. However, accurately modeling host−guest interactions involving OMSs remains challenging for classical force fields (FFs) based on the 12−6 Lennard−Jones (LJ) potential, as the polarization effect of the guest molecule induced by the positively charged OMS is not considered. Here, we introduce an FF based on the 12−6−4 LJ potential, which incorporates charge−induced dipole interactions and is parametrized against a diverse set of host−guest potential energy surfaces (PESs) obtained from density functional theory (DFT). The resulting FF, trained on a generic trimetallic cluster, performs well in both host−guest binding energetics and gas adsorption isotherms across different OMS-containing MOFs, including MOF-74 series and Cu-BTC. These results highlight the excellent transferability of our approach and its potential to enhance the accuracy and robustness of high-throughput MOF discovery workflows, particularly for gas adsorption and separation in large and diverse MOF databases.

36 MATERIALS SCIENCE↗
DOE OSTI2024

ReaxFF Reactive Force Field for Exploring Electronically Switchable Polarization in Zn 1– x Mg x O Ferroelectric Semiconductors

Cation misfit in traditional ferroelectric crystals offers a new material platform that can drive electronic components toward structural miniaturization and high-density integration, enabling deviation from von-Neumann architectures. Here, we explore ferroelectricity in Zn 1–x Mg x O, a nontraditional ferroelectric material with tunable properties. Using data from density-functional theory calculations, we have developed a ReaxFF reactive force field to explore the ferroelectric properties of Zn 1–x Mg x O and reveal the hysteresis behavior. We discover that ferroelectric switching can be observed at a critical thickness of 10 nm with a residual polarization of ~100 μC/cm 2 . Our analysis indicates that an increase in Mg-substitution correlates with a decrease in the coercive field. We also observe a strong temperature dependence of the coercive field in Zn 1–x Mg x O, with values decreasing as the temperature increases. Additionally, we find that the distribution of Mg atoms significantly impacts the coercive field, with a clustered distribution leading to a substantial increase. In particular, a decrease in coercive field values is observed when Mg atoms are randomly distributed, compared to uniform distribution. Furthermore, leveraging tunable hysteresis behavior offered by varied percentages and distribution of Mg-substitution provides valuable insights into the design of next-generation functional devices and will inspire further investigations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Investigation of the Effect of Framework Flexibility on Adsorption in SIFSIX-3-Cu using a Machine-Learned Force Field

Metal-organic frameworks (MOFs) are a promising class of adsorbents. The performance of MOF sorbents relies on high selectivity and low regeneration energy. This work focuses on the use of machine learned force fields (MLFFs) to model adsorption in a flexible MOF, SIFSIX-3-Cu. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies, forces, and stresses, using an iterative sampling scheme combining sampling based on molecular dynamics, Monte Carlo, and geometry optimization to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included in this model. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations using the MLFF predicted adsorption isotherms in good agreement with experimental data for a range of pressures (40 Pa – 104 Pa) in contrast to rigid models, which overpredict CO2 adsorption at low pressures. The improvement was the result of a description of the variability of fluorine-fluorine diagonal distances at adsorption sites. This detailed description of flexibility afforded by the MLFF resulted in more accurate predictions adsorption isotherms when compared to the experimentally measured values. These results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low pressure applications and provide a robust procedure for training MLFF models for MOFs.

Atomistic Simulation↗
DOE OSTI2024

phosaa14SB and phosaa19SB: Updated Amber Force Field Parameters for Phosphorylated Amino Acids

Phosphorylated amino acids are involved in many cell regulatory networks; proteins containing these post-translational modifications are widely studied both experimentally and computationally. Simulations are used to investigate a wide range of structural and dynamic properties of biomolecules, such as ligand binding, enzyme-reaction mechanisms, and protein folding. However, the development of force field parameters for the simulation of proteins containing phosphorylated amino acids using the Amber program has not kept pace with the development of parameters for standard amino acids, and it is challenging to model these modified amino acids with accuracy comparable to proteins containing only standard amino acids. In particular, the popular ff14SB and ff19SB models do not contain parameters for phosphorylated amino acids. Here, the dihedral parameters for the side chains of the most common phosphorylated amino acids are trained against reference data from QM calculations adopting the ff14SB approach, followed by validation against experimental data. Finally, library files and corresponding parameter files are provided, with versions that are compatible with both ff14SB and ff19SB.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Machine learning force field model for kinetic Monte Carlo simulations of itinerant Ising magnets

Here, we present a scalable machine learning (ML) framework for large-scale kinetic Monte Carlo (kMC) simulations of itinerant electron Ising systems. As the effective interactions between Ising spins in such itinerant magnets are mediated by conducting electrons, the calculation of energy change due to a local spin update requires solving an electronic structure problem. Such repeated electronic structure calculations could be overwhelmingly prohibitive for large systems. Assuming the locality principle, a convolutional neural network (CNN) model is developed to directly predict the effective local field and the corresponding energy change associated with a given spin update based on Ising configuration in a finite neighborhood. As the kernel size of the CNN is fixed at a constant, the model can be directly scalable to kMC simulations of large lattices. Our approach is reminiscent of the ML force field models widely used in first-principles molecular dynamics simulations. Applying our ML framework to a square-lattice double-exchange Ising model, we uncover unusual coarsening of ferromagnetic domains at low temperatures. Our work highlights the potential of ML methods for large-scale modeling of similar itinerant systems with discrete dynamical variables.

machine learning↗
DOE OSTI2024

Machine Learned Force Field Modeling of Metal Organic Frameworks for CO2 Direct Air Capture

Metal organic frameworks (MOFs) are a large class of porous materials and have garnered significant interest due to their large surface areas and their tunable physical and chemical properties. Numerous prior studies have been performed to screen large databases of this material class for promising DAC sorbent materials. These studies have often relied on classical model potentials. While density functional theory (DFT) calculations have been shown to be very accurate for modeling the interaction of CO2 with MOFs, such calculations are too computationally demanding for statistically significant adsorption predictions. To overcome this barrier, we developed methods for training models to achieve DFT-level accuracy for the forces and energies associated with MOF flexibility and CO2 adsorption using machine learned force fields (MLFFs). These methods were parametrized based on DFT calculations of CO2 in a flexible MOF and used to predict MOF structural properties as well as CO2 adsorption in several MOFs.

Findley, John↗
DOE OSTI2025

Molecular Dynamics Simulations of Supercritical Carbon Dioxide and Water using TraPPE and SWM4-NDP Force Fields

The increased levels of carbon dioxide (CO 2 ) emissions due to the combustion of fossil fuels and the consequential impact on global climate change have made CO 2 capture, storage, and utilization a significant area of focus for current research. In most electrochemical CO 2 applications, water is used as a proton donor due to its high availability and mobility and use as a polar solvent. Additionally, supercritical CO 2 is a promising avenue for electrochemical applications due to its unique chemical and physical properties. Consequently, understanding the interactions between water and supercritical CO 2 is of great importance for future electrochemical applications. Molecular dynamics (MD) simulation is a powerful tool that enables atomistic-resolution dynamics of molecular systems, which can complement and guide future experimental investigations. This study employed atomistic MD to study the cosolubilities, codiffusivities, and structure of supercritical CO 2 and water systems, with a polarizable water model (SWM4-NDP) and a nonpolarizable CO 2 model (TraPPE). Additionally, ab initio MD simulations were used to better understand how atomistic polarizable/nonpolarizable models compare to explicit modeling of electron densities. The polarizable water model exhibited substantial improvement in water-associated properties. In conclusion, we anticipate the development of a compatible polarizable CO 2 model to yield similar improvement, providing a pathway for realizing novel high-pressure electrochemical systems.

25 ENERGY STORAGE↗