Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic simulation”

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 217 records · Page 12

Parameter uncertainties for imperfect surrogate models in the low-noise regime

Abstract Bayesian regression determines model parameters by minimizing the expected loss, an upper bound to the true generalization error. However, this loss ignores model form error, or misspecification, meaning parameter uncertainties are significantly underestimated and vanish in the large data limit. As misspecification is the main source of uncertainty for surrogate models of low-noise calculations, such as those arising in atomistic simulation, predictive uncertainties are systematically underestimated. We analyze the true generalization error of misspecified, near-deterministic surrogate models, a regime of broad relevance in science and engineering. We show that posterior parameter distributions must cover every training point to avoid a divergence in the generalization error and design a compatible ansatz which incurs minimal overhead for linear models. The approach is demonstrated on model problems before application to thousand-dimensional datasets in atomistic machine learning. Our efficient misspecification-aware scheme gives accurate prediction and bounding of test errors in terms of parameter uncertainties, allowing this important source of uncertainty to be incorporated in multi-scale computational workflows.

Swinburne, Thomas D. (ORCID:0000000232554257)↗

Developing Machine Learning Interatomic Potential for Fe-Cr-Ni Alloys

Accurate prediction of creep and fatigue behavior of stainless steel at elevated temperatures in hydrogen environment requires fundamental understanding of alloy-hydrogen interaction at cross-scale including bulk lattice and key defects such as vacancies, grain boundaries, surfaces, stacking faults, dislocations, and precipitates. This project aims to predict creep behavior of 347H stainless steel with H using machine learning interatomic potentials based on first-principles density functional theory simulations. The Moment Tensor Potentials platform is adopted for this work since it demonstrates a fine balance between model accuracy and computational efficiency. The potential is well trained based on large amount of high-fidelity density functional theory calculations. The validation is carried out by comparing various important properties including short range order, coefficient of thermal expansion, elastic properties, stacking fault energy, grain boundary energy, and surface energy. This work lays the foundation for reliable atomistic simulation of high temperature hydrogen attack of stainless steel.

density functional theory (DFT)↗

The Role of Local Inhomogeneities on Dendrite Growth in LLZO-Based Solid Electrolytes

The majority of the ceramic solid electrolytes (LLZO, LATP) demonstrate polycrystalline grain/grain-boundary (G/GB) microstructure. Higher lithium (Li) concentration and lower mechanical stiffness result in current focusing at the GBs. Growth of Li dendrites through local inhomogeneities and subsequent short circuit of the cell is a major concern. Recent studies have revealed that bulk Li metal is a viscoplastic material that has low (~0.3 MPa) and high (~1.0 MPa) yield strength during deformation at smaller and larger rates of strain, respectively. It has been argued that during deposition at smaller current densities, due to its lower yield strength, Li metal should demonstrate plastic flow against stiff ceramic electrolytes, and Li dendrites will be prevented from penetrating through solid electrolytes. In this manuscript, a multiscale modeling framework has been developed for predicting properties of GBs and the bulk of ceramic electrolytes using atomistic calculations for input to mesoscale models. Using the parameters obtained from the atomistic simulations, the mesoscale model reveals that, given enough time, even at low charge rates, lithium dendrites can grow through the GBs of LLZO. The present multiscale model results also provide information regarding the dendrite growth velocity through LLZO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fingerprinting shock-induced deformations via diffraction

Abstract During the various stages of shock loading, many transient modes of deformation can activate and deactivate to affect the final state of a material. In order to fundamentally understand and optimize a shock response, researchers seek the ability to probe these modes in real-time and measure the microstructural evolutions with nanoscale resolution. Neither post-mortem analysis on recovered samples nor continuum-based methods during shock testing meet both requirements. High-speed diffraction offers a solution, but the interpretation of diffractograms suffers numerous debates and uncertainties. By atomistically simulating the shock, X-ray diffraction, and electron diffraction of three representative BCC and FCC metallic systems, we systematically isolated the characteristic fingerprints of salient deformation modes, such as dislocation slip (stacking faults), deformation twinning, and phase transformation as observed in experimental diffractograms. This study demonstrates how to use simulated diffractograms to connect the contributions from concurrent deformation modes to the evolutions of both 1D line profiles and 2D patterns for diffractograms from single crystals. Harnessing these fingerprints alongside information on local pressures and plasticity contributions facilitate the interpretation of shock experiments with cutting-edge resolution in both space and time.

36 MATERIALS SCIENCE↗

Irradiation-Induced Structural Disorder and Its Influence on the Mechanical Response of Polycrystalline MoS2

Molybdenum disulfide (MoS2) thin films are widely used as dry-film lubricants and protective coatings in aerospace and other radiation-exposed environments. Conventional synthesis routes produce polycrystalline films whose grain boundaries and other native defects cause their mechanical and tribological behavior to differ substantially from that of ideal single crystals. Under irradiation, these films progressively evolve from polycrystalline structures, composed of layered MoS2 grains, into highly disordered and eventually amorphous structures, altering both their tribological performance and mechanical integrity. Here, we employ reactive atomistic simulations to investigate irradiation-driven structural evolution in bulk polycrystalline MoS2. Using controlled primary knock-on atom (PKA) events, we characterize the progressive transition from a polycrystalline microstructure to an amorphous network by tracking defect accumulation and structural disorder. We then establish how this transition modifies the dominant deformation mechanisms and the temperature-dependent tensile response. Specifically, irradiation suppresses interlayer sliding and delamination, mechanisms which facilitate the deformation of the pristine polycrystal, resulting in defect-induced hardening. Broadly, our results establish direct process–structure–property relationships linking irradiation-induced defect accumulation, microstructural evolution, deformation mechanisms, and mechanical behavior, providing an atomistic framework for understanding the structural integrity and long-term reliability of irradiated MoS2 coatings.

Moore, Daniel [Sandia National Laboratories (SNL)]↗

Insights into Copper Sulfide Formation from Cu and S K edge XAS and DFT studies

An understanding of the fundamentals of the reaction between CuO with trace amounts of H 2 S to form CuS products is critical for the optimal utilization of this process in sulfur removal applications. Unfortunately, CuS is a complex material, featuring various Cu 2-x S compounds (with 0 ≤ x ≤ 1), distorted crystal phases, and varying electronic structures and coordination environments of Cu and S ions. In this work, we combine ex situ and in situ X-ray absorption spectroscopy (XAS) at S and Cu K edges, fixed bed sorption experiments, DFT simulations, and other characterization techniques to speciate the CuS products formed at different temperatures (298–383 K) and from CuO sorbents with different crystallite sizes (2.8–40 nm). The results of our analysis identify the formation of a distorted CuS layer at the surface of CuO crystals with disulfide groups with shorter Cu–S bonds and higher delocalization of the positive charge of the Cu center into (S 1– )2. This distorted CuS layer dominates the XAS signal at lower temperatures (298–323 K) and at the initial stages of sulfidation at higher temperatures (353 and 383 K) where conversion is low (<40%). First-principles atomistic simulations confirm the thermodynamic favorability of the formation of surface (S 1– )2 on both CuO (111) and ($\bar{1}11$) surfaces, providing further support for our experimental observations. Furthermore, these simulations reveal that the presence of disulfide bonds stabilized surface hydroxyl groups, leading to lower Gibbs Free Energies of their surface migration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Explicit solvent machine-learned coarse-grained model of sodium polystyrene sulfonate to capture polymer structure and dynamics

In this article, strongly charged polyelectrolytes (PEs) demonstrate complex solution behavior as a function of chain length, concentrations, and ionic strength. The viscosity behavior is important to understand and is a core quantity for many applications, but aspects remain a challenge. Molecular dynamics simulations using implicit solvent coarse-grained (CG) models successfully reproduce structure, but are often inappropriate for calculating viscosities. To address the need for CG models which reproduce viscoelastic properties of one of the most studied PEs, sodium polystyrene sulfonate (NaPSS), we report our recent efforts in using Bayesian optimization to develop CG models of NaPSS which capture both polymer structure and dynamics in aqueous solutions with explicit solvent. We demonstrate that our explicit solvent CG NaPSS model with the ML-BOP water model [Chan et al. Nat Commun 10, 379 (2019)] quantitatively reproduces NaPSS chain statistics and solution structure. The new explicit solvent CG model is benchmarked against diffusivities from atomistic simulations and experimental specific viscosities for short chains. We also show that our Bayesian-optimized CG model is transferable to larger chain lengths across a range of concentrations. Overall, this work provides a machine-learned model to probe the structural, dynamic, and rheological properties of polyelectrolytes such as NaPSS and aids in the design of novel, strongly charged polymers with tunable structural and viscoelastic properties

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Grain boundary relaxation in doped nano-grained aluminum

Here, simulation studies are done to understand the role of dopants that segregate preferentially to grain boundaries on the stability of nanocrystalline aluminum. A dopant design framework based on thermodynamic principles, is used to identify the specific dopant type with the highest potential to segregate to grain boundaries in nanocrystalline aluminum. Various elements are evaluated as potential dopants and magnesium is identified to have the highest tendency to segregate to grain boundaries and release the excess free energyleading to the relaxation of the grain boundaries. A systematic combination of atomic structure analysis is then done to correlate grain boundary relaxation and the mechanical response of the magnesium-doped nanocrystalline aluminum at ambient temperature. The atomistic simulations reveal that the preferential partitioning of magnesium dopants to the grain boundaries reduces the excess volume within this region which stabilizes the nanostructure. At low contents, the magnesium dopants are observed initially partition to the grain boundaries, but once saturation of the grain boundaries is reached, excess dopants are accommodated in the crystalline interiors. It is found that the addition of the magnesium dopants even in the dilute limit, enhances the strength of the nanocrystalline aluminum. The formation of large, disordered GBs in doped nanocrystalline aluminum under tensile load allowed it to accommodate the deformation and prohibit crack growth.

36 MATERIALS SCIENCE↗

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗

Kinetics of Carbon Condensation in Detonation of High Explosives: First-Order Phase Transition Theory Perspective

We report the kinetics of carbon condensation, or carbon clustering, in detonation of carbon-rich high explosives is modeled by solving a system of rate equations for concentrations of carbon particles. Unlike previous efforts, the rate equations account not only for the aggregation of particles but also for their fragmentation in a thermodynamically consistent manner. Numerical simulations are performed, yielding the distribution of particle concentrations as a function of time. In addition to that, analytical expressions are obtained for all the distinct steps and regimes of the condensation kinetics, which facilitates the analysis of the numerical results and allows one to study the sensitivity of the kinetic behavior to the variation of system parameters. The latter is important because the numerical values of many parameters are not reliably known at present. The theory of the kinetics of first-order phase transitions is found adequate to describe the general kinetic trends of carbon condensation, as described by the rate equations. Such physical phenomena and processes as the coagulation, nucleation, growth, and Ostwald ripening are observed, and their dependence on various system parameters is studied and reported. It is believed that the present work will become useful when analyzing the present and future results for the kinetics of carbon condensation, obtained from experiments or atomistic simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Recent advances and applications of deep learning methods in materials science

Deep learning (DL) is one of the fastest-growing topics in materials data science, with rapidly emerging applications spanning atomistic, image-based, spectral, and textual data modalities. DL allows analysis of unstructured data and automated identification of features. The recent development of large materials databases has fueled the application of DL methods in atomistic prediction in particular. In contrast, advances in image and spectral data have largely leveraged synthetic data enabled by high-quality forward models as well as by generative unsupervised DL methods. In this article, we present a high-level overview of deep learning methods followed by a detailed discussion of recent developments of deep learning in atomistic simulation, materials imaging, spectral analysis, and natural language processing. For each modality we discuss applications involving both theoretical and experimental data, typical modeling approaches with their strengths and limitations, and relevant publicly available software and datasets. We conclude the review with a discussion of recent cross-cutting work related to uncertainty quantification in this field and a brief perspective on limitations, challenges, and potential growth areas for DL methods in materials science.

36 MATERIALS SCIENCE↗

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↗

Visualization and validation of twin nucleation and early-stage growth in magnesium

The abrupt occurrence of twinning when Mg is deformed leads to a highly anisotropic response, making it too unreliable for structural use and too unpredictable for observation. Here, we describe an in-situ transmission electron microscopy experiment on Mg crystals with strategically designed geometries for visualization of a long-proposed but unverified twinning mechanism. Combining with atomistic simulations and topological analysis, we conclude that twin nucleation occurs through a pure-shuffle mechanism that requires prismatic-basal transformations. Also, we verified a crystal geometry dependent twin growth mechanism, that is the early-stage growth associated with instability of plasticity flow, which can be dominated either by slower movement of prismatic-basal boundary steps, or by faster glide-shuffle along the twinning plane. The fundamental understanding of twinning provides a pathway to understand deformation from a scientific standpoint and the microstructure design principles to engineer metals with enhanced behavior from a technological standpoint.

36 MATERIALS SCIENCE↗

The role of stacking fault tetrahedra on void swelling in irradiated copper

Abstract A long-standing and critical issue in the field of irradiated structural materials is that void swelling is significantly higher in face-centered cubic-structured (fcc) materials (1% dpa −1 ) as compared to that of body-centered cubic-structured (bcc) materials (0.2% dpa −1 ). Despite extensive research in this area, the underlying mechanism of the difference in swelling resistance between these two types of materials is not yet fully understood. Here, by combining atomistic simulations and STEM imaging, we find stacking fault tetrahedra (SFTs) are the primary cause of the high swelling rate in pure fcc copper. We reveal that SFTs in fcc copper are not neutral sinks, different from the conventional knowledge. On the contrary, they are highly biased compared to other types of sinks because of the SFT-point defect interaction mechanism. SFTs show strong absorption of mobile self-interstitial atoms (SIAs) from the faces and vertices, and weak absorption of mobile vacancies from the edges. We compare the predicted swelling rates with experimental findings under varying conditions, demonstrating the distinct contributions of each type of sink. These findings will contribute to understanding the swelling of irradiated structural materials, which may facilitate the design of materials with high swelling resistance.

36 MATERIALS SCIENCE↗

Stabilizing Nickel‐Rich Cathodes in Aqueous Process through Nanocellulose as Water Barrier

Nickel-rich LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NMC 811) cathode offers high voltage and high specific capacity, making it promising for high energy density batteries. However, large-scale manufacturing of aqueous-processed NMC 811 electrodes remains challenging due to proton exchange causing material decomposition and capacity loss. This work addresses this issue by constructing an in situ nanocellulose protective layer for NMC 811 particles via electrostatic interactions during the slurry preparation. For the first time, the interatomic spacing between inter-chains of nanocellulose is measured through wide-angle X-ray scattering and demonstrate the ability to effectively confine interlayer water using atomistic simulations. Moreover, this nanocellulose coverage simultaneously minimizes Li + surface segregation and mitigates water infiltration. Owing to less material decomposition during the aqueous processing, nanocellulose-protected NMC electrodes exhibit higher initial coulombic efficiency (83% vs 62% at 0.1C) and capacity (133 vs 59 mAh g −1 at 6C) than unprotected electrodes. Additionally, optimized aqueous-processed NMC electrodes offer comparable or even superior electrochemical properties compared to the electrodes fabricated using the conventional toxic organic solvent, N-methyl-2-pyrrolidone. Consequently, the developed approach enables affordable, sustainable aqueous processing for Nickel-rich NMC 811 cathodes with excellent electrochemical performances.

25 ENERGY STORAGE↗

Freestanding Ferroelectric Bubble Domains

Bubble-like domains, typically a precursor to the electrical skyrmions, arise in ultrathin complex oxide ferroelectric-dielectric-ferroelectric heterostructures epitaxially clamped with flat substrates. Here, it is reported that these specially ordered electric dipoles can also be retained in a freestanding state despite the presence of inhomogeneously distributed structural ripples. By probing local piezo and capacitive responses and using atomistic simulations, this study analyzes these ripples, sheds light on how the bubbles are stabilized in the modified electromechanical energy landscape, and discusses the difference in morphology between bubbles in freestanding and as-grown states. Furthermore, these results are anticipated to be the starting point of a new paradigm for the exploration of electric skyrmions with arbitrary boundaries and physically flexible topological orders in ferroelectric curvilinear space.

36 MATERIALS SCIENCE↗

Large Exchange Coupling Between Localized Spins and Topological Bands in MnBi 2 Te 4

Magnetism in topological materials creates phases exhibiting quantized transport phenomena with potential technological applications. The emergence of such phases relies on strong interaction between localized spins and the topological bands, and the consequent formation of an exchange gap. However, this remains experimentally unquantified in intrinsic magnetic topological materials. Here, this interaction is quantified in MnBi 2 Te 4 , a topological insulator with intrinsic antiferromagnetism. This is achieved by optically exciting Bi-Te p states comprising the bulk topological bands and interrogating the consequent Mn 3d spin dynamics, using a multimodal ultrafast approach. Ultrafast electron scattering and magneto-optic measurements show that the p states demagnetize via electron-phonon scattering at picosecond timescales. Despite being energetically decoupled from the optical excitation, the Mn 3d spins, probed by resonant X-ray scattering, are observed to disorder concurrently with the p spins. Together with atomistic simulations, this reveals that the exchange coupling between localized spins and the topological bands is at least 100 times larger than the superexchange interaction, implying an optimal exchange gap of at least 25 meV in the surface states. Here, by quantifying this exchange coupling, this study validates the materials-by-design strategy of utilizing localized magnetic order to manipulate topological phases, spanning static to ultrafast timescales.

36 MATERIALS SCIENCE↗

On the Structure–Property Relationship of Semi‐Coherent FeCr 2 O 4 /Cr 2 O 3 Spinel/Corundum Interfaces

Abstract Oxide heterointerfaces are extremely common in both natural and artificial composite structures, including corroded structural materials. Often, key properties such as segregation and atomic transport are dictated by the structure of these interfaces. However, despite this critical link, very few heterointerfaces have been studied in any detail at the atomic scale. Here, one important oxide heterointerface is examined, between spinel and corundum, using the chemical system FeCr 2 O 4 /Cr 2 O 3 as a representative and technologically important case. Using atomistic simulation techniques, it is found that the structure, particularly the local chemistry, of the interface depends on the crystal chemistry at the interface. This atomic and chemical structure further impacts important properties such as defect segregation and mass transport. It is found that defects can nucleate at some regions of these interfaces and migrate back and forth across the corundum layer, suggesting high atomic mobility that may be important for the evolution of spinel/corundum composite structures in extreme conditions.

36 MATERIALS SCIENCE↗