Evidence for a Correlation of Melt Fragility Index With Topological Phases of Multicomponent Glasses

Scores of glass compositions in the equimolar Ge x As x Se 100 − 2x ternary system are synthesized across the 0 < x < 26% range, and their homogeneity tracked by ex-situ Raman proﬁling. In synthesis, we alloyed the starting materials to homogenize until the variance, <1 × > Ge , in Ge content “x,” across a 1.5-g sized batch composition of < 0.01% was realized. We undertook Modulated-DSC, Raman scattering

<0.01% was realized.We undertook Modulated-DSC, Raman scattering, and molar volume experiments as a function of composition.Trends in T g (x) increase monotonically with x over the examined range, but the non-reversing enthalpy of relaxation at T g , H nr (x), displays a sharp square-well like reversibility window over the 9.5(2)% < x < 17.0(2)% range even in rejuvenated glasses.Trends in melt fragility index, m(x), established by measuring the T-dependence of the enthalpy relaxation time across T g , show a Gaussian-like variation with m(x) < 20 in the 9.5(2)% < x < 17.0(2)% range, and m > 20 outside that range, thus establishing a fragility window.The close correlation between the variations of m(x) and H nr (x) underscores that super-strong melts formed in the fragility window form Intermediate Phase (IP) glasses in the reversibility window, while fragile melts formed at non-IP compositions give rise to either flexible or stressed-rigid glasses.Molar volumes of glasses reveal a global reduction for IP glass compositions, thereby underscoring the compacted nature of the isostatically rigid networks formed in that phase.Evidence for specific dynamic features of the different phases is provided by molecular dynamics simulations, which indicate that the diffusivity is minuscule for x > 17% and increases substantially in the flexible phase.These results show that liquid dynamics and fragility encode glass topological phases in this ternary chalcogenide.The optimization of the glass-forming tendency for IP compositions, the stress-free nature of these networks, and the qualitative suppression of aging, are features that all point to their ideal nature.These ideas of IP glasses as ideal glasses are at variance with the notion that ideal glasses possess a low configurational entropy and form at T g values that approach the Kauzmann temperature.IP networks display, adaptability, high glass-forming tendency, stress-free nature, and high configurational entropy of networks.

INTRODUCTION
Chalcogenide glasses have served as useful test laboratory to examine the predictions of Topological Constraint Theory (TCT) since the theory was first introduced (Phillips, 1979) in 1979 and its vibrational consequences elucidated by M.F.Thorpe (Thorpe, 1983) in 1983.Chalcogenides have also emerged as the materials of choice for 3D-X-point memory devices(Intel 1 ; Malventano 2 ; Tang et al., 2009), phase-change memory (Raoux and Wuttig, 2009), Continuum emission in the IR (Goncalves et al., 2018;Tremblay et al., 2018), optical fibers, and waveguides (Tang et al., 2015;Chen et al., 2019).The interplay between these applications and underlying basic science including topological phases (TP) now offers new prospects to tune material functionality with applications.Amongst the chalcogenides, and particularly within the context of TCT, binary Ge x Se 100−x glasses have acquired much notoriety starting from the mid-1980s.Inelastic neutron scattering experiments revealed a low frequency mode near 5 meV (Boolchand et al., 1990;Kamitakahara et al., 1991) and one found that the fractional scattering strength of the mode normalized to the total vibrational density of states acquires a value of nearly 1/3 in pure Se as x→ 0. The result gives confirmation of the floppy mode concept that Thorpe introduced (Thorpe, 1983) to establish the percolation of rigidity near <r> = 2.40.This indicates that the elastic phase transition is driven by a control parameter <r>, whereas the density of low frequency modes represents the order parameter of the transition.The presence of these low frequency excitations led the 119 Sn Mössbauer Debye-Waller factors in Ge x Se 100−x glasses to display a giant softening as x→ 0. This follows because the meansquare displacement (<r 2 >) relates to the first inverse moment of the vibrational density of states.Measurements of 119 Sn meansquare displacement, <r 2 > as a function of Ge content x, show a vibrational threshold near x = 20.0%, or <r> of 2.40, in harmony with the notion of percolation of rigidity (Boolchand et al., 1990).And the connection of such floppy modes has been linked to the energy landscape in glasses (Naumis, 2005).
A significant step forward in the field was the observation of a reversibility window in binary Si x Se 100−x glasses (Selvanathan et al., 1999).In MDSC one observes the enthalpy of relaxation at T g to show a square-well like minimum (reversibility window), with the walls of the window.Such walls have been identified in a variety of theoretical models as a rigidity and a stress transition (Thorpe et al., 2000;Micoulaut and Phillips, 2003;Barré et al., 2005;Chubynsky et al., 2006).The first transition near x = 20% (<r> = 2.40), the rigidity transition, is followed by the second transition near x = 27% (<r> = 2.54), the stress transition (Thorpe et al., 2000).And it was not until 2010 (Micoulaut, 2010) that one recognized that the near vanishing of the enthalpy of relaxation in the reversibility window constituted a signature of the isostatically rigid nature of glassy networks (n c = 3) formed between the two elastic phase transitions.Such reversibility windows have now been documented in several scores of glasses (Yildirim et al., 2016b), unequivocally demonstrating the generality of the three topological phases, flexible-, intermediate, and stressed-rigid.
As larger batch sizes (2-20 g) of the chalcogenides were synthesized, it emerged in 2011 that such melts undergo delayed homogenization (Bhosle et al., 2012b;Chakraborty and Boolchand, 2014;Ravindren et al., 2014).The behavior can now be traced to the super-strong behavior of Intermediate Phase (IP) melts formed in these binary glass systems.The present equimolar Ge-As-Se ternary drew our attention for two reasons.Although MDSC and Raman scattering results did reveal (Qu et al., 2002) evidence of a reversibility window, the glass samples synthesized in 2002, were not homogeneous, a recognition that emerged only later in 2011.The second reason came from the recognition that the glass transition temperature T g (x) displayed an almost linear variation with x, even when x exceeded the chemical threshold near x c = 18.18%.Such is not the case in several binary Selenide and Sulfide glasses such as Ge-Se (or S), and As-Se (or S) where a T g maximum occurs near the chemical thresholds.The observation suggested (Boolchand et al., 2002) that fully polymerized glassy networks must apparently exist in the Ge x As x Se 100−2x ternary over a rather wide range of stoichiometry, 0 < x < 26%.That feature of the ternary was appealing (Qu et al., 2002;Qu, 2004), as it suggested the absence of any NSPS at compositions exceeding the chemical threshold.Thus, probing the nature of the topological phases, and especially the underlying rigidity-and stress-elastic phase transitions that separate these phases would be immune from effects due to NSPS.
The finding of delayed homogenization of binary Ge-Se melts/glasses was a strong impetus to re-examine the TPs in the present Ge x As x Se 100−2x ternary, and to clarify the role of glass synthesis in relation to the sharpness of the rigidity-and stresselastic phase transitions that border the Intermediate Phase (IP).The second goal of the present work was to also measure the variation of fragility index in such homogenized melts and to explore if trends in melt dynamics have a bearing on the nature of the glass TPs.
In the present work we report observing a rather abrupt thermally reversing window that is square-well-like in nature even in freshly quenched glasses after rejuvenation (T g cycling) of the especially homogenized Ge x As x Se 100−2x ternary glasses.Furthermore, we also examined the variation of fragility index (m(x)), and find evidence of a compositional range across which m < 20, which defines the "fragility window."Remarkably, the compositional width of the "fragility window" observed not only coincides with that of the "reversibility window," but both are somewhat wider in range, 9.5% < x < 17.0%, than the previous report (Qu et al., 2002;Qu, 2004).These experimentally measured anomalies are correlated with dynamic anomalies calculated from molecular dynamics simulations and a huge increase in diffusivity is detected once the glass-forming liquids become flexible.
Why is the homogeneity of a melt/glass so critical to examining topological phases?Previous molecular dynamics simulations have shown that the reversibility window arises from the fact that the flexible-to-rigid and unstressed-to-stressed transitions occur at different compositional thresholds, thereby giving birth to an intermediate phase that is rigid but free of stress (Bauchy, 2019).These two elastic transitions have been shown to be percolative in character (Jacobs and Thorpe, 1995;Thorpe et al., 2000) and, hence, one expects these transitions, in principle, to be intrinsically sharp in the composition space.In binary Ge x Se 100−x glasses, for example, the network connectivity is determined locally by the Ge fraction x.If the variance in Ge stoichiometry "< ×> Ge " across a glass specimen were to exceed 1%, one can expect the two elastic phase transitions to be smeared.That is exactly what we observed in our experiments on the Ge x As x Se 100−2x ternary in 2003.If the variance in Ge stoichiometry exceeded 3%, the reversibility window could be triangular (Bhosle et al., 2012b) or, worse yet, the window may not be observed.On the other hand, if the variance in Ge stoichiometry is <0.01%, as noted in the present work, one expects and, indeed, observes rather sharp rigidity-and stress-transitions.And, ideally, one would like to establish the Ge variance on an atomic scale.As a practical matter it should be at least on a scale of 50 microns, if not lower.The observation of abrupt elastic phase transitions in a topologically disordered network is not only novel and exciting but also constitutes evidence of the percolative nature of these elastic phase transitions.

Ft-Raman Profiling
The starting materials for the synthesis included small pieces (∼ 3 mm in diameter) of 99.999% pure Se nodules from Alfa Aesar, 99.999% pure Ge lumps from Alfa Aesar and stoichiometric As 2 Se 3 from Noah Technologies.We handled these materials in a model HE-493/MO-5 Glove box from Vacuum Atmospheres, which is designed to be flushed with dry N 2 gas continuously and with a relative humidity of <0.01%.We sealed small pieces of the starting materials of about 3 mm in size or less, in the desired weight ratio, in evacuated quartz tubes of 5 mm ID and 1 mm wall thickness.Prior to their use, we dried quartz tubes in a vacuum oven at 80 • C for at least 24 h.Batch size were kept at 1.5 g, and the starting materials encapsulated in evacuated quartz tubes at a pressure of 2 × 10 −7 Torr using a liquid Nitrogen trapped diffusion pump system.We heated quartz tubes to 950 • C for typically 10 days, keeping them vertical in a box furnace.Water impurity doping of chalcogenide glasses, unwittingly assists homogenization of melts, increases the enthalpy of relaxation of glasses, decreases T g , and increases the fragility index of melts (Bhosle et al., 2012a,b).These extrinsic features are not representative of the dry melts/glasses.
A Thermo-Nicolet Nexus 870 based FT-Raman system, using 1,064 nm radiation from a Nd-YAG laser, focused to a 50 µm spot size (micro-setting), was used to excite the scattering.The collection optics also has provision for a macro-setting with a laser spot size of 1 mm.Glass specimen wetting 5 mm ID evacuated quartz tubes were mounted on a digitally controlled x-y-z stage to focus the laser beam entering along the quartz tube diameter (x-axis adjustment) and onto the quartz tube glass specimen interface (z-axis adjustment) to optimize the signal.The y-axis adjustment of the sample stage permitted to access specific locations along the glass column for recording Raman scattering.Spectra were acquired at typically six locations about 2 mm apart along the glass specimen column of typically 20 mm in length (y-axis adjustment).The software permits comparing observed lineshapes normalized to the highest scattering strength mode, and that facilitated directly variations in glass specimen stoichiometry by comparing the scattering strength of the corner sharing (CS) of GeSe 4 tetrahedra (near 200 cm −1 ) with the broad band of a triad of modes (Dash et al., 2017) associated with Se n chains-uncorrelated, Se n chain-correlated, and Se 8 rings (near 250 cm −1 ).Periodically, quartz tubes were water quenched after lowering their T to 650 • C and Raman spectra taken at 6 locations along the 20 mm long melt column.By using a laser beam spot size of 50 µm or less, one probes glass homogeneity to a finer scale, which is essential to get a true representation of glass/melt homogeneity.It is important to emphasize that in these homogeneity studies the full batch composition was profiled.This of course is necessary since the stoichiometry of a small part of the batch is tied to the stoichiometry of the rest of the batch (Gunasekera et al., 2013).Figure 1 shows a partial summary of the observed spectra; and one finds how the starting melt, which is rich in Ge(Se) atoms at the bottom(top) of the vertical melt column, steadily homogenizes as Ge atoms diffuse up while Se atoms diffuse down the melt column.After 2 days of alloying the variance in Ge stoichiometry, < ×> Ge is found to be about 1.04%, but it systematically decreased by almost three orders of magnitude to nearly 0.01 % or lower after 10.5 days of alloying (Figure 2).In these profiling experiments, one could successfully decode the heterogeneity of the batch to a level of 0.01% range in the variance of Ge by merely sampling only a tiny volume fraction of 2% of the batch column.The existence of a gravitationally induced segregation that sets up in the melt column at the beginning of the alloying process as quartz tubes are held vertical, aids in tracking how Ge-and Se-atoms, respectively, diffuse the melt column (Gunasekera et al., 2013) to homogenize with no rocking.Had we rocked the melt column it < ×> Ge in Ge concentration "x" across the batch decreased from 1.04(5) to 0.28(5) to < 0.05% upon progressive alloying as revealed by the spread in the normalized scattering strength of the Corner-Sharing mode near 200 cm −1 in insets.The Omnic software normalizes each spectrum to the most intense peak near 250 cm −1 coming from polymeric Se n chains.The < ×> Ge variance is being probed on a 50 micron scale set by the laser-spot size.
FIGURE 2 | Observed reduction in the Ge variance < ×> Ge of the Ge x As x Se 100−2x glasses as a function of alloying time t r in days at indicated compositions x when probed at a 50 µm resolution determined by the laser spot size.Glass homogeneity steadily increases as < ×> Ge decreases by almost two orders of magnitude upon alloying the starting materials for typically 10-13 days.We show that delayed homogenization is tied to the Super-strong character of melts in the 9.5% < x < 17% range, the Intermediate Phase.
would have been difficult to detect the isolated super-strong melt inclusions spread across the melt column unless one probed more like 10-30% or more of the batch volume, which surely would have been an enormous task in itself.
There are other aspects of the homogenization studies that we mention in passing.The scale on which a batch is examined in the Raman profiling experiments matters because TPs are determined by atomic scale network connectivity or mean coordination number <r>, fixed by the local Ge, As, and Se stoichiometry.Furthermore, the scale on which a batch is examined in the Raman profiling experiments matters.A sample may be homogeneous after 2 days of alloying on a spatial scale of 1 mm (laser spot size) but not on a finer scale of 50 µm (laser spot size) as we have noted (Bhosle et al., 2012a).Finally, melt homogenization kinetics are broadly determined by the slowest diffusing species or the species with the highest viscosity, which occurs in the IP.Most compositions in the Ge x As x Se 100−2x ternary, at x > 5%, take typically 7-11 days to homogenize in our chosen geometry.The IP compositions such as x = 10 and 13.3% take a few days longer to homogenize than non-IP ones.The fragility index of IP melts apparently controls the kinetics of melt homogenization, as we will comment later in discussion of results.In very select cases where a stoichiometric glass is being synthesized at a composition where there exists a distinct well-defined underlying simple crystalline phase, delayed homogenization could be circumvented.In the present ternary, we are unaware of any crystalline phase or phases that occur and that may expedite the homogenization process.
Finally, in our present efforts we aimed to synthesize the most homogeneous melts to explore the physics of TPs in network glasses.Once these fundamental issues are understood, new strategies can be implemented to reduce the alloying time of glass synthesis by at least an order of magnitude to achieve homogeneous melts/glasses.

Modulated-DSC as a Probe of the Enthalpy of Relaxation of Glasses at T g and the Fragility Index of Melts
Modulated Differential Scanning Calorimetry (MDSC) has proven to be a reliable (Thomas, 2005;Bhosle et al., 2012b) and rather direct means to independently establish the glass transition temperature T g , and the modulation frequency corrected enthalpy of relaxation of a glass at T g ( H nr ) with minimal scanning-rate-related kinetic upshifts.The procedure for obtaining the frequency corrected enthalpy of relaxation was developed as early as in 1995 by Len Thomas (Thomas, 2005), and more recently confirmed in numerical simulations (Guo et al., 2012).The measured H nr term is a pivotal calorimetric observable of a glass, that characterizes not merely its non-ergodic state but is crucially linked to the flexibility and rigidity of networks (Micoulaut, 2010) provided these are dry and homogenous.Furthermore, MDSC also permits measuring the fragility index "m, " of a melt.We illustrate the principles of both these methods by showing results of ternary Ge x As x Se 100−2x glasses at select compositions.We used a Q2000 TA Instruments MDSC with glass samples in hermetically sealed T zero Al pans (Thomas, 2005) to achieve the best baselines in the calorimetric scans.
Enthalpy of Relaxation of a Glass at T g A significant advance in understanding the nature of T g in network glasses has emerged from the use of MDSC over the traditional DSC (Bhosle et al., 2012a,b;Gunasekera et al., 2013).In this approach, one programs a sinusoidal temperature variation on a linear T-ramp, to extract the part of the total heat flow that tracks the sinusoidal oscillations, called the reversing heat flow, and then obtain from them the difference signal between the total heat flow and the reversing heat flow, the non-reversing heat flow.Experiments reveal that the reversing heat flow shows a rounded step (shown in blue), while the nonreversing heat flow, shows a Gaussian-like peak (shown in red) that relates to the overshoot observed in the total heat flow (shown in green) as illustrated in Figure 3 for the case of a ternary glass at x = 11%.The red shaded area in the endotherm (curve 3, Figure 3) gives the enthalpy of relaxation upon heating the glass, and is used to calculate the frequency corrected enthalpy of relaxation (Thomas, 2005) of H nr (x = 11%) = 0.00(3) cal/g.The H nr term is closely tied to the physics of the glassy state.The deconvolution of the total endothermic heat flow into the reversing and non-reversing heat flow is of fundamental interest in glass science.The reversing heat flow permits establishing the T g from the inflection point of the step while the step height C p yields the change in the specific heat between the melt and the glass.The non-reversing heat flow provides a measure of non-ergodicity of the T g transition and it's vanishing for the case of a ternary glass at x = 11%, constitutes the signature of an isostatically-rigid glass, i.e., n c = 3 (Micoulaut, 2010).The reversing heat flow provides a measure of vibrational entropy change between the glass and the melt, while the non-reversing heat flow provides a direct measure of the configurational entropy difference between the glass and the melt.The vanishing of the H nr (x) term suggests that the configurational entropy of the glass is melt-like, in other words the excess entropy of the glass remains high as in the corresponding melt.Excess entropy is defined as the entropy difference between the super cooled liquid at T and the state of the glass at which the entropy can no longer be lost (Bestul and Chang, 1964).

Fragility Index of a Melt
There is broad recognition (Carpentier et al., 2003;Thomas, 2005) that MDSC experiments can be undertaken to analyze the glass transition endotherm either in terms of the reversing and non-reversing heat flow components or in the complex C p formalism.In the latter approach, the in-phase-part of the complex specific heat, C p ip , shows a step-like feature akin to the reversing heat flow, while the out-of-phase part of the complex specific heat, C p op , shows a Gaussian-like peak akin to the nonreversing enthalpy of relaxation.Starting from a T above T g , as one cools a melt across T g , one observes a glass transition exotherm.One analyses the thermal event in terms of C p ip and C p op as a function of modulation frequency and obtains the melt fragility index.In these experiments, one tracks the shift of C p op peak as the modulation frequency is systematically changed.
To illustrate the method, we show below (Figure 4) the complex C p results for a Ge x As x Se 100−2x melt at x = 14%.These results were obtained as a function of modulation time period t m in the range 60 s < t m < 120 s as the melt T is lowered from 220 to 140 • C. One finds that the step in C p ip steadily shifts to higher T as the modulation frequency ω = 2π/t m increases.On the other hand, C p op displays a peak that steadily shifts to higher T as ω increases.At the peak location the glass enthalpy of relaxation tracks the programmed modulated T-profile, i.e., ωτ = 1.By establishing the peak location of the relaxation peak, one fixes the enthalpy relaxation time τ given the modulation frequency (ω) used, and from the abscissa of the C p op , the T at which that enthalpy relaxation time τ occurs.By plotting log of τ as a function of T g /T (Figure 4 right) one obtains the fragility index.
The following relation then permits deducing the activation energy, E a , for enthalpy relaxation, The typical operating parameters are as follows; the cooling rate used was 0.5 • C/min, modulation time varied from 120 to 60 s, modulation amplitude was kept at ±1 • C. We calibrated the Q2000 MDSC calorimeter for temperature by performing an Indium melting scan.We calibrated the head for specific heat by using a sapphire standard.These are slow experiments typically requiring 48 h of scanning to examine one glass composition.
Analogies between dielectric relaxation studied in electrical measurements and enthalpy relaxation examined in a calorimetric measurement (MDSC) in several glass forming systems, both organic as well as inorganic, show a rather appealing feature (Carpentier et al., 2003;Descamps and Decroix, 2014).Dielectric relaxation involves variation of polarization with a time varying electric field, while enthalpy relaxation relates to the variation of heat flow with a time varying temperature.Dielectric relaxation proceeds at a higher frequency (0.05 Hz−150 KHz) than enthalpy relaxation (0.01-0.1 Hz), and thus occurs at a lower T regime.And if one plots how these relaxation times vary with T, one observes a smooth extrapolation of the dielectric relaxation times to the longer enthalpy relaxation times as T approaches T g (Carpentier et al., 2003).These results reassure us that melt dynamics probed by two different probes in different relaxation time regimes, show universality of behavior.These observations lend credence to the melt fragility index results on glasses deduced from MDSC.
The increase of viscosity (η) of a glass forming melt as its T is lowered to T g is often described (Vogel, 1921;Fulcher, 1925;Tammann and Hesse, 1926;Williams et al., 1955;Angell, 1995;Hodge, 1996) by plotting log η against T g /T.Only select melts display an Arrhenius variation of η on such a plot, such as fused SiO 2 , and these are termed strong, but a large number of melts display a super-Arrhenius behavior in which the activation energy E a of viscosity (η = η o exp E a /kT) steadily increases as T approaches T g .One describes (Angell, 1995;Hodge, 1996) such melts as fragile.And, formally, one defines the fragility index of a melt, m, as the dimensionless slope of the log η vs. T g /T, as T approaches T g (Equation 3).Within the fragile-strong classification, m varies from about a minimum value near 17 for super-strong melts to a maximum value of about 145 for the fragile ones (Mauro et al., 2009;Gunasekera et al., 2013;Chakraborty and Boolchand, 2014).
The Maxwell relation (Equation 4) connects viscosity to the shear stress relaxation time τ with G ∞ representing the infinite frequency shear modulus.And, since G ∞ is T-independent, one can write the fragility index as shown in Equation 1.
In the calorimetric measurements one measures the enthalpy relaxation time τ H near T g to deduce the fragility index (Carpentier et al., 2003).
A wide variety of glass forming melts, both organic as well as inorganic, display a fragile behavior (Angell, 1995;Debenedetti and Stillinger, 2001).And one expects that at high T, close to twice T g , the viscosity of a superstrong melt (m = 17) will exceed that of a fragile melt (m = 35) by one order in magnitude.Most remarkably, for the especially homogenized IP melt compositions examined presently and elsewhere (Gunasekera et al., 2013;Chakraborty and Boolchand, 2014;Mohanty, 2018) the superstrong nature of such melts with m < 20 (Figure 5C), leads understandably to a higher viscosity of such melts than those of the flexible phase or stressed-rigid melts.This key finding rests at the base of delayed homogenization of glass forming melts encompassing an IP composition.

Raman Scattering
A Thermo-Nicolet FTIR Nexus 870 research bench with a Raman module was used to record Raman scattering from the bulk glasses encapsulated in evacuated quartz tubes.The bulk glass specimen was positioned on an xyz stage of the instrument and the scattering was excited using the 1,064 nm radiation from a Nd-YAG laser tightly focused to a 50 µm spot size.In a typical run 126 mW of laser power was used, the radiation focused at the quartz/ bulk glass interface and an acquisition of 100 cycles took 7 min of acquisition time.In section Quantitative Mapping of melt homogeneity by FT Raman Profiling below, we provide details of the quantitative mapping of glass/melt homogeneity from these Raman scattering data that profiled the bulk glass.

Quantitative Mapping of Melt Homogeneity by FT Raman Profiling
The variance in Ge stoichiometry < ×> Ge of a ternary glass at a composition x = 9%, was deduced from the spectra of Figure 1 as follows.The scattering strength ratio of the CS mode near 200 cm −1 to the broad band observed near 250 cm −1 , which has contribution from predominantly Se n polymeric chains, As-centered pyramids (PYR) and -quasi-tetrahedra (QT) was deduced by analyzing Raman scattering.An empirical correlation of the Raman scattering strength ratios as a function of weighed glass composition in FT-Raman profiling experiments is generated self-consistently by having Raman results (Figure 1) from all the fully homogenized glass compositions as shown in Figure 6 by the red line joining the red filled circles.Next the profiled spectra of the glass specimen at x = 9% after 2 days of alloying was deduced and these data projected in Figure 6 as the black square data points.These data place the variance in Ge concentration across the batch after 2 days of alloying near 1.04% and provides a visual measure of glass heterogeneity.Upon alloying the glass for 7 days, the variance in Ge content in the batch decreased to 0.28%, and after 10.5 days of alloying became minuscule (<0.01%) when the batch was considered to be homogeneous.
The rigidity-and stress-elastic phase transitions in these disordered molecular networks are intrinsically sharp because of their percolative (Thorpe, 1985) nature, and glasses must be homogenized (Boolchand, 1985) at least to a scale of 50 µm or less to establish the intrinsic nature of these transitions.

Molar Volumes
Densities of our glass samples were determined by using Archimedes principle with 200 proof alcohol being the immersion fluid.We established the density of alcohol by using a Si crystal wafer of a known density (ρ si = 2.323 g/cm 3 ).We mounted the Si crystal on a quartz fiber that was suspended from a Mettler-Toledo model B154 digital microbalance, with a resolution of 0.1 mg, and weighed it in air.The quartz fiber is then submerged into the fluid and the balance tared to negate the buoyant force on the fiber.The Si crystal was then placed on the quartz fiber hook and submerged into the fluid and the weight of the Si crystal in alcohol recorded.Having the weight of the Si crystal in air and in alcohol then enabled one to calculate the density of the alcohol.  1 are normalized to the highest scattering strength band near 250 cm −1 , and the spread in the CS mode strength after 2 days of alloying provides a visual measure of glass heterogeneity.These data are plotted as the black squares in the plot above.Upon alloying the glass for 7 days, the variance in Ge content of the glass decreased to 0.28%, and then became minuscule (<0.01%) after 10.5 days of alloying.(B) An alternative representation of melt homogenization can be displayed by plotting the variance < ×> in Ge stoichiometry deduced from FT-Raman profiling of a 1.5-g batch composition weighed at x = 9% and its composition tracked as a function of t r , the alloying time, at 950 • C in days.Note that the variance decreases by 2 orders of magnitude in going from 2d (1.04%) to 10.5d (0.01%) and yields the expected weighed composition at x = 9%.
To establish the accuracy and reproducibility of our method, the density of a Ge single crystal was independently measured, and a value of 5.324(6) g/cm 3 was found, which is comparable to the published standard value of 5.323 g/cm 3 .
We emphasize that the glass samples used in density experiments are to be no <100 mg in weight.This allows the density to be measured to an accuracy of 0.25% or less.Each sample's density measurement was repeated 5 times.
We provide, as an example, the calculations used to establish the density of the alcohol, ρ alc , and the density of a Ge x As x Se 100−2x glass specimen at x = 14%, ρ sample .The molar volumes of the bulk glass were then deduced from the mass density.

First Principles Molecular Dynamics Simulations
We performed a series of First principles molecular dynamics simulations (Car and Parrinello, 1985), using the canonical (NVT) ensemble, for the computation of the mean square displacement and the diffusion constant, on the same Ge x As x Se 100−2x systems containing N = 249 atoms at 11 different compositions (2, 4, 5, 8, 9, 10, 15, 17, 18, 19, 21%) with the number of As, Ge and Se atoms fulfilling the desired stoichiometry.A periodically repeated cubic cell was used, whose size changes according to the number density of the glasses (see section Molar Volumes).For the simulation of the present liquids, we used density functional theory (DFT) in combination with plane wave basis sets.The electronic scheme has been chosen after a series of methodological investigations on different liquid chalcogenide glasses and liquids (Micoulaut et al., 2009;Massobrio et al., 2010;Bouzid et al., 2015).These showed that a generalized gradient approximation (GGA) for the exchange correlation energy improves substantially the description of both short and intermediate-range order as compared to the local density approximation (LDA) (Massobrio et al., 1999).In addition, it was demonstrated that the exchangecorrelation functional reducing the metallic character of the bonding (a Becke, Lee, Yang and Parr (BLYP) functional) led to systems with a reduced number of miscoordinated atoms and metallic character (Micoulaut et al., 2009), and improved the reproduction of structure functions obtained from neutron diffusion in the liquid and amorphous state (Bauchy et al., 2013;Micoulaut et al., 2013;Yildirim et al., 2016b).Using these recipes, the electronic structure of the As-Ge-Se liquids was described within DFT and evolved self-consistently during the motion using a GGA for the exchange and correlation parts of the total energy, according to Becke (B) and LYP, respectively (Becke, 1988;Lee et al., 1988).Valence electrons were treated explicitly, in conjunction with norm-conserving pseudopotentials to account for core-valence interactions.The wave functions were expanded at the γ-point of the supercell and the energy cutoff was set at 20 Ry.Starting configurations were taken from binary Ge-Se liquids (Yildirim et al., 2016b) and As atoms were appropriately inserted in order to meet the desired stoichiometry.Loss of the memory of the initial configurations has been achieved through preliminary runs at 2,000 K over 50 ps.For all investigated compositions and temperatures (1,500 and 1,200 K), statistical averages in the liquid state were obtained over 20-25 ps of trajectory with a time step of t = 0.12 fs and a fictitious mass of 1,000 a.u.

Glass Transition Temperature, Glass Enthalpy of Relaxation, and Melt Fragility Index
The calorimetric results on homogeneous Ge x As x Se 100−2x glasses show that T g (x) increases monotonically with x over a wide range 0 < x < 25% of compositions (Figure 5A).These results show that Ge-and As-atoms serve to cross-link (Micoulaut and Naumis, 1999) polymeric Se n chains as GeSe 2 tetrahedra both CS and Edge-Sharing (ES), and AsSe( 1/2 ) 3 pyramids form in the chemical composition range 8% <×< x c = 18.2%.At x > x c homopolar Ge-Ge and As-As bonds must form part of the same backbone as the one prevailing at x < x c, for T g (x) to continue to increase at x > x c .
The variation of the non-reversing enthalpy of relaxation at T g ( H nr ), in the glasses shows a sharply defined square-well like minimum (Figure 5B), and it is highlighted by the light blue panel, i.e., the thermally reversing window.The vanishing of the frequency corrected enthalpy of relaxation at T g , H nr, in the 9.5% < x < 17.0.% range constitutes evidence that it represents the region of phase space (Micoulaut, 2010;Bauchy and Micoulaut, 2015;Mantisi et al., 2015) where glass transitions become thermally-reversing in character.The H nr is minuscule in the IP, which is a clear evidence that the samples are dry and do not contain water impurities (Vaills et al., 2005).The dynamic reversibility of glass transitions corresponds to the overshoot in the endotherm vanishing, and is identified with the formation of isostatically-rigid local structures (Bauchy and Micoulaut, 2015), with the count of bonding constraints per atom, n c = 3, and satisfy the Maxwell-Phillips rigidity criteria (Maxwell, 1864;Phillips, 1979).Such networks are rigid but unstressed and possess a liquid-like entropy.
The variation of the melt fragility index, m(x), displaying a Gaussian like minimum (Figure 5C) with m < 20 for IP compositions, and serves to define a fragility window or, at least, a fragility minimum, in which melts are viewed to be super-strong, in the strong-fragile classification (Angell, 1995) of melts.Note that the dashed line drawn in Figure 5C marks the fragility index of SiO 2 glass taken as m = 20, and it divides the m(x) plot into three regions; (i) a low x region (x < 9.5%) and (ii) a high x region (x > 17%) and (iii) an intermediate x -region.Note that fragility window coincides with the reversibility window in x-space and the three regions mentioned above represent, respectively, the flexible phase (x < 9.5%), stressed-rigid phase (x > 17.0%), and the intermediate phase (9.5% < x < 17.0%).
Intermediate phase melts, with m < 20, are quite special; these are super-strong and their viscosity at 950 • C where we alloyed the starting materials are significantly greater than those of their non-IP composition counterparts.Finally, the close correlation in composition between melt-and glass-properties leads to the central finding of this work, viz, that melt fragility index encode glass TPs.Parallel findings have recently been noted in binary Ge-Se (or S) (Bhosle et al., 2012a,b;Gunasekera et al., 2013;Chakraborty and Boolchand, 2014) glasses and are suggestive that the melt-glass correlation could well be a general feature of glass science.

Compacted Nature of Intermediate Phase Networks
The variation in molar volumes of the present ternary glasses as a function of composition is plotted in Figure 5D.The average behavior shows V m to decrease from 18.5 to 17.5 cm 3 /mol, reflecting the increased cross-linking of the networks.The plot also shows an additional reduction of V m in the reversibility window shown by the hashed mark region.These data show that glass compositions are compacted (Boolchand et al., 2013), a feature observed in the IP of other glasses as well (Bhosle et al., 2012a,b;Chakraborty and Boolchand, 2014).The network compaction is broadly viewed as the consequence of the stressfree nature of the glass backbone where the effect of longer range forces such as van der Waals and Coulombic come into play, largely because the strong short-range covalent forces are balanced by the 3D network dimensionality.The average reduction in V m (x), with increasing x, is a consequence of the increased mean coordination number, <r> = 2 + 3x, that increases linearly with x leading to space filling, i.e., more atoms being packed per unit volume.Here the mean coordination number follows if one assumes that Ge, As, and Se atoms are, respectively, 4, 3-, and 2-fold coordinated conforming to the 8-N bonding rule.

Glass Structure Evolution With Composition From Raman Scattering
Raman scattering has been particularly useful in elucidating the glass structure evolution of the present chalcogenides.The observed Raman line shapes as a function of increasing As and Ge alloying content "x" are displayed in Figure 7.We are greatly assisted in mode assignments because of the availability (Jackson and Pederson, 1990;Pederson andJackson, 1990, 1991;Quong et al., 1993;Porezag and Pederson, 1996;Porezag, 1997;Briley and Pederson, 1998;Jackson et al., 1999) of the expected normal mode frequencies of the various As-centered and Ge-Centered local structures from NRLMOL cluster calculations (Table 1).Furthermore, Raman scattering in the Ge-Se binary glasses, intensely investigated in previous work (Bhosle et al., 2011(Bhosle et al., , 2012a)), have confirmed the Ge-centered local structure mode frequencies.The case of the As-Se binary glasses has also been investigated (Chen et al., 2010;Ravindren et al., 2014), although some ambiguities in mode assignments have persisted, issues that we will address here.
At low x (< 5%), the spectra are dominated by the broad band near 250 cm −1 coming from polymeric Se n chains and some contribution from Se 8 rings and compacted Se n chains.These display a triad of modes (234, 259, and 260 cm −1 ) as shown in Table 1 and have been confirmed in experiments on very dilute Ge doped Se glasses (Yannopoulos and Andrikopoulos, 2004;Dash et al., 2017).Detailed results in binary Ge x Se 100−x glasses (Bhosle et al., 2012a) reveal that the vibrational modes near 197 cm −1 comes from Corner-sharing and the one near 217 cm −1 from edge sharing GeSe 2 tetrahedra.With increasing x, and particularly near x = 15% the square shaped band in the 225-275 cm −1 range has contributions from modes of Ascentered pyramidal (235 cm −1 , 265 cm −1 ) and quasi-tetrahedral (198, 271 cm −1 ) units.
Perhaps the most striking result to emerge from these spectra on homogeneous glasses in the present work is the mode near 350 cm −1 (shown in the inset of Figure 8) that steadily increases Ethane like Ge 2 Se 6 units are included in the top panel, followed by those of the As-centered local structures base on pyramidal and quasi-tetrahedral local structure in the middle panel, and those of the triad of modes related to polymeric Se n Chains in the lowest panel (Jackson et al., 1999).

Ge x Se 100-x cluster
ω Th (cm −1 ) I Ram (Å 4 /amu) with x.Although the increase is slight in QT 1 , by comparing it to its bond bending vibrational mode, QT 2 , there is undoubtedly an increase in the concentration of the quasi-tetrahedra with increasing x.In a previous study (Qu et al., 2002;Qu, 2004), this particular mode was not observed in glasses probably because these were not especially homogenized.This particular mode is predicted near 364 cm −1 (Table 1) and involves stretch of the double bonded Se = As(Se 1/2 ) 3 in a quasi-tetrahedral local structure.With increasing x, the QT mode near 350 cm −1 sits at the edge of the more strongly excited mode of Ge ES units near 320 cm −1 and is thus eventually masked in the inset of Figure 8.The strongest QT mode occurs at 198 cm −1 and it overlaps with the mode of CS tetrahedral units, making a quantitative determination of the QT local structures at best difficult.

Dynamic Anomalies in the Liquid State
In the present contribution, we also focus on the dynamics of the obtained simulations.Structural properties will be considered and discussed elsewhere.The calculation of the mean-square displacement (msd) <r 2 (t)> of the three tagged atoms (As, Ge, Se) leads to a usual profile in a log-log plot representing <r 2 (t)> with time, i.e., a ballistic regime at short time where the behavior scales as t 2 and a diffusive regime that sets in for most compositions (1,500 K) at around 1-5 ps.Note that, because of the chosen high temperatures considered herein, no cage-like motion is obtained.The diffusion constant D i (i = As, Ge, Se) can be obtained using the Einstein relation limit:  1.In the high frequency range spanning 300-600 cm −1 , we illustrate in the inset evidence of the 350 cm −1 mode coming from the As-centered QT units.The red hash-marked mode represents second order scattering from the Se n chain mode near 250 cm −1 .See text.
and both D Se and D As are plotted as a function of content for the isotherms 1,500 and 1,200 K (Figure 9).Experimental diffusivity values for such system are unavailable but one can check that the values of D i are globally compatible with similar systems.For instance, it has been found (Yildirim et al., 2016a) that Ge diffusion is of about 5.10 −5 cm 2 .s−1 in Ge 20 Se 80 liquids at 1,373 K, i.e., rather close to the value obtained for a similar Se content in the present case (6.4.10 −5 cm 2 .s−1 for Ge 10 As 10 Se 80 at 1,500 K).Such values have been found to be compatible with viscosity results (Stølen et al., 2002;Gueguen et al., 2011) at high temperature once a Stokes-Einstein equation is used to estimate the viscosity (Sipp and Richet, 2002).An inspection of the behavior with composition clearly indicates that the diffusivity results contain features of the underlying TPs.D i is very small for highly cross-linked melts (x > 17%, e.g., D Se = 2.5 10 −5 cm 2 .s−1 for x = 21% and 1,500 K), increases slightly in the IP, and then becomes very large (about 10 times larger at x = 4%), with the substantial increase of D manifesting in the flexible phase.In contrast to the Ge-Se, wherein IP melts were found to display a minimum in diffusivity (Yildirim et al., 2016b), in the present Ge-As-Se, melts in the IP (9-17%) remain nearly constant in diffusivity.However, it should be noted that, for such elevated temperatures, extrema can be barely visible, as noted earlier for high temperature Ge-Se liquids (Yildirim et al., 2016a).
The nature of the dynamics is also very strongly dependent on composition and the distance of motion for the atoms (e.g., Se atoms in the inset of Figure 9) increases substantially in liquids that are flexible.For stressed-rigid phase compositions, the msd is rather small, e.g., <r 2 (t)>=19Å 2 for x = 21% at 1,200 K, corresponding to a distance of several interatomic distances only (4.5 Å) and larger distances can hardly be explored due to the large cross-link density.As the Se content increases, IP compositions seem to lead to a slight increase <r 2 (t)>, barely visible from the figure.The most noticeable change is obtained at the intermediate to flexible boundary and at both investigated temperatures, <r 2 (t)> starts to dramatically increase at 8-9% to reach nearly 200 Å 2 for x = 2%.

Correlating Melt Dynamics With Glass Topological Phases
Trends in Fragility Index and in Projected Melt Viscosity In the observed trend of m(x) (Figure 5C), we had noted before, that there is a global Gaussian-like variation of m(x) centered around the fragility window; with m < 20 in the fragility window, 9.5% < x < 17.0%, and m > 20 for compositions outside that window.The implications are that the viscosity of melts in the fragility window would be much greater than those of melt compositions outside that window.
We have used the Mauro, Yue, Allison, Gupta and Alan (MYEGA) model to calculate the viscosity of Ge x As x Se 100−2x melts using the following relation, as a function of T g /T given the measured fragility index, m, the alloying temperature of T = 950 • C, and the measured T g (x) of the glasses.In the Equation ( 5) above η ∞ is taken as log 10 η ∞ = −4 at m = 60 .We project the variation of the viscosity η(x) with melt composition in Figure 10C as the green data points.We find η(x) increases by almost two orders of magnitude in going from the fragile melt composition x = 9.0% to the super-strong melt composition at x = 15.0 % in the fragility window.In fragile melts at x > 17% the increase of T g is responsible for the increased viscosity from Equation ( 5).The broad implications of these projections are clear, viscosity of melts in the fragility window increases substantially leading melt diffusivity to reduce in that window in relation to diffusivity of fragile melt compositions at x < 9.5%.We believe these considerations are responsible for the delayed homogenization of all chalcogenide melts (Figure 2) synthesized presently.As an

Trends in Enthalpy of Relaxation at T g and Configurational Entropy Difference S config Between a Glass and Corresponding Melt
The present calorimetric results (Figure 10) unambiguously show the existence of two elastic phase transitions (Boolchand et al., 2001;Qu et al., 2002) rather than a single one, as was suspected several years ago (Thorpe, 1983;Tatsumisago et al., 1990;Böhmer and Angell, 1992), corresponding to the percolation of rigidity near <r> = 2.40.The observation of two elastic phase transitions in binary Si x Se 100−x glasses (Selvanathan et al., 1999) lead as early as 1998 to the recognition that, in real glasses, the rigidity and stress transitions will, in general, occur at two distinct and separate network connectivity.And it led to the birth of the IP as composed of rigid but unstressed networks formed between these two transitions.In the present ternary these two transitions occur near x = 9.5% or <r> = 2.285, and near x = 17.0% or <r> = 2.51.And we identify these, respectively, with the onset of rigidity and onset of stress in these glasses (Thorpe et al., 2000).The identification is consistent with the vanishing of the frequency corrected enthalpy of relaxation, which has been shown to be the signature (Micoulaut, 2010) of networks being isostatically rigid, n c = 3. Micoulaut (2010) and Bauchy and Micoulaut (2015) were one of the first to link enthalpy of relaxation at T g in calorimetric measurements with the optimally-constrained nature, n c = 3, of glassy networks.The two transitions are also observed in C p (x), which reveals a small step of 0.01 cal/g/ • C at both the rigidity-and the stress-transitions (Figure 10B).The two transitions are also observed in melt fragility index variation m(x), which show a fragile (m > 20) to superstrong (m < 20) transition near x = 9.5%, and then a superstrong (m < 20) to fragile (m > 20) transition near x = 17.0%, the stresstransition.

Correlating Melt Dynamics With Their Glass Topological Phases
The compositional variation of H nr (x) and m(x) observed (Figures 10A,C) in the present glass system demonstrates that superstrong melts in the fragility window upon cooling to T g will invariably form IP glasses in the reversibility window.On the other hand, fragile melts formed at non-IP compositions give rise to either flexible or stressed-rigid glasses.Furthermore, one can uniquely distinguish between the two kinds of fragile melts by recognizing that for glasses in the flexible phase the fragility index, m, decreases with an increase of T g , while the reverse would be the case for glasses in the stressed-rigid phase (see Figures 5A,C).The correlation observed between H nr (x) and m(x) here in the case of the present ternary has, indeed, close parallels to the case of the Ge-Se (or S) binary systems (Gunasekera et al., 2013;Chakraborty and Boolchand, 2014).Taken together, these results on especially homogenized melts/glasses are suggestive that melt dynamics encode the glass topological phases.The result that may well-represent a general principle of glass science.

Correlating Melt Fragility Index From Flexure-and Viscosity-With Present Modulated-DSC Experiments
Three decades ago, ternary Ge x As x Se 100−2x melts were investigated in terms of viscosity (Tatsumisago et al., 1990) and flexure measurements (Böhmer and Angell, 1992).Both sets of experiments made use of the same stock of glasses.These were synthesized by alloying the pure elements in evacuated quartz tubes at a T range in the 600-1,150 • C range for 24 h.The fragility index results as a function of glass composition from both sets of experiments (Böhmer and Angell, 1992) yielded remarkably similar results (Figure 11).And, for ease of comparison, we have included the present results (Figure 5C) in Figure 11 as well.It is comforting to see that in all three data sets, fragility index of melts within a narrow range of compositions near x = 13.3% or <r> = 2.40 display a strong character, while those at x < 9.5% and x > 17% become increasingly fragile.
On the other hand, the actual values of the fragility index, m, reported in the viscosity and flexure measurements are significantly greater from those measured in the present MDSC work.The difference is as large as 60 units in m near pure Se to about 30 units near <r> = 2.40 or x = 13.3 mole% and about eight units in m near the composition x = 20 in the stressed-rigid phase.We have also measured the fragility index of the present Ge-As-Se glasses using the Moynihan approach (Moynihan et al., 1974).If the Tg's in the DSC measurement are established from the inflection point of the heat flow, trends in m(x) (Figure 11) are quite similar to those found using MDSC (Figure 5C).In particular in the IP range m approach a value close to 20.
We can think of two possible factors contributing to this large disparity in melt fragility index on the same glass ternary system.The first factor could be the role of the experimental method.In viscosity and flexure, one probes the shear stress relaxation time, while in the calorimetric measurement the enthalpy relaxation time-and it has been suggested that different properties of glasses and melts may exhibit distinct relaxation times (Potuzak et al., 2011;Yu et al., 2015Yu et al., , 2017)).On the other hand what is of relevance here is the variation in the T dependence of the relaxation time, and not the relaxation time magnitude itself.For example, fragility index of melts measured by Dielectric spectroscopy (Carpentier et al., 2003), Dynamic Light scattering (Fabian and Sidebottom, 2009), and MDSC (Mohanty et al., 2019), involve measuring different kinds of relaxation times, but when these are plotted against inverse Temperature close to T g , yield nearly the same slope, i.e., fragility index.The second factor could be the sample make up.For small (1.5-g) batch sizes, as shown in Figures 1, 2, we noted at least 1.5 mole% spread in Ge stoichiometry near x = 9% for melts alloyed for 24 h at 950 • C. In the viscosity (Tatsumisago et al., 1990) and flexure (Böhmer and Angell, 1992) measurements, which require much larger samples than the ones synthesized presently, the synthesis of homogeneous bulk glasses is unlikely to be achieved with the alloying treatment described (Tatsumisago et al., 1990;Böhmer and Angell, 1992).Melt heterogeneity would invariably lead to increased fragility index as noted in recent work as well (Ravindren et al., 2014).The large differences in the magnitude of m between the present work and these earlier reports at the same composition, x, probably result from glass sample make up.Heterogeneity of melts would invariably lead the fragile component to rapidly diffuse through a melt leaving the least mobile super-strong component to be largely localized in small pockets.The delayed homogenization of chalcogenide melts, such as those reported here in the present ternary, may well be a generic factor contributing to the large spread in fragility index at the same melt/glass composition.

Aging of the Reversibility Window in Fully-With Partially-Homogeneous Ge x As x Se 100-2X Glasses
The challenges in the synthesis of homogeneous chalcogenide glasses came to the fore when FT-Raman profiling measurements (Bhosle et al., 2012a,b) were first brought to bear on understanding how melt/glass batches actually homogenize in real time upon alloying of the starting materials.The delayed nature of the homogenization over 7-10 days of alloying rather than 1 or 2 days was subsequently found to be a fairly pervasive feature of all chalcogenides examined since.And it was traced to the super-strong nature of melts in the IP with fragility index m < 20 that slow down the diffusive mixing of the batch composition (Gunasekera et al., 2013).In 2003, we had synthesized Ge x As x Se 100−2x glasses by alloying the starting materials at 960 • C for 48 h, which we now recognize only lead to partially homogeneous glasses-i.e., the variance, < ×>, of Ge stoichiometry across the 1.5 g batches are estimated near 1.5 mole% from the results of Figure 2. The reversibility window in these glasses displayed edges that were smeared when examined after 1 month of aging (Qu, 2004) for that reason.But upon long term aging (6 months, 3 years) at room temperature the reversibility window evolved and showed that, for glasses at x < 10%, the H nr (x) term increased and the reversibility window edge near x = 10% sharpened up (Figure 10A).The aging induced increase of H nr (x) at x > 17% was much less pronounced, understandably so because the T g values of these stressed-rigid glasses are measurably greater than those of the flexible ones at x < 10% (Figure 5A), leading to substantially reduced kinetics of aging as determined by the T/T g ratio.In spite of that limitation, one finds that both the rigidity-and stresstransitions became better defined upon aging, to closely approach the transitions observed in the especially homogenized glasses.
The aging behavior of the homogeneous glasses examined in the present work stands in sharp contrast to the findings on the partially homogenized glasses studied earlier.We could clearly observe sharp rigid-and stress-elastic phase transitions in the fresh rejuvenated glasses (Figure 5B).It helps one understand how homogeneity of melts/glasses is intimately tied to the observation of the IP in these glasses.Melt/glass homogenization is precursive to observation of a sharp squarewell like reversibility window permitting one to fix the three TPs (Boolchand et al., 2018).
MDSC experiments also show that the non-ergodic enthalpy of relaxation H nr (x), builds up significantly upon aging of non-IP glass compositions to about 2.0 cal/g after 3 years of room temperature aging (Qu, 2004) (Figure 10A).Thus, in fully relaxed glasses the underlying configurational entropy change associated with T g , S conf exceeds 0.08 cal/g/ • C. In the estimate, we take the width of T g to be 25 • C. Finally, the heat capacity change at T g , C p (x) associated with the reversing heat flow of 0.015 cal/g/ • C, suggests that the underlying change in vibrational entropy, S vib , is nearly half that of C p (x), a value of 0.008 cal/g/ • C. Thus S conf (= 80 mcal/g/ • C) is an order of magnitude greater than S vib (= 8 mcal/g/ • C).A correlation of strong (fragile)-melts with a small (large) magnitude of C p across T g has been suggested (Angell, 2000).That correlation is also observed in the present findings (Figure 10B).But the present results also show a far more striking correlation of melt fragility with the enthalpy of relaxation, H nr (x) of the glass.
We note that H nr (x) is vanishingly small in glasses resulting from superstrong melts, but it increases by several orders of magnitude in glasses that result from fragile melts.The H nr (x) term results from the non-reversing heat flow and it contributes an entropy change ( S conf ) that is configurational, which far exceeds S vib resulting from the jump in C p across T g , C p , which comes from the reversing heat flow.We view the presence of the rather high value of S conf term in the IP to suggest that the configurational entropy of IP melts decreases minimally ( H nr (x) = 0) for IP glasses, underscoring that networks characterizing the IP have a rather high entropy, i.e., a liquid-like entropy.On the other hand, the much higher value of H nr (x) in the relaxed non-IP glass compositions constitutes a signature that their phases have undergone substantial reconstruction upon aging and possess a lower configurational entropy than their high-T melt counterparts (Figure 5C).
In the Flexible phase (x < 10%), the large increase in H nr upon aging (Figure 5B), and in the projected melt configurational entropy S c (Figure 5C), we trace to the polymeric Se n chains.As x decreases in the 0% < x < 10% range, the Se n chains grow in length with n ∼ 250 in pure Se.The uncorrelated Se n chains characteristic of the fresh glass, particularly if they are long (n > 15) will reconstruct with each other and get correlated (Dash et al., 2017), while some of the long chains will transform to Se 8 rings, as documented (Dash et al., 2017) in the Se-rich (x < 8%) phase of binary Ge x Se 100−x glasses from Raman scattering.Specifically, the broad band near 250 cm −1 in Raman scattering of these Se-rich glasses is composed of a triad of modes identified respectively with (i) correlated Se n chains, (ii) uncorrelated Se n chains and (iii) Se 8 rings, and aging induced changes in the scattering strength of these modes have assisted in identifying the microscopic origin of aging (Dash et al., 2017) and would understandably apply to the flexible phase of the present ternary as well.
In the Stressed-rigid phase of the present ternary glasses (x > 17.0%), the glass backbone is composed of high mean coordination number <r> local structures.These include CS-(2.40) and ES-GeSe 2 (2.67) tetrahedra and Pyramidal As(Se 1/2 ) 3 (2.40).The backbone compacts as revealed by the steadily decreasing Molar volumes (Figure 5D).And aging will further lead to a compaction of the network resulting in the configurational entropy to decrease or H nr to increase (Figure 5B).On the other hand, in the IP, quite remarkably, the presence of an isostatically rigid and stress-free backbone, endows glasses with a new functionality to adapt and selforganize.And aging is qualitatively suppressed in this phase.

Glass Molecular Structure and Intermediate Phases
One of the more challenging issues in glass science is to predict the compositional windows across which networks can selforganize to display IPs and the underlying molecular structure manifestations.Some important clues on the microscopic origin of self-organization effects can be sought by examining the molecular structure of glasses.The IPs in the Ge-Se and As-Se binary systems, and the present Ge x As x Se 100−2x ternary are established.It is useful to inquire if these windows display any correlations, and, if so, what one might learn on the molecular structure of the IP from such correlations.We comment on some of these issues next.
Reversibility Windows of As-Se, Ge-Se Binary Exhaust the Window of the Present Ge x As x Se 100−2x Ternary Upon alloying traces of As-and Ge-in a base Se glass one expects the chains of polymeric Se n chains to be crosslinked by As-centered quasi-tetrahedra, Se=As(Se 1/2 ) 3 (<r> = 2.28) and As-centered pyramids, As(Se 1/2 ) 3 , (<r> = 2.40), Ge centered corner-sharing (CS)-GeSe 4 tetrahedra (<r> = 2.40), and Ge-centered edge sharing (ES)-GeSe 2 tetrahedra (<r>= 2.67) progressively.A count of bond-stretching and bondbending constraints shows that the four local structures, viz., Se = As(Se 1/2 ) 3 , As(Se 1/2 ) 3 , GeSe 4 , ES-GeSe2 are isostatically rigid (n c = 3).The presence of both As-centered and Gecentered local structures leads to chemical disorder, increases the configurational entropy of the backbone and thus assists (Yan, 2018), in formation of the IP.The presence of chemical disorder suppresses NSPS effects and thereby a T g maximum near the chemical threshold.The latter feature is observed in respective binary glasses as eluded earlier.And since the mean coordination numbers of these local structures steadily increase in the 2.28 < <r> < 2.67 range, one expects As-centered local structures to predominantly form in the 2.28 < <r> < 2.40 range and, at higher x, the Ge-centered local structures to form in the 2.40 < <r> < 2.51 range.In a qualitative fashion, one may thus view the superposition of both binary glass reversibility windows to broadly account for the wider reversibility window in the present ternary.What do we know about the IP in Ge-Se binary and As-Se binary glass systems?Bhosle et al. (2012b), reported upon the reversibility window in homogeneous binary Ge x Se 100−x glasses.They found the window to be abrupt and square-well like with the onset of rigidity near x r = 2.39 and onset of stress, x s = 2.52.That result was the culmination of an effort spread over 15 years, in which the shape of the reversibility window changed from triangular-like, to Gaussian-like to becoming Square-well like as melt/glass specimen were steadily homogenized (Bhosle et al., 2012b), and the variance in Ge-stoichiometry was systematically reduced from about 1.0-0.05%as the alloying time of melts at 950 • C was increased from 2 days to about 7 days.
In the first MDSC study of As-Se glasses (Georgiev et al., 2001) where the existence of a reversibility window was observed (Figure 12) in the 29% < x < 37% range, binary melts were alloyed near 700 • C for 3 days and then water quenched.The glasses were never T g cycled after the quench but merely aged at room temperature, thereby offsetting the segregation effect of the eutectic.And the reversibility window observed is characteristic of glasses composed of As-centered pyramidal and quasi-tetrahedral local structures.Here, one must recall that the Eutectic temperature, T E (= 147 • C) is close to the T g (= 104 • C) of glasses near x = 20% of As.It is for this reason that in aging studies of the reversibility window in subsequent studies (Golovchak et al., 2008;Ravindren et al., 2014) the integrity of the window was not preserved.In the more detailed study of Ravindren et al. (2014) where glasses were especially homogenized by alloying melts for extended periods, they were T g cycled after a water quench.In that study, one also found the reversibility window to steadily smear, particularly near x = 20%, the eutectic composition (Figure 13), as glasses aged and steadily segregated (Ravindren et al., 2014).In a glass forming system wherein a eutectic is present close to the reversibility window, one expects thermodynamic considerations based on segregation to override self-organization effects, particularly suppression of aging of the enthalpy of relaxation due to segregation effects.
Although the existence of the Eutectic in As-Se binary, as well as P-Se binary, is well-recognized (Okamoto, 2000), the role of the eutectic on self-organization effects has received recognition only recently.And it may be helpful to add a few comments on that topic for completeness.The eutectic is signature of segregation of respective melts into Se-rich and As-rich regions.We believe that the As-rich regions consist of ethylene-like (As) 2 (Se 1/2 ) 4 local structures containing homopolar As-As bonds in a Serich glassy matrix.The evidence for these local structures have emerged recently (Ravindren, 2015) from close analogies in Raman scattering of the P-Se and As-Se binary systems.They have suggested that ethylene like (As) 2 (Se 1/2 ) 4 local structure in analogy to (P) 2 (Se 1/2 ) 4 local structure exist, for which independent confirmation is given by MD simulations (Bauchy and Micoulaut, 2013) and 31 P NMR (Eckert, 1989).
To summarize, we illustrate in Figure 13, Topological phases of the Ge x Se 100−x and As x Se 100−x binary systems and the present Ge x As x Se 100−2x ternary on the Ge-As-Se composition triangular plot (Qu et al., 2002;Wang et al., 2009).The IP in As-Se binary (green, label A) and the Ge-Se binary (purple, label B) when combined add up to the IP of the present ternary as illustrated by the vertical lines delineating the <r> markers at 2.28, 2.40, and 2.51.Noteworthy in the plot is the onset of NSPS in the binary glass systems near the stoichiometric chemical compositions (Boolchand et al., 2002) GeSe 2 and As 2 Se 3 .In these NSPS regions it is unlikely that the network backbone connectivity has a unique value, often estimated by the usual 8-N bonding rule, taking the coordination number of Ge, As, and Se to be 4, 3, and 2, respectively.
Molecular Structure Evolution of Stressed-Rigid Networks in Binary Ge-Se and Ternary Ge x As x Se 100-2x Glasses Compared and Contrasted The unusual nature of glass structure evolution in the binary Ge x Se 100−x displaying molecular clustering in the Se-rich range (x < 33.3%) emerged from Mössbauer spectroscopy experiments using 129m Te tracer (Bresser et al., 1986).Both Te and Se atoms belong to the VI th column of the periodic table, are isovalent and take on a 2-fold coordination here.The larger covalent radius of Te (of 138 pm) than Se (of 120 pm) has profound consequences; Te selectively replaces available Se sites possessing a larger free volume, such as surface sites of characteristic molecular clusters rather than cluster interior sites (Bresser et al., 1986).
In these Mössbauer experiments, the 129 I daughter A-site represents the chemically ordered site, formed from a parent 129m Te atom bonded to 2 Ge near neighbors (nns), while the 129 I daughter B-site results from a 129m Te parent atom bonded to either two Se nns or a Se and a Ge nn (Bresser et al., 1986).In a continuous random network (CRN) description of these glasses, one expects the site intensity ratio N B /N A (x) as a function of Ge content "x" to show a smooth broken-line curve starting from an infinite value at x = 0, to approaching 0 at x = 33.3%for a chemically ordered stoichiometric GeSe 2 glass (Figure 14B).At x = 0, all Te atoms would have two Se nns (i.e., only B sites are possible), while at x = 33.3% each Te atom would have  (Georgiev et al., 2000) and Ge x As x Se 100−2x ternary projected on the Ge-As-Se compositional triangle.The glass forming region is taken from ref. Z.U.Borisova (Borisova, 1981).Starting from pure Se, as we move along the equimolar composition x, one spans the Flexible phase (0 < < r> < 2.285), the Intermediate Phase (2.285 < <r> < 2.51) followed by the stressed-rigid phase (<r> > 2.51), and with the red data points showing the 20 glass compositions examined in the present study.Note the black line joining the As 2 Se 3 glass composition to the GeSe 4 one, is the locus of glass compositions with <r> = 2.40, and it splits the IP of the equimolar ternary into almost nearly equal parts, labeled A (hashed green) and B (hashed purple), with the segment A (B) nearly overlapping with IPs observed in As x Se 100−x (Ge x Se 100−x ) binary systems, respectively.The broken line of purple color, joining As 2 Se 3 to GeSe 2 composition, represents the join of chemical threshold compositions (x chem ), with NSPS setting in as x approaches x chem in both binary glass systems but not the present equimolar ternary as sketched.Feature E centered near x = 20% of As in the As x Se 100−x join is the Eutectic in that binary system.See text.
only Ge nns if the GeSe 2 glass is chemically ordered, so that only A sites are possible.In the Mössbauer effect experiments, on the other hand, trends in I B /I A (x) mimic the expectation of a CRN model only at low x < 15% of Ge, but, as x increases (> 15%), the site intensity ratio I B /I A (x) displays a local maximum near x = x c = 23%, corresponding to the centroid of the IP, where molecular clustering has evolved, and the I B /I A (x) ratio remains finite at 1.60(4) as x = 33.3%.The results underscore that the stoichiometric glass is not chemically ordered (Bresser et al., 1981(Bresser et al., , 1986;;Boolchand et al., 1982;Boolchand and Bresser, 2000;Salmon and Petri, 2003).It is composed of large Se-rich clusters based on the high-T form of c-GeSe 2 (α-GeSe 2 ) (Bridenbaugh et al., 1979;Boolchand, 1985) whose surfaces are intrinsically reconstructed in the glass to form Se-Se dimers, which may be viewed as Topological defects.The late Dr. K. Murase correctly identified that the Raman stretch mode of Se-Se dimers on the surface of reconstructed clusters is near 246 cm −1 (Boolchand, 2000a).
Presence of these Se-Se bonds, on grounds of stoichiometry, requires Ge-Ge bonds to exist in stoichiometric GeSe 2 glass.Evidence of such bonds in small Ge-rich clusters based type of a variation at low x (<10%).But at higher x (>15%) a Molecular Clustered Network (MCN) prevails as the observed site intensity ratios I B /I A (x) far exceed the expected N B /N A (x) site concentration ratio based on a CRN (Bresser et al., 1986).See text.Also note that as x increases to 33.3% of Ge, the finite value of the I B /I A ratio of 1.66 (10) shows the stoichiometric GeSe 2 glass has an intrinsically broken chemical order (Bresser et al., 1986).
on ethane-like Ge 2 Se 6 local molecular units, emerged from Raman scattering, and independently from 119 Sn Mössbauer effect experiments.Since Sn is isovalent to Ge, 119 Sn Mössbauer Effect experiments probe the chemically specific Ge-local structures (Boolchand and Bresser, 2000) formed in the glass, which includes both the tetrahedrally coordinated Ge sites as well as the non-tetrahedrally coordinated Ge sites as in the ethane-like Ge-rich cluster possessing Ge-Ge bonds.
The molecular clustering picture of binary Ge x Se 100−x glasses developed above is elegantly corroborated in the measured V m (x) results that show a qualitative increase in the stressed-rigid regime at x > 26%.That behavior in V m (x) in Ge-Se binary differs qualitatively from the one observed in the present ternary chalcogenide.It raises the questions what is so special about the structure of these stressed-rigid binary Ge-Se glass compositions that their network structure leads to packing effects that are so different?The issue was at the center of controversy in the mid-1980s (Bridenbaugh et al., 1979;Bresser et al., 1981Bresser et al., , 1986;;Murase et al., 1983).In the equilibrium phase diagram of the Ge-Se binary there exist two crystalline phases of relevance here (i) a low T metastable structure of Ge 4 Se 9 stoichiometry (x = 30.8%) of V m = 17.62 cm 3 (Fjellvåg et al., 2001), and a high-T form α-GeSe 2 (x = 33.33%) of V m = 17.69 cm 3 (Burgeat et al., 1975).The α-GeSe 2 phase consist of a stacking of 2D sheets, with each sheet composed of chains of CS GeSe 2 tetrahedra that are laterally coupled by a pair of ES units (Phillips et al., 1980).The structure of Ge 4 Se 9 is strikingly similar to that of α-GeSe 2 , except the pair of ES tetrahedra linking the CS chains are altered in a subtle way: these become CS tetrahedral units linked by an additional Se-Se dimer bond.Returning to the binary Ge x Se 100−x glasses, once x > 31.5%, one observes ethanelike Ge-rich cluster to appear in 119 Sn Mössbauer spectroscopy (Boolchand et al., 1982;Boolchand and Bresser, 2000), as well as in Raman scattering (Bhosle et al., 2012a) where evidence of a mode near 180 cm −1 associated with these Ge-rich cluster is observed.Clearly, as the Ge content of glasses x exceeds 31.5%,much of the Ge additive now goes into forming the Ge-rich clusters that decouple from the Se-rich backbone, thus arresting the growth of T g (x), and one observes the slope dT g /dx to display a maximum at that composition (Bhosle et al., 2012b).The onset of backbone demixing as x > 31.5% leads the non-bonding van der Waals interaction between the clusters to grow.And we view the rapid increase of V m (x) of binary Ge x Se 100−x glasses in the 2.62 < <r> < 2.67 range (Figure 14A) as a direct consequence of demixing of the backbone into small Ge-rich and large Se-rich clusters.And as x = 33.3%,dT g /dx = 0, leading to a maximum of T g , which constitutes direct evidence for the stoichiometric glass to be segregated or NSPS (Boolchand, 2000b;Boolchand et al., 2002).Such segregation is also observed in binary As-Se, Ge-S, and As-S glasses near the chemical threshold, but it notably absent in binary Si-Se glasses (Boolchand, 2000b).
In the present ternary, upon alloying As an additive in the Ge-Se binary, the NSPS effects are also completely suppressed.Trends in V m (x) are qualitatively reversed.In the present ternary V m (x) decreases with x (Figure 14A), as more connected local structures based on As-centered -QT units and -PYR units, Gecentered -CS and -ES tetrahedral units emerge, and chemically disorder the network backbone.The presence of this chemical disorder precludes the growth of the 2D layer form of α-GeSe 2 fragments or the 2D form As 2 Se 3 like 6 membered rings (Kim and Boolchand, 1979;Phillips et al., 1980;Boolchand et al., 1986) to nucleate and form part of the backbone.It is for this reason the glassy backbone evolves in a topology that is Continuous Random Network (CRN) like steadily compacting with increasing x or <r>.And remarkably no evidence of NSPS is evidenced as suggested by the absence of a T g maximum near the chemical threshold (x = 18.18% or <r> = 2.55).Furthermore, an additional reduction of V m in the IP range (2.28 < <r> < 2.53) (Figure 14A) from the average value of V m (x) outside the IP is observed.It emerges because the strong covalent bonding forces (bond-stretching and bond-bending) due to the known local structures satisfy the Maxwell-Phillips criteria (n c = 3) for rigidity, permitting the weaker and long range van der Waals forces to come into play and compact (Boolchand et al., 2013) glassy backbones formed in that isostatically rigid reversibility window (Figure 14A).

Intermediate Phase Glasses as Ideal Glasses
Adam-Gibbs (Adam and Gibbs, 1965) suggested that the viscosity of a melt increases as its excess configurational entropy (S c ) decreases as it is cooled to T g described by the relation, Since the viscosity of a melt is viewed to be a thermallyactivated process, and from the two relations above, one must then require that the inverse of the activation energy for shear relaxation time provides a measure of the excess entropy loss (S c ) upon cooling of a melt, From Equation ( 9), one recognizes that the inverse of the activation energy for shear relaxation time provides a measure of the excess configurational entropy of the melt as it cools.
Since the fragility index deduced from calorimetric measurements are often correlated with those inferred from viscosity or shear relaxation times, one then expects the inverse activation energy, 1/E a , from calorimetric measurements to provide a measure of configurational entropy loss as a melt is cooled.
The inverse activation energy for enthalpy relaxation (1/E a ) in the present Ge-As-Se chalcogenide melts as a function of composition, Figure 5C, displays a maximum near the center of the IP.The result unambiguously shows that as an IP melt cools from a high T to T = T g , the configurational entropy S c remains high, and indeed, that is also the reason the resulting glass has a near vanishing non-reversing enthalpy, H nr (x).Such glass compositions are viewed as possessing liquid-like structures with a high S c and for that reason show minimal aging (Boolchand and Goodman, 2017).The low propensity for aging of isostatic glasses also echoes the fact that isostatic systems have been shown to exhibit minimum creep deformations under sustained load (Pignatelli et al., 2016;Bauchy et al., 2017).On the other hand, the non-IP glass compositions have their structure evolve with either polymeric Se n chains (as in t-Se, entropy sinks) forming at x < 9.5% or CS and ES-GeSe 2 tetrahedra and pyramidal As(Se 1/2 ) 3 units forming at x > 17%, and leading the H nr (x) term to increase, and thus lower the configurational entropy of the melt/glass (Figures 5B,C).IP glass compositions display a high glass-forming tendency, show minimal aging, and are thus Ideal glasses.These ideas of IP glasses as ideal glasses are at variance with the notion that ideal network glasses possess low configurational entropy and form at T close to the Kauzmann temperature (Angell, 1995).On the other hand, the high configurational entropy of IP glass compositions and their functionality to adapt are much closer to the recent ideas of a competition between entropy and energy to stabilize the intermediate phase (Yan, 2018).

CONCLUSIONS
The specially homogenized nature of the bulk Ge x As x Se 100−2x glasses synthesized in the present work have lead the frequency corrected Enthalpy of relaxation at T g , H nr (x), to display a square-well-like reversibility window even in the fresh (rejuvenated) state.Upon room temperature aging for 8 months, the window gets substantially deeper.The reversibility window is bordered by an abrupt rigidity-transition, near x r = 9.5(2)% and an equally abrupt stress-transition near x s = 17.0(3) %.
The specially homogenized nature of the bulk glasses synthesized in the present work has also led molar volumes (V m ) across the IP composition range, 9.5% < x < 17%, to show an additional small but systematic reduction of 2.5% in V m in the IP from the average behavior of V m across the range of compositions examined.The result demonstrates the compacted nature (Boolchand et al., 2013) of the isostatic networks formed in that special phase.
The melt fragility index, m(x), deduced from MDSC experiments show a Gaussian-like global minimum with m(x) < 20 for IP compositions, defining the fragility window.The fragility window coinciding with the reversibility window, underscores that super-strong melts give rise to IP glasses, while fragile melts to either flexible or stressed-rigid glasses, leading to the conclusion that melt dynamics encode glass Topological phases.
Trends in the projected melt viscosity, η(x) using the measured fragility index, m(x), and the MYEGA formalism, show almost three orders of magnitude increase in η(x) in the 0 < x < 20% range.Trends in the inverse of the projected viscosity, 1/η(x), broadly parallel the melt diffusivities D calculated from ab-initio MD simulation.These results are strongly suggestive that the delayed homogenization of melts during synthesis is due to the low diffusivities or high viscosity of melts at 950 • C in the 13% < x < 25% composition range.
Finally, the inverse of the activation energy for enthalpy relaxation time shows a maximum near x = 14%, the center of the IP melt compositions, strongly suggesting that such melts possess a high liquid-like configurational entropy, S c , in the Adams-Gibbs picture.The optimization of the glass forming tendency for IP compositions, the stress-free nature of these structures in which three types of isostatically-rigid local structures proliferate, and the qualitative suppression of aging in that phase, the high configurational entropy of IP networks, each point to the ideal nature of glasses formed in the IP.Some of these exceptional functionalities of IP networks observed in experiments, appear to find support in

FIGURE 1 |
FIGURE 1 | FT-Raman spectra of a Ge 9 As 9 Se 82 glass composition profiled after being alloyed at 950 • C in (A) 2d, in (B) 7d, and in (C) 10.5 d.The variance

FIGURE 3 |
FIGURE 3 | MDSC scan of a Ge x As x Se 100−2x glass at x = 11%, showing the total, reversing and non-reversing heat flow as curves labeled 1, 2, and 3 in the heating cycle.The corresponding heat flow terms in the cooling cycle are labeled as curves 4, 5, and 6.See text for the operating parameters.The mean value of T g of 151.3 • C in the heating cycle and of 143.7 • C in the cooling cycle yields the scanning rate independent T g .The frequency corrected enthalpy of relaxation, H nr (x) is obtained by subtracting the exothermic contribution shaded in red (curve 4) from the endothermic one (in curve 3).The typical operating conditions are as follows; linear T ramp at 3 • C/min, modulation time of 100 s, and modulation amplitude at ±1 • C.

FIGURE 4 |
FIGURE 4 | (A) Shows the variation of the C p ip and C p op contributions as a function of modulation frequency.As the modulation frequency increases in going from a t m of 120-60 s, the C p op peaks corresponding to ω•τ = 1, steadily shifts to higher temperature underscoring the melt has to relax faster (τ decreases) to track the increased (ω) modulated heat flow.(B) Shows a plot of log of τ as a function of T g /T permitting one to deduce the fragility index from the slope (Equation1).

FIGURE 5 |
FIGURE 5 | Displays the observed compositional trends in (A) glass transition temperature, T g (x), (B) the non-reversing enthalpy of relaxation at T g , H nr (x), (C) melt fragility index, m(x) and (D) Molar Volumes, V m (x), of especially homogenized Ge x As x Se 100−2x ternary glasses.The green striped region represents the volumetric window where the average behavior deviates from the mean behavior showing a compaction in the intermediate phase.The blue panel highlights the IP-compositions spanning the range 9.5% < x < 17.0%.
FIGURE 6 | (A)From the fully homogenized glasses we established the red, empirical correlation of the Raman scattering strength ratios with weighed glass composition once that composition has been homogenized in the FT-Raman profiling experiments, describing the variation in the scattering strength ratio of the Corner-Sharing mode near 200 cm −1 to the broad Se n chain band near 250 cm −1 .The six spectra plotted in Figure1are normalized to the highest scattering strength band near 250 cm −1 , and the spread in the CS mode strength after 2 days of alloying provides a visual measure of glass heterogeneity.These data are plotted as the black squares in the plot above.Upon alloying the glass for 7 days, the variance in Ge content of the glass decreased to 0.28%, and then became minuscule (<0.01%) after 10.5 days of alloying.(B) An alternative representation of melt homogenization can be displayed by plotting the variance < ×> in Ge stoichiometry deduced from FT-Raman profiling of a 1.5-g batch composition weighed at x = 9% and its composition tracked as a function of t r , the alloying time, at 950 • C in days.Note that the variance decreases by 2 orders of magnitude in going from 2d (1.04%) to 10.5d (0.01%) and yields the expected weighed composition at x = 9%.

FIGURE 7 |
FIGURE 7 | Observed FT-Raman scattering results on Ge x As x Se 100−2x ternary glasses as a function of increasing Ge and As content "x."The expected vibrational modes are labeled.Jackson et al. have calculated these normal modes by NRLMOL code summarized in Table 1.In the 300 cm −1 to the 600 cm −1 range we see weakly excited modes that are presented and analyzed in Figure 8.The assigned modes are Eth, Ethane; CS, Corner Sharing; ES, Edge Sharing; Se n , Selenium Chain; PYR, Pyramidal; Se 8 , Selenium 8 rings.

FIGURE 8 |
FIGURE 8 | Raman line shape deconvolution at several glass compositions in the low x regime of the present ternary.These modes are in harmony with the expected vibrational modes of the characteristic local structures as shown in Table1.In the high frequency range spanning 300-600 cm −1 , we illustrate in the inset evidence of the 350 cm −1 mode coming from the As-centered QT units.The red hash-marked mode represents second order scattering from the Se n chain mode near 250 cm −1 .See text.

FIGURE 9 |
FIGURE 9 | Average self-diffusion constant (broken lines) and species dependent (As, red; Se, green) self-diffusion constant D(x) for two isotherms (1,500 and 1,200 K) as a function of content x in Ge x As x Se 100−2x glass-forming liquids.The inset represents the Se average mean-square displacement after 18 ps simulation time.Gray zones represent the location of the intermediate phase.

FIGURE 10 |
FIGURE 10 | Observed variation of (A) the Non-reversing heat flow, H nr (x) in present glasses (•) compared to previous work of Tao and Boolchand ( , •), (B) the specific change C p (x), and (C) the fragility index, m, in the especially homogenized ternary Ge x As x Se 100−2x glasses in the present work.Also projected is the variation of melt viscosity, η at 950 • C expected using the measured m and the MYEGA equation.

FIGURE 11 |
FIGURE 11 | Variations in the fragility index of equimolar Ge x As x Se 100−2x ternary melts reported by Tatsumisago et al. ( ) and by Böhmer and Angell ( ) compared to those reported in the present work using MDSC (•) and DSC ( ).

FIGURE 12 |
FIGURE 12 | Reversibility window in Ge x As x Se 100−2x after 8 months of aging compared to the windows in As-Se and Ge-Se binary glasses taken from Georgiev et al. (2000) and Bhosle et al. (2012b).

FIGURE 13 |
FIGURE 13 | Topological phases observed in Ge x Se 100−x binary(Bhosle et al., 2012b), As x Se 100−x binary(Georgiev et al., 2000) and Ge x As x Se 100−2x ternary projected on the Ge-As-Se compositional triangle.The glass forming region is taken from ref. Z.U.Borisova(Borisova, 1981).Starting from pure Se, as we move along the equimolar composition x, one spans the Flexible phase (0 < < r> < 2.285), the Intermediate Phase (2.285 < <r> < 2.51) followed by the stressed-rigid phase (<r> > 2.51), and with the red data points showing the 20 glass compositions examined in the present study.Note the black line joining the As 2 Se 3 glass composition to the GeSe 4 one, is the locus of glass compositions with <r> = 2.40, and it splits the IP of the equimolar ternary into almost nearly equal parts, labeled A (hashed green) and B (hashed purple), with the segment A (B) nearly overlapping with IPs observed in As x Se 100−x (Ge x Se 100−x ) binary systems, respectively.The broken line of purple color, joining As 2 Se 3 to GeSe 2 composition, represents the join of chemical threshold compositions (x chem ), with NSPS setting in as x approaches x chem in both binary glass systems but not the present equimolar ternary as sketched.Feature E centered near x = 20% of As in the As x Se 100−x join is the Eutectic in that binary system.See text.

FIGURE 14
FIGURE 14 | (A) Trends in molar volumes, V m (x) of binary Ge x Se 100−x glasses compared to those of the present Ge x As x Se 100−2x ternary glasses.(B) 129 I Mössbauer site intensity ratios I B /I A (x), deduced from doping 129m Te tracer in Ge x Se 100−x glasses reveal a Continuous Random network (CRN)type of a variation at low x (<10%).But at higher x (>15%) a Molecular Clustered Network (MCN) prevails as the observed site intensity ratios I B /I A (x) far exceed the expected N B /N A (x) site concentration ratio based on a CRN(Bresser et al., 1986).See text.Also note that as x increases to 33.3% of Ge, the finite value of the I B /I A ratio of 1.66 (10) shows the stoichiometric GeSe 2 glass has an intrinsically broken chemical order(Bresser et al., 1986).
Molar Volume of Ge x As x Se 100-2x x = 14%

TABLE 1 |
Vibrational modes of the Ge-centered CS and ES tetrahedra and