Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomistic modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Understanding (La,Sr)(Co,Fe)O$_{3-δ}$ Phase Instability within SOECs Using a Combined Experimental and Atomistic Modeling Approach

Understanding the onset of degradation in the air electrode within solid oxide electrolysis cells (SOECs), and the subsequent impact on cell performance, is a critical step in mitigating the performance losses and stability issues of SOECs. In an effort to identify early onset degradation phenomena, SOECs were characterized as fabricated and after testing potentiostatically at 1.3 V for 1000 h at 750 °C. SOEC air electrodes composed of a 1:1 composite of La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3–δ (6428-LSCF) and Gd 0.1 Ce 0.9 O 1.95 (GDC) were studied using synchrotron X-ray diffraction (XRD), scanning transmission electron microscopy coupled with energy dispersive X-ray spectroscopy (STEM-EDS), and X-ray absorption near-edge spectroscopy (XANES) to evaluate the changes in the air electrode structurally and chemically. These techniques show the migration of Sr species from the air electrode through pores in the GDC barrier layer, progressing to the electrolyte boundary, where it accumulates and reacts with (Zr 0.84 Y 0.16 )O 2–δ (YSZ) to form SrZrO 3 . Microscopy results are paired with atomistic simulations to better understand the relationship between the thermodynamic instability of 6428-LSCF and cell fabrication/testing conditions. First-principles calculations reveal that LSCF-6428 is not stable during cell manufacturing and testing conditions, which supports the experimental identification of secondary phases in both as-fabricated and tested cells. Together, these results demonstrate that the challenging environments encountered by SOECs during cell manufacturing and operation lead to instabilities of the target 6428-LSCF anode material and underscore the need for more durable, high-performing SOEC components.

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

Atomistic modeling of interface strengthening in Al-Si eutectic alloys

Al-Si cast alloys are usually composed of α-Al and Al-Si eutectic. Si flakes and Al matrix generally hold cube-on-cube orientation relationship with the primary interface (111) Al ∥(111) Si . Extensive experimental studies demonstrated that Si flakes cannot significantly improve mechanical properties of Al-Si cast alloys. We hypothesize that the weak strengthening effect associated with Si flakes might be attributed to thermomechanical properties of Al-Si interfaces besides their morphologies. To characterize Al-Si interfaces with a large lattice mismatch (> 30%), we proposed the quasi-coincident site lattice (Q-CSL) as reference lattice, and demonstrated that the Q-CSL Al-Si coherent interface has three characteristic coherent structures, one stable and low energy structure and two metastable and high energy structures. The translation vectors for the same type of coherent Q-CSL structures are consistent with three displacement shift complete (DSC) vectors. The two metastable structures can be obtained by shifting the low energy structure with three partial DSC vectors. Semi-coherent interface is composed of the low energy Q-CSL patches and three sets of interface misfit dislocations with Burgers vectors same as the DSC vectors. Atomistic simulations revealed that Al-Si interface exhibits low shear resistance. Ideal shear strength of the Q-CSL coherent interface is 110 MPa and semi-coherent interface is 20 MPa. The low shear resistance is attributed to the glide of interface misfit dislocations. Al-Si interface also exhibits low formation and migration energies of point defects. Owing to low shear strength and low formation and migration energies of point defects, interface sliding or shear readily happen under mechanical loading or during dislocation-interface interactions. Lattice dislocations can cross slip onto or climb along Al-Si interfaces. These reactions decrease the number of accumulated dislocation loops around Si flakes and promote nucleation and emission of lattice dislocations from Al-Si interfaces to matrix, consequently reduce the repulsive force on approaching dislocations and weaken Si flakes strengthening effect. In situ tension and compression tests in a scanning electron microscope reveal relatively weak strengthening effect due to Si flakes, consistent with the computed dislocation interaction with interfaces and shear behavior of interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2026

From atomistic models to machine learning: Predictive design of nanocarbons under extreme conditions

The formation of technologically valuable nanocarbon structures under extreme conditions, such as those produced during high-explosive detonations, remains poorly understood but holds significant potential for the development of controlled synthesis pathways. While detonation shockwaves provide the high-pressure, high-temperature environment required for nanodiamond formation, subsequent cooling and decompression dictate whether the diamond phase is preserved or transformed into other nanocarbon structures. Here, in this study, we employ GPU-accelerated reactive molecular dynamics (ReaxFF) simulations to investigate the graphitization and structural remodeling of detonation nanodiamond under nonlinear quench and pressure-release trajectories. We further investigate how the initial nanodiamond morphology; cuboctahedral, octahedral, or hexagonal prism influences the resulting transformation products. Evolution of nanostructure, allotrope (via simulated x-ray diffraction), carbon hybridization, and ring statistics are tracked during a two-stage quench from 5000 K to 60 GPa. Rapid cooling combined with slow decompression optimizes cubic diamond retention, whereas slow cooling with rapid pressure release promotes surface-to-core graphitization, producing concentric sp 2 -hybridized layers and hollowed inner shells. Octahedral nanodiamonds evolve into carbon nano-onions, initially forming bucky diamonds that progressively transform into fully sp 2 -hybridized structures, while hexagonal prisms preferentially form parallel-stacked graphite layers resembling carbon dots. Transient hexagonal diamond (lonsdaleite) emerges as an interfacial phase, suggesting potential reversibility in the shock-induced graphite-to-diamond transformation pathway transformation route. To extend predictive capabilities, we trained machine learning (ML) regressors on over 10 5 node-hours of molecular dynamics (MD) trajectories. A multilayer perceptron (MLP) model reliably predicts the number of graphitized layers from temperature–pressure trajectories with a coefficient of determination (R 2 ) exceeding 0.90. This high predictive fidelity enables efficient, high-throughput mapping of the synthesis parameter space for optimized graphitization outcomes. Collectively, morphological control combined with optimized quench–decompression conditions promote the selective synthesis of nanocarbon allotropes. This work establishes a data-driven framework for the rational, a priori design of carbon nanomaterials for applications in energy storage, sensing, and biomedicine.

Detonation nanodiamond remodeling↗
DOE OSTI2015

Atomistic modeling of crystal structure of Ca{sub 1.67}SiH{sub x}

The atomic structure of calcium-silicate-hydrate (C{sub 1.67}-S-H{sub x}) has been investigated by theoretical methods in order to establish a better insight into its structure. Three models for C-S-H all derived from tobermorite are proposed and a large number of structures were created within each model by making a random distribution of silica oligomers of different size within each structure. These structures were subjected to structural relaxation by geometry optimization and molecular dynamics steps. That resulted in a set of energies within each model. Despite an energy distribution between individual structures within each model, significant energy differences are observed between the three models. The C-S-H model related to the lowest energy is considered as the most probable. It turns out to be characterized by the distribution of dimeric and pentameric silicates and the absence of monomers. This model has mass density which is closest to the experimental one.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Sensitivity of Dislocation-GB interactions to simulation setups in atomistic models

Dislocation-grain boundary (GB) interactions play an integral role in strengthening of crystalline materials, and are dependent on external loading conditions and atomic arrangements at the grain boundary. While molecular dynamics (MD) simulations can provide critical insights into relationships between these parameters and the exact interaction, the outcomes are sensitive to the computational setup. Here, in this work, we explore the effect of computational setup (system size, GB orientation and boundary conditions) on stress-controlled dislocation-grain boundary (DGB) reactions using MD simulations for a well-studied copper symmetric $\langle$1 1 2$\rangle$ tilt boundary. The study demonstrates that the DGB reactions (pinning, absorption, transmission, etc.) are sensitive to the system size and the orientation of the GB to the applied stress. Additionally, the current work shows that there is a critical system size for the dislocation-grain boundary reaction stresses to converge, one which was found to be much larger than the setups in comparable MD studies. This was attributed to the specific stress states and the influence of surfaces, fixed or free, as these computational setups are modified. Finally, a stress-controlled setup with minimal model artifacts is examined, and the effect of coupled Schmid and non-Schmid stress components on the dislocation-GB reactions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mechanistic insights into nitrogen activation on atomic Ru clusters in self-pillared pentasil using operando atomistic models and experimental kinetics

Alternative catalysts to the industrial Haber Bosch process have been of significant interest in the field of heterogeneous catalysis, yet realizing ammonia synthesis under mild conditions (e.g., 300 °C and 10 bar) is challenging due to the low per-pass conversion. One strategy is to promote the associative ammonia synthesis mechanism which eschews direct N-N bond cleavage. Here, in this work, we use self-pillared pentasil, a self-pillared hierarchical zeolite built by thin MFI zeolite nanosheets, as a support for subnanometric Ru clusters to synthesize ammonia. We show that Ru remains well-dispersed during reaction and further demonstrate that ammonia synthesis rates are higher than Cs-Ru/MgO. Reaction kinetics show a positive order in H 2 providing evidence for the associative mechanism, which then becomes negative in H 2 if Ru is allowed to aggregate into nanoparticles. Operando Density Functional Theory models for Ru speciation in SPP, free energy diagrams, and microkinetic modeling were then applied to develop a reaction mechanism that involves sequential hydrogenation of N 2 from metallic Ru clusters. For this hydrogenation to occur, there are site requirements for N 2 to adopt a bridge-bound configuration that facilitates sequential hydrogenation on single sites and metal clusters. These site requirements in turn inform the design of improved zeolite-supported ammonia synthesis catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Anomalous small-angle X-ray scattering of nanoporous two-phase atomistic models for amorphous silicon–germanium alloys

The present work deals with a detailed analysis of the anomalous small-angle X-ray scattering in amorphous silicon–germanium alloy using the simulation technique. We envisage the nanoporous two-phase alloy model consisting in a mixture of Ge-rich and Ge-poor domains and voids at the nanoscale. By substituting Ge atoms for Si atoms in nanoporous amorphous silicon network, compositionally heterogeneous alloys are generated with various composition-contrasts between the two phases. After relaxing the as-generated structure, we compute its radial distribution function, and then we deduce by the Fourier transform technique its anomalous X-ray scattering pattern. Using a smoothing procedure, the computed X-ray scattering patterns are corrected for the termination errors due to the finite size of the model, allowing so a rigorous quantitative analysis of the anomalous small-angle scattering. Our simulation shows that, as expected, the anomalous small-angle X-ray scattering technique is a tool of choice for characterizing compositional heterogeneities coexisting with structural inhomogeneities in an amorphous alloy. Furthermore, the sizes of the compositional nanoheterogeneities, as measured by anomalous small-angle X-ray scattering technique, are X-ray energy independent. A quantitative analysis of the separated reduced anomalous small-angle X-ray scattering, as defined in this work, provided a good estimate of their size.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nanopolysaccharide Builder: A User-Friendly Tool for Atomistic Models of Polysaccharide-Based Nanostructures

Here, we introduce Nanopolysaccharide Builder (NPB), a user-friendly software tool designed to construct polysaccharide nanostructures─mainly those based on cellulose, chitin, and chitosan─using experimental data or user-defined parameters. NPB enables the generation of cellulose and chitin allomorphs with customizable biochemical topologies and also facilitates the construction of large bundles that replicate nanostructures found in biological support systems, including plant cell walls and arthropod cuticles. The software outputs atomic Cartesian coordinates in Protein Data Bank (PDB) format and also provides atom connectivity files in PSF and PARM formats, ensuring seamless integration with major molecular dynamics (MD) engines such as NAMD, CHARMM, GROMACS, AMBER, OpenMM, and LAMMPS. Built on an interactive visualization framework, NPB features a graphical user interface (GUI) and supports both macOS and Linux operating systems. By enabling detailed atomic-scale studies of polysaccharide evolution in extracellular matrices and cell walls of algae, bacteria, fungi, and plants, NPB is poised to advance AI-guided research in sustainable chemical development and biomass utilization.

Wan, Zhangmin [Univ. of British Columbia, Vancouve↗
DOE OSTI2020

Graphite Lithiation under Fast Charging Conditions: Atomistic Modeling Insights

The charging of lithium ion batteries in a fast and safe manner is critical for promoting the mass adoption of electric vehicles. Li intercalation in graphite electrodes is known to be one of the bottlenecks during the fast charging process. The mechanism of Li diffusion in highly polarized graphite anode at high current rates remains, however, not well understood. In this work, Density Functional Theory (DFT) calculations are used to gain insights into the Li diffusion process in graphite when it is far from equilibrium under fast charging conditions. The effect of uncompensated charges on Li mobility is determined in the highly polarized regions of the anode close to the interfaces. The extra charge was found to increase the interlayer spacing in the diffusion layer and adjacent channels, increasing the diffusivity and promoting the formation of Li clusters. A concerted diffusion mechanism at the edge of high-concentration Li domains is proposed to enhance the diffusion of Li.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗