Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic models”

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.

931 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

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↗
DOE OSTI2024

A transferable classical force field to describe glyme based lithium solvate ionic liquids

A non-polarizable force field for lithium (Li + ) and bis(trifluoromethanesulfonyl)imide (TFSI – ) ions solvated in diglyme at around 0.2 mol fraction salt concentration was developed based on ab initio molecular dynamics (AIMD) simulations and a modified polymer consistent force field model. A force–torque matching based scheme, in conjunction with a genetic algorithm, was used to determine the Lennard-Jones (LJ) parameters of the ion–ion and ion–solvent interactions. This force field includes a partial charge scaling factor and a scaling factor for the 1–4 interactions. The resulting force field successfully reproduces the radial distribution function of the AIMD simulations and shows better agreement compared to the unmodified force field. The new force field was then used to simulate salt solutions with glymes of increasing chain lengths and different salt concentrations. The comparison of the MD simulations, using the new force field, with experimental data at different salt concentrations and AIMD simulations on equimolar concentrations of the triglyme system demonstrates the transferability of the force field parameters to longer glymes and higher salt concentrations. Furthermore, the force field appears to reproduce the features of the experimental x-ray structure factors, suggesting accuracy beyond the first solvation shell, for equimolar salt solutions using both triglyme and tetraglyme as the solvent. Altogether, the new force field was found to accurately reproduce the molecular descriptions of LiTFSI-glyme systems not only at various salt concentrations but also with glymes of different chain lengths. Thus, the new force field provides a useful and accurate tool to perform in silico studies of this family of systems at the atomistic level.

25 ENERGY STORAGE↗
DOE OSTI2026

Revealing Local Structures through Machine-Learning-Fused Multimodal Spectroscopy

Atomistic structures of materials offer valuable insights into their functionality. Determining these structures remains a fundamental challenge in materials science, especially for systems with defects. While both experimental and computational methods exist, each has limitations in resolving nanoscale structures. Core-level spectroscopies, such as X-ray absorption (XAS) or electron energy-loss spectroscopies (EELS), have been used to determine the local bonding environment and structure of materials. Recently, machine learning (ML) methods have been applied to extract structural and bonding information from XAS/EELS data. However, frameworks relying solely on a single data stream, defined as characterization data derived from a single element using one technique, are often insufficient because multiple local environments can yield similar spectral features, making it challenging to differentiate between competing structural hypotheses. Here, in this work, we address this challenge by integrating multimodal ab initio simulations, experimental data acquisition, and ML techniques for structure characterization. Our goal is to determine local structures and properties using EELS and XAS data from multiple elements and edges. To showcase our approach, we use various lithium nickel manganese cobalt (NMC) oxide compounds which are used for lithium ion batteries, including those with oxygen vacancies and antisite defects, as the sample material system. We successfully inferred local element content, ranging from lithium to transition metals, with quantitative agreement with experimental data. Beyond local element inference, we find that ML model based on multimodal spectroscopic data is able to determine whether local defects such as oxygen vacancy and antisites are present, a task which is impossible for single mode spectra or other experimental techniques. Furthermore, our framework is able to provide physical interpretability, bridging spectroscopy with the local atomic and electronic structures.

battery↗
DOE OSTI2020

A Quantum-Based Approach to Predict Primary Radiation Damage in Polymeric Networks

Initial atomistic-level radiation damage in chemically reactive materials is thought to induce reaction cascades that can result in undesirable degradation of macroscale properties. Ensembles of quantum-based molecular dynamics (QMD) simulations can accurately predict these cascades, but extracting chemical insights from the many underlying trajectories is a labor-intensive process that can require substantial a priori intuition. We develop here a general and automated graph-based approach to extract all chemically distinct structures sampled in QMD simulations and apply our approach to predict primary radiation damage of polydimethylsiloxane (PDMS), the main constituent of silicones. A postprocessing protocol is developed to identify underlying polymer backbone structures as connected components in QMD trajectories. Furthermore, these backbones form a repository of radiation-damaged structures. A scheme for extracting and updating a library of isomorphically distinct structures is proposed to identify the spanning set and aid chemical interpretation of the repository. The analyses are applied to ensembles of cascade QMD simulations in which the four element types in PDMS are selectively excited in primary knock-on atom events. Our approach reveals a much higher degree of combinatorial complexity in this system than was inferred through radiolysis experiments. Probabilities are extracted for radiation-induced network changes including formation of branch points, carbon linkages, cycles, bond scissions, and carbon uptake into the Si–O siloxane backbone network. The general analysis framework presented here is readily extendable to modeling chemical degradation of other polymers and molecular materials and provides a basis for future quantum-informed multiscale modeling of radiation damage.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Combining Reactive Quantum-Mechanical Molecular-Dynamics Simulations with Mutagenesis, Crystallography, and Enzyme Kinetics to Reveal Plausible Steps of Isocyanide Hydratase Catalysis

A complete understanding of enzyme mechanisms requires atomistic details of chemical reactions. Quantum-based molecular dynamics simulations (QMD) are a potential source of this information, but trade-offs between accuracy and computational cost have limited their use. We previously developed extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) methods that leverage a negligible compromise in accuracy to substantially decrease the cost of QMD simulations. Here, we develop a reactive QMD approach using the latest XL-BOMD formulation, which enables efficient simulations of highly reactive systems, and use it to investigate mechanisms of intermediate formation in isocyanide hydratase (ICH) catalysis. In QMD simulations, molecular analogs of ICH active site residues reacted with para-nitrophenyl isocyanide, forming a thioimidate. Analysis of simulated atomic configurational and charge dynamics revealed a pathway where protonation of the isocyanide carbon occurs prior to thioimidate formation and suggested a possible role of Asp17 as a proton donor in the early phase of ICH catalysis. To test whether the pathway seen using the reactive QMD approach might be relevant to ICH catalysis, we performed X-ray crystallography and pre-steady-state enzyme kinetics studies of wild-type and D17N mutant ICH. Both the structure and kinetics are sensitive to the D17N mutation in a manner that is consistent with the order of the reaction steps seen in the simulations. Mobile protons play essential roles in many enzymes, yet they are difficult to observe experimentally, making the ordering of proton-dependent steps ambiguous in many enzyme mechanisms. The ability to directly simulate model reactions for the design of experiments that provide information about enzyme mechanisms involving mobile protons demonstrates the significance of our reactive QMD approach and motivates further biological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Challenges of conventional iterative all-atom and coarse-grained multiscale molecular dynamics

In this work, we evaluate the biomolecular dynamics behaviors when conventionally iterating between all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations over multiple cycles. We implemented the workflow to iterate between AA and CG in OpenMM, namely the iterative multiscale MD (iMMD) simulation workflow. In particular, we aim to identify practical applications for iterating between AA and CG simulations in a conventional manner without any constraints or model modifications. We evaluate the iMMD workflow on four representative systems, spanning folding of two soluble proteins and protein-protein as well as protein-lipid interactions of two membrane proteins. We observe that iteration between AA and CG representations could help the soluble proteins exit undesirable metastable states to fold, resulting from random protein structural distortions due to cycling. Consequently, the most reliable use of iterative AA and CG simulations appears to be to accelerating complex lipid mixing for membrane-bound protein systems rather than sampling protein conformational space. Our work explores the practical usages and limitations for iterative AA and CG simulations using readily available AA and CG force fields. The evaluated iMMD workflow in OpenMM is made available at https://github.com/lanl/iMMD.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRSDate not supplied

Fracture Prediction of Epoxy Resin using Morse Bond Potential Embedded in GAFF using Molecular Dynamics Simulations

A generic method is developed to investigate the fracture behavior of an epoxy resin using atomistic molecular (MD) dynamics simulations. The epoxy system consisted of the stoichiometric mixture of tetra-/tri-functional epoxies and di-functional hardener molecules and was cured with a step-growth crosslinking algorithm. A hybrid force field (hFF) in which Morse bond function were added in second generation-general Amber Force Field (GAFF2) was used to capture bond breaking of the crosslinked epoxy system under uniaxial deformation. The Morse bond parameters for covalent bonds in the crosslinked backbone were fit to dissociation curves computed by CASPT2/6-311+G**. After the systems in various crosslink densities were obtained and equilibrated, the full length of stress-strain (𝜎-𝜀) curve was produced including initial elastic regime, yielding, plastic flow, strain hardening, and progressive failure. We found that the hFF was effective in studying the dissociation of crosslinked polymers and understanding the failure mechanism. In the results, the hFF delivered the same elastic property obtained from the unmodified force field GAFF. The stress-strain curve, however, began to be distinct in the plastic regime due to broken covalent bonds capable of altering the tensile behavior. Once the stress reached maximum, it progressively decreased until the ultimate failure (𝜎=0). Since hFF employed the existing GAFF functions with including a minimal change of its bond potential, this approach is easier to perform and computationally more efficient for unveiling the fracture behavior of polymer materials at the molecular level than other approaches with reactive force fields. We expect this approach will be utilized to accelerate the material-by-design process for thermosets by incorporate data from molecular models.

Changwoon Jang↗
DOE OSTI2026

Melatonin-Induced Modulation of Cholesterol-Enriched Model Neuronal Membranes

Melatonin, a hormone primarily produced by the brain’s pineal gland, not only regulates circadian rhythms, but also influences the structural and biophysical properties of neuronal membranes. Its amphiphilic nature enables direct incorporation into lipid bilayers and preferential interactions with cholesterol-rich lipid rafts, critical hubs for cellular signaling and membrane organization. Despite increasing recognition of its membrane activity, the molecular basis of melatonin’s interactions with coexisting liquid-ordered (L o ) and liquid-disordered (L d ) phases remains unclear. Here, in this study, we combine small-angle neutron scattering (SANS) and all-atom molecular dynamics simulations to examine model neuronal membranes composed of DSPC, DOPC, POPC, and cholesterol. Our results show that melatonin preserves domain morphology while adopting distinct orientations within the bilayer and at the membrane interface, allowing both lateral and transmembrane bridging across lipid phases. These findings establish the molecular underpinnings of melatonin’s modulation of membrane heterogeneity and provide strong support for its receptor-independent actions.

atomistic simulations↗
DOE OSTI2021

Grain size dependencies of intergranular solute segregation in nanocrystalline materials

Grain boundary segregation is of increasing importance for stabilizing nanocrystalline alloys, and while segregation is typically modeled as a scale-independent phenomenon, at the finest nanocrystalline grain sizes there are several sources of grain size dependence. Here we provide quantitative detail on two grain size dependencies not previously elaborated, both of which arise from the spectrality of intergranular segregation sites. The first of these arises because a change in grain size affects which portions of the segregation spectrum are occupied by the solutes, with no change in the shape or nature of the spectrum itself. The second of these arises from the presence of higher-order grain junctions, which have unique segregation behaviors and an increasing relative presence at finer grain sizes. This paper uses molecular statics to quantify these size dependencies in detail for the case of an Al(Mg) binary alloy. Together, these two effects combined have a significant effect on the effective McLean segregation energy of more than 6 kJ/mol (from approximately -18.5 kJ/mol at 5 nm to -12.5 kJ/mol at more than 40 nm), with the majority of grain size dependence arising from the spectrality of segregation sites.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Controlling Deformation in Al/Ti: How Interface Roughness and Orientation Drive Bimetal Mechanics

The microstructural characteristics and morphology of interfaces in metals can be crucial in governing the initiation of plasticity and early deformation mechanisms under extreme conditions. During high-strain-rate deformation, these interfaces significantly affect dislocation nucleation, twinning, and other mechanisms that directly impact material strength and failure. Despite their importance, a substantial knowledge gap remains between the observed macroscopic material behavior and the underlying role of bimetal interfaces in plasticity initiation. Here, to address this gap, large-scale molecular dynamics simulations are performed on Al/Ti bimetal structures to examine the effect of interface characteristics on the onset of plasticity under uniaxial compression. Specifically, this study investigates how interface roughness (flat vs waveform interfaces) modifies the initiation of plastic events in the microstructure. Atomistic simulations indicate that interface roughness (a microscopic behavior) reduces the stress required for dislocation nucleation, thereby reducing the peak stress relative to a flat interface. For the square interface, the step height strongly influences plasticity initiation by setting the separation of locally flat regions. Varying the interface rotation relative to the loading direction (macroscopic behavior) reveals that the peak stress for both flat and waveform interfaces initially decreases and then increases with rotation, accompanied by a shift in dominant mechanisms—from Al twinning (flat) to interface sliding and Ti-dominated phase transformation and twinning at lower angles. The change in the stress-strain slope during initial compression reflects how rotation alters the resolved shear stress and activates different slip systems. Overall, the study provides valuable insights into the role of bimetal interfaces in controlling plasticity initiation and early deformation pathways for Al/Ti under extreme loading conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hall–Petch strengthening limit through partially active segregation in nanocrystalline Ag-Cu alloys

The breakdown from grain-size strengthening to softening mechanisms is generally well understood for high-purity nanocrystalline materials when the mean grain size decreases to the nanometer range. In nanocrystalline alloys, however, the stabilization of nanosized grains by grain-boundary solute segregation complicates the above mechanisms. Moreover, current segregation models have little predictive power for determining the optimal solute content that maximizes Hall-Petch strengthening effects. In this article, using large-scale hybrid Monte-Carlo/molecular dynamic simulations, we present a systematic study of the Hall–Petch breakdown in Cu-segregated Ag alloys with grain sizes ranging from 8 nm to 59 nm, where three concentration-dependent regimes of plasticity are described: (1) Classical segregation strengthening behavior at low solute contents, (2) shear band-induced softening at high solute contents, and (3) a previously unknown, but extended plateau of maximum strengths for intermediate solute contents from 4 to 15 at.%, which we term as nanocrystalline Sterling silver. We find that flow strengths in nanocrystalline Sterling alloys naturally exhibit a zero-slope limit at the smallest grain sizes that is well below the ideal Hall–Petch strengthening trend. This phenomenon results from partially active grain-boundary segregation that acts to influence interfacial plasticity in some, but not all, grain boundary regions. Our findings amplify the atomic nature of solute segregation and interaction at grain boundaries and its complex roles on grain boundary-mediated plasticity mechanisms in nanocrystalline alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Competing conduction mechanisms in high performance carbon nanotube fibers

The performance of carbon nanotube (CNT) cables, a contender for copper-wire replacement, is tied to its metallic and semi-conducting-like conductivity responses with temperature; the origin of the semi-conducting-like response however is an underappreciated incongruity in literature. With controlled aspect-ratio and doping-degree, over 61 unique cryogenic experiments including anisotropy and Hall measurements, CNT cable performance is explored at extreme temperatures (65 mK) and magnetic-fields (60 T). A semi-conducting-like conductivity response with temperature becomes temperature-independent approaching absolute-zero, uniquely demonstrating for the first time the necessity of heterogeneous fluctuation induced tunneling; complete de-doping leads to localized hopping, contrasting graphite's pure metallic-like response. High-field magneto-resistance (including novel +22 % longitudinal magneto-resistance near room-temperature) is analyzed with hopping and classical two-band models, both yielding a similar parameter useful for conductor development. Varying field-orientation angle uncovers significant two- and four-fold symmetries that are shown to be from Aharonov-Bohm-like corrections to the curvature-induced bandgap, a first for macroscale CNT fibers. Tight-binding calculations using Green's Function formalism model the largest, coherent transport to-date in commensurate CNT bundles in magnetic-field, revealing non-uniform transmission across bundle cross-sections with doping restoring uniformity; independent of doping, transport in bundle-junction-bundle systems are predominantly from CNTs adjacent to the other bundle—demonstrating that smaller bundles are more efficient for electronic transport. The final impact is predicting the ultimate conductivity of heterogeneous CNT cables using temperature and field-dependent transport, surpassing conductivity of traditional metals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Minimal Information Set To Enable Verifiable Theoretical Battery Research

Batteries are an enabling technology for addressing sustainability through the electrification of various forms of transportation (1) and grid storage. (2) Batteries are truly multi-scale, multi-physics devices, and accordingly various theoretical descriptions exist to understand their behavior (3-5) ranging from atomistic details to techno-economic trends. As we explore advanced battery chemistries (6,7) or previously inaccessible aspects of existing ones, (8-10) new theories are required to drive decisions. (11-13) The decisions are influenced by the limitations of the underlying theory. Advanced theories used to understand battery phenomena are complicated and require substantial effort to reproduce. However, such constraints should not limit the insights from these theories. We can strive to make the theoretical research verifiable such that any battery stakeholder can assess the veracity of new theories, sophisticated simulations or elaborate analyses. We distinguish verifiability, which amounts to “Can I trust the results, conclusions and insights and identify the context where they are relevant?”, from reproducibility, which ensures “Would I get the same results if I followed the same steps?” With this motivation, we propose a checklist to guide future reports of theoretical battery research in Table 1. We hereafter discuss our thoughts leading to this and how it helps to consistently document necessary details while allowing complete freedom for creativity of individual researchers. Given the differences between experimental and theoretical studies, the proposed checklist differs from its experimental counterparts. (14,15) This checklist covers all flavors of theoretical battery research, ranging from atomic/molecular calculations (16-19) to mesoscale (20,21) and continuum-scale interactions, (9,22) and techno-economic analysis. (23,24) Finally, as more and more experimental studies analyze raw data, (25) we feel this checklist would be broadly relevant.

25 ENERGY STORAGE↗