ISSN: 2641-6921
Zakarian D* and Khachatrian A
Received: October 26, 2020; Published: November 19, 2020
*Corresponding author: Zakarian D, Frantzevich Institute for Problems of Materials Science NAS, Ukraine
DOI: 10.32474/MAMS.2020.03.000165
A method has been developed for calculating the melting temperature of materials from first principles based on the temperature dependence of the theoretical strength. An analytical regularity of the temperature dependence of the theoretical strength of materials (metals, multielement metal alloys, borides, diamond-like covalent and ion-covalent crystals, quasi-binary eutectic systems) is obtained. The resulting ratio allows one to estimate the theoretical strength of materials in the temperature range from zero to the melting point.
Keywords: Theoretical strength; temperature dependence; melting point; pseudopotential method; quasi-harmonic approximation
There is a scatter in the data on the melting temperature of refractory materials even in the case of sufficiently pure materials [1]. Often tabular data differ from each other by 300K or more (for example, in the case of TiB_{2}). Therefore, it is relevant to clarify the melting point from first principles, without using of experimental data. In this work, the method of a priori pseudopotentials [2-7] is used as the most suitable for the conditions of ab initio and computational speed. To calculate the physical characteristics, it is necessary to have the values of the energy of the electron-ion system of the material depending on the parameter of the crystal lattice or its volume. The total energy of the material can be represented as the sum U_{0} of the energies of the electron-ion system at a temperature T=0 and U_{T}, the energy of thermal vibrations of ions at a temperature T ≠ 0
U =U_{0} +U_{T} (1)
U_{0} is calculated from first principles in the second order of perturbation theory, and the energy of thermal vibrations can be considered using one of the approximate methods - Debye or Einstein. Debye's model works better at low temperatures, while Einstein's model describes well the heat capacity of crystals at room and higher temperatures [2]. For problems concerning the problem of superconductivity, only Debye's method is suitable. When studying the temperature dependence of the mechanical characteristics of crystals, it is more convenient to use the Einstein model, especially when it comes to systems with a complex structure (multielement metal alloys and transition metal borides). According to Einstein's model, all atoms vibrate at the same frequency, which is proportional to the elastic characteristics of the material. Elastic characteristics are determined using the pseudopotential method. When calculating the energy of thermal vibrations (T ≠ 0), we use the quasi-harmonic approximation in combination with the method of a priori pseudopotential; as a result, we obtain the dependence of the energy of the electron-ion system on temperature through the volume of the unit cell
U =U (Ω (T) [4]
The developed model in the quasi-harmonic approximation allows, within the framework of the harmonic approximation, to calculate the parameters of the crystal lattice during thermal expansion, and then the corresponding physical and mechanical characteristics of materials. Strength and elastic characteristics represent the derivatives of the total energy with respect to the lattice parameter. To calculate the melting point of materials, we solve the inverse problem. Using the example of materials with a verified melting point, we determine the dependence of the theoretical strength on temperature using the pseudopotential method and the quasi-harmonic approximation [4]. The results of such a calculation are presented in tabular and graphical (Figure 1) form. Table 1 shows the results of calculations for representatives of the class of metals - metal alloys, borides, diamond-like materials, and eutectic systems. As the results of calculations show, for all crystalline materials (borides, carbides, metals, alloys, multielement alloys), the theoretical strength has a maximum value at temperature T = 0 K and decreases with increasing temperature. At maximum temperatures for all studied materials, the strength is 84.59 - 85.4% of the maximum strength.
Table 1: Temperature dependence of the theoretical strength of materials σ(T), calculated from first principles, tabular values of the melting temperature T_{max}, K.
In [1], it was experimentally proved that the high strength of boride and cermet eutectic systems is preserved up to a temperature of 0.8 T_{max} (T_{max} is the melting point), which was theoretically confirmed in our calculations from first principles [2,3]. (Figure 2) based on the calculation results (Table 1), shows the dependence of relative strength σ(T)/σ_{0} on the parameter T/T_{max} for different materials for which the experimental values of the melting temperature (T_{max}) are known. According to mathematical analysis, the dependence of σ(T)/σ_{0} on T/T_{max} can be approximated by a function of the type:
σ (T) / σ _{0} = 1 − 0,123T / T_{max} − 0, 027 (T / T_{max})^{2}
Here σ_{0} is the strength at T = 0 K, σ (T) is the strength at the current temperature. At the melting point (T = T_{max}), the strength is σ(T_{max}) ≈ 0,85 σ_{0}.
Note that the T_{max} value is not used in the process of calculating the temperature dependence of strength from first principles. This means that with an unknown value of the melting point of the material, it is necessary to find the temperature dependence of the strength in the range (σ_{0} – 0,85∙ σ_{0}). The corresponding temperature Т_{х}, at which σ (Т_{х}) ≈ 0,85 σ_{0} will be the sought melting temperature Т_{х} ≈ Т_{max}. According to the data presented in table. for TiB_{2} σ (Т3500) = 0.8492 σ_{0}, for the eutectic system 2500 < Т_{max} <2750 K, etc. As a result, formula (2) was obtained, which describes the temperature dependence of the theoretical strength of materials. The dependence of the relative theoretical strength σ(T)/σ_{0 }on T/T_{max} has the same form for different classes of materials. Knowing the melting point and the value of the theoretical strength at T = 0, using formula (2), makes it possible to calculate the strength of the material at any temperature. As for the theoretical strength σ_{0}, it can be calculated from first principles [5-7].
Bio chemistry
University of Texas Medical Branch, USADepartment of Criminal Justice
Liberty University, USADepartment of Psychiatry
University of Kentucky, USADepartment of Medicine
Gally International Biomedical Research & Consulting LLC, USADepartment of Urbanisation and Agricultural
Montreal university, USAOral & Maxillofacial Pathology
New York University, USAGastroenterology and Hepatology
University of Alabama, UKDepartment of Medicine
Universities of Bradford, UKOncology
Circulogene Theranostics, EnglandRadiation Chemistry
National University of Mexico, USAAnalytical Chemistry
Wentworth Institute of Technology, USAMinimally Invasive Surgery
Mercer University school of Medicine, USAPediatric Dentistry
University of Athens , GreeceThe annual scholar awards from Lupine Publishers honor a selected number Read More...
The annual scholar awards from Lupine Publishers honor a selected number read more...