Diffusion Mechanisms in Bcc-Zr:A Molecular Dynamics Approach

1990 ◽  
Vol 209 ◽  
Author(s):  
F. Willaime ◽  
C. Massobrio

ABSTRACTBasing our calculations on a realistic N-body interatomic potential for Zr, we study the vacancy migration mechanism and determine the related diffusion coefficient in the bcc phase. The form of the potential energy along the nearest-neighborjump migration path is single-peaked. The vacancy jump rate determined by molecular dynamics simulations has a perfectly Arrhenian behavior and its activation energy is very close to the static value of the vacancy migration energy, both being very low (≈ 0.3 eV) . The diffusion coefficient is in very satisfactory agreement with experiments.

2000 ◽  
Vol 646 ◽  
Author(s):  
D. Farkas ◽  
B. Soulé de Bas

ABSTRACTMolecular dynamics simulations of the diffusion process in ordered B2 NiAl at high temperature were performed using an embedded atom interatomic potential. Diffusion occurs through a variety of cyclic mechanisms that accomplish the motion of the vacancy through nearest neighbor jumps restoring order to the alloy at the end of the cycle. The traditionally postulated 6-jump cycle is only one of the various cycles observed and some of these are quite complex. A detailed sequential analysis of the observed 6-jump cycles was performed and the results are analyzed in terms of the activation energies for individual jumps calculated using molecular statics simulations.


Materials ◽  
2021 ◽  
Vol 14 (15) ◽  
pp. 4326
Author(s):  
Pawel Goj ◽  
Aleksandra Wajda ◽  
Pawel Stoch

Iron-phosphate glasses, due to their properties, have many potential applications. One of the most promising seems to be nuclear waste immobilization. Radioactive 90Sr isotope is the main short-lived product of fission and, due to its high solubility, it can enter groundwater and pose a threat to the environment. On the other hand, Sr is an important element in hard tissue metabolic processes, and phosphate glasses containing Sr are considered bioactive. This study investigated the effect of SrO addition on a glass structure of nominal 30Fe2O3-70P2O5 chemical composition using classical molecular dynamics simulations. To describe the interaction between Sr-O ion pairs, new interatomic potential parameters of the Buckingham-type were developed and tested for crystalline compounds. The short-range structure of the simulated glasses is presented and is in agreement with previous experimental and theoretical studies. The simulations showed that an increase in SrO content in the glass led to phosphate network depolymerization. Analysis demonstrated that the non-network oxygen did not take part in the phosphate network depolymerization. Furthermore, strontium aggregation in the glass structure was observed to lead to the non-homogeneity of the glass network. It was demonstrated that Sr ions prefer to locate near to Fe(II), which may induce crystallization of strontium phosphates with divalent iron.


Membranes ◽  
2021 ◽  
Vol 11 (5) ◽  
pp. 355
Author(s):  
Tamar Zelovich ◽  
Mark E. Tuckerman

Fuel cell-based anion-exchange membranes (AEMs) and proton exchange membranes (PEMs) are considered to have great potential as cost-effective, clean energy conversion devices. However, a fundamental atomistic understanding of the hydroxide and hydronium diffusion mechanisms in the AEM and PEM environment is an ongoing challenge. In this work, we aim to identify the fundamental atomistic steps governing hydroxide and hydronium transport phenomena. The motivation of this work lies in the fact that elucidating the key design differences between the hydroxide and hydronium diffusion mechanisms will play an important role in the discovery and determination of key design principles for the synthesis of new membrane materials with high ion conductivity for use in emerging fuel cell technologies. To this end, ab initio molecular dynamics simulations are presented to explore hydroxide and hydronium ion solvation complexes and diffusion mechanisms in the model AEM and PEM systems at low hydration in confined environments. We find that hydroxide diffusion in AEMs is mostly vehicular, while hydronium diffusion in model PEMs is structural. Furthermore, we find that the region between each pair of cations in AEMs creates a bottleneck for hydroxide diffusion, leading to a suppression of diffusivity, while the anions in PEMs become active participants in the hydronium diffusion, suggesting that the presence of the anions in model PEMs could potentially promote hydronium diffusion.


2019 ◽  
Vol 97 (11) ◽  
pp. 795-804 ◽  
Author(s):  
Dong Xiang ◽  
Weihua Zhu

The density functional tight-binding molecular dynamics approach was used to study the mechanisms and kinetics of initial pyrolysis and combustion reactions of isolated and multi-molecular FOX-7. Based on the thermal cleavage of bridge bonds, the pyrolysis process of FOX-7 can be divided into three stages. However, the combustion process can be divided into five decomposition stages, which is much more complex than the pyrolysis reactions. The vibrations in the mean temperature contain nodes signifying the formation of new products and thereby the transitions between the various stages in the pyrolysis and combustion processes. Activation energy and pre-exponential factor for the pyrolysis and combustion reactions of FOX-7 were obtained from the kinetic analysis. It is found that the activation energy of its pyrolysis and combustion reactions are very low, making both take place fast. Our simulations provide the first atomic-level look at the full dynamics of the complicated pyrolysis and combustion process of FOX-7.


10.29007/6kp3 ◽  
2020 ◽  
Author(s):  
Renji Mukuno ◽  
Manabu Ishimaru

The structural changes of amorphous silicon (a-Si) under compressive pressure were examined by molecular-dynamics simulations using the Tersoff interatomic potential. a-Si prepared by melt-quenching methods was pressurized up to 30 GPa under different temperatures (300K and 500K). The density of a-Si increased from 2.26 to 3.24 g/cm3 with pressure, suggesting the occurrence of the low-density to high-density amorphous phase transformation. This phase transformation occurred at the lower pressure with increasing the temperature because the activation barrier for amorphous-to-amorphous phase transformation could be exceeded by thermal energy. The coordination number increased with pressure and time, and it was saturated at different values depending on the pressure. This suggested the existence of different metastable atomic configurations in a-Si. Atomic pair-distribution functions and bond-angle distribution functions suggested that the short-range ordered structure of high-density a-Si is similar to the structure of the high-pressure phase of crystalline Si (β-tin and Imma structures).


Author(s):  
Keivan Esfarjani ◽  
Gang Chen ◽  
Asegun Henry

Based on first-principles density-functional calculations, we have developed and tested a force-field for silicon, which can be used for molecular dynamics simulations and the calculation of its thermal properties. This force field uses the exact Taylor expansion of the total energy about the equilibrium positions up to 4th order. In this sense, it becomes systematically exact for small enough displacements, and can reproduce the thermodynamic properties of Si with high fidelity. Having the harmonic force constants, one can easily calculate the phonon spectrum of this system. The cubic force constants, on the other hand, will allow us to compute phonon lifetimes and scattering rates. Results on equilibrium Green-Kubo molecular dynamics simulations of thermal conductivity as well as an alternative calculation of the latter based on the relaxation-time approximation will be reported. The accuracy and ease of computation of the lattice thermal conductivity using these methods will be compared. This approach paves the way for the construction of accurate bulk interatomic potentials database, from which lattice dynamics and thermal properties can be calculated and used in larger scale simulation methods such as Monte Carlo.


Author(s):  
Andrey I. Kartamyshev ◽  
Alexey Gennadievich Lipnitskii ◽  
Anton Boev ◽  
Ivan Nelasov ◽  
Vyacheslav N. Maksimenko ◽  
...  

2020 ◽  
Vol 6 (3) ◽  
pp. 50 ◽  
Author(s):  
Steven Best ◽  
Jake B. Wasley ◽  
Carla de Tomas ◽  
Alireza Aghajamali ◽  
Irene Suarez-Martinez ◽  
...  

Amorphous carbons are disordered carbons with densities of circa 1.9–3.1 g/cc and a mixture of sp2 and sp3 hybridization. Using molecular dynamics simulations, we simulate diffusion in amorphous carbons at different densities and temperatures to investigate the transition between amorphous carbon and the liquid state. Arrhenius plots of the self-diffusion coefficient clearly demonstrate that there is a glass transition rather than a melting point. We consider five common carbon potentials (Tersoff, REBO-II, AIREBO, ReaxFF and EDIP) and all exhibit a glass transition. Although the glass-transition temperature (Tg) is not significantly affected by density, the choice of potential can vary Tg by up to 40%. Our results suggest that amorphous carbon should be interpreted as a glass rather than a solid.


Sign in / Sign up

Export Citation Format

Share Document