Publications

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.

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.

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.

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.

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.

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.

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.

Austenitic stainless steels (Fe-Cr-Ni) are resistant to hydrogen embrittlement but have not been studied using molecular dynamics simulations due to the lack of an Fe-Cr-Ni-H interatomic potential. Herein we describe our recent progress towards molecular dynamics studies of hydrogen effects in Fe-Cr-Ni stainless steels. We first describe our Fe-Cr-Ni-H interatomic potential and demonstrate its characteristics relevant to mechanical properties. We then demonstrate that our potential can be used in molecular dynamics simulations to derive Arrhenius equation of hydrogen diffusion and to reveal twinning and phase transformation deformation mechanisms in stainless steels.

Abstract not provided.

Abstract not provided.