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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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 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 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↗
DOE OSTI2025

Heavy boron doping effects on biaxially tensile strained germanium (>1.5%) investigated via structural characterization, effective lifetime assessment and atomistic modeling

Highly tensile strained germanium (ε-Ge) represents an essential material system for emerging electronic and photonics applications. Moreover, adjusting the doping levels to moderate or high concentrations can effectively tailor the properties of ε-Ge for specific applications. This article combines experimental characterization with a theoretical framework to examine the effects of heavy elemental boron (B) doping on pseudomorphic sub-50 nm ε-Ge. High resolution X-ray diffractometry is used to validate tensile strain levels of 1.53% and 1.68% in Ge epilayers, surpassing the indirect-to-direct band gap crossover point at ∼1.5% biaxial tensile strain. Cross-sectional transmission electron microscopy revealed visual evidence of stacking faults and surface roughening in 1.68% ε-Ge, although a coherent and abrupt Ge/III–V heterointerface is observed, devoid of interfacial misfit dislocations. Effective lifetime measurements demonstrated approximately twofold enhancement in 1.53% B-doped ε-Ge (N B ∼7 × 10 19 cm −3 ) compared to its unstrained B-doped counterpart, while no such improvement was observed in 1.68% B-doped ε-Ge. This lack of enhancement is attributed to the presence of stacking faults and surface roughness within the ε-Ge epilayer. Through density functional theory calculations, we independently demonstrate that substitutional B atoms induce local deformation of Ge–Ge bonds in both unstrained Ge and ε-Ge epilayers, resulting in an additive tensile strain. This phenomenon could potentially lead to dynamic reduction and overcoming of the critical layer thickness for the system, facilitating the nucleation and subsequent glide of 90° leading Shockley partial dislocations, thereby generating stacking faults. In essence, these findings establish an upper limit on the B-doping concentration that can be achieved in highly ε-Ge epilayers, and collectively, offer valuable insights into the significance of heavy doping in Ge-based heterostructures. As such, this study delineates a fundamental constraint for integrating heavily doped ε-Ge in high-performance optoelectronic systems, necessitating precise strain-doping co-optimization to avoid performance degradation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Atomistic models of In and Ga diffusion in Cu(In,Ga)Se 2

To improve the performance of Cu(In,Ga)Se 2 thin-film photovoltaic devices, a robust understanding of the dominant diffusion pathways of the alloy species In and Ga is needed. Here, the most probable defect complexes and mechanisms for In and Ga diffusion are identified with the aid of density functional theory. The binding energies and migration barriers for these complexes are calculated in bulk CuInSe 2 and CuGaSe 2 . Analytic models and kinetic lattice Monte Carlo simulations are employed to predict the diffusivity of In and Ga under variations in composition and temperature. Here, we find that a model based on coulombic interactions between group III antisites and vacancies on the Cu-sublattice produces results that match well with experiment.

14 SOLAR ENERGY↗
DOE OSTI2024

Exploring dielectric properties in atomistic models of amorphous boron nitride

Abstract We report a theoretical study of dielectric properties of models of amorphous Boron Nitride, using interatomic potentials generated by machine learning. We first perform first-principles simulations on small (about 100 atoms in the periodic cell) sample sizes to explore the emergence of mid-gap states and its correlation with structural features. Next, by using a simplified tight-binding electronic model, we analyse the dielectric functions for complex three dimensional models (containing about 10.000 atoms) embedding varying concentrations of sp 1 , sp 2 and sp 3 bonds between B and N atoms. Within the limits of these methodologies, the resulting value of the zero-frequency dielectric constant is shown to be influenced by the population density of such mid-gap states and their localization characteristics. We observe nontrivial correlations between the structure-induced electronic fluctuations and the resulting dielectric constant values. Our findings are however just a first step in the quest of accessing fully accurate dielectric properties of as-grown amorphous BN of relevance for interconnect technologies and beyond.

Materials Science↗
DOE OSTI2026

Atomistic modeling of lanthanide diffusion in refractory body-centered cubic molybdenum

Lanthanide fission products can strongly interact with candidate cladding alloys, but their transport properties in refractory metals remain poorly understood. Here, in this work, we investigate the atomic-scale diffusion behavior of La, Ce, Pr, and Nd in body-centered cubic (bcc) molybdenum, a potential candidate for advanced nuclear cladding. Self-consistent mean-field transport modeling is performed to evaluate the fission product transport and vacancy mobility, informed by first-principles and nudged elastic band calculations of vacancy formation energies, migration barriers, and solute–vacancy binding characteristics. Compared with bcc Fe, lanthanide solutes in bcc Mo exhibit slower tracer diffusion due to higher vacancy formation and migration energies. Furthermore, the calculations reveal that the influence of fission products on migration barriers in bcc Mo are not as extensive in range compared to bcc Fe. Among the studied lanthanides, La exhibits the strongest vacancy binding while also being the fastest diffuser in Mo. These findings highlight how refractory bcc alloys can reduce fission product infiltration, offering valuable insight into the development of durable cladding systems for advanced reactors.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Atomistic modeling of LiF microstructure ionic conductivity and its influence on nucleation and plating

We report the formation and degradation of the solid electrolyte interphase (SEI) and its underlying transport properties play an essential role in the overall performance of lithium-ion batteries. This paper presents classical molecular dynamics studies on polycrystalline inorganic lithium fluoride (LiF) layers to model and predict the SEI transport properties. The ionic conductivity is obtained from the lithium-ion diffusivity in polycrystalline structures of LiF using the Nernst-Einstein relation. The predicted molecular dynamics data are used in a continuum scale phase-field model to evaluate the plating kinetics under fast charging conditions. The analysis emphasizes that the SEI ionic conductivity properties impact the plating dynamics, where SEI's low ion conductivity value is prone to large plating and subsequent capacity degradation. The combination of atomic and continuum scale studies shown herein lays a foundation to tune in SEI transport properties to decrease the amount of lithium plating and improve the performance of fast-charging batteries.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗