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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

The Development of Directional Decohesion Finite Elements for Multiscale Failure Analysis of Metallic Polycrystals

Atomistic simulations of intergranular fracture have indicated that grain-scale crack growth in polycrystalline metals can be direction dependent. At these material length scales, the atomic environment greatly influences the nature of intergranular crack propagation, through either brittle or ductile mechanisms, that are a function of adjacent grain orientation and direction of crack propagation. Methods have been developed to obtain cohesive zone models (CZM) directly from molecular dynamics simulations. These CZMs may be incorporated into decohesion finite element formulations to simulate fracture at larger length scales. A new directional decohesion element is presented that calculates the direction of Mode I opening and incorporates a material criterion for dislocation emission based on the local crystallographic environment to automatically select the CZM that best represents crack growth. The simulation of fracture in 2-D and 3-D aluminum polycrystals is used to illustrate the effect of parameterized CZMs and the effectiveness of directional decohesion finite elements.

Saether, Erik↗
DOE OSTI2021

Disconnection-Mediated Transition in Segregation Structures at Twin Boundaries

Twin boundaries play an important role in the thermodynamics, stability, and mechanical properties of nanocrystalline metals. Understanding their structure and chemistry at the atomic scale is key to guide strategies for fabricating nanocrystalline materials with improved properties. In this work, we report an unusual segregation phenomenon at gold-doped platinum twin boundaries, which is arbitrated by the presence of disconnections, a type of interfacial line defect. By using atomistic simulations, we show that disconnections containing a stacking fault can induce an unexpected transition in the interfacial-segregation structure at the atomic scale, from a bilayer, alternating-segregation structure to a trilayer, segregation-only structure. This behavior is found for faulted disconnections of varying step heights and dislocation characters. Supported by a structural analysis and the classical Langmuir–McLean segregation model, we reveal that this phenomenon is driven by a structurally induced drop of the local pressure across the faulted disconnection accompanied by an increase in the segregation volume.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1978

Equation of state fits to the lower mantle and outer core

The lower mantle and outer core are subjected to tests for homogeneity and adiabaticity. An earth model is used which is based on the inversion of body waves and Q-corrected normal-mode data. Homogeneous regions are found at radii between 5125 and 4825 km, 4600 and 3850 km, and 3200 and 2200 km. The lower mantle and outer core are inhomogeneous on the whole and are only homogeneous in the above local regions. Finite-strain and atomistic equations of state are fit to the homogeneous regions. The apparent convergence of the finite-strain relations is examined to judge their applicability to a given region. In some cases the observed pressure derivatives of the elastic moduli are used as additional constraints. The effect of minor deviations from adiabaticity on the extrapolations is also considered. An ensemble of zero-pressure values of the density and seismic velocities are found for these regions. The range of extrapolated values from these several approaches provides a measure of uncertainties involved.

Butler, R.↗
DOE OSTI2023

Atomistic simulations and machine learning of solute grain boundary segregation in Mg alloys at finite temperatures

Understanding solute segregation thermodynamics is the first step in investigating grain boundary (GB) properties, such as strong yttrium (Y) effects on grain growth and texture evolution in micro-scale polycrystalline magnesium (Mg) alloys. To estimate the average GB segregation behavior in low-solute-concentration Mg alloys (e.g., 2 at.% Y), a state-of-the-art spectral approach is applied based on a per-site segregation energy spectrum for Y solute atoms at zero K obtained from molecular statistics (MS) simulations of ~10 4 GB sites in Mg symmetric tilt GBs (STGBs). Although selected MS simulation results are consistent with verification by density functional theory (DFT) calculations, estimates of average segregation tendency based on the zero-K energy spectrum deviate from experimental observations. To resolve this problem, thermodynamic integration (TI) methods based on molecular dynamics (MD) simulations are used to determine the per-site segregation free energies of Y at representative GB sites, which show contributions beyond harmonic approximations can be important for certain GB sites at high temperatures. A surrogate model of per-site segregation free energy is constructed from a small set of TI data points using stacking cross-validation regressors and physics-informed descriptors. This model is applied to predict the Y segregation free energy spectra for thousands of GB sites in Mg STGBs with uncertainty quantification. Finally, the average segregation tendency predicted by the spectral approach based on the free energy spectra agrees well (within the uncertainty range) with experimental observations for micro-scale polycrystalline Mg alloys at typical thermomechanical processing temperatures (500 ~ 800 K), where thermodynamic equilibrium states are likely to be achieved due to fast diffusion.

Atomistic simulations↗
DOE OSTI2021

Decoupling between Shockley partials and stacking faults strengthens multiprincipal element alloys

We report that mechanical properties are fundamental to structural materials, where dislocations play a decisive role in describing their mechanical behavior. Although the high-yield stresses of multiprincipal element alloys (MPEAs) have received extensive attention in the last decade, the relation between their mechanistic origins remains elusive. Our multiscale study of density functional theory, atomistic simulations, and high-resolution microscopy shows that the excellent mechanical properties of MPEAs have diverse origins. The strengthening effects through Shockley partials and stacking faults can be decoupled in MPEAs, breaking the conventional wisdom that low stacking fault energies are coupled with wide partial dislocations. This study clarifies the mechanistic origins for the strengthening effects, laying the foundation for physics-informed predictive models for materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine Learning Thermodynamics And Kinetics of Defects For Accelerated Materials Discovery

Atomistic defects play a pivotal role in functional and structural materials’ performance across a myriad of technology applications. Quantitative prediction of the thermodynamics and kinetics of defect formation and migration, respectively, typically requires accurate but expensive first-principles approaches, such as density functional theory (DFT). Their computational expense limits the throughput needed to perform high-throughput materials discovery/screening exercises or to perform materials modeling tasks relying on extensive sampling techniques. Therefore, in this Sandia National Laboratories Laboratory Directed Research and Development (LDRD) project (Project #229366), we developed a variety of machine learning techniques, trained on density functional theory calculations, to accelerate the discovery and modeling of materials in which vacancy and interstitial defects primarily dictate material performance. These include applications such as metal oxides for water-splitting or mixed ionic-electronic conduction, metal hydrides for hydrogen storage, and transition metal dichalcogenides for electronics, and the approaches developed herein can further be applied to many other domains that similarly depend on materials’ thermodynamic and kinetic defect properties for their desired functionality.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of Nanoconfinement on NMR Relaxation of Heptane in Kerogen from Molecular Simulations and Measurements

Kerogen-rich shale reservoirs will play a key role during the energy transition, yet the effects of nanoconfinement on the NMR relaxation of hydrocarbons in kerogen are poorly understood. Here, we use atomistic MD simulations to investigate the effects of nanoconfinement on the 1H NMR relaxation times T 1 and T 2 of heptane in kerogen. In the case of T 1 , we discover the important role of confinement in reducing T 1 by ~3 orders of magnitude from that of bulk heptane, in agreement with measurements of heptane dissolved in kerogen from the Kimmeridge Shale, without any models or free parameters. In the case of T 2 , we discover that confinement breaks spatial isotropy and gives rise to residual dipolar coupling which reduces T 2 by ~5 orders of magnitude from the value for bulk heptane. We use the simulated T 2 to calibrate the surface relaxivity and thence predict the pore-size distribution of the organic nanopores in kerogen, without additional experimental data.

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

Three-dimensional atomic structure and local chemical order of medium- and high-entropy nanoalloys

Medium- and high-entropy alloys (M/HEAs) mix several principal elements with near-equiatomic composition and represent a model-shift strategy for designing previously unknown materials in metallurgy, catalysis and other fields. One of the core hypotheses of M/HEAs is lattice distortion, which has been investigated by different numerical and experimental techniques. However, determining the three-dimensional (3D) lattice distortion in M/HEAs remains a challenge. Moreover, the presumed random elemental mixing in M/HEAs has been questioned by X-ray and neutron studies, atomistic simulations, energy dispersive spectroscopy and electron diffraction, which suggest the existence of local chemical order in M/HEAs. However, direct experimental observation of the 3D local chemical order has been difficult because energy dispersive spectroscopy integrates the composition of atomic columns along the zone axes and diffuse electron reflections may originate from planar defects instead of local chemical order. Here, in this work, we determine the 3D atomic positions of M/HEA nanoparticles using atomic electron tomography and quantitatively characterize the local lattice distortion, strain tensor, twin boundaries, dislocation cores and chemical short-range order (CSRO). We find that the high-entropy alloys have larger local lattice distortion and more heterogeneous strain than the medium-entropy alloys and that strain is correlated to CSRO. We also observe CSRO-mediated twinning in the medium-entropy alloys, that is, twinning occurs in energetically unfavoured CSRO regions but not in energetically favoured CSRO ones, which represents, to our knowledge, the first experimental observation of correlating local chemical order with structural defects in any material. We expect that this work will not only expand our fundamental understanding of this important class of materials but also provide the foundation for tailoring M/HEA properties through engineering lattice distortion and local chemical order.

36 MATERIALS SCIENCE↗
DOE OSTI2026

ChemGraph as an agentic framework for computational chemistry workflows

Atomistic simulations are essential in chemistry and materials science but remain challenging to run due to the expert knowledge required for the setup, execution, and validation stages of these calculations. We present ChemGraph, an agentic framework powered by artificial intelligence and state-of-the-art simulation tools to streamline and automate computational chemistry and materials science workflows. ChemGraph leverages graph neural network-based foundation models for accurate yet computationally efficient calculations and large language models (LLMs) for natural language understanding, task planning, and scientific reasoning to provide an intuitive and interactive interface. We evaluate ChemGraph across 13 benchmark tasks and demonstrate that smaller LLMs (GPT-4o-mini, Claude-3.5-haiku, Qwen-2.5-14B) perform well on simple workflows, while more complex tasks benefit from using larger models. Importantly, we show that decomposing complex tasks into smaller subtasks through a multi-agent framework enables GPT-4o to reach perfect accuracy and smaller LLMs to match or exceed single-agent GPT-4o's performance in these benchmarks.

Computational chemistry↗
DOE OSTI2025

Multiscale Explanation of the Missing Gallium Vacancy in Gallium Arsenide

Irradiation of gallium arsenide (GaAs) produces immobile vacancies and mobile interstitials. Yet, after decades of experimental investigation, the immobile Ga vacancy continues to evade detection, raising the question: where is the Ga vacancy? Static first-principles calculations predict a Ga vacancy should be readily observed. We find that short-time dynamical evolution of primary defects is the key to explaining this conundrum. Using a dynamical multiscale atomistically informed device engineering (AIDE) method, we discover that during the initial displacement damage, the Ga vacancy (3-/2-) defect level pins the Fermi level near the midgap, producing oppositely charged vacancies and interstitials. Driven by Coulomb attraction, fast As interstitials preferentially annihilate Ga vacancies. The Ga vacancy population plummets below detectable limits—and the now unpinned Fermi level recovers—before being experimentally observed. This dynamical model solves the mystery of the missing Ga vacancy and reveals the importance of a multiscale approach to explore the dynamical chemical behavior in experimentally inaccessible short-time regimes.

Diaz, Leopoldo [Sandia National Laboratories (SNL-↗
DOE OSTI2025

Establishing the Role of Metal, Interface, and Vacancy Sites in Pt/TiO 2 -Catalyzed Acetic Acid Hydrodeoxygenation

Catalytic hydrodeoxygenation (HDO) following catalytic fast pyrolysis (CFP) offers an approach to convert the vapor-phase product of biomass pyrolysis to a stable bio-oil product by reducing the oxygen content. Fundamental insights into the HDO of carboxylic acids, which are a corrosive and acidic CFP product, on promising catalyst materials, such as Pt/TiO 2 , are needed to inform the design of multifunctional HDO catalysts with improved carbon efficiency. In this contribution, density functional theory (DFT) calculations were used to assess the role of Pt-metal and Pt-TiO 2 -interface sites on acetic acid HDO (AA-HDO), and to determine the effect of interfacial oxygen vacancies at the Pt-TiO 2 interface, by calculating the reaction energetics for key AA-HDO surface intermediates and elementary steps on each site type. Pt-metal sites, modeled via Pt(111), preferred to form undesired decarboxylation products (CH 4 and CO 2 ), whereas Pt-TiO 2 -interface sites, modeled via an anatase-supported Pt nanowire, favored the formation of desired deoxygenation products (acetaldehyde and ethane). Interfacial-vacancy sites lowered the activation energy barrier for the first C-O bond-scission step in AA-HDO, predicted to be the rate-limiting step for AA-HDO at the Pt-TiO 2 interface in the absence of a vacancy. These atomistic insights reveal the importance of metal-metal oxide interface sites in AA-HDO selectivity and can be used to inform the rational design of improved HDO catalysts.

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

An atomistic survey of shear coupling in asymmetric tilt grain boundaries and interpretation using the disconnections framework

Grain Boundaries (GB) play an important role in determining the behavior of polycrystalline materials. While the mechanisms of motion and associated shear response for symmetric tilt grain boundaries (STGBs) are well studied, the same is not true for asymmetric tilt grain boundaries (ATGBs) despite their greater prevalence in polycrystals. Here, this study aims to investigate the shear response of a large collection of asymmetric tilt grain boundaries (ATGBs) using molecular dynamics (MD) simulations and interpret the data using a discrete disconnections model that works remarkably well for STGBs. MD simulations of shear-driven ATGBs show that the plastic shear (shear coupling factor) in the region swept by a GB exhibits a complex dependence on the inclination angle, and this dependence changes with the misorientation of the GB. In addition, the shear response was observed to be highly sensitive to the applied shear rate and temperature. Recognizing the spatial and temporal scale limitations of MD simulations, we extended the discrete disconnections mesoscale model of Khateret al. (2012) to calculate the nucleation barriers of disconnection modes and predict the effective shear coupling of an ATGB. We observed that the mesoscale model’s predictions of the shear coupling factor of ATGBs do not agree with those observed in MD simulations. Finally, we examine the hypotheses of our mesoscale model that contribute to disagreements between MD simulations and the mesoscale model and propose improvements to the mesoscale model for future work.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Generalized quantum master equations can improve the accuracy of semiclassical predictions of multitime correlation functions

Multitime quantum correlation functions are central objects in physical science, offering a direct link between the experimental observables and the dynamics of an underlying model. While experiments such as 2D spectroscopy and quantum control can now measure such quantities, the accurate simulation of such responses remains computationally expensive and sometimes impossible, depending on the system’s complexity. A natural tool to employ is the generalized quantum master equation (GQME), which can offer computational savings by extending reference dynamics at a comparatively trivial cost. However, dynamical methods that can tackle chemical systems with atomistic resolution, such as those in the semiclassical hierarchy, often suffer from poor accuracy, limiting the credence one might lend to their results. By combining work on the accuracy-boosting formulation of semiclassical memory kernels with recent work on the multitime GQME, here we show for the first time that one can exploit a multitime semiclassical GQME to dramatically improve both the accuracy of coarse mean-field Ehrenfest dynamics and obtain orders of magnitude efficiency gains.

Chemistry↗
DOE OSTI2024

Machine learning for accuracy in density functional approximations

Machine learning techniques have found their way into computational chemistry as indispensable tools to accelerate atomistic simulations and materials design. In addition, machine learning approaches hold the potential to boost the predictive power of computationally efficient electronic structure methods, such as density functional theory, to chemical accuracy and to correct for fundamental errors in density functional approaches. In this paper, recent progress in applying machine learning to improve the accuracy of density functional and related approximations is reviewed. Promises and challenges in devising machine learning models transferable between different chemistries and materials classes are discussed with the help of examples applying promising models to systems far outside their training sets.

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

Predicting the location of shear band initiation in a metallic glass

We report atomistic simulation results which indicate that the location of shear banding in a metallic glass (MG) can be ascertained with reasonably high accuracy solely from the undeformed static structure. Correlation is observed between the location of the initiation of shear bands in a simulated MG and the initial distribution of the density of fertile sites (DFS) for stress-driven shear transformations identified a priori based on a deep learning model devised in our recent work. In addition, we demonstrated that one can judge whether a glass is brittle or ductile solely based upon its initial DFS distribution. These validate that shear bands in MG arise from nonlinear instabilities and that the as-quenched glass structure contains inhomogeneities that influence these instabilities. In this study, we also reveal an important subtlety regarding the nondeterministic nature of athermal quasistatic shear simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Defect Characterization of Monocrystalline Silicon Solar Cells with Polysilicon Passivated Contact Using Electrically-Detected Magnetic Resonance (EDMR) Spectroscopy

As the c-Si based solar cell efficiencies are approaching over 26%, it is becoming critical to characterize the low concentrations of the defects – as low as 10^10-10^11 cm-3 (for e.g., iron contamination in high-lifetime Ga-doped wafers3 and n-type wafers), and further reduce them. Also, atomistic level understanding of the mechanisms of the low concentration process-induced-defects and reliability limiting defects (such as light and elevated temperature induced degradation, surface passivation degradation) is needed to design the mitigation strategies. The conventional characterization techniques are limited due to their detection limitations. Some of the techniques based on lifetime spectroscopies can still be used for low concentration characterization however, they are based on estimations and theoretical models and hence, indirect and cannot fully reveal information about the microscopic mechanism of the defects. Thus, we present the application of an ultrasensitive magnetic resonance-based technique for the direct spectroscopic detection of the defects in Si PV - electrically detected magnetic resonance (EDMR). In this work, we aim to focus on establishing a process flow for fabrication of minicells with (miniature replica of the larger-area cells) and setting up the routine for EDMR measurements on them with the EDMR instrumentation capability at NREL. For the EDMR measurements, sample size is limited by the dimensions of sample holder tube (width less than 3.2 mm, active area - 20 mm). Thus, we have designed c-Si based minicells with polysilicon (poly-Si) passivated contacts same as the larger-area cells that we fabricate in our group at NREL. We also modified our minicell process flow for fabricating the textured minicells for preserving the texture during processing and taking care of the laser-ablation edge damage which can significantly affect the performance of such small devices. We have achieved comparable performance on these newly fabricated minicells as that of our 4 cm2 devices with same structure (comparable VOC, JSC, FF). We also conducted EDMR measurements on the minicells and observed a distinct EDMR signal at g-value ~2.005 at temperatures 30K and above, as shown in Fig. 2. We associate this signal to the presence of silicon dangling bonds based on the g-value. We also observed an EDMR signal at g-value ~1.998 at temperature ~5K. The origin of this signal is still being investigated. Thus, we show the proof of concept of minicells and EDMR measurements with which we now aim to study some of the unknown defects in silicon solar cell devices.

EDMR↗
DOE OSTI2021

Anisotropic grain boundary area and energy distributions in tungsten

Describing microstructure evolution in tungsten requires a quantitative description of the anisotropic grain boundary energy. Here we present a grain boundary energy function for tungsten that specifies the energy of an arbitrary boundary given its five macroscopic crystallographic parameters. A comparison of measured grain boundary areas and the grain boundary energies given by the function at the Σ11, Σ17b, and Σ33a misorientations, which are problematic to determine by measurement or atomistic calculations, reveals inverse correlations that are similar to what have been observed in other metals.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Polymers in Deep Eutectic Solvents

The project investigated the behavior of polymers in ionic liquids such as deep eutectic solvents using an array of techniques. These included the development of atomistic force fields, coarse graining these force fields to the united atom level, large scale molecular dynamics (MD) simulations using new thermostat algorithms, and machine learning (ML) methods for the phase behavior. The project demonstrated the feasibility and accuracy of first principles force fields for ionic liquids, deep-eutectic solvents, and urea-water mixtures. A novel hierarchical coarse graining method was then used to develop accurate and efficient united-atom models, using which microsecond simulations were performed for polymers in ionic liquids. These simulations were in quantitative agreement with experiment, thus resolving previous controversies. Methods were also developed to obtain the potential of mean force between complex ions in solution. Finally, supervised ML methods were developed for the phase behavior of polymers in ionic liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗