Atomistic Simulation Activities at Sandia-Hydrogen Materials Consortium (H-Mat)
Abstract not provided.
Abstract not provided.
Materials Today Chemistry
Our ability to shape and finish a component by combined methods of fabrication including (but not limited to) subtractive, additive, and/or no theoretical mass-loss/addition during the fabrication is now popularly known as solid freeform fabrication (SFF). Fabrication of a telescope mirror is a typical example where grinding and polishing processes are first applied to shape the mirror, and thereafter, an optical coating is usually applied to enhance its optical performance. The area of nanomanufacturing cannot grow without a deep knowledge of the fundamentals of materials and consequently, the use of computer simulations is now becoming ubiquitous. This article is intended to highlight the most recent advances in the computation benefit specific to the area of precision SFF as these systems are traversing through the journey of digitalization and Industry-4.0. Specifically, this article demonstrates that the application of the latest materials modelling approaches, based on techniques such as molecular dynamics, are enabling breakthroughs in applied precision manufacturing techniques.
Materials Theory
The fundamental interactions between an edge dislocation and a random solid solution are studied by analyzing dislocation line roughness profiles obtained from molecular dynamics simulations of Fe0.70Ni0.11Cr0.19 over a range of stresses and temperatures. These roughness profiles reveal the hallmark features of a depinning transition. Namely, below a temperature-dependent critical stress, the dislocation line exhibits roughness in two different length scale regimes which are divided by a so-called correlation length. This correlation length increases with applied stress and at the critical stress (depinning transition or yield stress) formally goes to infinity. Above the critical stress, the line roughness profile converges to that of a random noise field. Motivated by these results, a physical model is developed based on the notion of coherent line bowing over all length scales below the correlation length. Above the correlation length, the solute field prohibits such coherent line bow outs. Using this model, we identify potential gaps in existing theories of solid solution strengthening and show that recent observations of length-dependent dislocation mobilities can be rationalized.
Pd readily absorbs hydrogen and its isotopes, and can be used to purify gas mixtures involving tritium. Tritium decays to He, forming He bubbles. Bubbles causes possible PCT effects swelling, He release, all leading to failures. Radioactive decay experiments take many years. Molecular dynamics (MD) studies can be quickly done. No previous MD methods can simulate He bubble nucleation and growth.
Pd readily absorbs hydrogen and its isotopes, and can be used to purify gas mixtures involving tritium. Tritium decays to He, forming He bubbles. Bubbles causes possible PCT effects swelling, He release, all leading to failures. Radioactive decay experiments take many years. Molecular dynamics (MD) studies can be quickly done. No previous MD methods can simulate He bubble nucleation and growth.
RSC Advances
In the past, experimentally observed dislocations were often interpreted using an isolated dislocation assumption because the effect of background dislocation density was difficult to evaluate. Contrarily, dislocations caused by atomistic simulations under periodic boundary conditions can be better interpreted because linear elastic theory has been developed to address the effect of periodic dislocation array in the literature. However, this elastic theory has been developed only for perfect dislocations, but not for dissociated dislocations. The periodic boundary conditions may significantly change the dissociation energy of dislocations and stacking fault width, which in turn, change the deformation phenomena observed in simulations. To enable materials scientists to understand the dislocation behavior under the periodic boundary conditions, we use isotropic elastic theory to analyze the thermodynamics of dissociated periodic dislocations with an arbitrary dislocation character angle. Analytical expressions for force, stacking fault width, and energies are presented in the study. Results obtained from the periodic dislocation array were compared with those obtained from isolated dislocations to shed light on the interpretation of experimentally observed and simulated dislocations.
Abstract not provided.
JOM
Aluminum alloys are being explored as lightweight structural materials for use in hydrogen-containing environments.To understand hydrogen effects on deformation, we perform molecular statics studies of the hydrogen Cottrell atmosphere around edge dislocations in aluminum. First, we calculate the hydrogen binding energies at all interstitial sites in a periodic aluminum crystal containing an edge dislocation dipole. This allows us to use the Boltzmann equation to quantify the hydrogen Cottrell atmosphere. Based on these binding energies, we then construct a continuum model to study the kinetics of the hydrogen Cottrell atmosphere formation. Finally, we compare our results with existing theories and discuss the effects of hydrogen on deformation of aluminum-based alloys.
Computational Materials Science
POur experiments indicated that upon a post-processing anneal, an additively manufactured 316L stainless steel exhibits cubic grains rather than the conventional equiaxed grains. Here, we have used kinetic Monte Carlo simulations to explore the origin of these cubic grains. First, we implemented a new kinetic Monte Carlo model in parallel code SPPARKS to simulate grain growth and recrystallization under a residual energy distribution. Our model incorporates physical properties and real-time, as opposed to generic properties and relative time. We further validated that our SPPARKS simulations reproduced the expected kinetic behavior of single-grain evolution. We then used the validated approach to simulate the anneal of an additively manufactured material under the same conditions used in our experiments. We found that the cubic grains can origin from a periodically varying residual energy that may be present in additively manufactured materials.
Abstract not provided.
Abstract not provided.
Physical Chemistry Chemical Physics
TlBr can surpass CZT as the leading semiconductor for γ- A nd X-radiation detection. Unfortunately, the optimum properties of TlBr quickly decay when an operating electrical field is applied. Quantum mechanical studies indicated that if this property degradation comes from the conventional mechanism of ionic migration of vacancies, then an unrealistically high vacancy concentration is required to account for the rapid aging of TlBr seen in experiments. In this work, we have applied large scale molecular dynamics simulations to study the effects of dislocations on ionic migration of TlBr crystals under electrical fields. We found that electrical fields can drive the motion of edge dislocations in both slip and climb directions. These combined motions eject enormous vacancies in the dislocation trail. Both dislocation motion and a high vacancy concentration can account for the rapid aging of the TlBr detectors. These findings suggest that strengthening methods to pin dislocations should be explored to increase the lifetimes of TlBr crystals.
In order to study the effects of Ni oxidation barriers on H diffusion in Zr, a Ni-Zr-H potential was developed based on an existing Ni-Zr potential. Using this and existing binary potentials H diffusion characteristics were calculated and some limited findings for the performance of Ni on Zr coatings are made.
Abstract not provided.
Physical Chemistry Chemical Physics. PCCP
It has been widely believed that crystalline TlBr can surpass CdZnTe to become the leading semiconductor for γ- and X- radiation detection. The major hurdle to this transition is the rapid aging of TlBr under the operating electrical field. Here, while ionic migration of vacancies has been traditionally the root cause for property degradation, quantum mechanical calculations indicated that the vacancy concentration needed to cause the observed aging must be orders of magnitude higher than the usual theoretical estimate. Recent molecular dynamics simulations and X-ray diffract ion experiments have shown that electrical fields can drive the motion of edge dislocations in both slip and climb directions. Furthermore, these combined motions eject a large number of vacancies. Both dislocation mot ion and vacancy ejection can account for the rapid aging of the TlBr detectors. Based on these new discoveries, the present work applies molecular dynamics simulations to “develop” aging-resistant TlBr crystals by inhibiting dislocation motions.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
ChemPhysChem
Magnesium-based materials provide some of the highest capacities for solid-state hydrogen storage. However, efforts to improve their performance rely on a comprehensive understanding of thermodynamic and kinetic limitations at various stages of (de)hydrogenation. Part of the complexity arises from the fact that unlike interstitial metal hydrides that retain the same crystal structures of the underlying metals, MgH 2 and other magnesium-based hydrides typically undergo dehydrogenation reactions that are coupled to a structural phase transformation. As a first step towards enabling molecular dynamics studies of thermodynamics, kinetics, and (de)hydrogenation mechanisms of Mg-based solid-state hydrogen storage materials with changing crystal structures, we have developed an analytical bond order potential for Mg−H systems. We demonstrate that our potential accurately reproduces property trends of a variety of elemental and compound configurations with different coordinations, including small clusters and bulk lattices. More importantly, we show that our potential captures the relevant (de)hydrogenation chemical reactions 2H (gas)→H 2 (gas) and 2H (gas)+Mg (hcp)→MgH 2 (rutile) within molecular dynamics simulations. This verifies that our potential correctly prescribes the lowest Gibbs free energies to the equilibrium H 2 and MgH 2 phases as compared to other configurations. It also indicates that our molecular dynamics methods can directly reveal atomic processes of (de)hydrogenation of the Mg−H systems.
Abstract not provided.
Abstract not provided.