Search NASA⌕ Search

SEARCH · Search NASA

Results for “Atomistic simulations”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quantum bath augmented stochastic nonequilibrium atomistic simulations for molecular heat conduction

Classical molecular dynamics (MD) has been shown to be effective in simulating heat conduction in certain molecular junctions since it inherently takes into account some essential methodological components which are lacking in the quantum Landauer-type transport model, such as many-body full force-field interactions, anharmonicity effects and nonlinear responses for large temperature biases. However, the classical MD reaches its limit in the environments where the quantum effects are significant (e.g. with low-temperatures substrates, presence of extremely high frequency molecular modes). Here, we present an atomistic simulation methodology for molecular heat conduction that incorporates the quantum Bose–Einstein statistics into an “effective temperature” in the form of a modified Langevin equation. We show that the results from such a quasi-classical effective temperature MD method deviates drastically when the baths temperature approaches zero from classical MD simulations and the results converge to the classical ones when the bath approaches the high-temperature limit, which makes the method suitable for full temperature range. In addition, we show that our quasi-classical thermal transport method can be used to model the conducting substrate layout and molecular composition (e.g. anharmonicities, high-frequency modes). Anharmonic models are explicitly simulated via the Morse potential and compared to pure harmonic interactions to show the effects of anharmonicities under quantum colored bath setups. Finally, the chain length dependence of heat conduction is examined for one-dimensional polymer chains placed in between quantum augmented baths.

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

Multiscale Modeling of Grain-Boundary Fracture: Cohesive Zone Models Parameterized From Atomistic Simulations

A multiscale modeling strategy is developed to study grain boundary fracture in polycrystalline aluminum. Atomistic simulation is used to model fundamental nanoscale deformation and fracture mechanisms and to develop a constitutive relationship for separation along a grain boundary interface. The nanoscale constitutive relationship is then parameterized within a cohesive zone model to represent variations in grain boundary properties. These variations arise from the presence of vacancies, intersticies, and other defects in addition to deviations in grain boundary angle from the baseline configuration considered in the molecular dynamics simulation. The parameterized cohesive zone models are then used to model grain boundaries within finite element analyses of aluminum polycrystals.

Glaessgen, Edward H.↗
NASA NTRS1998

Efficient 3D 'Atomistic' Simulation Technique for Studying of Random Dopant Induced Threshold Voltage Lowering and Fluctuations in Decanano MOSFETs

A 3D 'atomistic' simulation technique to study random dopant induced threshold voltage lowering and fluctuations in sub 0.1 micron MOSFETs is presented. It allows statistical analysis of random impurity effects down to the individual impurity level. Efficient algorithms based on a single solution of Poisson's equation, followed by the solution of a simplified current continuity equation are used in the simulations.

Asenov, Asen↗
DOE OSTI2024

Atomistic simulations for investigation of substrate and salt effects on lipid in-source fragmentation in secondary ion mass spectrometry: A follow-up study

In-source fragmentation (ISF) poses a significant challenge in secondary ion mass spectrometry (SIMS). These fragment ions increase the spectral complexity and can lead to incorrect annotation of fragments as intact species. The presence of salt that is ubiquitous in biological samples can influence the fragmentation and ionization of analytes in a significant manner, but their influences on SIMS have not been well characterized. To elucidate the effect of substrates and salt on ISF in SIMS, we have employed experimental SIMS in combination with atomistic simulations of a sphingolipid on a gold surface with various NaCl concentrations as a model system. Our results revealed that a combination of bond dissociation energy and binding energy between N-palmitoyl-sphingomyelin and a gold surface is a good predictor of fragment ion intensities in the absence of salt. However, ion-fragment interactions play a significant role in determining fragment yields in the presence of salt. Additionally, the charge distribution on fragment species may be a major contributor to the varying effects of salt on fragmentation. This study demonstrates that atomistic modeling can help predict ionization potential when salts are present, providing insights for more accurate interpretations of complex biological spectra.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory↗
DOE OSTI2021

Manifold learning for coarse-graining atomistic simulations: Application to amorphous solids

In this work, we introduce a generalized machine learning framework to probabilistically parameterize upper-scale models in the form of nonlinear PDEs consistent with a continuum theory, based on coarse-grained atomistic simulation data of mechanical deformation and flow processes. The proposed framework utilizes a hypothesized coarse-graining methodology with manifold learning and surrogate-based optimization techniques. Coarse-grained high-dimensional data describing quantities of interest of the multiscale models are projected onto a nonlinear manifold whose geometric and topological structure is exploited for measuring behavioral discrepancies in the form of manifold distances. A surrogate model is constructed using Gaussian process regression to identify a mapping between stochastic parameters and distances. Derivative-free optimization is employed to adaptively identify a unique set of parameters of the upper-scale model capable of rapidly reproducing the system's behavior while maintaining consistency with coarse-grained atomic-level simulations. The proposed method is applied to learn the parameters of the shear transformation zone (STZ) theory of plasticity that describes plastic deformation in amorphous solids as well as coarse-graining parameters needed to translate between atomistic and continuum representations. We show that the methodology is able to successfully link coarse-grained microscale simulations to macroscale observables and achieve a high-level of parity between the models across scales.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomistic simulations of He bubbles in Beryllium

Formation of He bubbles can have a significant effect on the microstructural evolution and properties of irradiated materials. In this report, we use atomistic simulations based on machine learning potentials to investigate the fundamental behavior of He bubbles in Be, with a specific focus on the shape, stability, and diffusivity of bubbles. Stability of He bubbles is quantified in terms of formation energies, which are determined as a function of the ratio of He/V. We find that He bubbles become unstable with respect to plastic deformation through punch-out dislocations around the bubble when the He/V ratio is larger than ~1.25, and the punch-out process induces the change of the regular bubble shape. In general, the bubble shape of He in Be is found to be ellipsoid-like. It is also found that for a fixed He/V ratio, the bubble attracts vacancies to become larger in size. If the bubble size is constant, the bubble attracts additional He atoms until the punch-out reaction occurs. The dominant diffusion mechanism of He bubbles changes from surface diffusion to volume diffusion as the temperature is increased, with a crossover occurring at about 900 K.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Extended Rice–Thomson analysis and atomistic simulations revealing grain boundary effects on fracture in refractory high-entropy alloys

Significance This work serves to extend the fundamental ductile vs. brittle fracture theory, specifically the Rice–Thomson criterion, by introducing a grain boundary ahead of an initiating crack which propagates at an oblique angle to impinge the boundary. Atomistic fracture simulations on two refractory complex concentrated alloys, the brittle NbMoTaW and the ductile Nb 45 Ta 25 Ti 15 Hf 15 , demonstrate qualitative correspondence with the extended Rice–Thomson criterion and experimental observations. Abstract Understanding how grain boundaries mediate fracture remains a critical challenge in designing ductile, high-performance refractory alloys. Here, we extend the Rice–Thomson criterion to account for the angle between cracks and the impinging grain boundaries (GBs), capturing the competition between intergranular fracture and dislocation-mediated plasticity. Using machine learning interatomic potentials, we performed molecular statics simulations to probe fracture mechanisms in nanocrystalline NbMoTaW and Nb 45 Ta 25 Ti 15 Hf 15 , each with two different grain sizes, revealing trends consistent with experimental observations and the extended Rice model. Comparison with averaged R-curves for bulk samples demonstrates that GBs enhance ductility in Nb 45 Ta 25 Ti 15 Hf 15 in both grain sizes investigated. In contrast, GBs only locally improve fracture resistance in NbMoTaW when cracks are temporarily pinned at GBs inclined at high angles from the crack, but generally promote brittle intergranular fracture. These contrasting behaviors are attributed to differences in GB cohesion, reflecting clear alloying trends that align with ab-initio calculations and trends observed experimentally. Our results bridge classical fracture theory, atomistic simulations, and experimental observations, providing a comprehensive understanding of the fracture mechanisms in nanocrystalline refractory complex concentrated alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

From anti-Arrhenius to Arrhenius behavior in a dislocation-obstacle bypass: Atomistic simulations and theoretical investigation

Dislocations are the primary carriers of plasticity in metallic materials. Understanding the basic mechanisms for dislocation movement is paramount to predicting the material mechanical response. Relying on atomistic simulations, we observe a transition from non-Arrhenius to Arrhenius behavior in the rate of an edge dislocation overcoming the long-range elastic interaction with a prismatic loop in tungsten. Close to the critical resolved shear stress, the process shows a non-Arrhenius behavior at low temperatures. However, as the temperature increases, the activation entropy starts to dominate, leading to a traditional Arrhenius-like behavior. Here, we have computed the activation entropy analytically along the minimum energy path following Schoeck’s method, which captures the cross-over between anti-Arrhenius and Arrhenius domains. Also, the Projected Average Force Integrator (PAFI), another simulation method to compute free energies along an initial transition path, exhibits considerable concurrence with Schoeck’s formalism. We conclude that entropic effects need to be considered to understand processes involving dislocations bypassing elastic barriers close to the critical resolved shear stress. More work needs to be performed to fully understand the discrepancies between Schoeck’s and PAFI compared to molecular dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Diffusion and Deformation Mechanism Maps in Li Metal from Atomistic Simulations

The rate of Li transport in the Li metal anode is important for the operation of Li metal-solid state batteries. Here, transport rates due to diffusion and creep are predicted in Li using atomistic simulations. First, molecular dynamics is used to estimate the rate of Li diffusion along dislocations and in grain boundary triple junctions. By combining this data with that from a prior study of grain boundary diffusion the dominant mechanisms and rates of self-diffusion in Li polycrystals are predicted as a function of grain size, grain shape, dislocation density, and temperature. Second, the dominant creep mechanisms are predicted and used to estimate critical current densities and void annihilation times. Grain boundary sliding and coble creep are the dominant mechanisms for micron-sized grains. Lastly, a continuum model for interfacial contact loss reveals that high dislocation densities of ∼10 12 /cm 2 enable achieving battery performance targets for Li grain sizes of ∼10 μm.

Batteries↗
DOE OSTI2024

Nanocompression of 20 nm Silver Nanoparticles: In situ Aberration‐Corrected TEM and Atomistic Simulations

Single‐crystalline nanoparticles play an increasingly important role in a wide variety of fields including pharmaceuticals, advanced materials, catalysts for fuel cells, energy materials, as well as environmental detection and monitoring. Yet, the deformation mechanisms of very small nanoparticles are still poorly understood, in particular the role played by single dislocations and their interaction with surfaces. Here, in this work, silver nanoparticles with particularly small dimensions (≈20 nanometers in diameter) are compressed in situ in an aberration‐corrected transmission electron microscopy (TEM) and molecular dynamics (MD) simulations. During compression, the emergence of both dislocations and nanotwins are observed. However, these defects prove to be unstable and disappear upon removal of the indenter. Atomistic simulations confirm the role played by image stresses associated with the nearby surfaces and the reduction in dislocation line length as it approaches the free surface, thereby supporting the experimental observations. These results provide justification for the frequent observation of the absence of dislocations in nanoparticles of a few nanometers in size during in situ experiments, even after significant deformation. This phenomenon contributes to the self‐healing of samples through dislocation ejection toward the surfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The {10$\overline{1}$2} non-cozone twin-twin interactions in Mg: A stability and mobility study using 3-D atomistic simulations

Due to easy activation of tensile twinning in Mg, multiple {10$\overline{1}$2} twin variants can interact with each other and form twin-twin boundaries over the course of plastic deformation. Previous studies using 2-D settings provide only a partial understanding of these interactions, especially the non-cozone ones. Here, atomistic simulations are used to study the 3-D structural characteristic and evolution of the non-cozone {10$\overline{1}$2} twin-twin junctions. The study reveals the existence of new twin-twin boundaries (TTBs) such as TTB BP and TTB K2 , formed after the interaction between the basal prismatic and conjugate twin interfaces with the coherent twin boundary. Here, for both non-cozone twin-twin interactions, the {10$\overline{1}$2} TTB and its associated twin-twin junctions are found to play a major role in the {10$\overline{1}$2} twin's stability and mobility. Specifically, they promote the growth of the 3-D twin in both the normal and forward directions during the interaction and hinder the detwinning process upon unloading.

36 MATERIALS SCIENCE↗
NASA NTRS2011

Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 are candidate materials for applications in extreme environments because of their high melting point, good mechanical properties and reasonable oxidation resistance. Unlike many ceramics, these materials have high thermal conductivity which can be advantageous, for example, to reduce thermal shock. Recently, we developed Tersoff style interatomic potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current. Results at room temperature and at elevated temperatures will be reported.

Lawson, JOhn W.↗
DOE OSTI2021

Comparison of Biotite Elastic Properties Recovered by Spherical Nanoindentations and Atomistic Simulations — Influence of Nano-Scale Defects in Phyllosilicates

Phyllosilicate minerals, due to their sheets structure and morphology, are known to cause anisotropy in bulk rock properties and make the bulk rock more compliant. Accurately characterizing the micromechanical behavior of phyllosilicate minerals from laboratory observations, which eventually translates to the bulk rock behavior, is still challenging due to their fine-grained nature. Recent advances in atomistic simulations open the possibility of theoretically investigating such mineral mechanical behavior. We compare the elastic properties of biotites recovered by spherical nanoindentation with those predicted from density functional theory (DFT) simulations to investigate to what extent theoretical predictions reproduce actual phyllosilicate properties. Spherical nanoindentation was conducted using schist rocks from Poorman Formation, South Dakota, USA, to recover continuous indentation stress-strain curves. Loading in the layer-normal orientation shows an average indentation modulus (M) of about 35 GPa, while loading in the layer-parallel orientation gives a higher average of about 95 GPa. To facilitate comparison, the elastic stiffness constants (cij) determined from DFT were converted to indentation modulus (M) using solutions proposed in this study. The majority of the nanoindentation modulus results are below the values inferred from the simulation results representing ideal defect-free minerals. We suggest that crystal defects present at the nano-scale, potentially ripplocations, are the dominant cause of the lower indentation modulus recovered from nanoindentation compared to those inferred from DFT simulations. Results highlight the importance of acknowledging the defects that exist down to the nano-scale as it modifies the mechanical properties of phyllosilicates compared to its pure defect-free form.

58 GEOSCIENCES↗
DOE OSTI2024

Intricate short-range order in GeSn alloys revealed by atomistic simulations with highly accurate and efficient machine-learning potentials

GeSn alloys hold promise for silicon-compatible integrated applications in electronics, photonics, and topological quantum devices. However, understanding their intricate structures using density functional theory (DFT) calculations is hindered by spatiotemporal constraints. To overcome this limitation, we develop highly accurate and efficient machine-learning interatomic potentials based on a neuroevolution potential approach with farthest point sampling on a comprehensive DFT data set. The application of the developed machine-learning potential in large-scale atomistic simulations bridges the spatiotemporal gap between modeling and advanced characterizations, and facilitates the discovery of structural intricacies in GeSn alloys. Through extensive statistical sampling, we identify a type of short-range order (SRO) that is distinguished by both its structural signature and electronic band gap from the SRO structure previously predicted. Modeling based on a large simulation cell reveals the coexistence of nano SRO domains with various degrees of ordering, demonstrating a complex spatial heterogeneity of SRO structure. Finally, our study not only reinforces the significance of fine-level structural information in alloys, but it also constitutes an effective framework for exploring SRO in a broad range of complex alloys based on highly accurate and effective machine-learning potentials.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 have a number of properties that make them attractive for applications in extreme environments. One such property is their high thermal conductivity. Computational modeling of these materials will facilitate understanding of fundamental mechanisms, elucidate structure-property relationships, and ultimately accelerate the materials design cycle. Progress in computational modeling of UHTCs however has been limited in part due to the absence of suitable interatomic potentials. Recently, we developed Tersoff style parameterizations of such potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current, in contrast to the more typical monotonic decay seen in monoatomic materials such as Silicon, for example. Results at room temperature and at elevated temperatures will be reported.

Lawson, John W.↗
NASA NTRS2012

Lattice Thermal Conductivity from Atomistic Simulations: ZrB2 and HfB2

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 have a number of properties that make them attractive for applications in extreme environments. One such property is their high thermal conductivity. Computational modeling of these materials will facilitate understanding of fundamental mechanisms, elucidate structure-property relationships, and ultimately accelerate the materials design cycle. Progress in computational modeling of UHTCs however has been limited in part due to the absence of suitable interatomic potentials. Recently, we developed Tersoff style parameterizations of such potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current, in contrast to the more typical monotonic decay seen in monoatomic materials such as Silicon, for example. Results at room temperature and at elevated temperatures will be reported.

Lawson, John W.↗
DOE OSTI2022

Structure and glide of Lomer and Lomer-Cottrell dislocations: Atomistic simulations for model concentrated alloy solid solutions

Lomer (L) and Lomer-Cottrell (LC) dislocations have long been considered to be central to work hardening in face-centered cubic (FCC) metals and alloys. These dislocations act as barriers of motion for other dislocations, and can serve as sites for twin nucleation. Recent focus on multicomponent concentrated FCC solid solution alloys has resulted in many reported observations of LC dislocations. While these and L dislocations are expected to have a role in the mechanical behavior of these alloys, little is understood about how variations in composition and associated fault energies change the response of these dislocations under stress. Here we present atomistic simulations of L and LC dislocations in a model Cu-Ni system and find that changes in composition and applied stress conditions result in a wide variety of responses, including changes in core configuration and (100) glide. The results are compared to and extend previous literature related to the nature of L/LC core structures and how they vary with respect to intrinsic materials properties and stress states. This study also provides insights into mechanisms such as twin nucleation that could have important implications for work hardening in FCC solid-solution alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗