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.

At least 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Ethanol-Induced Condensation and Decondensation in DNA-Linked Nanoparticles: A Nucleosome-like Model for the Condensed State

Inspired by the conventional use of ethanol to induce DNA precipitation, ethanol condensation has been applied as a routine method to dynamically tune “bond” lengths (i.e., the surface-to-surface distances between adjacent nanoparticles that are linked by DNA) and thermal stabilities of colloidal crystals involving DNA-linked nanoparticles. However, the underlying mechanism of how the DNA bond that links gold nanoparticles changes in this class of colloidal crystals in response to ethanol remains unclear. Here, we conducted a series of all-atom molecular dynamic (MD) simulations to explore the free energy landscape for DNA condensation and decondensation. Our simulations confirm that DNA condensation is energetically much more favorable under 80% ethanol conditions than in pure water, as a result of ethanol’s role in enhancing electrostatic interactions between oppositely charged species. Moreover, the condensed DNA adopts B-form in pure water and A-form in 80% ethanol, which indicates that the higher-order transition does not affect DNA’s conformational preferences. We further propose a nucleosome-like supercoiled model for the DNA condensed state, and we show that the DNA end-to-end distance derived from this model matches the experimentally measured DNA bond length of about 3 nm in the fully condensed state for DNA where the measured length is 16 nm in water. Altogether, this study provides an atomistic understanding of the mechanism underlying ethanol-induced condensation and water-induced decondensation, while our proposed nucleosome-like model allows the design of new strategies for interpreting experimental studies of DNA condensation.

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

The confluence of machine learning and multiscale simulations

Multiscale modeling has a long history of use in structural biology, as computational biologists strive to overcome the time- and length-scale limits of atomistic molecular dynamics. Contemporary machine learning techniques, such as deep learning, have promoted advances in virtually every field of science and engineering and are revitalizing the traditional notions of multiscale modeling. Deep learning has found success in various approaches for distilling information from fine-scale models, such as building surrogate models and guiding the development of coarse-grained potentials. However, perhaps its most powerful use in multiscale modeling is in defining latent spaces that enable efficient exploration of conformational space. In conclusion, this confluence of machine learning and multiscale simulation with modern high-performance computing promises a new era of discovery and innovation in structural biology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Grain boundary segregation and chemical ordering in CoCrFeMnNi multi-principal element alloy

Owing to their far-from-dilute compositions, multi-principal element alloys (MPEAs) can exhibit unique combinations of engineering properties. As nearly all MPEAs are polycrystalline aggregates, it is necessary to understand the interactions of various elemental species with grain boundaries (GBs). This is of particular importance in extreme environments, such as radiation and elevated temperatures, where such interactions have implications on the properties of MPEAs. Herein, we employ atomistic simulations to generate a series of [001] asymmetric tilt GBs in a model CoCrFeMnNi MPEA and quantify solute interactions and segregation to these boundaries. We employ the Warren-Cowley order parameters to investigate the interplay between GB segregation and chemical short-range order (SRO). At temperatures above 800 K, simulation results reveal the segregation of Cr and Mn to CoCrFeMnNi GBs and show weak dependence of boundary solute excess on GB geometry, at least for the boundaries explored in this work. At temperatures in the range of 673–800 K, formation of domains rich in Cr is observed at GBs in agreement with experimental observations. Quantitative analysis shows that solute excess of various alloying elements decreases rapidly with the increase in temperature in the range of 1000–1200 K. Furthermore, we show that GB regions exhibit SRO characteristics that are distinct from the bulk crystals, leading to spatial variations in SRO. In broad terms, our study highlights the need to account for GB interactions with alloying elements when designing advanced MPEAs with novel chemistries.

Wang, Yitao [Lehigh University, Bethlehem, PA (Uni↗
DOE OSTI2021

Ion Pairing Mediates Molecular Organization Across Liquid/Liquid Interfaces

Liquid/liquid interfaces play a central role in scientific fields ranging from nanomaterial synthesis and soft matter electronics to nuclear waste remediation and chemical separations. This diversity of functions arises from an interface’s ability to respond to changing conditions in its neighboring bulk phases. Understanding what drives this interfacial flexibility can provide novel avenues for designing new functional interfaces. However, limiting this progress is an inadequate understanding of the subtle intermolecular and interphase interactions taking place at the molecular level. Here, we use surface-specific vibrational sum frequency generation spectroscopy combined with atomistic molecular dynamics simulations to investigate the self-assembly and structure of model ionic oligomers consisting of an oligodimethylsiloxane (ODMS) tail covalently attached to a positively charged methyl imidazolium (MIM + ) head group at buried oil/aqueous interfaces. We show how the presence of seemingly innocuous salts can impart dramatic changes to the ODMS tail conformations in the oil phase via specific ion effects and ion-pairing interactions taking place in the aqueous phase. These specific ion interactions are shown to drive enhanced amphiphile adsorption, induce morphological changes, and disrupt emergent hydrogen-bonding structures at the interface. Tuning these interactions allows for independent control over the oligomer structure in the oil phase versus interfacial population changes and represents key mechanistic insight that is needed to control chemical reactions at liquid/liquid interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Machine learning approaches for structural and thermodynamic properties of a Lennard-Jones fluid

Predicting the functional properties of many molecular systems relies on understanding how atomistic interactions give rise to macroscale observables. However, current attempts to develop predictive models for the structural and thermodynamic properties of condensed-phase systems often rely on extensive parameter fitting to empirically selected functional forms whose effectiveness is limited to a narrow range of physical conditions. Here, we illustrate how these traditional fitting paradigms can be superseded using machine learning. Specifically, we use the results of molecular dynamics simulations to train machine learning protocols that are able to produce the radial distribution function, pressure, and internal energy of a Lennard-Jones fluid with increased accuracy in comparison to previous theoretical methods. The radial distribution function is determined using a variant of the segmented linear regression with the multivariate function decomposition approach developed by Craven et al. [J. Phys. Chem. Lett. 11, 4372 (2020)]. The pressure and internal energy are determined using expressions containing the learned radial distribution function and also a kernel ridge regression process that is trained directly on thermodynamic properties measured in simulation. The presented results suggest that the structural and thermodynamic properties of fluids may be determined more accurately through machine learning than through human-guided functional forms.

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

Development of Martini 2.2 parameters for N -glycans: a case study of the HIV-1 Env glycoprotein dynamics

N-linked glycans are ubiquitous in nature and play key roles in biology. For example, glycosylation of pathogenic proteins is a common immune evasive mechanism, hampering the development of successful vaccines. Due to their chemical variability and complex dynamics, an accurate molecular understanding of glycans is still limited by the lack of effective resolution of current experimental approaches. Here, we have developed and implemented a reductive model based on the popular Martini 2.2 coarse-grained force field for the computational study of N-glycosylation. We used the HIV-1 Env as a direct applied example of a highly glycosylated protein. Finally, our results indicate that the model not only reproduces many observables in very good agreement with a fully atomistic force field but also can be extended to study large amount of glycosylation variants, a fundamental property that can aid in the development of drugs and vaccines.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Molecular Models of Atomically Dispersed Uranium at MoS 2 Surfaces Reveal Cooperative Mechanism of Water Reduction

Single atoms of uranium supported on molybdenum sulfide surfaces (U@MoS 2 ) have been recently demonstrated to facilitate the hydrogen evolution reaction (HER) through electrocatalysis. Theoretical calculations have predicted uranium hydroxide moieties bound to edge-sulfur atoms of MoS 2 as a proposed transition state involved in the HER process. However, the isolation of relevant intermediates involved in this process remains a challenge, rendering mechanistic hypotheses unverified. The present work describes the isolation and characterization of a uranium-hydroxide intermediate on molybdenum sulfide surfaces using [(Cp* 3 Mo 3 S 4 )UCp*], a molecular model of a reduced uranium center supported at MoS 2 . Mechanistic investigations highlight the metalloligand cooperativity with uranium involved in the water activation pathway. The corresponding uranium-oxo analogue, [(Cp* 3 Mo 3 S 4 )Cp*U(═O)], was also accessed from the hydroxide cluster via hydrogen atom transfer and from [(Cp* 3 Mo 3 S 4 )UCp*] through an alternative direct oxygen atom transfer. These results provide an atomistic perspective on the reactivity of low-valent uranium at molybdenum sulfide surfaces toward water, modeling key intermediates associated with the HER of U@MoS 2 catalysts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Agentic framework for programmatic crystal structure generation using a fine-tuned worker–supervisor large language model

Platinum group metals (PGMs) underpin many catalytic technologies but face severe supply constraints, motivating the search for alternative materials and computational methods to accelerate discovery. While atomistic simulation tools such as Pymatgen and ASE have streamlined structure manipulation, they require detailed inputs, limiting accessibility for experimentalists and slowing early-stage exploration. Here, in this study, we present an AI-driven agentic framework that orchestrates worker–supervisor large language models (LLMs). The worker translates natural-language prompts of varying abstraction into valid crystallographic structures using a compact LLM fine-tuned with low-rank adaptation on a curated text–code–CIF dataset, emphasizing energy-efficient training. Benchmarking against the baseline CodeGen-350M-mono model shows that fine-tuning reduces hallucination rates from 100% to as low as 5% and improves structural match accuracy to up to 82% for fully specified inputs. Accuracy declines with decreasing prompt detail but remains nontrivial even when only stoichiometry and space group are provided, underscoring the LLM’s capacity for crystallographic inference. The supervisor Claude LLM evaluates the outputs and triggers iterative refinement through the worker’s built-in structure manipulation capabilities (e.g., supercell scaling, strain, vacancy, and substitution operations). We further demonstrate use cases for technologically relevant catalysts, including IrO 2 , pyrochlore Pb 2 Ir 2 O 7 , Ni 2 FeO 4 , and Ni 3 Mo, where the framework generates physically consistent structures that can be refined via geometry optimization. This work introduces a low-energy, language-driven pathway for integrating human and machine intelligence in materials design, paving the way for AI-assisted synthesis planning and high-throughput screening of complex oxides.

AI agent↗
DOE OSTI2022

Anomalous thermal transport under high pressure in boron arsenide

High pressure represents extreme environments and provides opportunities for materials discovery. Thermal transport under high hydrostatic pressure has been investigated for more than 100 years and all measurements of crystals so far have indicated a monotonically increasing lattice thermal conductivity. Here we report in situ thermal transport measurements in the newly discovered semiconductor crystal boron arsenide, and observe an anomalous pressure dependence of the thermal conductivity. We use ultrafast optics, Raman spectroscopy and inelastic X-ray scattering measurements to examine the phonon bandstructure evolution of the optical and acoustic branches, as well as thermal conductivity under varied temperatures and pressures up to 32 gigapascals. Using atomistic theory, we attribute the anomalous high-pressure behaviour to competitive heat conduction channels from interactive high-order anharmonicity physics inherent to the unique phonon bandstructure. Our study verifies ab initio theory calculations and we show that the phonon dynamics—resulting from competing three-phonon and four-phonon scattering processes—are beyond those expected from classical models and seen in common materials. In conclusion, this work uses high-pressure spectroscopy combined with atomistic theory as a powerful approach to probe complex phonon physics and provide fundamental insights for understanding microscopic energy transport in materials of extreme properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

In-Situ Ion Irradiation to Add Irradiation Assisted Grain Growth to the Marmot Tool

In the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the MARMOT mesoscale fuel performance tool is used to inform the development of mechanistic materials models for the BISON fuel performance tool. The grain size of the fuel has a large impact on its performance, directly impacting heat conduction, fission gas release, creep, and fracture. Thus, atomistic and mesoscale MARMOT simulations have been used to investigate grain boundary migration and grain growth in UO 2 . However, the current grain growth models in MARMOT does not consider the effect of irradiation on causing grain growth. This research project experimentally studied and implemented irradiation effects on grain growth of UO 2 into MARMOT.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

DeePMD-kit v2: A software package for deep potential models

DeePMD-kit is a powerful open-source software package that facilitates molecular dynamics simulations using machine learning potentials known as Deep Potential (DP) models. This package, which was released in 2017, has been widely used in the fields of physics, chemistry, biology, and material science for studying atomistic systems. The current version of DeePMD-kit offers numerous advanced features, such as DeepPot-SE, attention-based and hybrid descriptors, the ability to fit tensile properties, type embedding, model deviation, DP-range correction, DP long range, graphics processing unit support for customized operators, model compression, non-von Neumann molecular dynamics, and improved usability, including documentation, compiled binary packages, graphical user interfaces, and application programming interfaces. This article presents an overview of the current major version of the DeePMD-kit package, highlighting its features and technical details. Additionally, this article presents a comprehensive procedure for conducting molecular dynamics as a representative application, benchmarks the accuracy and efficiency of different models, and discusses ongoing developments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Pyrite (001) Interface Chemistry is Controlled by a Sulfoxy Termination

Pyrite (FeS 2 ) is the most common sulfide mineral on Earth, forming through inorganic reactions in the crust and oceanic hydrothermal systems and via microbially driven processes in anaerobic sediments. The pyrite–water interface is the site of a wide range of adsorption and reaction processes in Earth systems including oxidation that dramatically affects the geochemistry of surface waters and influences global carbon and oxygen cycles. Mechanistic geochemical models of pyrite interfacial reactivity, however, are limited by the lack of experimentally derived atomistic structures of the reduced and reacting surfaces. Here, in this work, we reveal the atomic-scale structure of the pyrite (001)-water interface that forms at very low oxygen partial pressures, relevant to suboxic environments in Earth. The interface structure and surface speciation were obtained using the crystal truncation rod method supported by ambient-pressure photoelectron spectroscopy and density functional theory calculations. The surface is dominantly composed of disulfide groups bound to a single oxygen atom, forming a sulfoxy group that has no known molecular or bulk mineral analog. This surface is interpreted as the first step in the oxidative dissolution of pyrite. The sulfoxy group is readily protonated through surface acid–base reactions that alter the structure of interfacial water and the free energy of interfacial reactions. Surface iron sites are not oxidized. Surprisingly, this interface can likely develop in equilibrium with bulk pyrite in some reducing and acidic solutions. This termination is therefore likely representative of pyrite surfaces under a vast range of experimental, industrial and Earth conditions.

oxidation↗
DOE OSTI2022

Structural and Dynamic Disorder, Not Ionic Trapping, Controls Charge Transport in Highly Doped Conducting Polymers

Doped organic semiconductors are critical to emerging device applications, including thermoelectrics, bioelectronics, and neuromorphic computing devices. It is commonly assumed that low conductivities in these materials result primarily from charge trapping by the Coulomb potentials of the dopant counter-ions. Here, we present a combined experimental and theoretical study rebutting this belief. Using a newly developed doping technique, we find the conductivity of several classes of high-mobility conjugated polymers to be strongly correlated with paracrystalline disorder but poorly correlated with ionic size, suggesting that Coulomb traps do not limit transport. A general model for interacting electrons in highly doped polymers is proposed and carefully parameterized against atomistic calculations, enabling the calculation of electrical conductivity within the framework of transient localisation theory. Theoretical calculations are in excellent agreement with experimental data, providing insights into the disordered-limited nature of charge transport and suggesting new strategies to further improve conductivities.

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

Atomistic simulations for investigation of substrate effects on lipid in-source fragmentation in secondary ion mass spectrometry

In beam-based ionization methods, the substrate plays an important role on the desorption mechanism of molecules from surfaces. Both the specific orientation that a molecule adopts at a surface and the strength of the molecule-surface interaction can greatly influence desorption processes, which in turn will affect the ion yield and the degree of in-source fragmentation of a molecule. In the beam-based method of secondary ion mass spectrometry (SIMS), in-source fragmentation can be significant and molecule specific due to the hard ionization method of using a primary ion beam for molecule desorption. To investigate the role of the substrate on orientation and in-source fragmentation, we have used atomistic simulations—molecular dynamics in combination with density functional theory calculations—to explore the desorption of a sphingolipid (palmitoylsphingomyelin) from a model surface (gold). We then compare SIMS data from this model system to our modeling findings. Using this approach, we found that the combined adsorption and binding energy of certain bonds associated with the headgroup fragments (C3H8N+, C5H12N+, C5H14NO+, and C5H15PNO4+) was a good predictor for fragment intensities (as indicated by relative ion yields). This is the first example where atomistic simulations have been applied in beam-based ionization of lipids, and it presents a new approach to study biointerfacial lipid ordering effects on SIMS imaging.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Creation of an Fe 3 P Schreibersite Density Functional Tight Binding Model for Astrobiological Simulations

The mineral schreibersite, e.g., Fe 3 P, is commonly found in iron-rich meteorites and could have served as an abiotic phosphorus source for prebiotic chemistry. However, atomistic calculations of its degradation chemistry generally require quantum simulation approaches, which can be too computationally cumbersome to study sufficient time and length scales for this process. In this regard, we have created a computationally efficient semiempirical quantum density functional tight binding (DFTB) model for iron and phosphorus-containing materials by adopting an existing semiautomated workflow that represents many-body interactions by linear combinations of Chebyshev polynomials. Here, we have utilized a relatively small training set to optimize a DFTB model that is accurate for schreibersite physical and chemical properties, including its bulk properties, surface energies, and water absorption. We then show that our model shows strong transferability to several iron phosphide solids as well as multiple allotropes of iron metal. Our resulting DFTB parametrization will allow us to interrogate schreibersite aqueous decomposition at longer time and length scales than standard quantum approaches, providing for more detailed investigations of its role in prebiotic chemistry on early Earth.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mechanistic FCCI model for BISON

In this report, multi-scale simulations to develop a mechanistic model of fuel-cladding chemical interaction (FCCI) conducted under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY23 are described. The model assumes that cladding wastage is dominated by the intermetallic phase (Fe,Cr)$_{17}$Nd$_2$, and calculates production and diffusion of Nd in the fuel, diffusion into the cladding, and formation and growth of the intermetallic wastage layer. Atomistic simulations are used to calculate diffusivity of Nd in the (Fe,Cr) system. Mesoscale simulations, informed by previous atomistic calculations, are used to determine the effective diffusivity of Nd through the fuel, accounting for empty and sodium-filled porosity within the fuel. The BISON model incorporates this data and calculates the wastage layer thickness as a function of time, assuming parabolic growth kinetics based on the best available data. The results of the new mechanistic BISON model are compared to existing assessment cases for the X447 experiment in EBR-II, and show good agreement with experimental data and the existing BISON empirical model. To improve predictions of kinetics for future work, a phase-field model that includes fuel, wastage, and cladding phases is developed. Preliminary testing indicates that the kinetics may differ significantly from the assumed parabolic growth law, suggesting direction for future improvement of the model that will allow it to be applied to other fuel designs with greater confidence.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Nitrate Reduction Modeling under Acidic Conditions with Late Transition Metals

The electrochemical reduction of nitrate (NO 3 R) to ammonia is a bold yet conceivable way of producing ammonia using renewable electricity. However, serious challenges remain in finding optimal electrocatalysts for the process. An atomistic understanding of the surface energetics behind the NO 3 R is needed in order to design an efficient catalyst. Herein, we combine energetics from density functional theory and microkinetic modeling to demonstrate how surface descriptors can help simplify the search for efficient NO 3 R electrocatalysts. We illustrate the strong correlations between transition-state energetics and O* binding energies for adsorbed nitrate and nitrite on transition metals. For intermediates from NO* and beyond, we compare the benefits of using either the N* or H* binding energies to predict reduction onset potentials. These insights enable us to develop a simple microkinetic model that elucidates the surface coverages of intermediates and the product selectivity of NO 3 R across a range of potentials and transition metals. As a result, we show that the model adequately corroborates with quasi-steady-state rates observed experimentally.

ammonia↗
DOE OSTI2021

Advancements in modeling fuel pulverization and cladding behavior during a LOCA

During a loss of coolant accident (LOCA), there is the possibility of nuclear fuel rods to undergo a three step process known as fuel fragmentation, relocation, and dispersal (FFRD). The chance of FFRD occurring increases as the fuel burnup increases. To support the nuclear industry's desire to increase the discharge burnup of nuclear fuels in light water reactors (LWRs) it is imperative to understand the mechanisms driving the evolution of FFRD. In this work, a multiscale modeling modeling approach is used to garner insight into underlying mechanisms leading to the ne fragmentation (also known as pulverization) of nuclear fuel during a LOCA. This report includes a summary of the atomistic and phase-field studies to develop a new pulverization criterion for use in the engineering scale Bison fuel performance code. Details are also provided on cladding modeling improvements related to hydrogen/hydride embrittlement and damage, and anisotropic thermal creep. The new models are used on the existing integral and separate effects LOCA validation cases available in Bison.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗