Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forcing fields”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Exploring the Structural, Dynamic, and Functional Properties of Metal–Organic Frameworks through Molecular Modeling

This review spotlights the role of atomic-level modeling in research on metal-organic frameworks (MOFs), especially the key methodologies of density functional theory (DFT), Monte Carlo (MC) simulations, and molecular dynamics (MD) simulations. The discussion focuses on how periodic and cluster-based DFT calculations can provide novel insights into MOF properties, with a focus on predicting structural transformations, understanding thermodynamic properties and catalysis, and providing information or properties that are fed into classical simulations such as force field parameters or partial charges. Classical simulation methods, highlighting force field selection, databases of MOFs for high-throughput screening, and the synergistic nature of MC and MD simulations, are described. By predicting equilibrium thermodynamic and dynamic properties, these methods offer a wide perspective on MOF behavior and mechanisms. Additionally, the incorporation of machine learning (ML) techniques into quantum and classical simulations is discussed. These methods can enhance accuracy, expedite simulation setup, reduce computational costs, as well as predict key parameters, optimize geometries, and estimate MOF stability. By charting the growth and promise of computational research in the MOF field, the aim is to provide insights and recommendations to facilitate the incorporation of computational modeling more broadly into MOF research.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electronic Structure Calculations of the Fundamental Interactions in Solvent Extraction Desalination

Solvent extraction desalination (SED) is a promising membrane-free desalination technology, but new solvents are still needed that can improve the molecular-level water selectivity, while minimizing solvent contamination in the aqueous phase. Improvements in solvent design for SED require a fundamental understanding of the molecular interactions governing solvation behavior. Herein, quantum chemical calculations are used to delineate different fundamental interactions in water-solvent systems involving 38 different solvents, with a particular focus on the role of the electrostatic potential characteristics. Furthermore, these interactions are rigorously analyzed using the independent gradient model based on Hirshfeld partition (IGMH), Symmetry-Adapted Perturbation Theory (SAPT), as well as the energy decomposition analysis based on force field (EDA-FF) to compare against a standard force field model. The solvation free energies and partition coefficients are also calculated to estimate water extraction performance in a bulk solvent system.

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

Computational Exploration of High Entropy Alloys as Promising Materials for Future Beam Windows

With the ever-increasing demand for high beam power, the currently used beam-intercepting devices (BIDs) such as targets, and beam windows may not be able to handle the high power required for future accelerator complexes or the lifetime may be reduced drastically. As beam power increases, the damage incurred by BIDs, including thermal shock, fatigue, and irradiation damage, also rises. Therefore, it is imperative to design materials that can withstand high beam power for longer lifetimes. High entropy alloys (HEAs) have emerged as potential alternative materials for designing next-generation BIDs. In this study, we primarily focus on materials for developing beam windows for next-generation accelerator complexes. We propose an integrated approach that combines various computational techniques to study and design new materials. Specifically, we use CALPHAD, density functional theory (DFT), and molecular dynamics (MD) to comprehensively investigate the defect properties of suitable HEAs, offering potential alternatives for future beam windows. We begin by scanning the extensive phase space provided by Cr-Mn-V-Ti-Al-Co HEAs, selecting 8 compositions after evaluating approximately 120,000 unique compositions using CALPHAD. We, then employ DFT-informed machine learning techniques to develop force-field parameters. Finally, MD simulations using these developed force-field parameters will be used to study the effects of radiation damage on the defect and mechanical properties of the selected alloys. This research explains the use of the CALPHAD approach and shows how critical modeling (DFT and MD) is in developing novel material such as HEAs. It also highlights the promising role of machine learning in this field. The results from this study will greatly improve the novel materials development to be used in next-generation accelerator components, leading to higher beam power and longer operational times of BIDs.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

The Role of Dimer Formation in the Nucleation of Superlattice Transformations and Its Impact on Disorder

The formation of defect-free two-dimensional nanocrystal (NC) superstructures remains a challenge as persistent defects hinder charge delocalization and related device performance. Understanding defect formation is an important step toward developing strategies to mitigate their formation. However, specific mechanisms of defect formation are difficult to determine, as superlattice phase transformations that occur during fabrication are quite complex and there are a variety of factors influencing the disorder in the final structure. Here, we use Molecular Dynamics (MD) and electron microscopy in concert to investigate the nucleation of the epitaxial attachment of lead chalcogenide (PbX, where X = S, Se) NC assemblies. We use an updated implementation of an existing reactive force field in an MD framework to investigate how initial orientational (mis)alignment of the constituent building blocks impacts the final structure of the epitaxially connected superlattice. This Simple Molecular Reactive Force Field (SMRFF) captures both short-range covalent forces and long-range electrostatic forces and allows us to follow orientational and translational changes of NCs during superlattice transformation. Our simulations reveal how robust the oriented attachment is with regard to the initial configuration of the NCs, measuring its sensitivity to both in-plane and out-of-plane misorientation. We show that oriented attachment nucleates through the initial formation of dimers, which corroborate experimentally observed structures. We present high-resolution structural analysis of dimers at early stages of the superlattice transformation and rationalize their contribution to the formation of defects in the final superlattice. Collectively, the simulations and experiments presented in this paper provide insights into the nucleation of NC oriented attachment, the impact of the initial configuration of NCs on the structural fidelity of the final epitaxially connected superlattice, and the propensity to form commonly observed defects, such as missing bridges and atomic misalignment in the superlattice due to the formation of dimers. Here, we present potential strategies to mitigate the formation of superlattice defects.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Towards sharp error analysis of extended Lagrangian molecular dynamics

The extended Lagrangian molecular dynamics (XLMD) method provides a useful framework for reducing the computational cost of a class of molecular dynamics simulations with constrained latent variables. The XLMD method relaxes the constraints by introducing a fictitious mass ε for the latent variables and solving a set of singularly perturbed ordinary differential equations. While favorable numerical performance of XLMD has been demonstrated in several different contexts in the past decade, mathematical analysis of the method remains scarce. Here, we propose the first error analysis of the XLMD method in the context of a classical polarizable force field model. While the dynamics with respect to the atomic degrees of freedom are general and nonlinear, the key mathematical simplification of the polarizable force field model is that the constraints on the latent variables are given by a linear system of equations. We prove that when the initial value of the latent variables is compatible in a sense that we define, XLMD converges as the fictitious mass ε is made small with $\mathscr{O}$(ε) error for the atomic degrees of freedom and with $\mathscr{O}$($\sqrt{ε}$) error for the latent variables, when the dimension of the latent variable d' is 1. Furthermore, when the initial value of the latent variables is improved to be optimally compatible in a certain sense, we prove that the convergence rate can be improved to $\mathscr{O}$(ε) for the latent variables as well. Numerical results verify that both estimates are sharp not only for d'=1, but also for arbitrary d'. In the setting of general d', we do obtain convergence, but with the non-sharp rate of $\mathscr{O}$($\sqrt{ε}$) for both the atomic and latent variables.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Computational screening of fly ash zeolite sorbents for boric acid removal

In the United States, many impoundments at coal-fired power plants contain elevated contaminants like arsenic, boron, barium, and selenium. Zeolites synthesized from fly ash show promise as sorbents for these contaminants. However, optimizing sorption capacity is challenging due to numerous possible topologies, silicon to aluminum (Si/Al) ratios, and cation types. In this study, molecular simulations are used to design cationic zeolites for boric acid adsorption. Force field models based on quantum mechanical calculations (PBE + D2) for Na-, Ca-, Mn-, and Fe-exchanged chabazite and LTA are presented. The new D2FF force fields reproduce DFT energies with about half the error of UFF. Zeolite performance depends on Si/Al ratio and cation type, with low Si/Al ratio chabazite (CHA) and phillipsite (PHI) zeolite frameworks exchanged with Ca 2+ or Na + /Ca 2+ mixtures showing the highest adsorption. In conclusion, these findings suggest tailored fly ash-derived zeolites could provide effective boron removal from leachate ponds.

CCR impoundment↗
DOE OSTI2026

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
DOE OSTI2025

Iodine recombination in xenon solvent: Clusters in the gas to liquid-like state transition

Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computational investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.

Cage effect↗
DOE OSTI2024

Assessing Melting and Solid–Solid Transition Properties of Choline Chloride via Molecular Dynamics Simulations

Choline chloride (ChCl) is used extensively as a hydrogen bond donor in deep eutectic solvents (DESs). However, determining its melting properties experimentally is challenging due to decomposition upon melting, leading to widely varying literature values. Accurate melting properties are crucial for understanding the solid–liquid phase behavior of ChCl-containing DESs. Here, we employ molecular dynamics simulations to compute the phase transitions of ChCl, testing a variety of atomistic force fields. We find that the results are sensitive to the choice of force field, but a melting temperature of 627 K and a melting enthalpy of 7.8 kJ/mol seem most reasonable, in good agreement with some literature values. Furthermore, we suggest these as the likely melting properties of ChCl, though the results are tentative due to limited experimental data for the liquid ChCl phase.

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

Quantitative Simulations of Siloxane Adsorption in Metal–Organic Frameworks

We present a transferable force field (FF) for simulating the bulk properties of linear and cyclic siloxanes and the adsorption of these species in metal–organic frameworks (MOFs). Unlike previous FFs for siloxanes, our FF accurately reproduces the vapor–liquid equilibria of each species in the bulk phase. The quality of our FF combined with the Universal Force Field using standard Lorentz–Berthelot combining rules for MOF atoms was assessed in a wide range of MOFs without open metal sites, showing good agreement with dispersion-corrected density functional theory calculations. Predictions with this FF show good agreement with the limited experimental data for siloxane adsorption in MOFs that is available. As an example of using the FF to predict adsorption properties in MOFs, we present simulations examining entropy effects in binary linear and cyclic siloxane mixture coadsorption in the large-pore MOF with structure code FOTNIN.

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

Mapping TpPa-1 covalent organic framework (COF) molecular interactions in mixed solvents via atomistic modeling and experimental study

Complex solvent environments continue to limit the widespread adoption of organic solvent nanofiltration (OSN) in many chemical industry applications. In this paper we employ a commercially available covalent organic framework (COF), TpPa-1, and force field models to molecularly map separation performance of TpPa-1 membrane in mixed solvents. To minimize time and length scale mismatch between atomistic modeling and experiments, solvent permeance was normalized with water in modeling and experimental results to enable direct comparison. Model outputs, such as organic solvent permeance and solute rejection rate, matched well with filtration results. Since the atomistic models assume that all mass transfer is via through-pore transport, the discrepancies between modeling and experimental results provide insights on the effect of linear polymer defects, adsorption and interstitial mass transfer on polycrystalline COF membrane performance. In sum, force field models can serve as digital twins of COF membranes to simulate separation processes while capturing the effects of COF structure, chemistry, and crystallinity on membrane performance in complex organic solvent environments. Finally, this approach will provide insight into future COF design and synthesis for persisting separation challenges.

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

Active learning for robust, high-complexity reactive atomistic simulations

Machine learned reactive force fields based on polynomial expansions have been shown to be highly effective for describing simulations involving reactive materials. Nevertheless, the highly flexible nature of these models can give rise to a large number of candidate parameters for complicated systems. In these cases, reliable parameterization requires a well-formed training set, which can be difficult to achieve through standard iterative fitting methods. In this paper, we present an active learning approach based on cluster analysis and inspired by Shannon information theory to enable semi-automated generation of informative training sets and robust machine learned force fields. The use of this tool is demonstrated for development of a model based on linear combinations of Chebyshev polynomials explicitly describing up to four-body interactions, for a chemically and structurally diverse system of C/O under extreme conditions. We show that this flexible training database management approach enables development of models exhibiting excellent agreement with Kohn–Sham density functional theory in terms of structure, dynamics, and speciation.

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

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

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

Hybrid Quantum Mechanical, Molecular Mechanical, and Machine Learning Potential for Computing Aqueous-Phase Adsorption Free Energies on Metal Surfaces

Performing reliable computer simulations of elementary processes occurring at metal–water interfaces is pivotal for novel catalyst design in sustainable energy applications. Computational catalyst design hinges on the ability to reliably and efficiently compute the potential energy surface (PES) of the system. Here, due to the large system sizes needed for studying processes at liquid water–metal interfaces, these systems can currently not be described using density functional theory (DFT). In this work, we used a hybrid quantum mechanical, molecular mechanical, and machine learning potential for studying the adsorption behavior of phenol, atomic hydrogen, 2-butanol, and 2-butanone on the (0001) facet of Ru under reducing conditions when Ru is not oxidized. Specifically, we describe the adsorbate and the surrounding metal atoms at the DFT level of theory. Here, we also considered the electrostatic field effect of the water molecules on adsorbate–metal interactions. Next, for the water–water and water–adsorbate interactions, we used established classical force fields. Finally, for the water–Ru surface interaction, for which no reliable force fields have been published, we used Behler–Parrinello high-dimensional neural network potentials (HDNNPs). Employing this setup, we used our explicit solvation for metal surface (eSMS) approach to compute the aqueous-phase effect on the low-coverage adsorption of selected molecules and atoms on the (0001) facet of Ru. In agreement with previous experimental and computational studies of oxygenated molecules over transition metal facets, we found that liquid water destabilizes the tested adsorbates on Ru(0001). Interestingly, our findings indicate that adsorbates on Ru are less affected by the presence of an aqueous phase than on other transition metals (e.g., Pt), highlighting the necessity of experimental investigations of Ru-based catalytic systems in liquid water.

Adsorption↗
DOE OSTI2023

Persistent ocean anomalies as a response to Northern Hemisphere heating induced by biomass burning variability

Biomass burning aerosol (BBA) emissions in the Coupled Model Intercomparison Project Phase 6 (CMIP6) historical forcing fields have enhanced temporal variability during the years 1997–2014 compared to earlier periods. Recent studies document that the corresponding inhomogeneous shortwave forcing over this period can cause changes in clouds, permafrost, and soil moisture, which contribute to a net terrestrial Northern Hemisphere warming relative to earlier periods. Here, we investigate the ocean response to the hemispherically asymmetric warming, using a 100-member ensemble of the Community Earth System Model version 2 Large Ensemble forced by two different BBA emissions (CMIP6 default and temporally smoothed over 1990–2020). Differences between the two subensemble means show that ocean temperature anomalies occur during periods of high BBA variability and subsequently persist over multiple decades. In the North Atlantic, surface warming is efficiently compensated for by decreased northward oceanic heat transport due to a slowdown of the Atlantic Meridional Overturning Circulation. In the North Pacific, surface warming is compensated for by an anomalous cross-equatorial cell (CEC) that reduces northward oceanic heat transport. The heat that converges in the South Pacific through the anomalous CEC is shunted into the subsurface and contributes to formation of long-lasting ocean temperature anomalies. The anomalous CEC is maintained through latitude-dependent contributions from narrow western boundary currents and basin-wide near-surface Ekman transport. Furthermore, these results indicate that interannual variability in forcing fields may significantly change the background climate state over long timescales, presenting a potential uncertainty in CMIP6-class climate projections forced without interannual variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

When more data hurts: Optimizing data coverage while mitigating diversity-induced underfitting in an ultrafast machine-learned potential

Machine-learned interatomic potentials (MLIPs) are becoming an essential tool in materials modeling. However, optimizing the generation of training data used to parametrize the MLIPs remains a significant challenge. This is because MLIPs can fail when encountering local environments too different from those present in the training data. The difficulty of determining a priori the environments that will be encountered during molecular dynamics simulation necessitates diverse, high-quality training data. Here, this study investigates how training data diversity affects the performance of MLIPs using the Ultra-Fast force field (UF 3 ) to model amorphous silicon nitride. We employ expert and autonomously generated data to create the training data and fit four force field variants to subsets of the data. Our findings reveal a critical balance in training data diversity: insufficient diversity hinders generalization, while excessive diversity can exceed the MLIP's learning capacity, reducing simulation accuracy. Specifically, we found that the UF 3 variant trained on a subset of the training data, in which nitrogen-rich structures were removed, offered vastly better prediction and simulation accuracy than any other variant. By comparing these UF 3 variants, we highlight the nuanced requirements for creating accurate MLIPs, emphasizing the importance of application-specific training data to achieve optimal performance in modeling complex material behaviors.

ab initio molecular dynamics↗
DOE OSTI2022

Structure of water-in-salt and water-in-bisalt electrolytes

We report a systematic diffraction study of two “water-in-salt” electrolytes and a “water-in-bisalt” electrolyte combining high-energy X-ray diffraction (HEXRD) with polarized and unpolarized neutron diffraction (ND) on both H 2 O and D 2 O solutions. The measurements provide three independent combinations of correlations between the different pairs of atom types that reveal the short- and intermediate-range order in considerable detail. The ND interference functions show pronounced peaks around a scattering vector Q ~ 0.5 Å -1 that change dramatically with composition, indicating significant rearrangements of the water network on a length scale around 12 Å. The experimental results are compared with two sets of Molecular Dynamics (MD) simulations, one including polarization effects and the other based on a non-polarizable force field. The two simulations reproduce the general shapes of the experimental structure factors and their changes with concentration, but differ in many detailed respects, suggesting ways in which their force fields might be modified to better represent the actual systems.

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

Multiscale Modeling of Vinyl-Addition Polynorbornenes: The Effect of Stereochemistry

Vinyl-addition polynorbornenes are candidates for designing high-performance polymers due to unique characteristics, which include a high glass transition temperature associated with a rigid backbone. Recent studies have established that the processability and properties of these polymers can be fine-tuned by using targeted substitutions. However, synthesis with different catalysts results in materials with distinct properties, potentially due to the presence of various stereoisomers that are difficult to quantify experimentally. Herein, we develop all-atom models of polynorbornene oligomers based on classical force fields and density functional theory. To establish the relationship between chemical architecture, chain conformations, and melt structure, we perform detailed molecular dynamics simulations with the fine-tuned atomistic force field and propose simpler coarse-grained descriptions to address the high molecular weight limit. All-atom simulations of oligomers suggest high glass transition temperatures in the range of 550–600 K. In the melt state (800 K), meso chains form highly rigid extended coils (C∞≈11) with amorphous structural characteristics similar to the X-ray diffraction data observed in the literature. In contrast, simulations with racemo chains predict highly helical tubular chain conformations that could promote assembly into crystalline structures.

Polymer Science↗