The Activation Energy Temperature Dependence for Viscous Flow of Chalcogenides

: For some chalcogenide glasses, the temperature dependence of the activation energy E(T) of viscous flow in the glass transition region was calculated using the Williams–Landel–Ferry (WLF) equation. A method for determining the activation energy of viscous flow as a function of temperature is proposed using the Taylor expansion of the function E ( T ) using the example of chalcogenide glasses As-Se, Ge-Se, Sb-Ge-Se, P-Se, and AsSe-TlSe. The calculation results showed that the temperature dependence of the activation energy for the Ge-Se, As-Se, P-Se, AsSe-TlSe, and AsSe systems is satisfactorily described by a polynomial of the second degree, and for Sb-Ge-Se glass by a polynomial of the third degree. The purpose of this work is to compare the values of the coefficients obtained from the Taylor series expansion of E ( T ) with the characteristics of the E ( T ) versus ( T − T g ) curves obtained directly from the experimental temperature dependence of viscosity. The nature of the dependence E ( T ) is briefly discussed.


Introduction
Viscosity arises because of a transfer of momentum between material layers moving at different velocities and quantifies the resistance of the material to flow with the Frenkel equation commonly used for the temperature dependence of viscosity [1][2][3][4][5]: where η 0 is the pre-exponential coefficient, E is the activation energy of viscous flow, T is the temperature, and R is the universal gas constant.For glass-forming melts, the activation energy of viscosity is a well-defined constant only asymptotically at high temperatures at which the viscosities are low, i.e., log(η/Poise) < 3 and E = E L = 80-300 kJ/mol [2], whereas the pre-exponential coefficient η 0 was shown to be the high-temperature viscosity limit within logarithmic approximation [4].In the glass transition range of temperatures T~T g , where T g is the glass transition temperature, the activation energy of the viscous flow of glass-forming melts has a strong temperature dependence E = E(T) increasing on the decrease in temperature from its lowest value E L in the melts (T >> T g ) to its highest value in the glassy state (T < T g ) E H = 400-800 kJ/mol [2].Following [3][4][5][6], we observe that the viscous flow in glass-forming melts exhibits three temperature ranges, as illustrated by Figure 1.
Namely, the temperature regimes of viscosity are as follows: 1. Low temperatures T < Tg: isostructural Arrhenius-type law viscosity of the glass (E = EH); 2. Intermediate temperatures Tg < T < TA: non-Arrhenius-type law formally expressed with an exponent E with variable activation energy of viscosity E = E(T); 3. High temperatures T > TA: Arrhenius-type law viscosity (E = EL).
The crossover temperature TA is assumed to be close to the liquidus temperature Tliq, at which the material is completely melted [5].It was recently shown [7] that regardless of the type of glass-forming liquid, the crossover temperature is given by the universal Equation (2): where Tm is the melting temperature and k = 1.1 ± 0.15 (see for details Figure 3b of reference [7]).In addition, the TA of certain glass families, such as float and nuclear waste glasses, can be defined using a fixed viscosity value which is independent of composition [8]. Figure 1 also provides the universal viscosity equation [3], which is valid at all temperatures although it requires at least four parameters to be fitted to experimental data-the so-called Douglas-Doremus-Ojovan (DDO) double-exponential model (see Chapter 8 of Ref. [8]).This work is devoted to estimating the temperature dependence of the activation energy of viscous flow for several chalcogenide compositions within the temperature range (2): Tg < T < TA.Chalcogenide vitreous materials are used in the development of nonvolatile memory techniques both for selectors and storing data in memory cells.Various compositions are used in infrared and fiber optics, sensors and devices utilizing the properties of these glasses.Research of such materials has great prospects in electronics and optics.We have selected chalcogenide glasses as the main objects of the study because for these materials there are experimental data on viscosity available in a wide temperature range.For them, the expansion of the activation energy function E(T) into a Taylor series is acceptable.The method proposed can be expanded to include other glass-forming compositions.
Low temperatures T < T g : isostructural Arrhenius-type law viscosity of the glass (E = E H ); 2.
Intermediate temperatures T g < T < T A : non-Arrhenius-type law formally expressed with an exponent E with variable activation energy of viscosity E = E(T); 3.
High temperatures T > T A : Arrhenius-type law viscosity (E = E L ).
The crossover temperature T A is assumed to be close to the liquidus temperature T liq , at which the material is completely melted [5].It was recently shown [7] that regardless of the type of glass-forming liquid, the crossover temperature is given by the universal Equation (2): where T m is the melting temperature and k = 1.1 ± 0.15 (see for details Figure 3b of reference [7]).In addition, the T A of certain glass families, such as float and nuclear waste glasses, can be defined using a fixed viscosity value which is independent of composition [8].
Figure 1 also provides the universal viscosity equation [3], which is valid at all temperatures although it requires at least four parameters to be fitted to experimental data-the so-called Douglas-Doremus-Ojovan (DDO) double-exponential model (see Chapter 8 of Ref. [8]).This work is devoted to estimating the temperature dependence of the activation energy of viscous flow for several chalcogenide compositions within the temperature range (2): T g < T < T A .Chalcogenide vitreous materials are used in the development of nonvolatile memory techniques both for selectors and storing data in memory cells.Various compositions are used in infrared and fiber optics, sensors and devices utilizing the properties of these glasses.Research of such materials has great prospects in electronics and optics.We have selected chalcogenide glasses as the main objects of the study because for these materials there are experimental data on viscosity available in a wide temperature range.For them, the expansion of the activation energy function E(T) into a Taylor series is acceptable.The method proposed can be expanded to include other glass-forming compositions.

Modeling the Viscosity
Many useful models of the viscosity of amorphous materials have been proposed for the temperature range (2) (see the overviews [1][2][3][4][5]), among which we can note the two most developed and frequently used models, namely the Vogel-Fulcher-Tammann (VFT) equation [9][10][11], and the Williams-Landel-Ferry (WLF) equation [12,13]: Appl.Sci.2024, 14, 4319 where η 0 and B 0 are temperature-independent constants and material-dependent parameters, a T is the relative viscosity a T = η(T)/η T g , C 1 and C 2 are empirical temperatureindependent constants, T r is the reference temperature typically taken as the glass transition temperature T r = T g .The validity of these equations has been demonstrated in many works for various glass-forming systems, e.g., see Refs.[1][2][3][4][5][14][15][16][17][18][19][20].We note that the activation energy of viscous flow starts to decrease exactly at temperatures exceeding T g , see e.g., Figure 1 of Ref. [5] and Figures 1, 2, and 4 of Ref. [19].The WLF equation is known to correctly describe the molecular kinetic processes in the glass transition region and provide fairly exact temperature dependence of relaxation time τ(T) and the viscosity η(T) within the glass transition range of temperatures for amorphous materials of various natures such as polymers, organic, and inorganic glasses, and metallic amorphous alloys [14][15][16][17][18][19][20].It is known that the WLF equation can be derived in [16] without specifying the functional dependence η(T) based on the series expansion of function ln(η(T)) near the T g with respect to the small dimensionless parameter λ = (T − T g )/T g << 1.This led to following Equation (5) in terms of natural logarithms (ln) rather than decimal (log) logarithms while assuming T r = T g : From here, the WLF equation constants C 1 and C 2 acquire a physical meaning, expressed through the derivative parameters A and B: Additionally, we observe that the VFT and WLF equations are equivalent if we suppose that: Equating Formulas (1) and (3) and taking into account that η 0 = const, we obtain the following expression for determining the activation energy of viscous flow on temperature: Accounting for (7) we can rewrite it using (5) as: We conclude from (9) that in order to calculate the variable activation energy of viscous flow within the temperature range (2) it is necessary to know the values of three parameters: C 1 , C 2 , and T g .From Equation (9) it follows the following expression for determining the activation energy at the liquid-glass transition E g = E(T g ): The validity and usefulness of Equations ( 9) and (10) were checked for several chalcogenide materials, as described below.
The resulting Equations ( 9) and ( 10) can be effectively used to predict E g .The calculation results show a good agreement with data obtained using relaxation spectrometry methods.For example, it was found that the activation energy at the glass transition temperature of silicate glass is 241 kJ/mol according to Formula (10), and E g = 242 kJ/mol according to relaxation spectroscopy data [21].

Materials and Methods
We have analyzed the activation energies of the viscous flow within the temperature range (2) for the following three families of glasses: As-Se, Sb-Ge-Se, and Ge-Se.Figures 2-4 show the temperature dependence of the activation energy E(T).Theoretical dependences were calculated using Equation (9), whereas the points in the plots represent experimental data taken from [20][21][22][23][24][25].

Materials and Methods
We have analyzed the activation energies of the viscous flow within the temperature range (2) for the following three families of glasses: As-Se, Sb-Ge-Se, and Ge-Se.Figures 2-4 show the temperature dependence of the activation energy E(T).Theoretical dependences were calculated using Equation (9), whereas the points in the plots represent experimental data taken from [20][21][22][23][24][25].  [20], curve calculation using Formula (9), Tg = 316 K.

Materials and Methods
We have analyzed the activation energies of the viscous flow within the temperature range (2) for the following three families of glasses: As-Se, Sb-Ge-Se, and Ge-Se.Figures 2-4 show the temperature dependence of the activation energy E(T).Theoretical dependences were calculated using Equation (9), whereas the points in the plots represent experimental data taken from [20][21][22][23][24][25].  [20], curve calculation using Formula (9), Tg = 316 K.
The dots in these figures indicate the values obtained from the logarithms of viscosity log(η(T)) according to reference [19] and using the equation: Figures 2-4 reveal that the experimental points practically coincide with those calculated, which is confirmed by the low values of the Pearson chi-square divergence criterion.This confirms the good applicability of the WLF equation.Table 1 presents data for two empirical coefficients of the WLF equation and the glass transition temperature for the analyzed compositions, as well as the activation energy of the liquid-glass transition at T g , calculated using Formula (10).The viscosity of the chalcogenide systems under consider-ation at various temperatures was measured by the penetration method with error on a logarithmic scale ± 0.05 [20][21][22][23][24][25].Ge content-92 mol.%,Se-8 mol.%.Points and experimental data from [22], curve calculation using Formula (9), Tg = 352 K.
The dots in these figures indicate the values obtained from the logarithms of viscosity log(η(T)) according to reference [19] and using the equation: Figures 2-4 reveal that the experimental points practically coincide with those calculated, which is confirmed by the low values of the Pearson chi-square divergence criterion.This confirms the good applicability of the WLF equation.Table 1 presents data for two empirical coefficients of the WLF equation and the glass transition temperature for the analyzed compositions, as well as the activation energy of the liquid-glass transition at Tg, calculated using Formula (10).The viscosity of the chalcogenide systems under consideration at various temperatures was measured by the penetration method with error on a logarithmic scale ± 0.05 [20][21][22][23][24][25].Ge content-92 mol.%,Se-8 mol.%.Points and experimental data from [22], curve calculation using Formula (9), T g = 352 K.

Theoretical
Let us expand the temperature-dependent activation energy of viscous flow E(T) near the T g into a Taylor series in terms of a small temperature increment with a small parameter available in the system λ = (T − T g )/T g << 1, limiting the expansion to the first three terms of the series [24]: The small parameter in expansion ( 12) is the dimensionless quantity λ, but T g to the first and second powers is hidden in the coefficients D 1 and D 2 .From Equation (1), it follows that: According to Equation ( 5), for parameter A, we obtain the relation: where we accounted for the definition of D 1 from (12).The parameter A is given in Equations ( 5) and ( 6) by: Equations ( 14) and ( 15) in turn give the following equality: which considering E g = C 1 RT g transforms into the following equation: We find the formula for coefficient D 2 from the calculations as follows: 1 2 The coefficient B is given accordingly by: Further, in this relation for D 1 , we consider expression (16), resulting in: From ( 6) and ( 15), we have equality B = C 1 (T g /C 2 ) 2 , which, on substitution into (20), finally gives D 2 : If we go a slightly different route and use a joint transformation of expressions ( 12) and ( 1) with the concomitant finding of differentials A and B, then we can obtain equivalent data.Thus, the activation energy of viscous flow is represented in the form of a Taylor series in the temperature range (2) in the form: The following assumption holds: in a given series, it is possible to limit ourselves to only the linear term if it is much larger than the next one.This is possible if the following condition is met: Thus, based on three known quantities (C 1 , C 2 , and T g ), and using Equations ( 12), (18), and (22), it is possible to calculate the dependence of the activation energy of viscous flow on the temperature in the temperature range (2).Table 2 contains the coefficients D 1 and D 2 for the chalcogenide glasses analyzed, which are compared with the characteristics of the activation energy E(T) versus (T − T g ) plots obtained directly from the experimentally known temperature dependence of viscosity [20][21][22][23][24][25].The parameters of the function E(T) versus (T − T g ) expressed as a polynomial of the second degree were calculated using the Microsoft Excel spreadsheet software, with results shown in Figures 5-7, where the dots indicate the dependence of the activation energy on the temperature difference, obtained according to Formula (1) using experimental data from [20][21][22][23][24][25].

AsSe
The parameters of the function E(T) versus (T − Tg) expressed as a polynomial of the second degree were calculated using the Microsoft Excel spreadsheet software, with results shown in Figures 5-7, where the dots indicate the dependence of the activation energy on the temperature difference, obtained according to Formula (1) using experimental data from [20][21][22][23][24][25].
The parameters of the function E(T) versus (T − Tg) expressed as a polynomial of the second degree were calculated using the Microsoft Excel spreadsheet software, with results shown in Figures 5-7, where the dots indicate the dependence of the activation energy on the temperature difference, obtained according to Formula (1) using experimental data from [20][21][22][23][24][25].In a sufficiently large temperature range, however, the contribution of the third-degree term in the Taylor series may increase, which affects the estimate of the effective value of D 2 .In this regard, in the calculations of the required energy, a polynomial was used to the appropriate extent.Nevertheless, taking into account the application of this assumption, the data obtained turned out to conform well with the experiment.
Thus, we conclude that the temperature dependence of the activation energy in a wide temperature range for chalcogenide glasses Ge-Se, As-Se, P-Se, As-Se-Tl-Se, and As-Se In a sufficiently large temperature range, however, the contribution of the third-degree term in the Taylor series may increase, which affects the estimate of the effective value of D2.In this regard, in the calculations of the required energy, a polynomial was used to the appropriate extent.Nevertheless, taking into account the application of this assumption, the data obtained turned out to conform well with the experiment.
Thus, we conclude that the temperature dependence of the activation energy in a wide temperature range for chalcogenide glasses Ge-Se, As-Se, P-Se, As-Se-Tl-Se, and As-Se is satisfactorily described by a second-degree polynomial, and for Sb-Ge-Se glass by a third-degree polynomial dependence.

Discussion
As research shows, many substances with disordered structures and initially different natures, such as different classes of glasses, polymers, metal alloys, etc., have similar mechanisms of transitioning from a liquid to a glassy state.This fact is also justified by the discovery in the region of the liquid-glass transition of general rules suitable for different classes of substances [1][2][3][4][14][15][16][17][18][19][26][27][28][29].A significant increase in the activation energy of viscous flow with decreasing temperature is a consequence of a change in the configuration of the structure of the melt, which passes into the glassy state [3,30,31].
It is known [3] that the activation energy of a viscous flow consists of a pair of components due to the processes of breaking bonds (formation of defects according to Mott [32]) and due to the movement of broken bonds, called configurons [33].This division of the above energy into a pair of components was first proposed by Ya.I. Frenkel [1].It was then further developed in the works of Duglas [34], Angell and Rao [28], Filipovich [30], Nemilov [31], Doremus [35], and later by the Sheffield School of Glass [36], also known as the DDO model [8].In a more rigorous approach, E(T) should be understood as the free energy of activation, as is customary in the classical works of Eyring [37].
The developed model of delocalized atoms [37][38][39][40] assumes that for melts transforming into a glassy state, the activation energy E(T) in Equation ( 1) is also decomposed into two components:

Discussion
As research shows, many substances with disordered structures and initially different natures, such as different classes of glasses, polymers, metal alloys, etc., have similar mechanisms of transitioning from a liquid to a glassy state.This fact is also justified by the discovery in the region of the liquid-glass transition of general rules suitable for different classes of substances [1][2][3][4][14][15][16][17][18][19][26][27][28][29].A significant increase in the activation energy of viscous flow with decreasing temperature is a consequence of a change in the configuration of the structure of the melt, which passes into the glassy state [3,30,31].
It is known [3] that the activation energy of a viscous flow consists of a pair of components due to the processes of breaking bonds (formation of defects according to Mott [32]) and due to the movement of broken bonds, called configurons [33].This division of the above energy into a pair of components was first proposed by Ya.I. Frenkel [1].It was then further developed in the works of Duglas [34], Angell and Rao [28], Filipovich [30], Nemilov [31], Doremus [35], and later by the Sheffield School of Glass [36], also known as the DDO model [8].In a more rigorous approach, E(T) should be understood as the free energy of activation, as is customary in the classical works of Eyring [37].
The developed model of delocalized atoms [37][38][39][40] assumes that for melts transforming into a glassy state, the activation energy E(T) in Equation ( 1) is also decomposed into two components: where E ∞ is the potential for the transition of a structural kinetic unit (atom, molecule) to a new local position, and E S (T) is the potential for a configurational change of structure, which is a function of temperature [28,38,39]: Here ∆ε e is the delocalization energy of the atom, characterizing the critical deviation of the kinetic unit from the point with the minimum energy.Such a deviation of a kinetic unit is possible with a spatial change in its immediate environment.This process can be considered as a transformation of the structural configuration at the local level.The energy of change in the structure configuration E S (T) tends to zero in on the case of increasing temperature when RT >> ∆ε e .Then, expression (1), taking into account Equations ( 24) and ( 25), takes the form of the Arrhenius equation with one potential E ∞ .
Thus, the component E∞ can have the value of the activation energy of the viscous flow at a temperature tending to infinity; in this case E ∞ = E H .When the temperature decreases to values around T g , the energy begins to decrease exponentially exp(−∆ε e /kT) [27,38].As a result of such a sharp decrease, ∆ε e becomes comparable to the energy of thermal vibrations of the lattice (~3 RT) and the relative number of delocalized atoms (N e /N).As a result of a sharp decrease in the delocalization energy, the structure of the glassy substance (melt) begins to compact, and rearrangement of the kinetic unit will require a local change in the configuration in its vicinity, which leads to a sharp increase in the potential for changing the configuration of the structure E S (T).In this way, we can explain the exponential increase in the activation energy of the melt in the glass transition region.The model noted in the work explains how a change in the relative position of the kinetic units of the structure affects the dynamics of the activation energy values when the temperature changes in the range close to T g .

Conclusions
The article describes the possibilities of determining the activation energy E(T) at various temperatures for the chalcogenide glasses under study.To analyze the dependence of E(T) at temperatures close to the glass transition temperature, the WLF equation can be used as a first approximation with sufficient reliability.A method is shown for expanding the above dependence into a Taylor series in the range near T g .The calculation results showed that the temperature dependence of the activation energy for the Ge-Se, As-Se, P-Se, AsSe-TlSe, and AsSe systems is satisfactorily described by a polynomial of the second degree, and for Sb-Ge-Se glass by a polynomial of the third degree.For Sb-Ge-Se systems, the values of coefficients D 1 and D 2 at corresponding temperatures are in the range from 0.135 to 0.575 kJ/mol•K and from 0.0007 to 0.0087 kJ/mol•K 2 , respectively.For other studied compositions, the values of D 1 and D 2 are approximately in the same ranges.Using the model described in this article, it is shown that the exponential increase in activation energy at T g follows from an increase in the potential for restructuring the structure during the transition from the liquid to the glassy state.

Figure 2 .
Figure 2. Temperature dependence of the activation energy of viscous flow E(T) of As-Se molten glass.The content of As is 3 mol.%,Se-97 mol.%.Points and experimental data from[20], curve calculation using Formula (9), Tg = 316 K.

Figure 2 .
Figure 2. Temperature dependence of the activation energy of viscous flow E(T) of As-Se molten glass.The content of As is 3 mol.%,Se-97 mol.%.Points and experimental data from[20], curve calculation using Formula (9), T g = 316 K.

Figure 2 .
Figure 2. Temperature dependence of the activation energy of viscous flow E(T) of As-Se molten glass.The content of As is 3 mol.%,Se-97 mol.%.Points and experimental data from[20], curve calculation using Formula (9), Tg = 316 K.

Figure 5 .
Figure 5. Results of processing the temperature dependence of the activation energy of viscous flow for As-Se molten glass.The content of As is 3 mol.%,Se-97 mol.%.Points and experimental data taken from [20], curve polynomials.

Figure 6 .Figure 5 .
Figure 6.Results of processing the temperature dependence of the activation energy of viscous flow for Ge-Se molten glass.The content of Ge is 8 mol.%,Se-92 mol.%.Points and experimental data from [22], curve polynomials.

Figure 5 .
Figure 5. Results of processing the temperature dependence of the activation energy of viscous flow for As-Se molten glass.The content of As is 3 mol.%,Se-97 mol.%.Points and experimental data taken from [20], curve polynomials.

Figure 6 .Figure 6 .
Figure 6.Results of processing the temperature dependence of the activation energy of viscous flow for Ge-Se molten glass.The content of Ge is 8 mol.%,Se-92 mol.%.Points and experimental data from [22], curve polynomials.

Figure 7 .
Figure 7. Results of processing the temperature dependence of the activation energy of viscous flow for Sb-Ge-Se molten glass.The content of Sb is 20 mol.%, Ge-10 mol.%,Se-70 mol.%.Points and experimental data from [24], curve polynomials.

Figure 7 .
Figure 7. Results of processing the temperature dependence of the activation energy of viscous flow for Sb-Ge-Se molten glass.The content of Sb is 20 mol.%, Ge-10 mol.%,Se-70 mol.%.Points and experimental data from [24], curve polynomials.