Soft anharmonic phonons and ultralow thermal conductivity in Mg 3 (Sb, Bi) 2 thermoelectrics
The atomistic origin of ultralow thermal conductivity in promising new thermoelectrics is revealed by experiments and simulations.
36 MATERIALS SCIENCE↗
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.
Neural Thermal Scattering (NeTS) Modules for Graphite & Beryllium [Slides]
This presentation discusses the motivations behind the project, which include facilitating compact formulation for TSL data, extending AI knowledge, and providing advanced reactor simulation framework. It also discusses how the project accounts for beryllium and graphite complexity, including atomistic, dynamical, and neural. Additionally, new material-informed neural thermal scattering (NeTS) Modules are examined. The implications, conclusions, and plan future work for the project are also discussed.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Energetics of Single Substitutional Impurities in NiTi
Shape-memory alloys are of considerable current interest, with applications ranging from stents to Mars rover components. In this work, we present results on the energetics of single substitutional impurities in B2 NiTi. Specifically, energies of Pd, Pt, Zr and Hf impurities at both Ni and Ti sites are computed. All energies are computed using the CASTEP ab initio code, and, for comparison, using the quantum approximate energy method of Bozzolo, Ferrante and Smith. Atomistic relaxation in the vicinity of the impurities is investigated via quantum approximate Monte Carlo simulation, and in cases where the relaxation is found to be important, the resulting relaxations are applied to the ab initio calculations. We compare our results with available experimental work.
Good, Brian S.↗
Crystallization of the P 3 Sn 4 Phase upon Cooling P 2 Sn 5 Liquid by Molecular Dynamics Simulation Using a Machine Learning Interatomic Potential
We performed molecular dynamics simulations to study the crystallization of the P3Sn4 phase from P 2 Sn 5 liquid using a machine learning (ML) interatomic potential with desirable efficiency and accuracy. Our results capture the liquid properties of P 2 Sn 5 at 1300 K, which is well above the melting temperature. The phase separation and crystallization are observed when P 2 Sn 5 liquid is cooled down below 832 and 505 K, respectively. The simulation results are in good agreement with the experimentally observed phase transformation behaviors and provide useful insights into the complex nucleation and crystallization process at the details of atomistic scale. Our work also demonstrated that ML interatomic potentials based on neural network deep learning are robust and capable of accurately describing the energetics and kinetics of complex materials through molecular dynamics simulations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Effects of Alloying Elements on Twinning in Ni-Based Superalloys
Micro-twinning is the major creep deformation mechanism in Ni-based superalloys at temperatures above 700 °C. Recent experiments suggest that superlattice stacking faults in y' phase may serve as the precursors to twin formation. Segregation of alloying elements to these precursors may have a significant effect on formation and extension of micro-twins. Using atomistic modeling we investigate and explain the effects of Nb and Cr alloying additions on these processes. The simulation shows that Nb increases the creep resistance which is mostly associated with impeding the reordering of the high energy double complex stacking fault. Cr, on the other hand, promotes twin growth, degrading the high temperature creep properties. These results can help to understand the effects of elemental composition of the alloy on creep resistance.
Ni-based superalloys↗
Effects of alloying elements on twinning in Ni-based superalloys
Micro-twinning is the major creep deformation mechanism in Ni-based superalloys at temperatures above 700 °C. Recent experiments suggest that superlattice stacking faults in y' phase may serve as the precursors to twin formation. Segregation of alloying elements to these precursors may have a significant effect on formation and extension of micro-twins. Using atomistic modeling we investigate and explain the effects of Nb and Cr alloying additions on these processes. The simulation shows that Nb increases the creep resistance which is mostly associated with impeding the reordering of the high energy double complex stacking fault. Cr, on the other hand, promotes twin growth, degrading the high temperature creep properties. These results can help to understand the effects of elemental composition of the alloy on creep resistance.
Ni-based superalloys↗
Atomistic Modeling of Irradiation-Induced Defects and Clusters in Additively-Manufactured Austenitic Stainless Steel
Understanding the progression of irradiation-induced defects in real materials and alloys requires insight into the interaction between point defects, defect clusters, and the existing microstructure. This is particularly true when it comes to additively manufactured (AM) materials with complicated microstructures and residual stresses. In this work, we employ molecular dynamics (MD) to simulate the creation of point defects via primary radiation damage modeling and point defect interactions with extended defects such as voids and dislocations using simplified but representative alloy compositions. We analyze the effect of alloying elements on the formation energies of the defect clusters and measure interactions such as sink strengths and binding energies. This study provides insight into the unique behaviors of irradiation damage in AM materials and information necessary to bridge the gap in length and time scales to connect the atomistic behavior of defects to longer length scale simulations such as cluster dynamics.
36 MATERIALS SCIENCE↗
Hole formation effect on surface morphological response of plasma-facing tungsten
We report simulation results on the effect of helium (He) bubble bursting-mediated surface hole formation on the surface morphological response of tungsten plasma-facing components (PFCs) in nuclear fusion devices. Our analysis is based on an atomistically informed, continuum-scale model, which is capable of accessing the spatiotemporal scales relevant to the fuzz nanostructure formation process on the surface of PFC tungsten. Our simulations account, in an empirical fashion, for two types of subsurface bubble dynamical phenomena in the nanobubble region of PFC tungsten during He plasma irradiation, involving bubble bursting and surface crater formation. We demonstrate that the hole formation effect on the PFC tungsten surface accelerates the growth rate of nanotendrils and the onset of fuzz formation. As a result, the predicted incubation time for surface nanotendril growth is shortened, in agreement with experimental data of incubation fluence at comparable plasma exposure conditions. We also explore systematically the dependence of the PFC surface morphological response on the areal density of holes introduced at regular time intervals onto the He-implanted tungsten surface, a parameter in our analysis that serves as a proxy for the rate of He bubble bursting. More importantly, our simulations capture fine surface features in the PFC tungsten surface morphology and predict that the average spacing between nanotendrils is on the order of 100 nm, consistent with the experimental findings.
Chen, Chao-Shou (ORCID:0000000307910042)↗
Tracking the protein conformational motions driving HIV-1 membrane fusion
HIV-1 Env (trimeric gp120/gp41) is the surface protein responsible for membrane fusion. The Env binds to the receptor proteins, which induces gp120 shedding leading to conformational changes of gp41 from the pre-fusion to post-fusion state, allowing its fusion peptide to embed in the host cell membrane and bringing the viral and host cell membranes together. The gp41 refolding is a target of several peptide inhibitors. Yet, the molecular mechanism of this dynamic process is still not well understood. In this study, we successfully simulate the conformational change of gp41 from pre-fusion to post-fusion state in atomistic resolution using all-atom structure-based models. We reveal that maintaining the directionality of protomer interactions in both pre-fusion and post-fusion states is crucial for gp41 refolding. Additionally, we find that HR1 inherently extends as a three-helical bundle toward the host-cell membrane without any bias. Importantly, we identify native contacts in the pre-fusion state that are critical for the proper refolding of gp41 towards the post-fusion state. Lastly, by incorporating the membrane-fusion inhibitors, T20 and SFT, we identify the most vulnerable stage in the fusion pathway that exhibits the greatest sensitivity to these drugs, which could aid in a better understanding of drug resistance mechanisms.
59 BASIC BIOLOGICAL SCIENCES↗
Localization of Zn 2+ ions affects the structural folding and mechanics of Nereis virens Nvjp-1
Several biological organisms utilize metal-coordination bonds to produce remarkable materials, such as the jaw of the marine worm Nereis virens, where metal-coordination bonds yield remarkable hardness without mineralization. Though the structure of a major component of the jaw, the Nvjp-1 protein, has recently been resolved, a detailed nanostructural understanding of the role of metal ions on the structural and mechanical properties of the protein is missing, especially with respect to the localization of metal ions. In this work, atomistic replica exchange molecular dynamics with explicit water and Zn 2+ ions and steered molecular dynamics simulations were used to explore how the initial localization of the Zn 2+ ions impacts the structural folding and mechanical properties of Nvjp-1. We found that the initial distribution of metal ions for Nvjp-1, and likely for other proteins with high amounts of metal-coordination, has important effects on the resulting structure, with larger metal ion quantity resulting in a more compact structure. These structural compactness trends, however, are independent from the mechanical tensile strength of the protein, which increases with greater hydrogen bond content and uniform distribution of metal ions. Our results indicate that different physical principles underlie the structure or mechanics of Nvjp-1, with broader implications in the development optimized hardened bioinspired materials and the modeling of proteins with significant metal ion content.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Large-scale atomistic studies of sliding friction in polycrystalline aluminum interfaces
In this work, we discuss large-scale non-equilibrium molecular dynamics (NEMD) simulations of ductile metal sliding comprising up to 1.8 × 10 9 atoms over time scales of 100 ns. The results of these simulations have identified a variety of physical mechanisms that are important in determining the steady-state frictional force for a wide range of velocities at compressed metal–metal interfaces. These include grain growth and refinement, the evolution of large plastic strains and strain rates, material mixing, and melting. These phenomena can be included in a strain, strain rate, and grain size model that gives good agreement with the NEMD simulations and can be applied to macroscopic continua.
74 ATOMIC AND MOLECULAR PHYSICS↗
SNW 2000 Proceedings. Oxide Thickness Variation Induced Threshold Voltage Fluctuations in Decanano MOSFETs: a 3D Density Gradient Simulation Study
We use the density gradient (DG) simulation approach to study, in 3D, the effect of local oxide thickness fluctuations on the threshold voltage of decanano MOSFETs in a statistical manner. A description of the reconstruction procedure for the random 2D surfaces representing the 'atomistic' Si-SiO2 interface variations is presented. The procedure is based on power spectrum synthesis in the Fourier domain and can include either Gaussian or exponential spectra. The simulations show that threshold voltage variations induced by oxide thickness fluctuation become significant when the gate length of the devices become comparable to the correlation length of the fluctuations. The extent of quantum corrections in the simulations with respect to the classical case and the dependence of threshold variations on the oxide thickness are examined.
Asenov, Asen↗
New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation
A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.
Saether, E.↗
Defect engineering of porous aromatic frameworks via end capping improves dioxane removal from water
Amorphous porous organic polymers show promise for energy-efficient adsorptive separations, but it is difficult to understand or improve their performance through intentional structural modification. Herein, we report the synthesis of porous aromatic frameworks (PAFs) with pore structures tailored by the incorporation of a monofunctionalized end-capping monomer that disrupts framework topology. Combining experimental characterization with molecular simulations, we show that this defect engineering strategy yields less densely crosslinked networks, which leads to pore collapse and the presence of unique adsorption sites. These defect-engineered PAFs exhibit enhanced removal of 1,4-dioxane, an important environmental pollutant, from water. Increasing the concentration of end-capping monomers produces PAFs with narrower pore size distributions and improved 1,4-dioxane uptake. These results illustrate that defect engineering can effectively modulate polymer connectivity and porosity for applications in selective adsorptive separations. In conclusion, this technique avoids post-synthetic treatments and presents another approach to tailor amorphous polymeric adsorbents.
36 MATERIALS SCIENCE↗
Nanotwin assisted reversible formation of low angle grain boundary upon reciprocating shear load
Severe plastic deformation of metals is known to lead to superior properties that cannot be achieved by any traditional metallurgic process. Origin of the superior properties is perceived to be closely associated with grain refinement, a fundamental process during the severe plastic deformation, which is essentially the formation of new grain boundaries. However, the atomistic mechanism of grain boundary formation remains largely obscure. Here, by using in-situ transmission electron microscopy and molecular dynamic simulation, we reveal, for the first time at atomic level, the shear-induced low-angle grain boundary (LAGB) formation processes in Au nanocrystal. In this work, we discover the LAGB formation is accomplished through inward propagation of nanotwins accompanied by dislocations gliding on the twin boundaries, a nanotwin-mediated dislocation slip mechanism, which shows reversible characteristic under reciprocating shear load and is affected by the nanocrystal microstructure and orientation. Our result unveils unprecedented atomistic insights on shear driven grain refinement towards nanostructure of superior properties.
36 MATERIALS SCIENCE↗
Coarse-Grained and Atomistic Modeling of Polyimides
A coarse-grained model for a set of three polyimide isomers is developed. Each polyimide is comprised of BPDA (3,3,4,4' - biphenyltetracarboxylic dianhydride) and one of three APB isomers: 1,3-bis(4-aminophenoxy)benzene, 1,4-bis(4-aminophenoxy)benzene or 1,3-bis(3-aminophenoxy)benzene. The coarse-grained model is constructed as a series of linked vectors following the contour of the polymer backbone. Beads located at the midpoint of each vector define centers for long range interaction energy between monomer subunits. A bulk simulation of each coarse-grained polyimide model is performed with a dynamic Monte Carlo procedure. These coarsegrained models are then reverse-mapped to fully atomistic models. The coarse-grained models show the expected trends in decreasing chain dimensions with increasing meta linkage in the APB section of the repeat unit, although these differences were minor due to the relatively short chains simulated here. Considerable differences are seen among the dynamic Monte Carlo properties of the three polyimide isomers. Decreasing relaxation times are seen with increasing meta linkage in the APB section of the repeat unit.
Clancy, Thomas C.↗
A Fast Algorithm for Massively Parallel, Long-Term, Simulation of Complex Molecular Dynamics Systems
The advances in theory and computing technology over the last decade have led to enormous progress in applying atomistic molecular dynamics (MD) methods to the characterization, prediction, and design of chemical, biological, and material systems,.
Algorithm Molecular Dynamics Systems Parallel Plat↗
Vibronic and Environmental Effects in Simulations of Optical Spectroscopy
Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.
Chemistry↗