On Structural Rearrangements Near the Glass Transition Temperature in Amorphous Silica

The formation of clusters was analyzed in a topologically disordered network of bonds of amorphous silica (SiO2) based on the Angell model of broken bonds termed configurons. It was shown that a fractal-dimensional configuron phase was formed in the amorphous silica above the glass transition temperature Tg. The glass transition was described in terms of the concepts of configuron percolation theory (CPT) using the Kantor-Webman theorem, which states that the rigidity threshold of an elastic percolating network is identical to the percolation threshold. The account of configuron phase formation above Tg showed that (i) the glass transition was similar in nature to the second-order phase transformations within the Ehrenfest classification and that (ii) although being reversible, it occurred differently when heating through the glass–liquid transition to that when cooling down in the liquid phase via vitrification. In contrast to typical second-order transformations, such as the formation of ferromagnetic or superconducting phases when the more ordered phase is located below the transition threshold, the configuron phase was located above it.


Introduction
Silica is the most common oxide on the Earth, characterized by a mean abundance of about 37 wt%. Amorphous silica exists either as glass while the temperature is below the glass transition temperature, which is T g = 1480 K at normal pressure, or as a melt above this temperature. Amorphous silica has a three-dimensional interconnected structure of SiO 4 tetrahedra up to about 15 GPa, with an average Si-O-Si bond angle of 144-153 • , which decreases by about 1 • per 100 K with an increase of temperature above T g at the ambient pressure [1]. The silica melt is thermodynamically stable above the melting temperature, which is T m = 1986 K, while below this temperature, both melt and glass are metastable structures. The changes in the structure that occur at T g when a metastable liquid vitrifies into a glass are of general interest, and investigations in this area persist, including works that are related to understanding the role of a modified random network in glasses and the glass transition, where silica is a relatively simple system for analysis [2][3][4][5][6].
Glass is an amorphous substance with the mechanical properties of an isotropic solid; it is found in nature and is used in various fields of human activity. However, both the nature of the glassy state and glass transition at the atomic-molecular level are not well understood [7][8][9][10][11][12][13][14][15]. At present, there is no unified approach to the physical nature of the glass transition. Moreover, the glass transition is often not considered as a thermodynamic phase transition; instead, a rule of thumb is introduced that states that an amorphous material is considered to be glass if its viscosity is equal to or higher than 10 12 Pa·s (10 13 poise) [16], although the glass transition temperature T g is in practice determined from the characteristic kink of the temperature dependences of the specific volume or enthalpy at T g [7]. Although the glass transition manifests itself as a typical secondorder phase transformation and the derivative thermodynamic parameters of amorphous materials, such as the thermal expansion coefficient or the heat capacity, experience a characteristic jump exactly at T g , its description in terms of Landau's theory is not an easy task, e.g., there is no clarity about the order parameter describing this transition.
Significant efforts are continuing to reveal the thermodynamic origins of the glass transition in various materials, including oxides, metals, organics (see, e.g., [17][18][19][20][21][22][23][24][25]), and the newly discovered Metal-organic framework (MOF) glasses [26]. Molecular dynamics (MD) experiments show that below T g , percolation clusters are formed in amorphous materials made up of Voronoi coordination polyhedra of high-density atomic configurations [27,28]. Both the analytical approach and MD simulations also show that, during the glass transition, the geometry of an amorphous material changes due to the formation of macroscopic percolation clusters [27][28][29], and percolation clusters are formed in the liquid state made up of Voronoi coordination polyhedra of low-density atomic configurations, while in the solid (glassy) state, there are no such clusters [27]. In many works, structural analysis was shown to be a powerful tool for understanding the changes that occur in an amorphous material at the glass transition [30,31], e.g., analysis of X-ray diffraction (XRD) patterns of amorphous materials shows a step-wise change of the slope of the temperature dependence of the first sharp diffraction minimum exactly at T g [25]. We note that similar methods can be used to study the structural rearrangement of other amorphous materials at the glass transition.
Glasses are isotropic and homogeneous at the macroscopic scale; however, the structure of glasses comprises the following attributes: (i) short-range order (SRO) with moleculartype units, i.e., building blocks, such as tetrahedral structures in silicates at the atomic scale; (ii) medium-range order (MRO) at a larger size range, which extends from second-and third-neighbor environments to percolating and fractal structures; and (iii) a disordered state (DS), which is homogeneous and isotropic, as observed at macroscopic sizes [22,32,33]. A diverse variety of glasses exist with a wide range of structural SRO building blocks from which the DS is composed of. The connections between the building blocks of glasses occur via the chemical bonding system. Typical oxide glasses have SRO building blocks that are connected via covalent and heteropolar (mixed covalent-ionic) bonding, whereas chalcogenide glasses are built using covalently bonded two-, three-, and four-coordinated SRO building blocks. Metallic glasses consist mainly of SRO icosahedral building blocks that are interconnected via metallic bonding [28]. Organic glasses have cross-linked polymeric molecular chains as SRO building blocks and van der Waals bonds between them. Glassy water has tetrahedral units as SRO building blocks, which are connected via hydrogen bonds.
The silica glass in normal conditions forms a fully polymerized network structure of corner-sharing tetrahedral SiO 4 units and its structure is characterized by a random distribution of chemically ordered rings, where the silicon atoms are linked by bridging oxygen atoms [32]. As such, the amorphous SiO 2 looks like a relatively simple system with covalent bonds connecting the silicon with oxygens and thus allowing for more easily revealing the changes that occur due to the increase in temperature above T g when a glass transforms into a liquid.
A liquid phase can be characterized as a dynamically homogeneous medium since there is always a tendency in it to form microscopic metastable ordered structures. Moreover, with a decrease in temperature, the sizes and lifetimes of these structures increase. On large time and space scales, the liquid is homogeneous and disordered, but locally, at the size of metastable clusters and during the lifetime of metastable formations, the liquid can have an SRO as well. Local ordering in supercooled liquids is confirmed by Fisher clusters, which are associated with long-wave density fluctuations and are revealed in glass-forming liquids and polymers [34]. With correlation lengths up to 300 nm, density fluctuations turned out to have a fractal structure with a dimensionality less than 3 [35]. Experimental observations confirmed the presence of local ordering in liquids, e.g., ordered structures were found above the liquidus of silicate melts [36]. Furthermore, direct observation of the spatial distribution of microparticles of emulsions confirmed the dynamic heterogeneities of liquids [37]. Thus, the structures of disordered materials in the liquid and solid states are significantly different. Although macroscopically, the distribution of atoms or molecules is topologically disordered both in the glassy and in the molten state, in the low-temperature (non-ergodic) region, the structure of the amorphous material is rather solid-like and the geometric structure of the interatomic (intermolecular) bonds is three-dimensional with a small fraction of defects in the form of broken chemical bonds termed configurons. Above T g , i.e., in the ergodic region, the structure is characterized by the presence of fractal dynamic molecular structures up to and even above the liquidus temperature, as well as macroscopic clusters that are made up of configurons.
The aim of this work was to analyze the structural changes at the glass-liquid transition. We focused here on amorphous SiO 2 since we could trace the topological characteristics of the covalent bonding system. We showed that the glass transition led to the formation of a configuron phase above T g that had a fractal structure and grew with the increase in temperature. The glass transition described in terms of the configuron phase formation above T g corresponded to a second-order phase transformation within the Ehrenfest classification.

Configurons
The term configuron was coined by Angell and co-authors, who introduced the congruent bond lattice (CBL) model with the aim of replacing the set of atoms with a congruent structure of weakly interacting bonds that is easier to analyze [38][39][40][41]. The CBL model is applicable to the case of any isotropic material, including crystalline substances, to understand the related transition because any condensed phase can be represented by a three-dimensional array of the interactions between its particles [41]. In the CBL, the system of weakly interacting chemical bonds congruently represents the initial system of strongly interacting ions, which are Si 4+ and O 2− ions in the amorphous SiO 2 . Figure 1 illustrates the introduction of CBL for the glassy A 2 O 3 system. can have an SRO as well. Local ordering in supercooled liquids is confirmed by Fisher clusters, which are associated with long-wave density fluctuations and are revealed in glass-forming liquids and polymers [34]. With correlation lengths up to 300 nm, density fluctuations turned out to have a fractal structure with a dimensionality less than 3 [35]. Experimental observations confirmed the presence of local ordering in liquids, e.g., ordered structures were found above the liquidus of silicate melts [36]. Furthermore, direct observation of the spatial distribution of microparticles of emulsions confirmed the dynamic heterogeneities of liquids [37]. Thus, the structures of disordered materials in the liquid and solid states are significantly different. Although macroscopically, the distribution of atoms or molecules is topologically disordered both in the glassy and in the molten state, in the low-temperature (non-ergodic) region, the structure of the amorphous material is rather solid-like and the geometric structure of the interatomic (intermolecular) bonds is three-dimensional with a small fraction of defects in the form of broken chemical bonds termed configurons. Above Tg, i.e., in the ergodic region, the structure is characterized by the presence of fractal dynamic molecular structures up to and even above the liquidus temperature, as well as macroscopic clusters that are made up of configurons. The aim of this work was to analyze the structural changes at the glass-liquid transition. We focused here on amorphous SiO2 since we could trace the topological characteristics of the covalent bonding system. We showed that the glass transition led to the formation of a configuron phase above Tg that had a fractal structure and grew with the increase in temperature. The glass transition described in terms of the configuron phase formation above Tg corresponded to a second-order phase transformation within the Ehrenfest classification.

Configurons
The term configuron was coined by Angell and co-authors, who introduced the congruent bond lattice (CBL) model with the aim of replacing the set of atoms with a congruent structure of weakly interacting bonds that is easier to analyze [38][39][40][41]. The CBL model is applicable to the case of any isotropic material, including crystalline substances, to understand the related transition because any condensed phase can be represented by a three-dimensional array of the interactions between its particles [41]. In the CBL, the system of weakly interacting chemical bonds congruently represents the initial system of strongly interacting ions, which are Si 4+ and O 2− ions in the amorphous SiO2. Figure 1 illustrates the introduction of CBL for the glassy A2O3 system. At finite temperatures T > 0 some of the chemical bonds are broken; moreover, the higher the T, the more bonds are broken. Each broken bond, along with the associated strain-releasing local adjustment of centers of atomic vibrations, is treated following Angell as an elementary excitation termed a configuron. Figure 2 illustrates the formation of a configuron in an A2O3 glass. At finite temperatures T > 0 some of the chemical bonds are broken; moreover, the higher the T, the more bonds are broken. Each broken bond, along with the associated strain-releasing local adjustment of centers of atomic vibrations, is treated following Angell as an elementary excitation termed a configuron. Figure 2 illustrates the formation of a configuron in an A 2 O 3 glass. The sizes dc of newly formed configurons, which are assumed to be spherical, are not necessarily equal to the initial bond lengths d, which is 1.62 Å in silica. They are typically slightly larger dc > d; however, they can be of the same size as the initial bond, as well as slightly smaller dc < d. It is important that the notion of a configuron is not limited to only well-localized bonds, such as those of ionic or covalent type, and can be used for metallically bonded materials, i.e., for metallic glasses as well (Table 1). Indeed, when atomic vibrational motions become anharmonic enough, the probability is high that groups of particles in an amorphous material will have simultaneously moved far enough from their original positions that it becomes energetically unfavorable for them to return, which corresponds to crossing of a local potential barrier. This motion affecting neighboring particles as concomitant shifts, which occur in the centers of vibration of other particles in the vicinity such that adjustments to the primary motion spread out through the material. As the new configuration will generally differ in energy from the initial one, the difference corresponds to a configuron that is an excited configurational state for the system. Phonon momentum but not energy will be conserved in such processes [38]. In the case of metallic glasses, the configuron as an elementary excitation is based on Egami's ideas of changes of the local atomic connectivity by losing or gaining one nearest neighbor [42], see Figure 3. The sizes d c of newly formed configurons, which are assumed to be spherical, are not necessarily equal to the initial bond lengths d, which is 1.62 Å in silica. They are typically slightly larger d c > d; however, they can be of the same size as the initial bond, as well as slightly smaller d c < d. It is important that the notion of a configuron is not limited to only well-localized bonds, such as those of ionic or covalent type, and can be used for metallically bonded materials, i.e., for metallic glasses as well (Table 1). Indeed, when atomic vibrational motions become anharmonic enough, the probability is high that groups of particles in an amorphous material will have simultaneously moved far enough from their original positions that it becomes energetically unfavorable for them to return, which corresponds to crossing of a local potential barrier. This motion affecting neighboring particles as concomitant shifts, which occur in the centers of vibration of other particles in the vicinity such that adjustments to the primary motion spread out through the material. As the new configuration will generally differ in energy from the initial one, the difference corresponds to a configuron that is an excited configurational state for the system. Phonon momentum but not energy will be conserved in such processes [38]. In the case of metallic glasses, the configuron as an elementary excitation is based on Egami's ideas of changes of the local atomic connectivity by losing or gaining one nearest neighbor [42], see Figure 3.
The configuron that is formed does not support particles that are bound to each other; therefore, instead of an elastic bond, the rigid connection between particles disappears and the previously bound particles can move away. With an increase in temperature, more and more bonds are broken until percolation via broken bonds occurs. As shown by Kantor and Webman [43], the rigidity threshold of an elastic percolating network is identical to the percolation threshold. Thus, the formation of the percolation cluster made of configurons results in a kind of disintegration loss of rigidity or, in other words, to a transformation of a solid to a liquid; hence, the glass transition occurs and the glass, which is a rigid elastic solid, transforms to a plastic liquid. The transformation of a glass into a liquid, i.e., the glass transition, is described based on the CBL model, which was specifically designed to treat the glass transition phenomena. This is done using the general methodology of percolation theory, i.e., the concept of configuron percolation theory (CPT). The configuron that is formed does not support particles that are bound to each other; therefore, instead of an elastic bond, the rigid connection between particles disappears and the previously bound particles can move away. With an increase in temperature, more and more bonds are broken until percolation via broken bonds occurs. As shown by Kantor and Webman [43], the rigidity threshold of an elastic percolating network is identical to the percolation threshold. Thus, the formation of the percolation cluster made of configurons results in a kind of disintegration loss of rigidity or, in other words, to a transformation of a solid to a liquid; hence, the glass transition occurs and the glass, which is a rigid elastic solid, transforms to a plastic liquid. The transformation of a glass into a liquid, i.e., the glass transition, is described based on the CBL model, which was specifically designed to treat the glass transition phenomena. This is done using the general methodology of percolation theory, i.e., the concept of configuron percolation theory (CPT).

Glass transition in Silica
To consider the amorphous SiO2 structure, we used the CBL where instead of the vitrification of molten material, we started with its solid state, i.e., the glass, and considered its melting, tracing the changes in the spatial distribution of configurons due to a gradual increase in temperature from zero to temperatures above Tg. The bond system forms a topologically disordered lattice that is congruent to the disordered lattice of the amorphous SiO2. Assuming only two allowed states for bonds, namely, binding (ground) and broken (excited), allowed us to take advantage of the statistics of two-level systems and simplified further formalism [44,45]. The system of N strongly interacting ions was replaced by a system of N' = NZ weakly interacting bonds, where broken bonds were viewed as antibonds with energy equal and opposite to those of bonds, where Z is the coordination number (Z = 4 for SiO2). Note that for amorphous materials in which cations are not bonded via bridging atoms, such as bridging oxygens in silica or chlorine in ZnCl2, e.g., for amorphous Fe, the number N' = NZ/2. The distance between the cation and the oxygen anion in the SiO4 coordination tetrahedron of amorphous SiO2 is d = 1.62 Å; thus, the volume V0 = π d 3 /6 could be assigned to each bond. The configuron volume V is typically larger than V0; thus, the change in molar volume of material Vm(T) at finite temperatures is typically positive and is proportional to the number of configurons formed: where Vmv(T) is the molar volume of the material when considering only vibrations [46,47] and Vc(T) is the molar volume of configurons, which is given by: where, Nc(Т) is the number of configurons that are formed at temperature T, and the added volume due to the configuron formation is: Figure 3. Configuron formation in amorphous nickel due to a change in the local atomic connectivity by losing one nearest neighbor in the first coordination shell. Reprinted with permission from [25]. Copyright 2020 American Chemical Society.

Glass Transition in Silica
To consider the amorphous SiO 2 structure, we used the CBL where instead of the vitrification of molten material, we started with its solid state, i.e., the glass, and considered its melting, tracing the changes in the spatial distribution of configurons due to a gradual increase in temperature from zero to temperatures above T g . The bond system forms a topologically disordered lattice that is congruent to the disordered lattice of the amorphous SiO 2 . Assuming only two allowed states for bonds, namely, binding (ground) and broken (excited), allowed us to take advantage of the statistics of two-level systems and simplified further formalism [44,45]. The system of N strongly interacting ions was replaced by a system of N' = NZ weakly interacting bonds, where broken bonds were viewed as antibonds with energy equal and opposite to those of bonds, where Z is the coordination number (Z = 4 for SiO 2 ). Note that for amorphous materials in which cations are not bonded via bridging atoms, such as bridging oxygens in silica or chlorine in ZnCl 2 , e.g., for amorphous Fe, the number N' = NZ/2. The distance between the cation and the oxygen anion in the SiO 4 coordination tetrahedron of amorphous SiO 2 is d = 1.62 Å; thus, the volume V 0 = π d 3 /6 could be assigned to each bond. The configuron volume V is typically larger than V 0 ; thus, the change in molar volume of material V m (T) at finite temperatures is typically positive and is proportional to the number of configurons formed: where V mv (T) is the molar volume of the material when considering only vibrations [46,47] and V c (T) is the molar volume of configurons, which is given by: where, N c (T) is the number of configurons that are formed at temperature T, and the added volume due to the configuron formation is: Typically, the added volume is positive ∆V 0 > 0; however, this is not an a priori condition and the materials that have the volumes of configurons V = π d c 3 /6 smaller than the initial bond volumes, i.e., ∆V 0 < 0, will exhibit a negative temperature expansion coefficient when the negative added volume is considerable.
The fraction of configurons that are formed or the degree of breakage of the material network as a function of temperature is determined according to statistical thermodynamics using the expression [29,44,45]: where G c = H c − TS c is the free Gibbs energy of the configurons, H c is the enthalpy of formation, S c is the entropy of formation of the configurons, and R is the universal gas constant. The number of configurons is thus: Configurons that are generated, annihilated, and moving at relatively higher concentrations associate with each other, which is a generic feature of any collection of independent particles that move, reproduce, and die. Indeed, such collections may undergo wild fluctuations at the local and global scales, inducing characteristic patchiness in the spatial distribution of the individuals observed in the spread of epidemics [48][49][50], the growth of bacteria on Petri dishes [51,52], the dynamics of ecological communities [53], the mutation propagation of genes [54], susceptible-infected-removed spreading processes [55], and in the distribution of neutrons in nuclear reactors, which was recently shown to be an effect of clustering [56]. As long as there is a tiny number of small clusters made up of configurons, they can be neglected; however, as the temperature increases, there are more configurons and configuron clusters, which grow in size, and when the threshold value determined by percolation theory [29,57] is reached, a macroscopic percolation cluster of configurons is formed in the system, which penetrates the entire volume of the material. The macroscopic cluster of configurons changes the behavior of the bond system since a path of facilitated motion of atoms appears. Based on Frenkel's [58] and Trachenko et al.'s [59][60][61][62] works, we conclude that when additional availability for atomic motions is ensured, the material shifts from solid-like to gas-like type behavior. In line with Benigni's statement on liquid-like (the B phase in the two-state model) [63], the state of atoms, which are included in the percolation cluster made up of configurons, is assimilated into a gas-like type with consequent contributions to the heat capacity of material and its mechanical properties. Thus, the glass-liquid transition temperature can be found from the condition of reaching the percolation threshold f c [29,44,45]: For the glass transition temperature the equation, this gives: The threshold f c in the idealized case can be taken to be equal to the critical packing density for the three-dimensional space given by the Scher and Zallen invariant [57,64]: The thermodynamic parameters of configurons in amorphous silica are known: H c = 237 kJ/mol and S c = 17.54 R [45], which gives the glass transition temperature of amorphous silica: T g = 1475 K. This temperature practically coincides with the experimentally known temperature T g,exp = 1480 K [1], especially when accounting for the experimental errors [32,65]. The universal critical threshold f c = ϑ c works well for silica and germania [44,45], as well as for many metallic systems [66]; however, for complex oxide systems, f c << 1, which indicates that the percolation threshold is reached at much lower concentrations of broken bonds and the effective configuron radii in these systems are larger [45,67] or there is patchiness of configuron clustering following the generic pattern of dynamic systems [48][49][50][51][52][53][54][55][56]. The formation of stable and ultrastable glasses leads to higher glass transition temperatures. Their liquid states are changed after heating them above T g . Consequently, the percolation threshold is increased. Glacial phases have still higher transition temperatures T g with higher percolation thresholds. These events occur for singular values of enthalpy [23,66], introducing more percolation thresholds. These threshold values must be singular to represent various organizations of elementary bricks. The percolation cluster does not exist below T g ; it appears in the system when the temperature reaches T g , growing further with an increase in temperature. The density of the percolation cluster of configurons is also a universal property that is determined by the percolation theory using: where p is a numerical coefficient close to one, Θ(x) is the Heaviside step function, and β = 0.41 is the critical index in the three-dimensional space [57]. Note that close to T g , we have approximately (f(T) − f c ) = (T − T g ) f c (1 − f c )) H c /RT g 2 [44]. The percolation cluster conventionally penetrates the whole volume of material, dividing it into cells with a size that is the so-called correlation length: where ν = 0.88 is the critical index in the three-dimensional space [57]. The ξ 0 is an elementary length that is approximately equal to the bond length d. As the percolation threshold is approached, we have f(T) → f c , and the correlation length diverges. This reflects the fact that the percolation cluster of broken bonds penetrates the entire volume of an amorphous material. At sizes that are significantly larger than ξ(T), the material can be considered as homogeneous, while at sizes that are similar to or smaller than ξ(T), the material has an inhomogeneous structure, which is described using fractal geometry. As the correlation length equalizes, the sample size L at the glass transition temperature is lower compared with T g of large size samples and is given by [44]: This result is in line with the known size dependence of the glass transition temperature (T g − T g (L)) ∝ 1/L [68] and enables an explanation of the size effects of glass transition in thin films [69]. It can also shed light on the effects that are observed in phase-separated glasses, where the changes in T g may be partly due to the finite size effects of the glass transition [70][71][72][73][74].
A percolation cluster is a fractal object that has the Hausdorff dimension [57] D = 2.55 ± 0.05. The density of a percolation cluster made of configurons (Equation (9)) changes from 0 to 1, with it being equal to 0 up to T g and then growing until it reaches 1 (Figure 4), and was proposed to be considered as an order parameter in the sense of the general Landau theory of phase transformations [44,75].  The number of configurons belonging to the percolation cluster formed above Tg is calculated according to percolation theory: The number of configurons belonging to the percolation cluster formed above T g is calculated according to percolation theory: The volume of a configuron in the percolation cluster V' is not the same as the volume of a separated configuron V; therefore, (V' − V) = 0 and there is an additional change of volume ∆V due to configurons that belong to the percolation cluster: This change is not necessarily positive for any material, although it is typically large ∆V > 0, e.g., the increase in volume of ZnCl 2 is ∆V = 6.9 mL/mol [40].
Thus, above T g , the equation for the temperature dependence of the volume change of materials (Equation (1)) needs to be modified to account for the percolation cluster formed in the liquid phase: The volume added by the percolation cluster is: which accounts for the volume of a configuron in the percolation cluster not being the same as the volume of a separate configuron. Hence, the dependence of the molar volume on temperature, including the configurational contributions (Equation (12)) does not experience an abrupt change and remains a continuous function of temperature, although it exhibits a kink at T g due to contributions of the term in Equation (13) that results from the formation of the percolation cluster made of configurons. Figure 5 depicts a schematic of the volume change of amorphous and crystalline materials due to an increase in temperature. Note that the antibond formation in the excited state does not modify this description of the volume change above T g .

Discussion
The configurational entropy of a material is determined by the number of configurons Nc that are accessible to it at the given temperature T, as follows [76]: (16) Figure 5. Schematic of specific volumes of stable glass (T < T g ), supercooled melt (T g < T < T m ), and stable melt (T > T m ). For silica, T g = 1980 K and T m = 1986 K. The molar volume of silica melt is 27.4 cm 3 /mol. Note that the slope of lines is significantly exaggerated for the amorphous silica to account for the very low expansion coefficient, which is not higher than 4.9 × 10 −6 K −1 for the glassy phase and not higher than 7 × 10 −6 K −1 for the melt [1].

Discussion
The configurational entropy of a material is determined by the number of configurons Nc that are accessible to it at the given temperature T, as follows [76]: where k B is the Boltzmann constant, W is the total number of distinct packing states that are available to a system, and W = N c (T) while T < T g . This term exists at any finite temperature T > 0 and cannot be arbitrarily set to zero as is most often done in the two-level models, such as that described in [63]. However, above T g , a new term needs to be added to the configurational entropy of liquid (Equation (16)), which is due to the formation of the percolation cluster that is made up of configurons, where atoms exploit their new degree of freedoms as more space becomes available for the cluster's motion: Figure 6 demonstrates the natural appearance above T g of the contribution of configurons that are associated with the percolation cluster for amorphous diopside, where the numerical data is taken from [32] (p 257). The existence of two variations of the entropy above T m is shown because two values of enthalpy exist up to the temperature T n+ : (i) The effective entropy S is decreased at T n+ (θ n+ = (T n+ − T m )/T m ): (ii) The effective entropy S is increased at T n+ : Consequently, the entropy variation at T n+ has two values ∆S (θ n+ ) = ±S m θ n+ /(1 + θ n+ ) V pc (T) can be positive or negative. It is negative for systems that are annealed between T m and T n+ or slowly cooled, favoring bond formation. It is positive for quenched systems, favoring antibonds when an enthalpy increase corresponds to a volume increase. The sign of this quantity depends on the heating and cooling rates and mainly on the thermal history above T m .
Vpc (T) can be positive or negative. It is negative for systems that are annealed between Tm and Tn+ or slowly cooled, favoring bond formation. It is positive for quenched systems, favoring antibonds when an enthalpy increase corresponds to a volume increase. The sign of this quantity depends on the heating and cooling rates and mainly on the thermal history above Tm. Figure 6. The temperature dependence of the entropy of amorphous and crystalline phases of diopside with the explicit contribution of a configuron percolation cluster above the Tg. ΔSf shows the entropy difference between the liquid and the crystalline phases at the melting temperature.
The account of configuron phase formation above Tg hence shows that the glass transition is similar to the second-order phase transformations within the Ehrenfest classification [77][78][79], which is based on the change of state variables at the transition temperature Ttr. The second-order transitions are continuous for the volume, the enthalpy, and the entropy, which show a kink at the Ttr, while the second derivative of the chemical potential Figure 6. The temperature dependence of the entropy of amorphous and crystalline phases of diopside with the explicit contribution of a configuron percolation cluster above the T g . ∆S f shows the entropy difference between the liquid and the crystalline phases at the melting temperature.
The account of configuron phase formation above T g hence shows that the glass transition is similar to the second-order phase transformations within the Ehrenfest classification [77][78][79], which is based on the change of state variables at the transition temperature T tr . The second-order transitions are continuous for the volume, the enthalpy, and the entropy, which show a kink at the T tr , while the second derivative of the chemical potential and heat capacity change discontinuously. In line with the Landau theory of second-order phase transformations, they occur when a new state of reduced symmetry develops continuously from the more disordered phase, which is typically at the higher temperature. He observed that near a phase transition, an approximate form for the free energy can be constructed with an order parameter that is zero before the transition and nonzero after the phase transition [79]. Examples of second-order transitions include the ferromagneticparamagnetic transition, where the magnetization is the order parameter and the structural phase transition from a cubic phase to a tetragonal phase occurs, where the order parameter is the parameter (c/a − 1), where c is the length of the long side of the tetragonal unit cell and a is the length of the short side of the tetragonal unit cell. For second-order phase transitions, the order parameter increases continuously from zero, starting at the T tr . As we see from Equation (9) and Figure 2, the transition of a glass to a liquid at T g exactly follows the universal case of second-order phase transformation, although it occurs in the metastable phases as both a glass and a supercooled liquid are not thermodynamically equilibrium states of matter. Due to the formation of the configuron phase, the specific heat capacity of amorphous materials at the transition follows the universal law: where C 0 is a coefficient that accounts for the contribution of the configuron phase. As (1 − β) = 0.59, the heat capacity contribution from the formation of configuron phase diverges at T g , which, in practice, is used to determine T g of glass-forming liquids. The same temperature dependence is true for the thermal expansion coefficient of amorphous materials. This type of behavior of heat capacity and thermal expansion allows for an exact identification of T g when increasing the temperature through the transition of glass into a liquid.
The configurons are quasiparticles that are specific to initially solid materials (glasses) and can only be properly identified when the system is not too far from the percolation threshold and therefore not significantly above the T g ; the system loses any individuality at T m and above it. The vitrification of liquids occurs when cooling down, which induces the formation of new bonds rather than configurons until the system finally achieves the rigidity threshold, which is the same as the percolation threshold following Kantor-Webman theorem [43]. The vitrification is not associated with discontinuities in the temperature dependence of heat capacity, i.e., when decreasing the temperature, because melt vitrification is not accompanied by configuron phase formation. The contribution from the annihilation of configurons, which is the same as bond formation when decreasing the melt temperature, has an opposite sign (exo-instead of endo-) to that of increasing the temperature. Hence, the CPT naturally explains the experimentally known hysteresis in the behavior of heat capacity near the glass transition when increasing or decreasing temperatures near T g . Notably, an account of the clusters formed in the liquids enables an explanation of liquid-liquid phase transformations and the Mpemba effect and its inverse, which exists in all melts. The complete annihilation of configurons occurs by heating above the temperature T n+ > T m toward the homogeneous liquid state. Each configuron has the energy of an antibond, which is equal in size but with an opposite sign to that of a bond that would appear by annealing the homogeneous liquid between T m and T n+ and would disappear by heating above T n+ [80]. Figure 4 reveals an atypical feature of the glass transition within CPT, namely, the location of the more ordered configuron phase above the T tr . The unusual location in the higher temperature range is due to the increase in configuron concentration with temperature. This makes the glass-liquid transition rather similar to the formation of condensed Rydberg matter at high levels of excitation of atoms [81,82], although the condensation of excited atoms into condensed Rydberg matter occurs via a first-order phase transition. Notably, the glass transition can be interpreted in terms of first-order transformations [12,17,20,23].
The configuron phase that is formed exactly at T g grows with an increase in temperature above T g . Set theory, which is a branch of mathematical logic that studies abstract sets as collections of objects [83], can be used to characterize the configuron phase, where the Hausdorff-Besicovitch dimension is defined as the limit: where N(ε) is the number of boxes of side length ε that are required to cover the set of configurons. It is concluded that the set has the Hausdorff-Besicovitch dimension D when N(ε) grows proportionally to 1/ε D as ε tends to zero. As shown above in Figure 5, the dimensionality of the set of configurons (configuron phase) changes at T g from 0 to D = 2.55 ± 0.05 above it [67]. This stepwise change of dimensionality of the set of configurons is due to the formation of a configuron phase above T g and has a kink in the first sharp diffraction minimum of scattered X-ray or neutrons as a consequence [15,25]. The presence of this kink gives evidence of structural differences between glasses and melts and is used in practice to identify T g of metallic systems [25]. Changes during the transition of glass into a liquid form the basis of the determination of T g by various techniques, such as differential scanning calorimetry (DSC), thermomechanical analysis (TMA), dynamic mechanical analysis (DMA), rheological methods, dielectric analysis, electron spin resonance (ESR), and nuclear magnetic resonance (NMR), where a combination of various techniques adds value to understanding the nature of glasses and the glass transition [84]. Notably, the structural changes of amorphous materials during the glass transition, in addition to X-ray scattering data, are also confirmed by atomic force microscopy [85,86] and the analysis of fifth-order susceptibility [87]. As Table 1 reveals, the main result of configuron formation in an amorphous material is the shift of atoms from the first coordination shell, regardless of the type of bonding. This means that the analysis of the maxima (PDF max ) and minima (PDF min ) of X-ray pair distribution functions is one of the most informative tools that allow for revealing the glass transition and structural changes that accompany this transformation [19,25,[88][89][90][91][92][93][94]. Figure 7 demonstrates the stepwise changes in the PDF min based on the data of amorphous Ti 2 Ni taken from [25], along with the changes in the Hausdorff-Besicovitch dimensionality of the set of configurons D, which are of future research interest for amorphous silica.
A slower cooling of a melt results in a lower T g (this is not shown in Figure 5). Below the melting point T m , any liquid tends to crystallize and, with exceptions, such as that of atactic polymers, crystallizes upon sufficiently slow cooling. Crystallization can be avoided kinetically by increasing the cooling rate, although, at any final cooling rate q, the supercooled liquid will contain a certain fraction, denoted here by x, of the crystallized material. The critical cooling rate for glass production is defined as the minimum cooling rate at which the degree of crystallinity of the frozen liquid does not exceed a certain critical value x c , where, for good glasses, x c~1 0 −6 − 10 −2 [95]. In an ideal amorphous material, for which x = 0, the percolation threshold is determined using the Scher and Zallen invariant f c = ϑ c ; however, in a real material, a smaller part of the volume (1 − x) is available for the development of clusters that are made up out of configurons. Similarly, the percolation threshold is similar to the system polydisperse particles [96] and is modified to be f c = ϑ c (1 − x). The fraction of the volume of the crystallized material can be found according to the Kolmogorov-Avrami theory [97,98], where at a constant nucleation rate I v , the growth rate of the crystalline phase u and temperature T depends on time t according to x ∼ = π I v u 3 t 4 /3. Thus, the glass transition temperature depends on the logarithm of the cooling rate T g = H c /{S c + R ln[[1 − f c + π I v (T m − T g ) 4 u 3 /3q 4 ]/f c ]} [44,46], which is known experimentally [7]. A proper account of relaxation effects in amorphous materials shows that T g cannot be decreased indefinitely by decreasing the cooling rate and the glass transition temperature interval is quite limited [24]. Kinetically controlled shifts of T g within a relatively narrow temperature interval by the decrease of cooling rate q are thus explained, revealing the logarithmic dependence of T g on q within CPT. The nucleation of bonds also occurs between T m and T n+ and not only below T m [36,80]. Consequently, clusters of bound colloids are built below T n+ , and when the number of new bonds attains the percolation threshold, the crystallization occurs at T m if the liquid is slowly cooled instead of being quenched.
the glass transition, in addition to X-ray scattering data, are also confirmed by atomic force microscopy [85,86] and the analysis of fifth-order susceptibility [87]. As Table 1 reveals, the main result of configuron formation in an amorphous material is the shift of atoms from the first coordination shell, regardless of the type of bonding. This means that the analysis of the maxima (PDFmax) and minima (PDFmin) of X-ray pair distribution functions is one of the most informative tools that allow for revealing the glass transition and structural changes that accompany this transformation [19,25,[88][89][90][91][92][93][94]. Figure 7 demonstrates the stepwise changes in the PDFmin based on the data of amorphous Ti2Ni taken from [25], along with the changes in the Hausdorff-Besicovitch dimensionality of the set of configurons D, which are of future research interest for amorphous silica. A slower cooling of a melt results in a lower Tg (this is not shown in Figure 5). Below the melting point Tm, any liquid tends to crystallize and, with exceptions, such as that of Figure 7. Schematic of the pair-distribution function PDF(R) = g(R) of amorphous Ti 2 Ni (a). Reprinted with permission from [25]. Copyright 2020 American Chemical Society. The prominent changes during the glass transition of PDF min (b,c) and the Hausdorff-Besicovitch dimension of configurons set (d).
The above analysis of structural changes of an amorphous material during the glass transition reveals that the transition occurs continuously and that it is similar to secondorder phase transformations. The fractal dimensional configuron phase formed above T g is responsible for the liquid-like behavior of the material. Although the transition is continuous and reversible, vitrification occurs differently when cooling the liquid down as no configuron clusters are available to form or disintegrate in the melt. At the same time, the CPT shows a continuous two-exponential function for the viscosity, which provides a rather exact description of viscosity-temperature relationships for both amorphous silica and other materials [24,44,75,99]. However, it should be emphasized that CPT with the equations used above operate with the volume of bonds (V 0 ) and configurons (V) and not with the so-called excess or free volume of the material, which is the specific volume of the material (the molecular weight divided by the density) minus the volume of molecules, which is typically assumed to be the volume of hard spheres representing atoms. In line with that, we recall that Doremus (characterizing the predictive ability of viscosity models) accentuated that the bonding between molecules, defects, and the structure are much more important than the free volume [100].

Conclusions
The role of structural changes during the glass transition is of primary importance to understand the reasons behind the drastic changes of amorphous materials' behavior due to an increase in temperature when crossing T g . The higher the temperature, the more chemical bonds are broken and, thus, the material's building blocks (either atoms or molecules) become increasingly free to move, gradually changing the solid-like behavior to a more gas-like type behavior. Based on the Kantor-Webman theorem, we conclude that when the percolation threshold for unbound blocks is achieved, then the rigidity threshold of the material is also achieved and, thus, the properties of materials qualitatively change from the solid-like to liquid-like. The crossing of the threshold is continuous; therefore, it is almost an undetectable structural change of the atomic system, i.e., the spatial distribution of building blocks of the material. However, the structural change can be readily revealed for the bonding system instead of the atomic one. The system of configurons, which are broken bonds with an account of relaxation effects near the bond broken, is a point-like system that resembles a gaseous system up to a temperature close to T g . The set theory characterizes this system as a set with a Hausdorff dimensionality of 0. Exactly at T g , the system of configurons forms a macroscopic percolation cluster for the first time, which is known to be a fractal that is characterized by the Hausdorff dimensionality 2.55 ± 0.05. Thus, the set of configurons changes its dimensionality in a stepwise manner from 0 to 2.55 ± 0.05 at T g , which can be interpreted as a kind of symmetry change characterizing the structural difference between a glass and a melt. The qualitative difference between glasses and melts, which could not be revealed for the atomic system, becomes obvious for the system of chemical bonds, which resembles viewing the real image as a positive photograph instead of a negative image on the photographic film.