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

CHGNet as a pretrained universal neural network potential for charge-informed atomistic modelling

Large-scale simulations with complex electron interactions remain one of the greatest challenges for atomistic modelling. Although classical force fields often fail to describe the coupling between electronic states and ionic rearrangements, the more accurate ab initio molecular dynamics suffers from computational complexity that prevents long-time and large-scale simulations, which are essential to study technologically relevant phenomena. Here we present the Crystal Hamiltonian Graph Neural Network (CHGNet), a graph neural network-based machine-learning interatomic potential (MLIP) that models the universal potential energy surface. CHGNet is pretrained on the energies, forces, stresses and magnetic moments from the Materials Project Trajectory Dataset, which consists of over 10 years of density functional theory calculations of more than 1.5 million inorganic structures. The explicit inclusion of magnetic moments enables CHGNet to learn and accurately represent the orbital occupancy of electrons, enhancing its capability to describe both atomic and electronic degrees of freedom. We demonstrate several applications of CHGNet in solid-state materials, including charge-informed molecular dynamics in Li x MnO 2 , the finite temperature phase diagram for Li x FePO 4 and Li diffusion in garnet conductors. We highlight the significance of charge information for capturing appropriate chemistry and provide insights into ionic systems with additional electronic degrees of freedom that cannot be observed by previous MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular Dynamics Study of the Effect of Grafting Density on Ion Diffusivity in a MARTINI Coarse-Grained Strong Polyelectrolyte Brush

Because surface-grafted polyelectrolyte brushes (PEBs) are responsive to external stimuli, such as electric fields and ionic strength, PEBs are attractive for applications ranging from drug delivery to separation technologies. Essential to PEB utilization is understanding how critical parameters like grafting density (σ) impact the PEB structure and the dynamics of the PEB and counterions. To study the effect of σ on PEB and the counterion structure and dynamics, we fine-tune a coarse-grained model that retains the chemical specificity of a strong polyelectrolyte, poly[(2-(methacryloyloxy)ethyl) trimethylammonium chloride] (PMETAC), using the MARTINI force field. Using “salt-free” conditions where the counterion concentration balances the charge on the brush, we build coarse-grained (CG) molecular dynamics simulations for MARTINI PMETAC brushes (N = 150 monomers; M W = 31.2 kg/mol) at experimentally relevant values of σ = 0.05, 0.10, 0.20, and 0.40 chains/nm 2 . Using 5 μs simulations, we investigate the effects of grafting density on the PEB structure, ion dissociation dynamics, polymer mobility, and counterion diffusivity. Results show that competition between electrostatic interactions, steric hindrance, and polymer mobility controls counterion diffusivity. Finally, the interplay of these factors leads to diffusivity that depends non-monotonically on σ, with counterion diffusivity peaking at an intermediate σ = 0.10 chains/nm 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling the Behavior of Complex Aqueous Electrolytes Using Machine Learning Interatomic Potentials: The Case of Sodium Sulfate

Understanding the structure and thermodynamics of solvated ions is essential for advancing applications in electrochemistry, water treatment, and energy storage. While ab initio molecular dynamics methods are highly accurate, they are limited by short accessible time and length scales whereas classical force fields struggle with accuracy. Herein, we explore the structure and thermodynamics of complex monovalent-divalent ion pairs using Na 2 SO 4 (aq) as a case study by applying a machine learning interatomic potential (MLIP) trained on density functional theory (DFT) data. Our MLIP-based approach reproduces key bulk properties such as density and radial distribution functions of water. We provide the hydration structure of the sodium and sulfate ions in the 0.1–2 M concentration range and the one-dimensional and two-dimensional potentials of mean force for the sodium–sulfate ion pairing at the low concentration limit (0.1 M), which are inaccessible to DFT. At low concentrations, the sulfate ion is strongly solvated, leading to the stabilization of solvent-separated ion pairs over contact ion pairs. Minimum energy pathway analysis revealed that coordinating two sodium ions with a sulfate ion is a multistep process whereby the sodium ions coordinate to the sulfate ion sequentially. Finally, we demonstrate that MLIPs allow the study of solvated ions beyond simple monovalent pairs with DFT-level accuracy in their low concentration limit (0.1 M) via statistically converged properties from ns-long simulations.

anions↗
DOE OSTI2020

Detecting Fractional Chern Insulators in Optical Lattices Through Quantized Displacement

The realization of interacting topological states of matter such as fractional Chern insulators (FCIs) in cold atom systems has recently come within experimental reach due to the engineering of optical lattices with synthetic gauge fields providing the required topological band structures. However, detecting their occurrence might prove difficult since transport measurements akin to those in solid state systems are challenging to perform in cold atom setups and alternatives have to be found. Here, we show that for a ν = 1/2 FCI state realized in the lowest band of a Harper-Hofstadter model of interacting bosons confined by a harmonic trapping potential, the fractionally quantized Hall conductivity σ xy can be accurately determined by the displacement of the atomic cloud under the action of a constant force which provides a suitable experimentally measurable signal for detecting the topological nature of the state. Using matrix-product state algorithms, we show that, in both cylinder and square geometries, the movement of the particle cloud in time under the application of a constant force field on top of the confining potential is proportional to σ xy for an extended range of field strengths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ~8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E int ) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)- based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E int values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ~±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

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

Stress field and interaction forces between dislocations and precipitate distributions

A computational method is developed for calculation of the stress field and interaction forces between dislocations and precipitates of arbitrary shape and distribution. The internal stress generated by precipitates due to coherency strain is implemented within the discrete dislocation dynamics (DDD) framework. The s-version finite element method (s-FEM), which models a precipitate of arbitrary shape using a local mesh is used to calculate coherency stress fields. The method facilitates meshing precipitate volumes of arbitrary geometry, and multiple local meshes can be superimposed at various positions of the global mesh. Accuracy and convergence conditions of the method are established. For a single precipiate, the method is shown to be 3.5 times faster than a standard FEM calculation for the same accuracy, and the gain in speed increases with the number of precipitates. The CRSS of spherical precipitates is found to be higher than disk-shaped precipitates of the same volume fraction.

42 ENGINEERING↗
DOE OSTI2024

Simple and Accurate One-Body Energy and Dipole Moment Surfaces for Water and Beyond

Water is often the testing ground for new, advanced force fields. While advanced functional forms for intermolecular interactions have been integral to the development of accurate water models, less attention has been paid to a transferable model for intramolecular valence terms. In this work, we present a one-body energy and dipole moment surface model, named 1B-UCB, that is simple yet accurate and can be feasibly adapted for both standard and advanced potentials. 1B-UCB for water is comparable in accuracy to those with much more complex functional forms, despite having drastically fewer parameters. The parametrization protocol has been implemented as part of the Q-Force automated workflow and requires only a quantum mechanical Hessian calculation as reference data, hence allowing it to be easily extended to a variety of molecular systems beyond water, which we demonstrate on a selection of small molecules with different symmetries.

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

Test samples and infrastructure for accelerator magnet developments

LHC still does not operate at full energy due to training features in magnets. Higher field practical superconductors (Nb 3 Sn) show much worse training behavior than NbTi. Most magnet performance issues relate to quenching, and it is this phenomenon and its dependencies that need to be understood for us to succeed addressing them in magnets. For pragmatic and scientific reasons underlying mechanisms are hard to investigate in real magnets. Available data suggest that emulating local conditions at quench location may be enough to arrive at a good understanding of what drives observed behavior. "Local" conditions imply emulation by small samples of cables/"stacks" powered within a facility with controllable conditions, including external magnetic and force fields. This paper argues about the necessity of this approach, while reminding current understanding and directions of quench studies, and pointing out inadequacy of available techniques. Improved facilities and change of perspective are required for consistent and affordable development.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

JARVIS-Leaderboard: a large scale benchmark of materials design methods

Abstract Lack of rigorous reproducibility and validation are significant hurdles for scientific development across many fields. Materials science, in particular, encompasses a variety of experimental and theoretical approaches that require careful benchmarking. Leaderboard efforts have been developed previously to mitigate these issues. However, a comprehensive comparison and benchmarking on an integrated platform with multiple data modalities with perfect and defect materials data is still lacking. This work introduces JARVIS-Leaderboard, an open-source and community-driven platform that facilitates benchmarking and enhances reproducibility. The platform allows users to set up benchmarks with custom tasks and enables contributions in the form of dataset, code, and meta-data submissions. We cover the following materials design categories: Artificial Intelligence (AI), Electronic Structure (ES), Force-fields (FF), Quantum Computation (QC), and Experiments (EXP). For AI, we cover several types of input data, including atomic structures, atomistic images, spectra, and text. For ES, we consider multiple ES approaches, software packages, pseudopotentials, materials, and properties, comparing results to experiment. For FF, we compare multiple approaches for material property predictions. For QC, we benchmark Hamiltonian simulations using various quantum algorithms and circuits. Finally, for experiments, we use the inter-laboratory approach to establish benchmarks. There are 1281 contributions to 274 benchmarks using 152 methods with more than 8 million data points, and the leaderboard is continuously expanding. The JARVIS-Leaderboard is available at the website: https://pages.nist.gov/jarvis_leaderboard/

36 MATERIALS SCIENCE↗
DOE OSTI2024

Geometry-complete perceptron networks for 3D molecular graphs

Abstract Motivation The field of geometric deep learning has recently had a profound impact on several scientific domains such as protein structure prediction and design, leading to methodological advancements within and outside of the realm of traditional machine learning. Within this spirit, in this work, we introduce GCPNet, a new chirality-aware SE(3)-equivariant graph neural network designed for representation learning of 3D biomolecular graphs. We show that GCPNet, unlike previous representation learning methods for 3D biomolecules, is widely applicable to a variety of invariant or equivariant node-level, edge-level, and graph-level tasks on biomolecular structures while being able to (1) learn important chiral properties of 3D molecules and (2) detect external force fields. Results Across four distinct molecular-geometric tasks, we demonstrate that GCPNet’s predictions (1) for protein–ligand binding affinity achieve a statistically significant correlation of 0.608, more than 5%, greater than current state-of-the-art methods; (2) for protein structure ranking achieve statistically significant target-local and dataset-global correlations of 0.616 and 0.871, respectively; (3) for Newtownian many-body systems modeling achieve a task-averaged mean squared error less than 0.01, more than 15% better than current methods; and (4) for molecular chirality recognition achieve a state-of-the-art prediction accuracy of 98.7%, better than any other machine learning method to date. Availability and implementation The source code, data, and instructions to train new models or reproduce our results are freely available at https://github.com/BioinfoMachineLearning/GCPNet.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Actinium–DOTA coordination in water from hybrid ML/MM: Structure, free energies, and water-exchange pathways

Quantitative simulation of trivalent ƒ-block chelates in water remains challenging because bonded and non-bonded force-field models make different approximations for coordination structure, exchange dynamics, and ion–ligand interactions in highly charged systems. Here, we develop a hybrid machine-learning/molecular-mechanics (ML/MM) framework for Ac 3+ –DOTA in explicit solvent by training an E(3)-equivariant neural network potential (MACELES) on mechanically embedded QM/MM data for Ac aquo and Ac–DOTA species and coupling it to NAMD 2.14 with particle-mesh Ewald electrostatics. Nanosecond ML/MM trajectories remain numerically stable and preserve chelate integrity, yielding a compact DOTA inner shell with an inner-sphere water coordination number of CN Ac,O w ≈ 1.7 arising from a dynamic equilibrium between one- and two-water states (37.5% and 59.9% of frames; three waters 2.5%). A 5 ns potential of mean force shows two low-lying basins at CN Ac,O w ≈ 1 and CN Ac,O w ≈ 2. DFT end-state free energies are consistent with the ML/MM profile, and DFT minimum-energy paths provide a qualitative electronic-structure reference for the observed basin connectivity. State-resolved kinetics reveal picosecond water-exchange pathways that couple hydration changes to transient DOTA arm fluctuations, and training-set comparisons show that temperature-matched Ac–DOTA data optimize energy/force accuracy while more diverse solvated data improve charge prediction. Overall, the present hybrid ML/MM model provides a practical description of Ac 3+ –DOTA hydration thermodynamics and short-time exchange behavior in explicit water at MD-like cost.

Actinium↗
DOE OSTI2023

What do far-infrared spectra of solitary water in “water-in-solvent” systems reveal about water’s solvation and dynamics?

Classical molecular dynamics simulations of water in ionic and dipolar solvents were used to interpret the far-infrared (FIR) rotation/libration spectra of “solitary water” in terms of water’s rotational dynamics and interactions with solvents. Seven solvents represented by nonpolarizable all-atom force fields and a series of idealized variable-charge solvents were used to span the range of solvent polarities (hydrogen bonding) studied experimentally. Simulated spectra capture the solvent dependence observed, as well as the relationship between the frequencies of water libration (νL) and OH stretching bands (νOH). In more strongly interacting solvents, simulated νL are ∼20% higher than those of experiment. In all solvents, the simulated spectra are composites of rotational motions about the two axes perpendicular to water’s dipole moment, and the different frequencies of these two motions are responsible for the breadth of the libration band and the bimodal shape observed in halide ionic liquids. Simulations overestimate the separation of these two components in most solvents. The character of water rotational motions changes markedly with solvent polarity, from quasi-free rotation in nonpolar and weakly polar solvents to highly constrained libration in strongly hydrogen bonding environments. The changeover to librational motions dominating the spectrum occurs between solvents such as benzene (νL ∼ 250 cm−1) and acetonitrile (νL ∼ 400 cm−1). For solvents in the latter category, the mean frequency of the experimental FIR band provides a direct measure of mean-squared torques and, therefore, force constants associated with interactions constraining water’s librational motion.

Chemistry↗
DOE OSTI2024

Direct Synthesis of Layer-Tunable and Transfer-Free Graphene on Device-Compatible Substrates Using Ion Implantation Toward Versatile Applications

Direct synthesis of layer-tunable and transfer-free graphene on technologically important substrates is highly valued for various electronics and device applications. State of the art in the field is currently a two-step process: a high-quality graphene layer synthesis on metal substrate through chemical vapor deposition (CVD) followed by delicate layer transfer onto device-relevant substrates. Here, we report a novel synthesis approach combining ion implantation for a precise graphene layer control and dual-metal smart Janus substrate for a diffusion-limiting graphene formation to directly synthesize large area, high quality, and layer-tunable graphene films on arbitrary substrates without the post-synthesis layer transfer process. Carbon (C) ion implantation was performed on Cu–Ni film deposited on a variety of device-relevant substrates. A well-controlled number of layers of graphene, primarily monolayer and bilayer, is precisely controlled by the equivalent fluence of the implanted C-atoms (1 monolayer ~4 × 10 15 C-atoms/cm 2 ). Upon thermal annealing to promote Cu-Ni alloying, the pre-implanted C-atoms in the Ni layer are pushed toward the Ni/substrate interface by the top Cu layer due to the poor C-solubility in Cu. As a result, the expelled C-atoms precipitate into a graphene structure at the interface facilitated by the Cu-like alloy catalysis. After removing the alloyed Cu-like surface layer, the layer-tunable graphene on the desired substrate is directly realized. The layer-selectivity, high quality, and uniformity of the graphene films are not only confirmed with detailed characterizations using a suite of surface analysis techniques but more importantly are successfully demonstrated by the excellent properties and performance of several devices directly fabricated from these graphene films. Molecular dynamics (MD) simulations using the reactive force field (ReaxFF) were performed to elucidate the graphene formation mechanisms in this novel synthesis approach. With the wide use of ion implantation technology in the microelectronics industry, this novel graphene synthesis approach with precise layer-tunability and transfer-free processing has the promise to advance efficient graphene-device manufacturing and expedite their versatile applications in many fields.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Revealing the Defect‐Driven Ferroelectric Mechanisms of Aluminum Nitride

Wurtzite III-nitride compounds are CMOS-compatible with widespread industrial interest to exercise ferroelectricity, despite their polar structure being highly resistant to polarization reversal. Here, we induce and tune ferroelectric properties in w-AlN via direct-write ion-beam processing, using nanoscale patterned defect engineering as a post-growth alternative to conventional cation substitution. Nanometric piezoresponse spectroscopy of the focused He + beam patterned defect concentrations in ferroelectric Al 0.92 B 0.08 N measures a localized 10x enhancement in effective piezoresponse and 40% reduction in switching barrier. The irradiation-induced point defects convert piezoelectric AlN into a ferroelectric system with site-saturated nucleation and raise the dielectric susceptibility, switched polarization, and effective piezoelectric coefficient. Enhanced defect-lattice interactions in AlN increase carrier conduction and phonon scattering loss but preserve long-range crystallinity. Here, based on atomistic analysis of nudged elastic band density functional theory calculations and reactive force field simulations, both nitrogen vacancies and defect complexes disrupt bond ordering, facilitating a line-by-line low-barrier switching of pristine AlN.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accelerated Discovery of CH 4 Uptake Capacity Metal–Organic Frameworks Using Bayesian Optimization

Abstract High‐throughput computational studies for discovery of metal–organic frameworks (MOFs) for separations and storage applications are often limited by the costs of computing thermodynamic quantities. Recent such studies at the time of writing may use ab initio results for a narrow selection of MOFs or empirical force‐field methods for larger selections. Here, a proof‐of‐concept study is conducted using Bayesian optimization on CH 4 uptake capacity of hypothetical MOFs for an existing dataset (Wilmer et al., Nature Chem. 2012, 4 , 83). It is shown that less than 0.1% of the database needs to be screened with the Bayesian optimization approach to recover the top candidate MOFs. This opens the possibility for efficient screening of MOF databases using accurate ab initio calculations for future adsorption studies on a minimal subset of MOFs. Furthermore, Bayesian optimization and the surrogate model presented here can offer interpretable material design insights and the framework will be applicable in the context of other target properties.

Taw, Eric↗
DOE OSTI2021

Proton Traffic Jam: Effect of Nanoconfinement and Acid Concentration on Proton Hopping Mechanism

Abstract The properties of the water network in concentrated HCl acid pools in nanometer‐sized reverse nonionic micelles were probed with TeraHertz absorption, dielectric relaxation spectroscopy, and reactive force field simulations capable of describing proton hopping mechanisms. We identify that only at a critical micelle size of W 0 =9 do solvated proton complexes form in the water pool, accompanied by a change in mechanism from Grotthuss forward shuttling to one that favors local oscillatory hopping. This is due to a preference for H + and Cl − ions to adsorb to the micelle interface, together with an acid concentration effect that causes a “traffic jam” in which the short‐circuiting of the hydrogen‐bonding motif of the hydronium ion decreases the forward hopping rate throughout the water interior even as the micelle size increases. These findings have implications for atmospheric chemistry, biochemical and biophysical environments, and energy materials, as transport of protons vital to these processes can be suppressed due to confinement, aggregation, and/or concentration.

Adams, Ellen M.↗
DOE OSTI2023

Dynamic Monkey Bar Mechanism of Superionic Li‐ion Transport in LiTaCl 6

Abstract The LiTaCl 6 solid electrolyte has the lowest activation energy of ionic conduction at ambient conditions (0.165 eV), with a record high ionic conductivity for a ternary compound (11 mS cm −1 ). However, the mechanism has been unclear. We train machine‐learning force fields (MLFF) on ab initio molecular dynamics (AIMD) data on‐the‐fly and perform MLFF MD simulations of AIMD quality up to the nanosecond scale at the experimental temperatures, which allows us to predict accurate activation energy for Li‐ion diffusion (at 0.164 eV). Detailed analyses of trajectories and vibrational density of states show that the large‐amplitude vibrations of Cl − ions in TaCl 6 − enable the fast Li‐ion transport by allowing dynamic breaking and reforming of Li−Cl bonds across the space in between the TaCl 6 − octahedra. We term this process the dynamic‐monkey‐bar mechanism of superionic Li + transport which could aid the development of new solid electrolytes for all‐solid‐state lithium batteries.

Lei, Ming↗
DOE OSTI2021

Proton Traffic Jam: Effect of Nanoconfinement and Acid Concentration on Proton Hopping Mechanism

Abstract The properties of the water network in concentrated HCl acid pools in nanometer‐sized reverse nonionic micelles were probed with TeraHertz absorption, dielectric relaxation spectroscopy, and reactive force field simulations capable of describing proton hopping mechanisms. We identify that only at a critical micelle size of W 0 =9 do solvated proton complexes form in the water pool, accompanied by a change in mechanism from Grotthuss forward shuttling to one that favors local oscillatory hopping. This is due to a preference for H + and Cl − ions to adsorb to the micelle interface, together with an acid concentration effect that causes a “traffic jam” in which the short‐circuiting of the hydrogen‐bonding motif of the hydronium ion decreases the forward hopping rate throughout the water interior even as the micelle size increases. These findings have implications for atmospheric chemistry, biochemical and biophysical environments, and energy materials, as transport of protons vital to these processes can be suppressed due to confinement, aggregation, and/or concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗