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 271 records · Page 15

Ion correlations explain kinetic selectivity in diffusion-limited solid-state synthesis reactions

Establishing viable solid-state synthesis pathways for novel inorganic materials remains a major challenge in materials science. Previous pathway design methods using pairwise reaction approaches have navigated the thermodynamic landscape with first-principles data but lack kinetic information, limiting their effectiveness. This gap leads to suboptimal precursor selection and predictions, especially for reactions forming competing phases with similar formation energies, where ion diffusion is a critical influence. Here we demonstrate an inorganic synthesis framework by incorporating machine learning-derived transport properties through ‘liquid-like’ product layers into a thermodynamic cellular reaction model. In the Ba–Ti–O system, known for its competitive polymorphism, we obtain accurate predictions of phase formation with varying BaO:TiO2 ratios as a function of time and temperature. We find that diffusion–thermodynamics interplay governs phase compositions, with cross-ion transport coefficients critical for predicting diffusion-limited selectivity. This work bridges length scales and timescales by integrating solid-state reaction kinetics with first-principles thermodynamics and spatial reactivity.

Atomistic models↗

Complexions at the iron-magnetite interface

Synthesizing distinct phases and controlling crystalline defects are key concepts in materials design. These approaches are often decoupled, with the former grounded in equilibrium thermodynamics and the latter in nonequilibrium kinetics. By unifying them through defect phase diagrams, we can apply phase equilibrium models to thermodynamically evaluate defects—including dislocations, grain boundaries, and phase boundaries—establishing a theoretical framework linking material imperfections to properties. Using scanning transmission electron microscopy (STEM) with differential phase contrast (DPC) imaging, we achieve the simultaneous imaging of heavy Fe and light O atoms, precisely mapping the atomic structure and chemical composition at the iron-magnetite (Fe/Fe 3 O 4 ) interface. We identify a well-ordered two-layer interface-stabilized phase state (referred to as complexion) at the Fe[001]/Fe 3 O 4 [001] interface. Using density-functional theory (DFT), we explain the observed complexion and map out various interface-stabilized phases as a function of the O chemical potential. The formation of complexions increases interface adhesion by 20% and alters charge transfer between adjacent materials, impacting transport properties. Our findings highlight the potential of tunable defect-stabilized phase states as a degree of freedom in materials design, enabling optimized corrosion protection, catalysis, and redox-driven phase transitions, with applications in materials sustainability, efficient energy conversion, and green steel production.

36 MATERIALS SCIENCE↗

Computational study of the effect of core–skin structure on the mechanical properties of carbon nanofibers

The effect of the core–skin structure on the mechanical properties of carbon nanofibers is investigated in large-scale molecular dynamics simulations of tensile deformation of carbon nanofibers with the core–skin and homogeneous structures. Contrary to an established notion of the deleterious effect of the skin layer on the strength of carbon fibers, the presence of a high-quality skin layer is found to increase both the Young’s modulus and tensile strength of the nanofiber. A detailed analysis of the fracture process indicates that the nanofiber strengthening is related to the ability of skin layer to suppress crack nucleation at the core–skin interface. Furthermore, the computational predictions suggest that the design of new approaches to carbon fiber manufacturing and processing leading to the generation of a high-quality skin layer while avoiding the introduction of structural defects at the core–skin interface may yield a significant enhancement of the mechanical properties of carbon fibers.

36 MATERIALS SCIENCE↗

A Review: Microstructural and Phase Evolution in Alloys during Extended Plastic Deformation

Materials in service and during processing are often subjected to plastic deformation. For multi-phase metallic alloys, simple geometric models and atomistic simulations suggest that two distinctive regimes in these materials’ evolution during deformation exist. At low strains, evolutions are often dominated by the kinetic roughening of interfaces, which results from the superdiffusive transport of matter in sheared crystals. Furthermore, at high strains and temperatures where thermal diffusion is sluggish, on the other hand, shearing- induced forced atomic mixing dominates these evolutions, resulting in significantly enhanced solubility. Distinguishing these two regimes is shown to provide a convenient framework for rationalizing and analyzing recent experiments and simulations on wear of layered structures in the low strain-regime and nonequilibrium phase co-existence and self-organization in highly immiscible or reactive alloy systems in the high-strain regime.

36 MATERIALS SCIENCE↗

Exploring NaCl-PuCl 3 molten salts with machine learning interatomic potentials and graph theory

Actinide molten salts are the basis of the liquid fuels used in molten salt reactors. Due to the inherent difficulties associated with high temperature and hazardous conditions, experimental investigations of fundamental properties of these materials are usually challenging. In this work, we describe the structure and transport of NaCl-PuCl 3 mixtures using computational techniques. Three compositions were considered (16, 25, and 36 mol% PuCl 3 ) over a temperature range (730 – 1257K) using ab initio molecular dynamics, which provided the necessary data sets for training machine learned interatomic potentials. Further, molecular dynamics simulations based on these potentials were then used to determine structure and transport properties. A substantial change was noted in the structure factor when increasing the PuCl 3 content from 25 to 36 mol%. This change is linked to the aggregation of larger Pu 3+ clusters. In addition, the similarity of the atomic environments of metal cations in molten salt systems to their solid states counterparts was investigated using an unsupervised learning technique. Finally, graph theory was employed to explore the structure and size of actinide networks. Consistent with the structure factor, a dense Pu 3+ intermolecular structure is observed within the 36 mol% PuCl 3 mixture. The structure of cation-cation inter-junctions is also discussed. In all cases, the diffusion of Pu 3+ is significantly lower than that of Na + and Cl - .

36 MATERIALS SCIENCE↗

Probing degradation at solid-state battery interfaces using machine-learning interatomic potential

Solid-state batteries featuring fast ion-conducting solid electrolytes are promising next-generation energy storage technologies, yet challenges remain for practical deployment due to electro-chemo-mechanical instabilities at solid-solid interfaces. These interfaces, which include homogeneous/internal interfaces such as grain boundaries (GBs) and heterogeneous/external interfaces between solid-electrolyte and electrode materials, can impede Li-ion transport, deteriorate performance, and eventually lead to cell failure. Here, in this study, we leverage large-scale molecular simulations, enabled by validated machine-learning interatomic potentials, to directly probe the onset of interfacial degradation at the garnet Li 7 La 3 Zr 2 O 12 (LLZO) solid-electrolyte/LiCoO 2 (LCO) cathode interface. By surveying different interfacial geometries and compositions, it is found that Li-deficient interfaces can lead to severe interfacial disordering with cation mixing and Co interdiffusion from LCO into LLZO. By contrast, Li-sufficient interfaces are less disordered, although elemental segregation with local ordering is observed. As a consequence of Co interdiffusion, Co-rich regions are formed at the GBs of LLZO due to cation segregation and trapping effects. This behavior is independent of the GB tilting axis, degree of disorder at the GBs, and Co concentration, which implies Co clustering at GBs is a general phenomenon in polycrystalline LLZO and can dictate its overall transport and mechanical properties. Our findings elucidate the underlying fundamental mechanisms that give rise to experimentally observed physicochemical properties and provide guidelines for interface design that can mitigate interfacial degradation and improve cycling performance.

25 ENERGY STORAGE↗

Influence of Ag metal dispersion on the thermal conversion of ethanol to butadiene over Ag-ZrO 2 /SiO 2 catalysts

Atomistic scale models were developed and coupled with experimental investigation to deliver a functional understanding of catalytic activity and selectivity in the conversion of ethanol to 1,3-butadiene over Ag/ZrO 2 /SiO 2 . A detailed evaluation of the structural and electronic properties of the resultant catalyst models led to the identification of critical active sites of the catalyst. More importantly, the extent of Ag dispersion on the SiO 2 support and relative proximity to ZrO 2 were found to vary with the oxidation state of Ag and local coordination environment (Ag-O SiO2 ) allowing for critical control of ethanol conversion towards butadiene or ethylene. Simulations revealed that less dispersed or clustered Ag contain predominantly Ag 0 charge state and promote conversion of ethanol to ethylene. The well-dispersed Ag/ZrO 2 /SiO 2 catalyst instead contain a larger fraction of cationic Ag δ+ , and predominantly promote ethanol dehydrogenation and subsequent production of butadiene. The theoretical insights drawn were validated and confirmed experimentally using TEM, XRD and reactivity measurements demonstrating the effect of Ag dispersion on the selectivity of ethanol conversion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomistic determination of Peierls barriers of dislocation glide in nickel

The Peierls barrier measures the lattice resistance to dislocation glide in crystalline solids. We use the nudged elastic band (NEB) method to calculate the Peierls barriers for screw and edge dislocation glide in a face-centered cubic (FCC) metal of Ni. The minimum energy paths (MEPs) across single or sequential Peierls barriers are determined under shear loading. The NEB results show the decreasing Peierls barrier with increasing shear stress, giving the Peierls stress at which the Peierls barrier vanishes. The effects of boundary condition and system size on Peierls barriers are studied by comparing strain- and stress-controlled NEB results. Furthermore, the free-end NEB methods are applied to determine MEPs with improved computational efficiency. The NEB results are also used to evaluate the energetic driving force of dislocation glide, which is consistent with that determined from the Peach-Koehler force. The accuracy of the present NEB results based on an empirical interatomic potential is assessed by comparison with a machine-learning potential. This work demonstrates the robust and efficient quantification of Peierls barriers to dislocation glide in an FCC metal, and it lays a solid foundation for the atomistic determination of Peierls barriers in compositionally complex alloys with the FCC structure in future studies.

42 ENGINEERING↗

Defect engineering of porous aromatic frameworks via end capping improves dioxane removal from water

Amorphous porous organic polymers show promise for energy-efficient adsorptive separations, but it is difficult to understand or improve their performance through intentional structural modification. Herein, we report the synthesis of porous aromatic frameworks (PAFs) with pore structures tailored by the incorporation of a monofunctionalized end-capping monomer that disrupts framework topology. Combining experimental characterization with molecular simulations, we show that this defect engineering strategy yields less densely crosslinked networks, which leads to pore collapse and the presence of unique adsorption sites. These defect-engineered PAFs exhibit enhanced removal of 1,4-dioxane, an important environmental pollutant, from water. Increasing the concentration of end-capping monomers produces PAFs with narrower pore size distributions and improved 1,4-dioxane uptake. These results illustrate that defect engineering can effectively modulate polymer connectivity and porosity for applications in selective adsorptive separations. In conclusion, this technique avoids post-synthetic treatments and presents another approach to tailor amorphous polymeric adsorbents.

36 MATERIALS SCIENCE↗

Local atomic environment analysis of short and long-range solute-solute interactions in a symmetric tilt grain boundary

We develop an efficient molecular-statics algorithm for predicting local solute-solute interaction effects on grain-boundary segregation. This algorithm is applied to the segregation and interaction of Ni and Cu solutes in a Σ9 (221) symmetric tilt grain boundary in Ag. It is found that the strength of interaction between two solutes directly correlates with the change of average segregation energy in the local atomic environment of a fixed grain-boundary solute. It is shown that local interaction effects at grain boundaries can be either attractive or repulsive wherein long-range attractive effects increase the propensity for solute cluster nucleation. This study underscores the atomic nature of grain-boundary solute segregation and provides a tool that could yield meaningful insight into heterogeneous solute segregation phenomena in polycrystalline alloys beyond their dilute limit.

36 MATERIALS SCIENCE↗

Crystal structure and shape selection in the growth of 3D metallic crystallites on layered materials: Fe on MoS 2

Nucleation and growth of supported 3D metal clusters or crystallites during deposition on MoS 2 , or on other weakly-adhering layered materials, can potentially produce diverse growth shapes, and even crystal structures differing from the bulk metal. For Fe deposition on MoS 2 , SEM and AFM observations reveal three distinct crystallite shapes. By comparison with atomistic structure models incorporating realistic Fe-MoS 2 interface structures, here we conclude that these are: triangular fcc(111) pyramids with sloped {100} side facets; bcc(110) A-frame tents with sloped {100} side facets; and bcc(110) mesas with vertical {100} and {110} side facets. The following picture is proposed for the competitive formation of clusters and crystallites with different structures: (i) small nanoclusters formed at the onset of deposition exhibit facile fluxional dynamics allowing sampling of different crystal structures and shapes; (ii) sufficient fluxionality implies a Boltzmann distribution of sampled structures, and thus coexistence of different structures follows from the demonstrated similar energies for those structures; (iii) growing clusters reach a threshold size above which the characteristic time scale for restructuring exceeds that for cluster growth. Thereafter, clusters are locked-in to a specific crystal structure and shape as revealed by imaging of larger crystallites. Despite a penalty for fcc(111) over bcc(111) pyramids based on bulk energetics, favorable surface and interface energies makes them preferable for smaller sizes.

36 MATERIALS SCIENCE↗

Computational Methods in Heterogeneous Catalysis

The unprecedented ability of computations to probe atomic-level details of catalytic systems holds immense promise for the fundamentals-based bottom-up design of novel heterogeneous catalysts, which are at the heart of the chemical and energy sectors of industry. Here, we critically analyze recent advances in computational heterogeneous catalysis. First, we will survey the progress in electronic structure methods and atomistic catalyst models employed, which have enabled the catalysis community to build increasingly intricate, realistic, and accurate models of the active sites of supported transition-metal catalysts. We then review developments in microkinetic modeling, specifically mean-field microkinetic models and kinetic Monte Carlo simulations, which bridge the gap between nanoscale computational insights and macroscale experimental kinetics data with increasing fidelity. Here, we finally review the advancements in theoretical methods for accelerating catalyst design and discovery. Throughout the review, we provide ample examples of applications, discuss remaining challenges, and provide our outlook for the near future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advances in Clayff Molecular Simulation of Layered and Nanoporous Materials and Their Aqueous Interfaces

As a general-purpose force field for molecular simulations of layered materials and their fluid interfaces, Clayff continues to see broad usage in atomistic computational modeling for numerous geoscience and materials science applications due to its (1) success in predicting properties of bulk nanoporous materials and their interfaces, (2) transferability to a range of layered and nanoporous materials, and (3) simple functional form which facilitates incorporation into a variety of simulation codes. Here, we review applications of Clayff to model bulk phases and interfaces not included in the original parameter set and recent modifications for modeling surface terminations such as hydroxylated nanoparticle edges. Overall, we conclude with a discussion of expectations for future developments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal Shape Stability of fcc Metal Nanocrystals Synthesized with Faceted Nonequilibrium Shapes

Highly refined capabilities of the shape-controlled solution-phase synthesis of metal nanocrystals (NCs) allow the generation of NCs with faceted nonequilibrium shapes, which optimize properties for target applications such as catalysis and plasmonics. Often, for such applications and also for TEM analysis, the NCs are removed from the solution-phase environment. We explore the postsynthesis evolution of these metastable NCs in a high-vacuum TEM environment. Specifically, here we analyze their reshaping toward the equilibrium Wulff shapes mediated by surface diffusion, where such reshaping degrades the above-mentioned optimized properties. Typical sizes for these NCs range from 5 to 30 nm or 10 3 –10 6 atoms, and reshaping often occurs on the time scale of minutes for temperatures around, say, 400 °C. We discuss the development of predictive stochastic atomistic-level models for NC evolution with a realistic description of surface diffusion. These models, in contrast to Molecular Dynamics, can naturally address the relevant time and length scales for these systems. KMC simulation results for the stochastic models are described, focusing on the reshaping of slightly elongated nanorods and of mildly truncated octahedra and nanocubes. In addition, we review appropriate theoretical formulations for reshaping, which involves the nucleation and growth on 2D islands or layers on outer facets of the NC. We note the limitations of classical nucleation theory in some scenarios and demonstrate the successes of a more fundamental and general master equation-based analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shape Stability of Truncated Octahedral fcc Metal Nanocrystals

Metallic nanocrystals (NCs) can be synthesized with tailored nonequilibrium shapes to enhance desired properties, e.g., octahedral fcc metal NCs optimize catalytic activity associated with {111} facets. However, maintenance of optimized properties requires stability against thermal reshaping. Thus, we analyze the reshaping of truncated fcc metal octahedra mediated by surface diffusion using a stochastic atomistic-level model with energetic input parameters for Pd. Furthermore, the model describes NC thermodynamics by an effective nearest-neighbor interaction and includes a realistic treatment of diffusive hopping for undercoordinated surface atoms. Kinetic Monte Carlo simulation reveals that the effective barrier, E eff , for the initial stage of reshaping is strongly tied to the degree of truncation of the vertices in the synthesized initial octahedral shapes. This feature is elucidated via exact analytic determination of the energy variation along the optimal reshaping pathway at low-temperature (T), which involves transfer of atoms from truncated {100} vertex facets to form new layers on {111} side facets. Deviations from predictions of the low-T analysis due to entropic effects are more prominent for higher T and larger NC sizes.

36 MATERIALS SCIENCE↗

Reaction mechanism and kinetics for CO 2 reduction on nickel single atom catalysts from quantum mechanics

Experiments have shown that graphene-supported Ni-single atom catalysts (Ni-SACs) provide a promising strategy for the electrochemical reduction of CO 2 to CO, but the nature of the Ni sites (Ni-N 2 C 2 , Ni-N 3 C 1 , Ni-N 4 ) in Ni-SACs has not been determined experimentally. Here, we apply the recently developed grand canonical potential kinetics (GCP-K) formulation of quantum mechanics to predict the kinetics as a function of applied potential (U) to determine faradic efficiency, turn over frequency, and Tafel slope for CO and H 2 production for all three sites. We predict an onset potential (at 10 mA cm –2 ) U onset = –0.84 V (vs. RHE) for Ni-N 2 C 2 site and U onset = –0.92 V for Ni-N 3 C 1 site in agreement with experiments, and U onset = –1.03 V for Ni-N 4 . We predict that the highest current is for Ni-N 4 , leading to 700 mA cm –2 at U = –1.12 V. To help determine the actual sites in the experiments, we predict the XPS binding energy shift and CO vibrational frequency for each site.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Directing reaction pathways via in situ control of active site geometries in PdAu single-atom alloy catalysts

Abstract The atomic scale structure of the active sites in heterogeneous catalysts is central to their reactivity and selectivity. Therefore, understanding active site stability and evolution under different reaction conditions is key to the design of efficient and robust catalysts. Herein we describe theoretical calculations which predict that carbon monoxide can be used to stabilize different active site geometries in bimetallic alloys and then demonstrate experimentally that the same PdAu bimetallic catalyst can be transitioned between a single-atom alloy and a Pd cluster phase. Each state of the catalyst exhibits distinct selectivity for the dehydrogenation of ethanol reaction with the single-atom alloy phase exhibiting high selectivity to acetaldehyde and hydrogen versus a range of products from Pd clusters. First-principles based Monte Carlo calculations explain the origin of this active site ensemble size tuning effect, and this work serves as a demonstration of what should be a general phenomenon that enables in situ control over catalyst selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic lattice distortions driven by surface trapping in semiconductor nanocrystals

Nonradiative processes limit optoelectronic functionality of nanocrystals and curb their device performance. Nevertheless, the dynamic structural origins of nonradiative relaxations in nanocrystals are not understood. Here, femtosecond electron diffraction measurements corroborated by atomistic simulations uncover transient lattice deformations accompanying radiationless electronic processes in semiconductor nanocrystals. Investigation of the excitation energy dependence shows that hot carriers created by a photon energy considerably larger than the bandgap induce structural distortions at nanocrystal surfaces on few picosecond timescales associated with the localization of trapped holes. On the other hand, carriers created by a photon energy close to the bandgap result in transient lattice heating that occurs on a much longer 200 ps timescale, governed by an Auger heating mechanism. Elucidation of the structural deformations associated with the surface trapping of hot holes provides atomic-scale insights into the mechanisms deteriorating optoelectronic performance and a pathway towards minimizing these losses in nanocrystal devices.

36 MATERIALS SCIENCE↗