Engineering PapersSearch

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 55 records · Page 3

Maximizing efficiency of dataset compression for machine learning potentials with information theory

Machine learning interatomic potentials (MLIPs) balance high accuracy and lower costs compared to density functional theory calculations, but their performance often depends on the size and diversity of training datasets. Large datasets improve model accuracy and generalization but are computationally expensive to produce and train on, while smaller datasets risk discarding rare but important atomic environments and compromising MLIP accuracy/reliability. Here, we develop an information-theoretical framework to quantify the efficiency of dataset compression methods and propose an algorithm that maximizes this efficiency. By framing atomistic dataset compression as an instance of the minimum set cover (MSC) problem over atom-centered environments, our method identifies the smallest subset of structures that contains as much information as possible from the original dataset while pruning redundant information. The approach is extensively demonstrated on the GAP-20 and TM23 datasets and validated on 64 varied datasets from the ColabFit repository. Across all cases, MSC consistently retains outliers, preserves dataset diversity, and reproduces the long-tail distributions of forces even at high compression rates, outperforming other subsampling methods. Furthermore, MLIPs trained on MSC-compressed datasets exhibit reduced error for out-of-distribution data even in low-data regimes. We explain these results using an outlier analysis and show that such quantitative conclusions could not be achieved with conventional dimensionality reduction methods. The algorithm is implemented in the open-source QUESTS package and can be used for several tasks in atomistic modeling, from data subsampling, outlier detection, and training improved MLIPs at a lower cost.

36 MATERIALS SCIENCE

Nonadiabatic Force Matching for Alchemical Free-Energy Estimation

We propose a method to compute free-energy differences from nonadiabatic alchemical transformations by using flow-based generative models. The method, nonadiabatic force matching, hinges on estimating the dissipation along an alchemical switching process in terms of a nonadiabatic force field that can be learned through stochastic flow matching. The learned field can be used in conjunction with short-time trajectory data to evaluate upper and lower bounds on the alchemical free energy that variationally converge to the exact value if the field is optimal. Applying the method to evaluate the alchemical free energy of atomistic models shows that it can substantially reduce the simulation cost of a free-energy estimate at a negligible loss of accuracy when compared with thermodynamic integration.

Computational chemistry

Predicting RNA structure and dynamics with deep learning and solution scattering

Advanced deep learning and statistical methods can predict structural models for RNA molecules. However, RNAs are flexible, and it remains difficult to describe their macromolecular conformations in solutions where varying conditions can induce conformational changes. Small-angle x-ray scattering (SAXS) in solution is an efficient technique to validate structural predictions by comparing the experimental SAXS profile with those calculated from predicted structures. There are two main challenges in comparing SAXS profiles to RNA structures: the absence of cations essential for stability and charge neutralization in predicted structures and the inadequacy of a single structure to represent RNA’s conformational plasticity. We introduce a solution conformation predictor for RNA (SCOPER) to address these challenges. This pipeline integrates kinematics-based conformational sampling with the innovative deep learning model, IonNet, designed for predicting Mg 2+ ion binding sites. Validated through benchmarking against 14 experimental data sets, SCOPER significantly improved the quality of SAXS profile fits by including Mg 2+ ions and sampling of conformational plasticity. We observe that an increased content of monovalent and bivalent ions leads to decreased RNA plasticity. Therefore, carefully adjusting the plasticity and ion density is crucial to avoid overfitting experimental SAXS data. SCOPER is an efficient tool for accurately validating the solution state of RNAs given an initial, sufficiently accurate structure and provides the corrected atomistic model, including ions.

59 BASIC BIOLOGICAL SCIENCES

Oscillatory redox behavior in oxides: Cyclic surface reconstruction and reactivity modulation via the Mars–van Krevelen mechanism

The breaking of translational symmetry at oxide surfaces gives rise to coordinatively unsaturated cations/anions and surface restructuring—key factors that govern surface reactivity. Using direct in situ environmental transmission electron microscopy (TEM) observations along with atomistic modeling, we report oscillatory redox behavior in CuO under H 2 , where cyclic surface reconstruction and reactivity modulation occur via the Mars–van Krevelen (MvK) mechanism. We observe self-switching between oxygen-rich and oxygen-deficient surface reconstructions, alternately activating and deactivating the surface for H 2 O formation. During periods of chemical inactivity, the oxygen-deficient surface undergoes slow reoxidation via lattice oxygen diffusing from subsurface and bulk reservoirs, restoring the active oxygen-rich surface termination. The inherent disparity in chemical activity among undercoordinated surface ions, along with sluggish subsurface-to-surface oxygen replenishment, drives this oscillatory redox cycle, modulating H2-induced loss of lattice oxygen at the surface and its delayed replenishment from the subsurface. This creates spatiotemporally separated redox steps at the oxide surface. The phenomena and atomistic insights presented here have significant implications for manipulating the surface reactivity of oxides by tuning the separation of these redox steps.

36 MATERIALS SCIENCE

Dynamic Interfacial Architectures: Cruciferin‐Stabilized Oil/Water Interfaces for Sustainable Emulsions

Stabilizing oil-water interfaces in emulsions by plant-based proteins provides sustainable and tunable ways for designing emulsions with specific properties, for food, healthcare, and pharmaceuticals. Cruciferin, a protein from rapeseed, has great potential as green emulsifier, but details about its structure and mobility at oil-water interfaces are largely unknown. Here, these properties are studied with small angle neutron and x-ray scattering, and neutron spin echo spectroscopy, analyzed by atomistic modelling of scattering curves and coarse-grained modelling, to gain insight into interface coverage, and molecular conformation and mobility at the interface. Cruciferin assumes trimeric conformations at the interface, as in solution, but with its protrusions from the central core of the subunits (“arms”) more compressed. Interfacial mobility is only marginally lower than in solution, indicating the arms still transiently extend and preserve a network, for the first time revealing the mechanism how cruciferin forms highly elastic 2d gel-like oil-water interfaces, as observed in macroscopic rheology. The high interfacial mobility may help in self-repairing non-stabilized interfacial fractions, reducing coalescence. These findings provide a deeper molecular level understanding of proteins at oil-water interfaces, which can stimulate development of new plant-based emulsion products, and contribute to the global protein transition.

cruciferin

Molecular‐Level Insight into the Chlorofluorocarbons Adsorption by Defective Covalent Organic Polymers

Abstract Halocarbons have important industrial applications, but because of their contribution to global warming and the fact that they can cause ozone depletion, they are considered highly toxic. Hence, the techniques that can capture and recover the used halocarbons with energy‐efficient methods have been recently received greater attention. In this contribution, we report the capture of dichlorodifluoromethane (R12), which has high global warming and ozone depletion potential, using covalent organic polymers (COPs). The defect‐engineered COPs were synthesized and demonstrated outstanding sorption capacities, ~226 wt % of R12 combined with linear‐shaped adsorption isotherms. We further identified the plausible microscopic adsorption mechanism of the investigated COPs via grand canonical Monte Carlo simulations applied to non‐defective and a collection of atomistic models of the defective COPs. The modeling work suggests that significant R12 adsorption performance is attributed to a gradual increment of porosities due to isolated/interconnected micro‐/meso‐pore channels and the change of the long‐range ordering of both COPs. The successive hierarchical‐pore‐filling mechanism promotes R12 molecular adsorption via moderate van der Waals adsorbate‐adsorbent interactions in the micropores of both COPs at low pressure followed by adsorbate‐adsorbate interactions in the extra‐voids created at moderate to high pressure ranges. This continuous pore‐filling mechanism makes defective COPs as promising sorbents for halocarbon adsorption.

Shen, Jian

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Molecular origin of viscoelasticity and influence of methylation in mesophase pitch

The viscoelastic and thermomechanical properties of pitches are responsible for their melt-spinning behavior, which is a critical step for manufacturing high-performance pitch-based carbon fibers. Here, we systematically explore the impact of methyl group modifications on the viscoelastic and thermal properties of mesophase pitches. We employ a range of atomistic modeling approaches, including Density Functional Theory (DFT), Density Functional Tight Binding (DFTB), and Classical Molecular Mechanics (MM), to provide detailed insights into the molecular interactions and structural changes. Further, our results revealed the molecular mechanisms that promote layered structures leading to the anisotropic nature of the viscoelastic behavior of mesophase pitch. Furthermore, we propose a modified molecular representation of naphthalene-based mesophase pitch based on the analysis of x-ray diffraction measurements. This study provides fundamental insights into the molecular structures of mesophase pitch and the role of methyl groups controlling its viscosity, which offer valuable insights into mesophase-based carbon fiber production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Systematic discovery of new nano-scale metastable intermetallic eutectic phases in laser rapid solidified Aluminum-Germanium alloy

Laser surface remelting of as-cast Al-Ge eutectic alloy is shown to produce ultrafine lamellar eutectic morphology with interlamellar spacing refined up to ∼60 nm and composed of FCC Al solid solution and unusual Al x Ge y intermetallic phases that do not form during near-equilibrium solidification. The microstructures are characterized and analyzed using a combination of selected area electron diffraction, high-resolution scanning transmission electron microscopy, energy dispersive X-ray spectroscopy to obtain high-resolution elemental maps, and atomistic modeling using density functional theory followed by atomic-scale image simulation. Depending on the local solidification conditions, the crystallography of the Al x Ge y intermetallic phases in the eutectic microstructure is either monoclinic (C 2/c) or monoclinic (P 2 1 ), with high densities of defects in both cases. This is in sharp contrast to the as-cast alloys that showed nominally pure Al and Ge phases with significant solute partitioning and equilibrium FCC and diamond cubic crystal structures, respectively. Corresponding kinetic phase diagrams are proposed to interpret the evolution of nano-lamellar eutectic morphologies with equilibrium Al and metastable Al x Ge y phases, and to explain increased solid solubility in the Al phases manifested by precipitation of ultrafine clusters of Ge. Furthermore, the reasons for the formation of these metastable eutectics under laser rapid solidification are discussed from the perspective of the competitive growth criterion.

Al-Ge eutectic

Changes in dislocation punching behavior due to hydrogen-seeded helium bubble growth in tungsten

The accumulation of gas atoms in tungsten is a topic of long-standing interest to the plasma-facing materials community due the metal's use as a divertor material in some tokamak fusion reactors. The nucleation and growth of He/H gas bubbles (along with their isotopes) can result from impinging fluxes of these gases which give rise to damage at the W divertor surface. The inclusion of He or H in W has been studied extensively by the community, finding that He bubbles modify the surface through periodic dislocation punching and bursting mechanisms while H bubbles impact the metal through plastic-strain induced material failure. However, the mechanisms which are present during the combined flux of both He and H is not well-studied atomistically. Motivated by this, an atomistic modeling study is conducted using molecular dynamics to assess the behavior of mixed concentration He:H bubbles in W. Here we find that the introduction of H into a growing He bubble results in a dramatic change in the nature and presence of dislocation loops which are typically generated via dislocation punching in over-pressurized He bubbles. Most notably, at high H concentrations, there is a switchover in energetic favorability from glissile 1/2<111> dislocations to sessile <100> dislocations. This thermodynamic crossover could imply a significant reduction in W surface morphology changes than with pure He bubbles and, additionally, we show this to have implications on the trapping of H in the bubbles and their associated dislocations.

36 MATERIALS SCIENCE

Atomic-Scale Visualization of Surface Segregation and Ordering of Pt in a Dilute Cu(Pt) Alloy under a Hydrogen Atmosphere

Surface segregation is a common phenomenon in alloys exposed to reactive atmospheres, yet the atomic mechanisms underlying surface structure and composition dynamics remain largely unexplored. Using a combination of environmental transmission electron microscopy observations and atomistic modeling, here we report the surface segregation process of Pt atoms in a dilute Pt(Cu) alloy and determine the distribution of Pt atoms at both atomically flat and stepped surfaces of the Pt(Cu) alloy at elevated temperatures and in a hydrogen gas atmosphere. Through directly probing Pt segregation, we find that Pt atoms segregated on the (100) surface exhibit a p(2 × 2) ordering, with ~25% Pt occupancy. In contrast, on the stepped (410) surface, hydrogen adsorption induces Pt segregation, initially occurring at the step edges, which then expands to the terrace sites upon increased hydrogen coverage, resulting in an ordered distribution of segregated Pt atoms with ~22% occupancy. Finally, these observations offer mechanistic insights into the structure and composition dynamics of the topmost atomic layer of the alloy in response to environmental stimuli and hold practical implications for the design and optimization of catalysts based on Pt group metals.

36 MATERIALS SCIENCE

Scale Invariance of Hot Spot Formation in TATB High Explosives

Shock-induced detonation of insensitive high explosives based on 1,3,5-triamino-2,4,6-trinitrobenzene starts with formation of hot spots at microstructural defects but has eluded atomistic modeling treatment at micron length scales. To this end, we performed multimicron scale all-atom molecular dynamics (MD) simulations of hot spots that form during the collapse of cylindrical pores with diameters between 10 and 300 nm. Our MD simulations show that hot spots formed at pores larger than 20 nm exhibit temperature fields with scale-invariant features for sizes up to at least 300 nm. Through a continuum-based grain-scale modeling framework, we span and extend beyond the size scales currently accessible to MD and find that hot spot scale invariance is a general feature that arises when the mechanical strength is insensitive to strain rate. Finally, our results demonstrate the applicability of all-atom MD to simulate the complicated dynamical evolution of micron-sized systems and bolster confidence in insights from MD simulations of materials that exhibit strength with negligible rate dependence over the relevant intervals.

36 MATERIALS SCIENCE

Surface Self-Diffusion Induced Sintering of Nanoparticles

Despite the critical role of sintering phenomena in constraining the long-term durability of nanosized particles, a clear understanding of nanoparticle sintering has remained elusive due to the challenges in atomically tracking the neck initiation and discerning different mechanisms. Through the integration of in situ transmission electron microscopy and atomistic modeling, this study uncovers the atomic dynamics governing the neck initiation of Pt–Fe nanoparticles via a surface self-diffusion process, allowing for coalescence without significant particle movement. Real-time imaging reveals that thermally activated surface morphology changes in individual nanoparticles induce significant surface self-diffusion. Further, the kinetic entrapment of self-diffusing atoms in the gaps between closely spaced nanoparticles leads to the nucleation and growth of atomic layers for neck formation. This surface self-diffusion-driven sintering process is activated at a relatively lower temperature compared to the classic Ostwald ripening and particle migration and coalescence processes. The fundamental insights have practical implications for manipulating the morphology, size distribution, and stability of nanostructures by leveraging surface self-diffusion processes.

36 MATERIALS SCIENCE

Atomic-Scale Dynamic Mechanisms of Embedded MoS 2 Wires

Nanowires composed of a 1:1 stoichiometry of transition metals and chalcogen ions can be fabricated from two-dimensional transition metal dichalcogenides (TMDs) by using electron beam irradiation. Wires fabricated through in situ experiments can be geometrically connected to TMD sheets in various ways, and their physical properties can vary accordingly. Understanding the structural transformation caused by electron beams is critical for designing wire-sheet structures for nanoelectronics. In this study, we report the behavior of nanowires formed inside a monolayer MoS 2 sheet by combining phase-contrast images and large-scale atomistic modeling. Here we investigate the effect of vacancies on the dynamic evolution of wires, such as rotations with different edge structures and breaking, by considering the interactions between MoS wires and MoS 2 nanosheets. The obtained insights can be applied to other monolayer TMDs to guide the behavior of TMD wires and fabricate favorable geometries for various applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

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

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

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