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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Active learning of reactive Bayesian force fields applied to heterogeneous catalysis dynamics of H/Pt

Abstract Atomistic modeling of chemically reactive systems has so far relied on either expensive ab initio methods or bond-order force fields requiring arduous parametrization. Here, we describe a Bayesian active learning framework for autonomous “on-the-fly” training of fast and accurate reactive many-body force fields during molecular dynamics simulations. At each time-step, predictive uncertainties of a sparse Gaussian process are evaluated to automatically determine whether additional ab initio training data are needed. We introduce a general method for mapping trained kernel models onto equivalent polynomial models whose prediction cost is much lower and independent of the training set size. As a demonstration, we perform direct two-phase simulations of heterogeneous H 2 turnover on the Pt(111) catalyst surface at chemical accuracy. The model trains itself in three days and performs at twice the speed of a ReaxFF model, while maintaining much higher fidelity to DFT and excellent agreement with experiment.

42 ENGINEERING↗
DOE OSTI2023

Incoherent nonadiabatic to coherent adiabatic transition of electron transfer in colloidal quantum dot molecules

Electron transfer is a fundamental process in chemistry, biology, and physics. One of the most intriguing questions concerns the realization of the transitions between nonadiabatic and adiabatic regimes of electron transfer. Using colloidal quantum dot molecules, we computationally demonstrate how the hybridization energy (electronic coupling) can be tuned by changing the neck dimensions and/or the quantum dot sizes. This provides a handle to tune the electron transfer from the incoherent nonadiabatic regime to the coherent adiabatic regime in a single system. We develop an atomistic model to account for several states and couplings to the lattice vibrations and utilize the mean-field mixed quantum-classical method to describe the charge transfer dynamics. Here, we show that charge transfer rates increase by several orders of magnitude as the system is driven to the coherent, adiabatic limit, even at elevated temperatures, and delineate the inter-dot and torsional acoustic modes that couple most strongly to the charge transfer dynamics.

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

Strain release by 3D atomic misfit in fivefold twinned icosahedral nanoparticles with amorphization and dislocations

Multiple twinning to form fivefold twinned nanoparticles in crystal growth is common and has attracted broad attention ranging from crystallography research to physical chemistry and materials science. Lattice-misfit strain and defects in multiple twinned nanoparticles (MTP) are key to understand and tailor their electronic properties. However, the structural defects and related strain distributions in MTPs are poorly understood in three dimensions (3D). Here, we show the 3D atomic misfit and strain relief mechanism in fivefold twinned icosahedral nanoparticles with amorphization and dislocations by using atomic resolution electron tomography. We discover a two-sided heterogeneity in variety of structural characteristics. A nearly ideal crystallographic fivefold face is always found opposite to a less ordered face, forming Janus-like icosahedral nanoparticles with two distinct hemispheres. The disordered amorphous domains release a large amount of strain. Molecular dynamics simulations further reveal the Janus-like icosahedral nanoparticles are prevalent in the MTPs formed in liquid-solid phase transition. This work provides insights on the atomistic models for the modelling of formation mechanisms of fivefold twinned structures and computational simulations of lattice distortions and defects. We anticipate it will inspire future studies on fundamental problems such as twin boundary migration and kinetics of structures in 3D at atomic level.

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

Direction-specific enhanced diffusion of CO 2 in chiral hexagonal boron nitride nanotubes

To meet performance requirements, the next generation of gas separation membranes will need both high gas permeability and selectivity, attainable if we could coax adsorbates to minimize random Brownian motion and produce direction-specific diffusion along a desired axis. In this atomistic modeling study, we detail how direction-specific diffusion of CO 2 can be achieved in chiral hexagonal boron nitride nanotubes (hBNNTs) by means of a non-Knudsen diffusion mechanism. Our findings detail how this mechanism of diffusion is driven by interactions with the tube walls and enables the CO 2 molecules to diffuse along the nanotube’s z-axis with minimized collisions and directional changes. hBNNTs with chiral indices exhibit CO 2 diffusion rates faster than non-chiral tubes of comparable and larger diameters. Of the hBNNTs studied, a (7,3) tube appears to be ideally sized (3.7 Å radius) exhibiting CO 2 diffusion that is 3.4 times faster than diatomic N 2 . Applying this mechanism of diffusion to hypothetical sheet membranes prepared with aligned chiral (7,3) hBNNTs results in membranes with a calculated CO 2 /N 2 permselectivity of 170 and a CO 2 permeability limit of nearly 1.35 ×10 7 Barrer, readily surpassing the Robeson upper bound for CO 2 /N 2 separations.

CO2↗
DOE OSTI2025

Machine learned potential for high-throughput phonon calculations of metal—organic frameworks

Metal–organic frameworks (MOFs) are highly porous and versatile materials studied extensively for applications such as carbon capture and water harvesting. However, computing phonon-mediated properties in MOFs, like thermal expansion and mechanical stability, remains challenging due to the large number of atoms per unit cell, making traditional Density Functional Theory (DFT) methods impractical for high-throughput screening. Recent advances in machine learning potentials have led to foundation atomistic models, such as MACE-MP-0, that accurately predict equilibrium structures but struggle with phonon properties of MOFs. In this work, we developed a workflow for computing phonons in MOFs within the quasi-harmonic approximation with a fine-tuned MACE model, MACE-MP-MOF0. The model was trained on a curated dataset of 127 representative and diverse MOFs. The fine-tuned MACE-MP-MOF0 improves the accuracy of phonon density of states and corrects the imaginary phonon modes of MACE-MP-0, enabling high-throughput phonon calculations with state-of-the-art precision. The model successfully predicts thermal expansion and bulk moduli in agreement with DFT and experimental data for several well-known MOFs. These results highlight the potential of MACE-MP-MOF0 in guiding MOF design for applications in energy storage and thermoelectrics.

Elena, Alin Marin↗
DOE OSTI2023

Geometrically frustrated interactions drive structural complexity in amorphous calcium carbonate

Amorphous calcium carbonate is an important precursor for biomineralization in marine organisms. Key outstanding problems include understanding the structure of amorphous calcium carbonate and rationalizing its metastability as an amorphous phase. Here we report high-quality atomistic models of amorphous calcium carbonate generated using state-of-the-art interatomic potentials to help guide fits to X-ray total scattering data. Exploiting a recently developed inversion approach, we extract from these models the effective Ca∙∙∙Ca interaction potential governing the structure. This potential contains minima at two competing distances, corresponding to the two different ways that carbonate ions bridge Ca 2+ -ion pairs. We reveal an unexpected mapping to the Lennard-Jones–Gauss model normally studied in the context of computational soft matter. The empirical model parameters for amorphous calcium carbonate take values known to promote structural complexity. We thus show that both the complex structure and its resilience to crystallization are actually encoded in the geometrically frustrated effective interactions between Ca 2+ ions.

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

The quest for safer nuclear fuels

Here, to increase the safety and performance of nuclear reactors, researchers are working on developing ‘accident tolerant’ fuels. Ideally, these fuel pellets will be less susceptible to swelling and be more resistant to extreme heat. This led chemists to develop uranium dioxide (UO 2 ) doped with divalent elements, E 2+ . These dopants can introduce vacancies into the lattice that lead to larger grain sizes, which create longer pathways for fission gasses to migrate out of the grains. In undoped fuel, fission gas build-up at grain boundaries causes the fuel pellet to swell. Therefore, having larger grains that slow fission gas migration into grain boundaries is highly desired. Divalent Cr 2+ ions are a commonly explored dopant for UO 2 fuels. These materials have larger grain sizes while maintaining the dissolution kinetics of undoped UO 2 , making their long-term storage and disposal paths similarly safe. Atomistic models have suggested that Mn 2+ is a conceivable alternative doping material, with the potential to increase the sizes of grains even more than Cr 2+ dopants. As such, experimental studies are required to understand the physical chemistry and reactivity of these potential fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Origins of the change in mechanical strength of silicon/gold nanocomposites during irradiation

Abstract Silicon-based layered nanocomposites, comprised of covalent-metal interfaces, have demonstrated elevated resistance to radiation. The amorphization of the crystalline silicon sublayer during irradiation and/or heating can provide an additional mechanism for accommodating irradiation-induced defects. In this study, we investigated the mechanical strength of irradiated Si-based nanocomposites using atomistic modeling. We first examined dose effects on the defect evolution mechanisms near silicon-gold crystalline and amorphous interfaces. Our simulations reveal the growth of an emergent amorphous interfacial layer with increasing dose, a dominant factor mitigating radiation damage. We then examined the effect of radiation on the mechanical strength of silicon-gold multilayers by constructing yield surfaces. These results demonstrate a rapid onset strength loss with dose. Nearly identical behavior is observed in bulk gold, a phenomenon that can be rooted to the formation of radiation-induced stacking fault tetrahedra which dominate the dislocation emission mechanism during mechanical loading. Taken together, these results advance our understanding of the interaction between radiation-induced point defects and metal-covalent interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structural water in amorphous carbonate minerals: ab initio molecular dynamics simulations of X-ray pair distribution experiments

Water is known to play a controlling role in directing mineralization pathways and stabilizing metastable amorphous intermediates in hydrous carbonate mineral MCO 3 ·nH 2 O systems, where M 2+ is a divalent metal cation. Despite this recognition, the nature of the controls on crystallization is poorly understood, largely owing to the difficulty in characterizing the dynamically disordered structures of amorphous intermediates at the atomic scale. Here, we present a series of atomistic models, derived from ab initio molecular dynamics simulation, across a range of experimentally relevant cations (M = Ca, Mg, Sr) and hydration levels (0 ≤ n ≤ 2). Theoretical simulations of the dependence of the X-ray pair distribution function on the hydration level n show good agreement with available experimental data and thus provide further evidence for a lack of significant nanoscale structure in amorphous carbonates. Upon dehydration, the metal coordination number does not change significantly, but the relative extent of water dissociation increases, indicating that a thermodynamic driving force exists for water dissociation to accompany dehydration. Mg strongly favors monodentate conformation of carbonate ligands and shows a marked preference to exchange monodentate carbonate O for water O upon hydration, whereas Ca and Sr exchange mono- and bidentate carbonate ligands with comparable frequency. Finally, water forms an extensive hydrogen bond network among both water and carbonate groups that exhibits frequent proton transfers for all three cations considered suggesting that proton mobility is likely predominantly due to water dissociation and proton transfer reactions rather than molecular water diffusion.

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

Size-dependent shape distributions of platinum nanoparticles

Transmission electron microscopy revealed size-dependent shape distributions in platinum nanoparticles, which were consistent with trends observed by applying Boltzmann statistics to the energy computed with atomistic models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Tuning martensitic transformations via coherent second phases in nanolaminates using free energy landscape engineering

In this work we explore the possibilities and limitations of using a coherent second phase to engineer the thermo-mechanical properties of a martensitic alloy by modifying the underlying free energy landscape that controls the transformation. We use molecular dynamics simulations of a model atomistic system where the properties of a coherent, nanoscale second phase can be varied systematically. With a base martensitic material that undergoes a temperature-induced transformation from a cubic austenite to a monoclinic martensite, simulations show significant ability to engineer the transformation temperatures, from a ~50% reduction to a ~200% increase, with 50 at.% of the cubic second phase. We establish correlations between the properties of the second phase, the transformation characteristics, and the microstructure via the free energy landscape of the twophase systems. Coherency stresses have a strong influence on the martensitic variants observed and can even cause the non-martensitic second phase to undergo a transformation. Reducing the stiffness of the second phase increases the transformation strain and modifies the martensitic microstructure, increasing the volume fraction of the transformed material. This increase in transformation strain is accompanied by a significant increase in Af and thermal hysteresis, while Ms remains unaltered. Our findings on the tunability of martensitic transformations can be used for informed searches of second phases to achieve desired material properties, such as achieving room temperature, lightweight shape memory alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simulations of nonradiative processes in semiconductor nanocrystals

The description of carrier dynamics in spatially confined semiconductor nanocrystals (NCs), which have enhanced electron–hole and exciton–phonon interactions, is a great challenge for modern computational science. These NCs typically contain thousands of atoms and tens of thousands of valence electrons with discrete spectra at low excitation energies, similar to atoms and molecules, that converge to the continuum bulk limit at higher energies. Computational methods developed for molecules are limited to very small nanoclusters, and methods for bulk systems with periodic boundary conditions are not suitable due to the lack of translational symmetry in NCs. This perspective focuses on our recent efforts in developing a unified atomistic model based on the semiempirical pseudopotential approach, which is parameterized by first-principle calculations and validated against experimental measurements, to describe two of the main nonradiative relaxation processes of quantum confined excitons: exciton cooling and Auger recombination. Here, we focus on the description of both electron–hole and exciton–phonon interactions in our approach and discuss the role of size, shape, and interfacing on the electronic properties and dynamics for II–VI and III–V semiconductor NCs.

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

Chemical order transitions within extended interfacial segregation zones in NbMoTaW

Interfacial segregation and chemical short-range ordering influence the behavior of grain boundaries in complex concentrated alloys. In this study, we use atomistic modeling of a NbMoTaW refractory complex concentrated alloy to provide insight into the interplay between these two phenomena. Hybrid Monte Carlo and molecular dynamics simulations are performed on columnar grain models to identify equilibrium grain boundary structures. Here our results reveal extended near-boundary segregation zones that are much larger than traditional segregation regions, which also exhibit chemical patterning that bridges the interfacial and grain interior regions. Furthermore, structural transitions pertaining to an A2-to-B2 transformation are observed within these extended segregation zones. Both grain size and temperature are found to significantly alter the widths of these regions. An analysis of chemical short-range order indicates that not all pairwise elemental interactions are affected by the presence of a grain boundary equally, as only a subset of elemental clustering types are more likely to reside near certain boundaries. The results emphasize the increased chemical complexity that is associated with near-boundary segregation zones and demonstrate the unique nature of interfacial segregation in complex concentrated alloys.

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

Inherent-state melting and the onset of glassy dynamics in two-dimensional supercooled liquids

Below the onset temperature T o , the equilibrium relaxation time of most glass-forming liquids exhibits glassy dynamics characterized by a super-Arrhenius temperature dependence. In this supercooled regime, the relaxation dynamics also proceeds through localized elastic excitations corresponding to hopping events between inherent states, i.e., potential-energy-minimizing configurations of the liquid. Despite its importance in distinguishing the supercooled regime from the high-temperature regime, the microscopic origin of T o is not yet known. Here, we construct a theory for the onset temperature in two dimensions and find that an inherent-state melting transition, described by the binding–unbinding transition of dipolar elastic excitations, delineates the supercooled regime from the high-temperature regime. The corresponding melting transition temperature is in good agreement with the onset temperature found in various two-dimensional (2D) atomistic models of glass formers and an experimental binary colloidal system confined to a water–air interface. Additionally, we find the predictions for the renormalized elastic moduli to agree with the experimentally observed values for the latter 2D colloidal system. We further discuss the predictions of our theory on the displacement and density correlations at supercooled conditions, which are consistent with observations of the Mermin–Wagner fluctuations in experiments and molecular simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Tuning the surface reactivity of oxides by peroxide species

The Mars–van Krevelen mechanism is the foundation for oxide-catalyzed oxidation reactions and relies on spatiotemporally separated redox steps. Herein, we demonstrate the tunability of this separation with peroxide species formed by excessively adsorbed oxygen, thereby modifying the catalytic activity and selectivity of the oxide. Using CuO as an example, we show that a surface layer of peroxide species acts as a promotor to significantly enhance CuO reducibility in favor of H 2 oxidation but conversely as an inhibitor to suppress CuO reduction against CO oxidation. Together with atomistic modeling, we identify that this opposite effect of the peroxide on the two oxidation reactions stems from its modification on coordinately unsaturated sites of the oxide surface. By differentiating the chemical functionality between lattice oxygen and peroxide, these results are closely relevant to a wide range of catalytic oxidation reactions using excessively adsorbed oxygen to activate lattice oxygen and tune the activity and selectivity of redox sites.

36 MATERIALS SCIENCE↗