In this work, the results of the evolutionary variable-composition search for binary compounds in the Ti–Si system are presented. The evolutionary algorithm did not find any new stable structures of ...silicides at 0 K and 0GPa. On the other hand, many low-energy metastable and unstable structures are predicted. The 33 predicted and 10 known from literature but previously unstudied structures of compounds with low formation energies with respect to the ground-state line are analyzed. The mechanical properties, electronic band structures, densities of states, and temperature dependencies of thermodynamic stability of 17 dynamically stable previously unstudied structures are calculated.
Display omitted
Display omitted
•Vibrational free energies of hydrogen solid solution and titanium hydrides are calculated.•Solubility of hydrogen in titanium is predicted.•The influence of hydrogen on concentration ...of H-vacancy complexes is estimated.
The saturation solubility of hydrogen with respect to titanium hydrides, and the concentration of H-vacancy complexes in α-Ti are predicted on the basis of a thermodynamic model and ab initio calculated free energies. The presence of hydrogen increases the vacancy complex concentration by several orders of magnitude, with most of these complexes containing three H atoms. We show also that the solubility of hydrogen within the titanium is strongly influenced by the temperature-dependent free energy terms, and that by including these terms in the thermodynamic model it is possible to obtain good agreement with experiment.
In this paper, we evaluated the influence of vibrational and electronic free energy on the thermodynamic phase stability of the Ti5Si3, Ti3Si, Ti6Si3, and other silicides in the binary Ti-Si system ...and the solubility of Si in α-Ti within the density functional theory and harmonic approximation. We showed, that vibrational and electronic free energy thermodynamically stabilize the Ti3Si silicide and destabilize the Ti6Si3 silicide with respect to the Ti5Si3 and α-Ti. According to our calculations, the transition from equilibrium of Si solid solution with the Ti5Si3 silicide to the equilibrium of Si solid solution with the Ti3Si silicide should be at ∼ 900 K. The account of both the vibrational and electronic free energy increase the calculated solubility limit of Si in α-Ti and make it closer to the experimental one.
Adoption of Li-ion batteries for stationary energy storage requires to prolong their operational lifetimes to several decades calling to inhibit all degradation mechanisms. One of such mechanisms is ...associated with gradual formation of Cu–Li alloys at copper current collector - a component traditionally assumed inert in a Li-ion cell. Recent studies shows that Li can penetrate inside a Cu upon electrochemical cycling. We demonstrate that this process may be strongly intensified due to the presence of grain boundaries (GB) in copper. Using density functional theory calculations, we discover that Li segregate at copper GBs, which facilitates its GB diffusion according to vacancy-assisted mechanism. The vacancy formation and migration energies at GBs are significantly decreased in comparison to the values in bulk copper causing strong accelerating of Li diffusion. The segregation is caused by atomic-size effect, while the decrease of migration barriers is explained by a more stable coordination of ionized Li at a saddle point. To prevent Li penetration and subsequent formation of Li–Cu alloys, the amount of grain boundaries in thin film copper current collectors should be minimized.
Display omitted
•Li weakly segregates at Cu grain boundaries.•Li segregation energy linearly depends on the site’s volume at Cu GB.•Li activation energy is smaller at grain boundary than that in bulk Cu.•Grain boundaries promotes Li penetration into copper.
In molecular dynamics simulations, there are dozens of potentials of different type, which take into account the angular dependence of the potential energy, for describing the interactions in carbon. ...Each of them describes simultaneously only a part of its physical properties. Meanwhile, a recently developed N-body approach provides the ability to set the required precision of this angular dependence. In particular, in this work we will show how accuracy in describing three-body interactions influences the accuracy of reproduction of physical properties of carbon on the example of one and two basis functions for these interactions. For both potentials, the identical optimization procedure, including fitting database and weights of each target value, was performed. We found that the last one reproduces the structural, elastic and defect properties of diamond, graphite and graphene phases of carbon in better agreement with experimental and theoretical data than the potential with only one basis function. For additional benchmark of this potential we assessed the melting point of graphite, which is predicted in an acceptable agreement with the experimental value.
Display omitted
•A new interatomic potential for carbon was developed.•An increase of basis functions improves the description of carbon structures.•The potential with one functions predicts graphite structures as unstable.•The potential with two functions predicts three main carbon structures as stable.
We present a new interatomic potential for atomistic modeling of α- and β-phases of zirconium. The potential was developed within the framework of the N-body approach, which takes into account pair, ...three-body and many-body interactions. The new potential predicts a wide range of zirconium properties in agreement with known experimental and DFT data including structural energies, elastic constants, energies of point and planar defects, thermodynamic properties, such as melting temperature and α−β transition temperature. The potential predicts the mechanical instability of the β-phase of zirconium at 0 K, as well as relaxation effects of internal degrees of freedom in elastic constants. The predicted interstitial stacking fault energy is in a good agreement with the DFT data. The formation energies of self-interstitial atoms in the α-phase are in agreement with the known DFT values, predicting the basal octahedral (BO) configuration as the most stable. The potential predicts a spontaneous β−α transition with decreasing temperature. The calculated α−β phase transition temperature is 20% higher than the experimental value. The potential accurately predicts the melting temperature of zirconium, due to which a good agreement with the experimental data was obtained on the temperature dependencies of the diffusion coefficients for the α- and β-phases of zirconium. For both phases, the main contribution to the diffusivity is due to vacancies, while the corresponding values for interstitials are by one or two orders of magnitude lower. The potential predicts the deviation from the Arrhenius law of the temperature dependence of the diffusion coefficient for the β-phase of zirconium, due to the nonlinear dependence of the Gibbs energy of a vacancy formation on temperature. The constructed potential can be used for molecular dynamics simulations of α- and β-phases of zirconium and further development of interatomic potentials within the framework of the N-body approach for alloys containing zirconium.
Display omitted
•A new interatomic potential for zirconium was developed.•The new potential accurately predicts defect and thermal properties of Zr.•The new potential predicts the spontaneous phase transition β–α.•Diffusivity in α and β phases of Zr is in agreement with the experimental data.
Tungsten, as the most refractory metal, is applied in fusion reactor in parts subjected to high temperatures and strong neutron irradiation. These factors lead to intense diffusion processes causing ...degradation of the material. Experimental investigations under such conditions are usually highly complicated and cannot provide a comprehensive understanding of the occurring phenomena. Therefore, their combination with theoretical approaches is required. One of the most robust approaches to simulate diffusion processes is molecular dynamics simulations based on classical interatomic potentials. It allows modeling relatively large samples consisting of several grains, grain boundaries, dislocations, and other types of defects for a reasonable computational time. The reliable simulations of the diffusion process require interatomic potentials satisfying the following criteria: prediction of melting point and thermal expansion as close as possible to the experimental values because the diffusion coefficient strongly depends on the homologous temperature and size factor. In the present paper, we present the new interatomic potential for tungsten, developed within the N-body approach, which reproduces the experimental value of melting temperature (3695 K) and thermal expansion at temperatures up to a melting point. The calculated diffusion coefficient demonstrates adequate agreement with experimental results. The constructed potential is applicable for simulation of processes involving diffusion, one of which is the irradiation damage.
Display omitted
•We constructed new interatomic potential for W.•New potential is faster by two orders of magnitude than the GAP one.•Predicted self-diffusion coefficient in W agrees with the experimental data.•The diffusion coefficient deviates from Arrhenius equation at high temperatures.
We present a new classical interatomic potential for molecular dynamics simulations of copper. The potential was developed within the N-body approach, which, in addition to pair and many-body ...interactions, can take into account three-body ones with the required precision. Our potential accurately predicts a wide range of properties of copper: structural, elastic, point and planar, and thermal such as melting point. The latter is essential for the proper reproduction of the homologous temperature and comparison of the calculated results with the corresponding experimental data, such as diffusion coefficients. The potential predicts the intrinsic stacking fault energy within the range of experimental measurements, which indicates its ability to correctly describe the deformation behavior of copper. In this work, we performed simulations of lattice and grain-boundary diffusion of copper with our potential. It shows a moderate overestimation of the lattice diffusivity with a deviation less than one order of magnitude. To simulate the grain boundary diffusion, we built a general-type grain boundary employing the original method. The calculated grain-boundary diffusivity is in remarkable agreement with the literature data and almost coincides with them at higher temperatures in the considered temperature range. The constructed potential can be effectively used for molecular dynamics simulations of defect properties and diffusion phenomena and composing the interatomic potentials within the N-body approach for alloys containing copper.
Display omitted
•A new interatomic potential for copper was developed.•The new potential accurately predicts defect and thermal properties of Cu.•The predicted Cu lattice diffusion is in adequate agreement with experiment.•The diffusivity in general-type grain boundaries (GTGB) was calculated.•The diffusivity in GTGB is in close agreement with experiment.
In this work, the interatomic potentials for modeling diffusion in the C15 Cr2Ta Laves phase were constructed within the N-body approach. The potential for Ta–Ta interactions reproduces the lattice ...parameter, cohesive energy, elastic constants, equation of state, thermal expansion, point defect energies, and phonon dispersion of body-centered cubic Ta in qualitative agreement with density functional theory (DFT) data and almost quantitative agreement with available experimental data in the temperature range of stability of C15 Cr2Ta Laves phase. The potential for Cr–Ta interactions reproduces the elastic constants, point defect properties, equilibrium volumes, and formation enthalpies of Cr–Ta structures in qualitative agreement with DFT data. Also, it reproduces the lattice stability and high-temperature formation enthalpy of C15 Cr2Ta Laves phase in very close agreement with the available experimental data. The calculations of diffusion coefficients with constructed potentials showed that diffusion in C15 Cr2Ta lattice is governed by Cr atoms which cannot move without creation of vacancies. The constructed potentials can be used in further investigations of diffusion processes in Cr–Ta phases at temperatures up to 2000 K, while the obtained results on diffusion coefficients in C15 Cr2Ta Laves phase would be useful in the rational design of Cr–Ta based alloys.
Display omitted
Display omitted
•We constructed an interatomic potential for titanium.•The potential correctly describe ω-Ti as the ground state at zero temperature.•The potential accurately predicts 0D and 2D ...defects in α-Ti.•Predicted melting point of β-Ti coincides with experiment within error of calculations.•Temperature range of mechanical stability of β-Ti agrees with the experimental data.
In this paper, we present an interatomic potential predicting thermal properties of phases of titanium in their temperature range of stability at zero pressure. The potential was developed within the approach proposed by A.G. Lipnitskii and V.N. Saveliev for atomic systems with metallic and covalent types of bonds, which exactly describes three-particle interactions. The parameters of the potential were optimized to the database of results of density-functional calculations and known experimental data for hexagonal close-packed phase of Ti. The developed potential correctly describe ω-Ti as the ground state at zero temperature. Prediction of the point-defect energies, stacking fault energies, surface energies, specific heat, and thermal expansion is in reasonable agreement with the experimental and density-functional data. This indicates the transferability of the potential to describe phases of titanium in a wide range of temperatures. Melting point of β-Ti predicted using the developed potential is in excellent agreement with the experimental value. For the developed potential the lowest temperature of the mechanical stability of the β phase equals to 1156 K, at which the spontaneous transition to the α phase occurs. The volume of the transformation from β to α phase is correctly predicted by the developed potential in agreement with the experimental value. MD simulations in combination with Gibbs-Helmholtz integration show that the α phase is mechanically and thermodynamically stable at any temperature up to melting temperature of the β-Ti. The constructed potential can be applied for modeling the β phase at its range of mechanical stability as well as α phase at temperatures up to melting point of the β-Ti.