CALCULATION OF THE PHASES COEXISTENCE SPACES IN THE SYSTEM Hg-Mn-Te-Se.
Introduction. Four-component solid solutions based on compounds of the A2B6 type are a good choice for creating optoelectronic devices, since they operate in a wide spectral range. But there is a possibility of loss of thermodynamic stability in solid phases in II-IV semiconductor four-component systems, which can lead to the appearance of unstable states of solid solutions. As a result, critical phenomena may exist in unstable solid phases, which lead to degradation of the properties of optoelectronic devices [1-5].
An important stage in the development of electronic technologies and solving problems related to the behavior of multicomponent systems was the introduction of modern methods of computer modeling into this field. This allows us to analyze thermodynamic models of multicomponent systems and predict their behavior.
In particular, computer modeling makes it possible to calculate and construct the phase diagrams necessary for the analysis of the processes of obtaining self-organized formed nanostructures, and also takes into account the possibility of the existence of bifurcation spaces and spaces of coexistence of phases of various orders.
Objective. The aim of the paper is to calculate, on the basis of the differential topological approach [6, 7], the higher free energy derivatives of the system from the second through the eighth inclusive in four concentration parameters for analyzing the four-component solid solution [Hg.sub.1-x][Mn.sub.x][Te.sub.1-y][Se.sub.y] for the existence of second-order phase coexistence spaces.
Thermodynamic modeling of critical phenomena in the system Hg - Mn - Te - Se.
When analyzing the four-component solid solution [Hg.sub.1-x][Mn.sub.x][Te.sub.1-y][Se.sub.y] as a system in which several phases are assumed, the condition for the existence of two phases that are in equilibrium is used [6, 7], by means of which the existence of second-order phase coexistence spaces is verified:
[d.sup.2]G/d[x.sup2] = [d.sup.3]G/d[x.sup.3] = 0; [d.sup.4]G/d[x.sup.4] > 0. (1)
For the analysis of a multicomponent system, a thermodynamic model was used in which the free Gibbs energy for a four-component solution was considered as a mixture of four hypothetical binary compounds in the form :
G = G([X.sub.AC] [X.sub.BC] [X.sub.AD] [X.sub.BD]), (2)
where G--free energy of the system, [X.sub.j]--the concentration of binary components.
The state functions of such a solid solution can be described using the indicated concentrations [X.sub.ij]. The basis of the regular solution model is the assumption of the random distribution of heterogeneous atoms that form a solid substitutional solution along the sites of the corresponding sublattices. Accordingly, the expression for the concentrations of the binary components takes the form:
[X.sub.AC] = (1 - x)(1 - y), [X.sub.BC] = x(1 - y), [X.sub.AD] = (1 - x)y, [X.sub.BD] = xy, (3)
where x and y--concentration parameters.
To assess the possibility of the emergence of critical spaces and coexistence spaces in the solid solution under consideration, analytical expressions for higher derivatives of the potential energy of the system are obtained and investigated.
The higher derivatives for the free Gibbs energy of the solid solution [Hg.sub.1-x][Mn.sub.x][Te.sub.1-y][Se.sub.y] by the concentration parameters are calculated by matrix-vector differentiation of multidimensional systems. To obtain matrices of higher derivatives, the method of direct sums was used .
The higher derivatives of the free energy of the system under consideration were obtained by the following algorithm. In the framework of the regular solution model [6, 7], the second derivatives of the free Gibbs free energy [[mu].sub.ij] of the solid solution [Hg.sub.1-x][Mn.sub.x][Te.sub.1-y][Se.sub.y] were obtained, of which an matrix of 4x4 dimensions was formed:
[mathematical expression not reproducible], (4)
Further, the matrix [A.sub.2] was differentiated in four corresponding concentrations [X.sub.ij]. As a result of differentiation, four fourth-order partial derivatives of the third order of 4x4 were obtained. Then, using the method of direct sums, the block-diagonal matrix of the total third derivative of the free energy of the system under investigation was formed:
[mathematical expression not reproducible], (5)
The determinant of the obtained matrix [A.sub.3] is defined as the sum of the determinants of the matrices of partial derivatives located on the main diagonal:
det [A.sub.3] = det [A.sub.13] + det [A.sub.23] + det [A.sub.33] + det [A.sub.43], (6)
[A.sub.12] = d[A.sub.2]/d[X.sub.AC]; [A.sub.23] = d[A.sub.2]/d[X.sub.AD]; [A.sub.33] = d[A.sub.2]/d[X.sub.BD]; [A.sub.43] = d[A.sub.2]/d[X.sub.BC]. (7)
Preparation of higher order system free energy derived from the fourth to the eighth, and the related determinants det [A.sub.i], where i is the order of the derivative was performed using the same algorithm. Analytical expressions of higher derivatives were obtained using Maxima computer mathematics .
After obtaining analytical expressions for the determinants of the matrices of the free energy derivatives of the system under consideration, the positions of the zero contours of the considered derivatives were calculated.
For this, the values of the concentration parameters x and y in the range from 0 to 1 in 0.1 increments were substituted into the expressions for the corresponding determinants, after which the condition det [A.sub.i]=0 was checked. Further, the areas in which the values of the corresponding determinants take positive and negative values were determined. Positive values were marked by light areas, negative values--darker.
Results. By calculating the positions of the zero contours, the analytical expressions for the total derivatives of the Gibbs free energy of the four-component solid solution [Hg.sub.1-x][Mn.sub.x][Te.sub.1-y][Se.sub.y], the phase diagrams of the studied system of the position of the areas of coexistence of the derivatives with allowance for (1) are numerically determined and plotted. Also their zero contours from the second to the eighth derivative inclusive for the temperature range of 800... 1000 K were determinated.
In the simulation, the values of the quasibinary interaction parameters in solid solutions and crystal-physical parameters of binary systems were used (Table):
Parameters of the interaction in the solid phase for quasibinary systems Solid solution Interaction Parameters [Hg.sub.1-x][Mn.sub.x] [[alpha].sub.13] = [[alpha].sub.HgSe-MnSe] [Te.sub.1-y][Se.sub.y] [[alpha].sub.24] = [[alpha].sub.HuTe-MnTe] [[alpha].sub.12] = [[alpha].sub.MnTe-MnSe] [[alpha].sub.34] = [[alpha].sub.HgTe-HiSe] Solid solution Parameter value [[alpha].sub.sj-lk], J/mole [Hg.sub.1-x][Mn.sub.x] 16000 [Te.sub.1-y][Se.sub.y] 14200 900 2933
In Fig. 1-4 are graphical representations of the positions of the zero contours of the derivatives of the free Gibbs energy in the [Hg.sub.1-x][Mn.sub.x][Te.sub.1-y][Se.sub.y] system in terms of the concentration parameters, starting with the second derivative and the eighth derivative inclusive for the 800 K temperature and the second-order phase coexistence spaces (Fig. 4, b) The area of positive values of the derivative is shown in dark color and, accordingly, the lighter one is the area of negative values.
Conclusion. The considered procedure for calculating of higher derivatives of free energy with respect to four concentration parameters makes it possible to correctly estimate the position of the second-order phase coexistence spaces in multicomponent solid solutions. The constructed thermodynamic model of the solid solution [Hg.sub.1-x][Mn.sub.x][Te.sub.1-y][Se.sub.y] can also be used to predict the formation and coexistence of three and four different phases in the solid solution considered.
[1.] Sorokin V.S., Sorokin S.V., Kaygorodov V.A., & Ivanov S.V. (2000). The Instability and Immiscibility Regions in MgxZn1-xSySe1-y alloys. J. Cryst. Growth., 214/215, 130-134.
[2.] Maksimov O., Guo S.P., & Tamargo M.C. (2002). Be-Chalcogenide alloys for improved R-G-B LEDs: Be[x]Zn[y]Cd[1-x-y]Se on InP. Phys. Stat. Sol., 229(b), 2, 1005-1009.
[3.] Hsieh C.H., Huang Y.S., & Ho, C.H. et al. (2004). Temperature dependence of the band-edge transitions of ZnCdBeSe. Jpn. J. Appl. Phys., 43, 2, 459-466.
[4.] Romcevic N., Romcevic M., & Golubovic A. et al. (2005). Far-infrared and Raman spectroscopy of [Cd.sub.1-x][Mn.sub.x][Te.sub.1-y][Se.sub.y]: Phonon properties. Journal of Alloys and Compounds, 397, 52-57.
[5.] Romcevic M., Kulbachinskii V.A., & Romcevic N. et al. (2005). Optical Properties of [Hg.sub.1-x][Mn.sub.x][Te.sub.1-y][Se.sub.y]. Infrared Physics & Technology, 46, 379-387.
[6.] Kazakov A., Kvatashidze L., & Shapovalov G. (2014). Computer modeling of critical coexistance spaces at phase diagrams of multicomponent solid solution. Informatics and Mathematical Methods in Simulation, 4, 4, 349-356.
[7.] Kazakov A., & Shapovalov G. (2016). Simulation of the regions of phase coexistence in Inx[Ga.sub.1-x][As.sub.y][P.sub.1-y] solid solutions using various thermodynamic models. Technical sciences and technologies, 5, 3, 96-103.
[8.] Kornn G., & Korn T. (1974). Handbook on mathematics. Moscow: Nauka.
[9.] Computer Algebra System. Maxima. Maxima, a Computer Algebra System. Retrieved from http://maxima.sourceforge.net.
[phrase omitted]; Kazakov Anatolii, ORCID: http://orcid.org/0000-0002-9402-6556
[phrase omitted]; Burtnyi Dmitro, ORCID: http://orcid.org/0000-0001-7157-5628
[phrase omitted]; Shapovalov Hennadii, ORCID: http://orcid.org/0000-0001-7157-5627
Received Desember 11, 2017
Accepted March 29, 2018
A. Kazakov, DSc, Prof., D. Burtnyi, G. Shapovalov
Odessa National Polytechnic University, 1 Shevchenko Ave., Odessa, Ukraine, 65044; e-mail: email@example.com
Caption: Fig. 1. Graphical representations of the positions of the zero contours of the second derivative (a) and the third derivative (b) of the free energy in the Hg-Mn-Te-Se system at the cross section for the existence of solid solutions of the phase diagram
Caption: Fig. 2. Graphical images of the positions of the zero contours of the fourth derivative (a) and the fifth derivative (b) of free energy in the Hg-Mn-Te-Se system at the cross section for the existence of solid solutions of the phase diagram
Caption: Fig. 3. Graphical representations of the positions of the zero contours of the sixth derivative (a) and the seventh derivative (b) of the free energy in the Hg-Mn-Te-Se system on the cross section for the existence of solid solutions of the phase diagram
Caption: Fig. 4. Graphical representations of the positions of the zero contours of the eighth derivative (a) offree energy in the Hg-Mn-Te-Se system on the cross section for the existence of solid solutions of the phase diagram and second-order phase coexistence spaces (b, shown in a lighter color)
[Please note: Some non-Latin characters were omitted from this article]
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||ELECTRONICS. RADIO ENGINEERING. TELECOMMUNICATION FACILITIES|
|Author:||Kazakov, A.; Burtnyi, D.; Shapovalov, G.|
|Publication:||Odes'kyi Politechnichnyi Universytet. Pratsi|
|Date:||Mar 1, 2018|
|Previous Article:||INFORMATION MODEL OF DISTANCE LEARNING SYSTEM IN TERMS OF DATA COMMUNICATION IN HETEROGENEOUS INTERNET NETWORKS.|
|Next Article:||THERMODYNAMIC OF OBTAINING OF MONODISPERSE PARTICLES SI[O.sub.2] BY TETRAETHOXYSILANE HYDROLYSIS IN THE SI-O-H-C-N SYSTEM.|