Metallic Glasses : A New Approach to the Understanding of the Defect Structure and Physical Properties

The work is devoted to a brief overview of the Interstitialcy Theory (IT) as applied to different relaxation phenomena occurring in metallic glasses upon structural relaxation and crystallization. The basic hypotheses of the IT and their experimental verification are shortly considered. The main focus is given on the interpretation of recent experiments on the heat effects, volume changes and their link with the shear modulus relaxation. The issues related to the development of the IT and its relationship with other models on defects in metallic glasses are discussed.


Introduction
Metallic glasses (MGs) constitute an amazing example of a man-made non-crystalline state, which is not observed in nature.These materials are very promising in the sense of technological applications that results in the growing interest to the investigation of their physical properties.At first glance, since MGs do not have directional chemical bonding, one can expect that the understanding of their structure and properties should be a simpler task than for other types of glasses.However, any commonly accepted theory describing their formation and main structural features has been absent thus far, and any general theory of non-crystalline substances is still lacking as well.This largely constrains the development of new type MGs with the physical properties predicted in advance.Since MGs are prepared by melt quenching while the melt is obtained by fusion of the crystalline state, one should naturally expect a relationship between these states, including a connection between the properties of the glass and those of the maternal crystal (whose melt is used for the glass production).These interdependencies should be taken into account by a theoretical model of glass.On the other hand, since it is commonly accepted that glass is a frozen liquid, this model should also imply certain relation to the melting mechanism and formation of the liquid state.By that, it should include as a major ingredient the notions on structural defects, which are intrinsically related to the whole glass prehistory, i.e., maternal crystal→melt→glass.To date, quite a few theoretical models describing the structural features, defects and different properties of MGs have been suggested [1][2][3][4][5][6][7][8][9][10][11][12].In our opinion, however, the above general requirements are satisfied by another approach-the Interstitialcy theory (IT) of condensed matter states suggested by Granato [13,14].It was demonstrated in recent years that the IT provides a powerful tool for the understanding and predicting different relaxation phenomena in MGs and unambiguously shows a genetic relationship of the glass with the maternal crystal.Therefore, this brief overview is firstly devoted to an analysis of the major hypotheses of the IT and main experiments related to its verification.Since certain issues related to the IT have already been discussed [15,16], this work is largely devoted to the new experiments and their interpretations while earlier information is mentioned briefly.Finally, a relation of the IT to other models of the metallic glass structure and its defects is discussed.

The Interstitialcy Theory
The Interstitialcy theory is essentially based on the experimental investigations of the effect of low-temperature (4 K), low-dose soft neutron irradiation on the elastic moduli of single crystal copper carried out by Granato's group in the 1970s [17,18].The irradiation leads to the formation of isolated Frenkel pairs (vacancy+interstitial) remaining mostly immobile in the liquid helium temperature range.A careful analysis of the experimental data led to the conclusion that the formation of interstitial defects results in a strong decrease of the shear modulus (diaelastic effect) C 44 according to ∆C 44 /C 44 = −B i c i , where c i is the interstitial concentration and B i ≈ 30 is the shear sensitivity.For the vacancies, the shear sensitivity is by an order of magnitude smaller, B v ≈ 2. The effect of point defects on the bulk modulus is also small, comparable to the diaelastic effect produced by the vacancies.At the same time, similar results were obtained upon electron irradiation of single crystal aluminum [19].It was argued that a strong diaelastic effect can be observed only if the interstitials do not make up an octahedral configuration, as considered before the 1970s (red circle in Figure 1a), but form a dumbbell (split) structure, which is characterized by the two atoms trying to occupy the same lattice cite (red circles in Figure 1b).To underline the difference between these defects, the latter defect is also called an interstitialcy [20].It is of major importance that the octahedral interistitial is a spherically symmetrical defect (similarly to the vacancy) while the dumbbell interstitial is clearly anisotropic and, thus, constitutes an elastic dipole strongly interacting with the external shear stress [19,21].It is this feature of the defect, which leads to a strong diaelastic effect.It is now generally accepted that dumbbell intertitials exist in all main lattices and constitute the defect state with the lowest formation enthalpy [22][23][24][25][26].An important consequence of the interstitial dumbbell structure consists in the appearance of the low-frequencies in the vibrational spectrum, which are by a few times smaller than the Debye frequency [28].This in turn leads to the high formation entropy of the defect (for copper, S f ≈ (5 − 10) k B according to the Granato's estimate [14]) and, respectively, to a decrease of the Gibbs formation free energy at high temperatures.The extrapolation of the elastic moduli of irradiated crystals towards high defect concentration in the experiments [17,18] showed that C 44 → 0 at c i ≈ 0.02 ÷ 0.03 providing a guess that such big defect concentration should lead to the crystal→liquid transition, because the liquid is characterized by a vanishing (or very small) shear modulus [29,30].This allowed Granato to suggest that melting of metals should be related to the rapid thermoactivated generation of dumbbell interstitials.Another important point realized by Granato was the understanding that the defect formation enthalpy is proportional to the unrelaxed shear modulus G, in line with earlier investigations [31,32].The above hypotheses led Granato to the formulation of the interstitialcy theory, which includes melting as an integral part [13,14] (see also a discussion given below).In fact, the mathematical formalism of the IT is based on the two equations, where H F i is the interstitial formation enthalpy, α is a dimensionless parameter close to unity, µ is the shear modulus of the defectless crystal, Ω is the volume per atom and B i = αβ with β being the dimensionless shear susceptibility.This quantity was estimated by Granato as β = −3C 4444 /C 44 ≈ 40, where C 44 is the shear modulus of the crystal and C 4444 is its fourth-rank (anharmonic) shear modulus.On the other hand, the shear susceptibility is related to the internal energy U of the crystal as 4 , where ε is the shear strain.The above equations show that the shear susceptibility constitutes a fundamental parameter of the material since it is proportional to the ratio of the fourth-rank shear modulus to the "usual" shear modulus.
The general approach of the IT to the crystal→melt→glass transformation and the relationship between these states consists in the following.Crystal melting is related to a rapid increase of the concentration of dumbbell interstitials, which remain identifiable structural units in the liquid state (as confirmed by later computer modeling [33]).Rapid melt quenching freezes the melt defect structure in the solid glass and different relaxation processes occurring in it upon heat treatment can be interpreted in terms of the changes of the defect concentration by using Equations ( 1) and (2).For the glassy state, the quantities G and µ in these equations correspond to the shear moduli of glass and maternal crystal, respectively.Despite the simplicity of these equations, it has been found that the IT, although originally derived for simple metals, provides explanations for many experiments on multicomponent MGs, as reviewed earlier [15,16] and discussed below.

Shear and Dilatation Contributions into the Defect Elastic Energy
The IT is actually built on the hypothesis that the shear component of the elastic energy created by interstitial defects is predominant while the dilatation contribution can be neglected.This agrees with later calculations by Dyre [34] who showed that it is the shear strain component, which produces the main contribution into the elastic energy far from a point defect in a solid.He derived a simple relation for the ratio of the dilatation U dil and shear U shear components of this energy and concluded that the former is much smaller, i.e., where B is the bulk modulus.This equation, however, was derived within the linear elasticity approach for a spherically symmetric defect and does not account for the energy of the defect nucleus.For further verification of the above Granato's hypothesis, a molecular-static modeling of interstitial defects in four FCC metals was performed [27].To compare the contributions U dil and U shear into the total elastic energy, the local relative change of the Voronoi polyhedra volume V i for each atom with respect to the Voronoi polyhedra volume of an atom in the ideal lattice V 0 was accepted as a measure of the volume change upon defect formation, i.e., . The shear component of the elastic energy was taken as a difference U shear = H f − U dil with H f being the interstitial formation enthalpy.The calculations showed that the ratio U dil /U shear for the interstitial dumbbell is nearly twice as that given by Equation (3) due to the accounting of the elastic energy of the defect nucleus.This correction allows using this equation for 63 elemental metals.The result shown in Figure 2 implies that this ratio does not exceed 0.15 for more than 90% of metals.Assuming that MGs contain similar interstitial-type defects in line with the IT, the same calculation procedure was applied to 189 metallic glass composition and nearly the same result was obtained (see Figure 2).Thus, in both cases, the contribution of the dilatation energy is indeed much smaller than that given by the shear energy, just as supposed by Granato [13,14].The same data for interstitial-type defects in 189 metallic glasses are also shown [27].With permission from JETP Letters, 2019.

Increase of the Interstitial Concentration before Melting
As mentioned above, the IT argues that melting of metals is related to a rapid increase of the concentration of dumbbell interstitials.This is a crucial statement of the theory.Since dumbbell interstitials produce tenfold bigger shear softening as compared with vacancies, this effect can be detected experimentally provided that there are no contributions coming from other defects (e.g., dislocations) in the crystal.Such experiments were recently performed on single crystal aluminum and coarse-crystalline indium [35,36].The main results of these measurements are shown in Figure 3.Despite the usual opinion that the equilibrium interstitial concentration c i is negligible at any temperature [37], it is seen that c i rapidly increases for both metals upon approaching the melting temperature T m .In aluminum, this concentration remains smaller than the vacancy concentration c v while for indium c i becomes even bigger than c v near T m .Besides that, the data on Al clearly demonstrate an increasing tendency of c i -growth at high temperatures T ≥ 926 K (0.99 T m ) while the formation Gibbs free energy start to rapidly decrease in this region [35].These features agree with the predictions of the IT.It is worth noting that increasing concentration of dumbbell interstitials can explain the non-linear growth of the heat capacity of simple metals near T m , whose nature remains hitherto unclear [38].It should also be mentioned that temperature dependence of the interstitial concentration obtained for aluminum [35] allowed an estimate of their formation entropy S F i , which was found to be ≈7 k B [39], in full agreement with Granato's value [14].[35,36].
With permission from Pleiades Publishing, LTD, 2019.

Identification of Interstitial-Type Defects in the Glassy State
The topological pattern of dumbbell interstitials in crystals is very clear-two atoms, trying to occupy the same lattice cite (Figure 1b).However, any similar topological picture in the glassy state is absent.In this case, one can try to identify these defects by searching structural regions, which display the properties similar to those of dumbbell interstitials.Thus far, such attempts have been performed using computer models of glassy copper and aluminum [40,41].It was found that certain nano-sized regions reveal large non-affine displacements and can be characterized by a strong sensitivity to the applied shear stress and distinctive local shear strain fields, which are described by the local shear susceptibility as well as by the diaelastic compliance and diaelastic polarizability tensors.Another feature consists in the characteristic low-and high-frequency modes (far below and above the Debye frequency, respectively) in the vibration spectra of atoms belonging to these regions.All these peculiarities are quite similar to those of dumbbell interstitials in crystals.Thus, interstitial-type defects indeed exist in non-crystalline simple metals.Two-component metallic structures should be analyzed in a similar way.Besides that, numerous simulations of MGs show the presence of the atoms characterized by < 0, 2, 8, 0 > Voronoi indexes (or close to them) [7], which constitute a characteristic feature of dumbbell interstitials in crystals (Figure 1b).One the other hand, an interstitial defect is characterized by the two atoms with these indexes.We are unaware of any attempts for searching two adjacent atoms with < 0, 2, 8, 0 > or close indexes in computer models of MGs.

Shear Susceptibility
The large magnitude of the shear susceptibility β (see Equation ( 2)), which determines the influence of defects on the shear modulus (diaelastic effect), constitutes a salient ingredient of the IT.This quantity is controlled by the non-linearity of the solid, specifically by the magnitude of the non-linear shear modulus C 4444 or, in an isotropic approximation, by the quantity γ 4 = 1 16 , where U el is the elastic energy and ε is the shear strain.To determine the value of the shear susceptibility β, the effect of elastic loading on the ultrasound velocity was studied on two (Zr-and Pd-based) MGs [42,43].The shear susceptibility derived in these experiments is smaller than that originally estimated by Granato (β ≈ 40) but has nonetheless the same order of magnitude.Moreover, this quantity is close to the estimates derived by other methods [44] (see Section 3.5 and Table 1).

Relation between the Shear Modulus and Heat Effects
A simple analysis of Equations ( 1) and (2) shows that the IT implies an intrinsic relation between shear modulus changes and heat effects in MGs.Indeed, integration of Equation ( 1) taking account of Equation ( 2) leads to an expression for the enthalpy increment ∆H upon insertion of interstitial-type defects, which can be expressed through the shear moduli of glass and maternal crystal, ∆H ≈ (µ − G)/βρ, where β is the shear susceptibility and ρ is the density.Then, one can arrive at an expression relating the heat flow with the change of the shear modulus, where Ṫ is the rate of temperature change.This relationship was repeatedly tested and found to give an excellent description of the heat flow on the basis on shear modulus relaxation data not only upon structural relaxation below the glass transition temperature T g and in the supercooled liquid state but upon crystallization as well [45,46].The latter fact is really surprising and actually suggests that the whole excess internal energy ∆U (with respect to the crystalline maternal state) is mostly related to the interstitial-type defect system frozen-in upon glass production.An analysis gives a simple expression for this quantity [47], It is seen that the excess internal energy is simply controlled by the shear moduli of glass and maternal crystal.Upon crystallization, the defect system disappears and this energy is released as heat, i.e., ∆U ≈ Q cryst , where Q cryst is the heat of crystallization.A specially designed experiment confirmed this idea [48].From a physical viewpoint, this result means that the whole heat content of glass (i.e., the excess enthalpy with respect to the crystalline maternal state) is mostly determined by the interstitial-type defect system frozen-in upon glass production.
On the other hand, considering the initial and relaxed states of MGs and calculating the corresponding differences for the heat flow and shear modulus, ∆W = W rel − W and ∆G = G rel − G, Equation (4) after simplification can be rewritten as [49] This relationship shows that the quantities ∆W and ∂∆G/∂T should be proportional to each other.Figure 4 shows these quantities derived from calorimetric and shear modulus data for bulk glassy Zr 65 Al 10 Ni 10 Cu 15 .It is clearly seen that they are indeed proportional.This allows determination of proportionality constant and calculation of the shear susceptibility β.Table 1 shows β-values thus derived for a few MGs.It is seen that these values belong to a relatively narrow range 15 ≤ β ≤ 22 that, in general, agrees with the IT.Other methods used for the determination of β give close results [49].Since the shear susceptibility determines the shear softening effect (via Equation ( 2)), the heat effects (according to Equations ( 4) and ( 5)) and also related to the anharmonicity of the interatomic potential, it appears to be a major integral parameter of the glassy structure.6) derived from calorimetric and shear modulus measurements [49].The data correspond to structural relaxation below the glass transition.With permission from Elsevier, 2019.

Refinement of the Parameters of the Interstitialcy Theory
As reviewed earlier [15,16] and discussed in the present work, the IT provides a good description of different aspects of MGs relaxation behavior.At the same time, some model parameters of the IT were introduced in a phenomenological way.It is therefore desirable to clarify their physical meaning and relationship with material parameters.In the initial model, an interstitial defect was considered by Granato as an elastic string.At the same time, a split interstitial can be treated as an elastic dipole [47,50].The corresponding "dipole" approach is based on the expansion of the energy into a series in powers of the elastic strain created by the dipoles [50].It was found that that the Granato and "dipole" approaches give practically identical expressions for the elastic energy and shear modulus [47].A comparison of these approaches leads to an expression for the parameter α introduced in the original version of the IT (see Equations ( 1) and ( 2) , where ε ij is the elastic strain field created by the interstitial.Thus, the parameter α characterizes the "strength" of the defect.The shear susceptibility within the "dipole" approach was calculated as β = −γ 4 Ω t /µ, where µ is the shear modulus of the maternal crystal and Ω t = 1.38 is a parameter characterizing the elastic anisotropy of the interstitial [50].This estimate is about two times smaller than that given by Granato.It should be noted that the "dipole" approach was found very useful upon further IT development, especially in the part that accounts for the effect of the concentration of interstitials on their interaction (see Section 6).

Reconstruction of Temperature Dependence of the Shear Modulus Using Calorimetric Data
Equation (4) for the heat flow W(T) can be used in the "opposite" way, i.e., for the calculation of the shear modulus relaxation using input calorimetric data .This leads to a relation [51] where the subscript "rt" refers to the room temperature.

Heat Absorption Occurring upon Heating of Relaxed Glass
Within the IT framework, the heat absorbed upon heating from room temperature T rt to a temperature T sql in the supercooled liquid region can be calculated as [52] Q where G T rt , G T sql , µ T rt and µ T sql are the shear moduli of glass and maternal crystal at temperatures T rt and T sql , respectively, other quantities are the same as above.Since the shear moduli G T sql and µ T sql in the supercooled liquid state do not depend on the thermal prehistory, the only quantity in Equation ( 8), which varies upon structural relaxation, is the room-temperature shear modulus.The moduli G T sql , µ T rt and µ T sql are only temperature dependent.The temperature T sql can be accepted as a constant and the quantities G T sql , µ T rt and µ T sql are then constants as well.Thus, the heat Q is dependent on a single variable G T rt .The latter can be changed by preliminary heat treatment.Figure 6 shows the experimental data on the heat absorbed upon warming up of Pd 40 Ni 40 P 20 glass into the supercooled liquid as a function of the shear modulus measured at 330 K.The data points correspond to different preannealing treatments as indicated.It is seen that Q(G T 330 )-data nicely fall onto a straight line, as implied by Equation (6).The slope of this line within a few percent error agrees with its theoretical value given by this equation as ∂Q/∂G 330 = 1 βρ [52].Similarly, one can describe [53] the widely known so-called "sub-T g enthalpy relaxation" effect, which consists in the growth of the heat absorption near the glass transition temperature T g in MGs subjected to prolonged preannealing well below T g [54,55].The data discussed above in Sections 3.5, 5.1 and 5.2 convincingly demonstrate a close relationship between shear modulus relaxation and heat effects in MGs, in full agreement with the IT predictions.The solid line gives the lest square fit [52].With permission from Elsevier, 2019.

Density Changes upon Structural Relaxation and Crystallization
The interpretation of volume changes using the IT is based on the expected change of the volume ∆V upon creation of a dumbbell interstitial defect, which can be represented as ∆V/Ω = −1 + r i , where r i is the so-called relaxation volume reflecting the relaxation of the structure after defect creation, and Ω is the volume per atom [24,56].Then, if a defect concentration c is created, the volume increases by ∆V and the relative volume change becomes ∆V/V = (r i − 1)c.Using the shear modulus given by Equation (2), one arrives at the relative change of the density upon isothermal structural relaxation at a particular temperature as ∆ρ(t) where ∆ρ(t) = ρ(t) − ρ 0 , the densities ρ(t) and ρ 0 correspond to the shear moduli G(t) and G 0 = G(t = 0), respectively, α and β are the same as in Equation ( 2).An example of relative density changes as a function of shear modulus change of glassy Pd 30 Cu 30 Ni 10 P 20 measured at room temperature after isothermal annealing is given in Figure 7 [57].It is seen that this dependence can be fitted by a straight line, in accordance with Equation ( 9).This equation implies the slope of this line equal to (r i − 1)

αβ
. With α = 1 (as usually assumed), the shear susceptibility for this glass β = 19 (Table 1) and the relaxation volume r i = 1.6 (as for FCC metals [24]), one arrives at the slope equal to 0.031, in a good agreement with the experimental slope of 0.037 (Figure 7).Equation ( 9) can be then modified as [57] With β and r i given above, one arrives at dlnG/dlnV = 32, rather close to the experimental value dlnG/dlnV = 25 reported in Ref. [58] for the same glass.
For the density change upon crystallization, the IT gives [57] ∆ρ where ∆ρ = ρ cryst − ρ with ρ and ρ cryst being the densities of glass and maternal crystal, respectively.It is seen that, if the relaxation volume r i > 1, the density change ∆ρ/ρ > 0 upon crystallization.It was shown that Equation ( 11) then provides a reasonable explanation of crystallization-induced density changes of Zr-based MGs [57].However, for loosely packed crystalline structures, one can expect that the relaxation volume r i is less than unity.In this case, Equation ( 11) predicts a decrease of the density upon crystallization.Indeed, the literature gives a few examples of crystallization-induced density decrease of about 1% [59,60] or even more [61,62].For warming up from room to the temperature of the full crystallization, the above reasoning leads to the relationship ∆ρ(T) where ∆ρ(T) is the density change upon heating.It was shown that this equation provides a good description of ∆ρ/ρ-changes occurring upon heating up to the temperature of the full crystallization [63].It can be concluded that changes of the density are controlled by the shear moduli of glass and maternal crystal, which in turn reflect the evolution of the interstitial-type defect system.

Relation between the Enthalpies of Relaxation, Crystallization and Melting
According to the general IT approach discussed above (Section 2), the interstitial-type defects in glass are inherited from the melt.Provided that the melt quenching rate is big enough, one can expect that the defect concentrations in the glass and melt are nearly equal, c glass ≈ c melt .The quantity c glass determines the whole excess heat content (enthalpy) of the glass.Subsequent structural relaxation below T g leads to a decrease of this concentration by ∆c rel that results in a growth of the shear modulus from G up to G rel = G exp(αβ∆c rel ) and corresponding release of the enthalpy ∆H rel .Upon crystallization, the remaining defect concentration c melt − ∆c rel ≈ c cryst drops down to zero (the defects disappear), the shear modulus increases up to its value µ in the crystalline state with simultaneous release of the enthalpy ∆H cryst .The whole heat release of the initial glass after crystallization is then given by Equation ( 5), where ∆U ≈ ∆H cryst , as discussed above (Section 3.5).Since the defect concentration quenched-in from the melt is c melt ≈ ∆c rel + c cryst , the heat absorbed upon melting should be approximately equal to the total heat release upon structural relaxation and crystallization , i.e., in terms of the corresponding enthalpy changes, The result of a specially designed experiment [64] aimed at the verification of this relationship is reproduced in Figure 8, which shows that the data taken on eleven Zr-, Pd-and La-based MGs can be approximated by a straight line with the unity slope.Thus, the relationship in Equation ( 13) is indeed valid within the experimental error (about 10%) confirming the idea on a connection between the defects occurring upon melting of the maternal crystal and those disappearing upon structural relaxation and crystallization of the glass.With permission from Elsevier, 2019.

Relation of the Boson Heat Capacity Peak to the Defect Structure
A peculiar universal feature of atomic dynamics in non-crystalline materials consists in the presence of excess (over the Debye contribution) low-frequency vibrational modes [65,66].These low-frequency modes are detected as a peak in the low temperature (5-15 K) heat capacity C plotted as C/T 3 vs.T.These features are usually referred to as the boson peak, which is known for metallic glasses as well [67].The nature of the boson peak constitutes a matter of intensive ongoing debates [68][69][70][71].Granato argued that the boson peak originates from low-frequency resonance vibration modes of interstitial-type defects frozen-in upon glass production [72].He showed that the boson peak height should be proportional to the defect concentration.A refined equation for the boson peak height has the form [73] where T B is the boson peak temperature, C d is the heat capacity related to the interstitial-type defect system, Θ D and ω D are the Debye temperature and Debye frequency of the maternal crystal, respectively, ω r is the characteristic frequency of interstitial resonance vibrations, f is the number of resonance modes per interstitial-type defect, R is the universal gas constant and other quantities are specified above.The defect concentration c can be monitored by measurements of the shear moduli of glass and maternal crystal as implied by Equation ( 2).Thus, the boson peak within the framework of the IT is considered to be a "fingerprint" of the defect glass structure.
An experiment aimed to check the prediction given by Equation ( 14) was carried out on glassy Zr 65 Cu 15 Ni 10 Al 10 [73].The main result of this experiment is shown in Figure 9, which gives the measured height of the boson peak as a function of the defect concentration c calculated with Equation (2) using room-temperature measurements of the shear modulus after different annealing treatments.The annealing protocol was designed to perform measurements on both fully amorphous and partially crystalline samples.It is seen that independent of the state of the samples (amorphous/partially crystalline), the boson peak height linearly increases with the defect concentration, in line with Granato's prediction [72].The derivative dh B /dc calculated from Figure 9 provides reasonable estimates for the resonant vibration frequencies of the defects assumed to be responsible for the boson peak [73].On the other hand, since the interstitial-type defect structure determines the excess enthalpy ∆H of glass, as discussed above in Section 5.3, one can expect that the boson peak height h B should also be related to ∆H.The calculation gives the relation between the boson peak height and excess enthalpy of glass ∆H as [74] where the quantity Γ is defined by Equation ( 14), ρ is the density, and µ, α and β are same as in Equation ( 2).The relationship in Equation ( 15) directly connects the boson peak height with the excess enthalpy of the glass.The latter may be considered as an independent variable, which can be changed by the annealing leading to either structural relaxation within the glassy state or partial crystallization.Using differential scanning calorimetry, the excess enthalpy can be determined as dT, where W is the heat flow measured by the calorimeter, Ṫ is the heating rate, T i can be accepted equal to room temperature and T f is the temperature leading to the full crystallization.Figure 10 shows the dependence of the boson peak height h B calculated using Equation ( 15) together with the experimental h B -data as function of the variable ln [µ/(µ − ρβ∆H)] [74].It is seen, first, that the calculated and experimental h B -points are quite close.The experimental dependence h B (∆H) nicely falls onto a straight line and the slope of this dependence equals to (5.14 ± 0.29) × 10 −4 J/(mol × K 4 ).
This agrees with this slope given by Equation (15) as Γ αβ = 4.9 × 10 −4 J/(mol × K 4 ).Thus, the IT-based approach reproduces the boson peak height and provides a good description of its height on the excess enthalpy of the glass.An interesting study of the boson peak was recently reported by Brink et al. [75].They performed a molecular dynamic simulation of an equiatomic CuNiCoFe alloy in the crystalline and amorphous states alternatively with chemical disorder (high-entropy state), structural disorder and reduced density.They found that the density reduction and fluctuations of the elastic constants cannot be responsible for the boson peak.However, they revealed that the boson peak in the crystal increases with the concentration of dumbbell interstitials while other defects (e.g., dislocations) do not contribute to it.Interstitial atoms even at a small concentration lead to a boson peak, which is close to that in the glass of the same composition.At that, the vibrational modes of interstitial defects in the crystal resemble those of glass.Finally, the authors of [75] concluded that the softened regions provided by interstitials resemble the "soft spots" discussed in the literature on MGs [70,76] and the boson peak is due to quasi-localized defect-related modes.

Relation between the Properties of Glass and Maternal Crystal
In general, one can expect that the physical properties of MGs should be somehow related with the properties of their maternal crystalline states, which were used for the production of these glasses by melt quenching.Indeed, for instance, the properties of PdNiCuP glasses and their relaxation upon annealing should necessarily be related to the properties of intermetallic and/or metal-phosphide crystalline phases.To our knowledge, however, this issue was not raised in the literature.
Meanwhile, the IT is intrinsically based on the crystal→glass relationship.It starts from the main IT equation for the shear modulus of the glass G, which is scaled by the shear modulus of the maternal crystal µ (defectless state) as G = µ exp(−αβc).That is why the shear moduli of glass and maternal crystal determine major thermodynamic parameters of the glass, including its excess internal energy (≈enthalpy) as given by Equation ( 5).As a result, the shear moduli G and µ explicitly enter all relations for the heat effects (Equations ( 4), ( 6), ( 8) and ( 13)), volume relaxations (Equations ( 9)-( 12)) and low temperature excess heat capacity (boson peak) (Equations ( 14) and ( 15)).By that, in most cases, one must know not only G and µ at particular temperatures but also their exact temperature dependences, otherwise any quantitative agreement with the experiment cannot be achieved.Temperature dependences of the shear moduli of glass and maternal crystal reflect the interstitial-type defect structure of glass and the physical origin of crystal→glass relationship within the framework of the IT is intrinsically related to this defect structure, which controls the fundamental properties of the glass [77].In a certain sense (not structural), one can accept Granato's statement that ". . .glasses and dense liquids are crystals containing a few percent of interstitials" [78].

Development of the Interstitialcy Theory
Despite the amazing matching of the IT predictions with a number of relaxation phenomena in MGs, there exist a few phenomena, which cannot be easily interpreted within the framework of the original Granato's theory.First, an increase of the apparent defect concentration above T g due to the thermal activation, which follows from the original IT version (e.g., see Figure 10 in Ref. [16]), faces certain difficulties related to the high formation enthalpies necessary for this process.The same issue applies to the understanding of so-called rejuvenation of relaxation properties of MGs by quenching (or even relatively slow cooling) from the supercooled liquid region (i.e., from above T g ) [79,80].However, this difficulty can be avoided by assuming that an interstitial-type defect in the glass has a few energy states and the transitions between them are accessible by thermal activation.In the original IT version, this possibility was not assumed.One can suggest that the occurrence of a spectrum of energy states is related to the high defect density in the glass and melt, which is estimated to be a few percent [13,16].Such high density of the defects will result in their strong interaction and this must be taken into account.In Granato's original approach, the defect interaction was considered only qualitatively.
It is long known from the physics of crystals that this interaction can lead to the formation of interstitial clusters consisting of N individual interstitials, from N = 2 up to N = 7 and even more [23,28,81].Clustering is energy profitable since the formation enthalpy per interstitial decreases with the number of interstitials [23,81].The interstitial cluster consisting of N = 7 split interstitials in the FCC lattice represents a perfect icosahedron with < 0, 0, 12, 0 > Voronoi indexes [81], as illustrated by Figure 11 for the FCC cell.Thus, if melting of metallic crystals is related to an increase of the concentration of split interstitials up to a few percent, as considered by the IT, then one can expect the formation of clusters consisting of N = 2 to N = 7 interstitial-type defects.It is to be emphasized in this relation the commonly accepted notion that icosahedral clusters constitute a major structural feature of metallic melts and their concentration increases upon supercooling defining thus the dynamic slowdown of the internal movements and eventual glass formation [7,82].One can expect, therefore, that the solid glass will contain, first, individual interstitial-type defects (adjacent atoms with < 0, 2, 8, 0 > Voronoi indexes and/or close to them) and their small clusters (N = 2, 3), which correspond to the defect part of the structure.Larger clusters define the icosahedral-type structural backbone.Thus, the metallic glass constitutes a heterogeneous structure.The properties of the clusters should be evidently dependent on the number N. Specifically, the shear sensitivity B i should decrease with N defining a reduction of the diaelastic effect produced by the clusters.By that, the vibrational entropy of interstitial-type defects in the clusters should decrease due to their interaction leading to the mutual damping of the low-frequency vibration modes, as was qualitatively noted in the original Granato's model [14].The quasi-equilibrium balance between all these clusters will be dependent on temperature and thermal prehistory defining the evolution of glass properties [83,84].Some qualitative estimates of clustering kinetics are given elsewhere [85,86].Further detailed work in this direction is challenging.

Comparison with Other Models
Quite a few models are suggested in the literature for a description of defects in metallic glasses.A common shortage in most of them is the lack of understanding of their nature.However, it turns out that many of these models are quite consistent with the IT, as sketched below.
Egami suggested that the main properties of MGs are determined by the defects, which create shear, hydrostatic compression and/or tension (so-called τ-, pand n-defects) [8].Meanwhile, the dumbbell interstitials are featured by the same combination of stresses.Figure 12 shows the changes of the Voronoi polyhedra volume (indicated by the numbers) for a dumbbell interstitial in a Cu crystal with respect to the ideal lattice.It is seen that the interstitial produces both positive and negative changes of the Voronoi polyhedra volume.In other words, the defect creates both hydrostatically compressed and hydrostatically tensioned regions.One can expect that the same should be applicable for an interstitial-type defect in the glass.Taking into account that shear stress field is one of its main characteristics, one can conclude that the defect demonstrates the properties similar to those of Egami's τ-, pand n-defects.
-10.53% +9.10%It is often mentioned that supercooled liquids contain "strings" ("string-like" solitons), which become frozen in the solid glass [12,[87][88][89][90][91].It is also sometimes said that these defects resemble the signatures of dumbbell interstitials in crystals [33,88].Meanwhile, the "string" idea is compatible with the IT since a split interstitial in the original Granato's model was considered as a string segment [13,14].Similar arguments can be applied for the defects viewed as "shear transformation zones" [3,92], "soft spots" [70,76], "soft zones" [92], "flow units" [10,93], "liquid-like regions" [11], "geometricallly unfavored motifs" [76,94] and "regions with large non-affine displacements" [3,9,76].Since the region of an interstitial-type defect is characterized by the big shear susceptibility, the shear deformation of the surrounding material upon action of the applied stress is bigger that of the matrix and contains a large non-affine component [41].Naturally, the above terms ("soft zones", "liquid-like regions", etc.) can be applied to this region.It should be noted that these defects also contribute to the low-frequency part of the vibration spectrum of a glassy structure [76,93] similar to what is assumed by the IT.
Another approach to the understanding of defects in MGs is the "free volume" model, which was suggested long ago [1,95,96] and subjected to numerous modifications since then (e.g., Ref. [97]).In this approach, the defects are considered as regions of the reduced density (vacancy-like "free volume"), which affect relaxation and deformation phenomena in MGs.The premise for its popularity consists, on the one hand, in a decreased density of MGs (frozen-in "free volume") with respect to their maternal crystals and, on the other hand, in the existence of numerous correlations of MGs' properties with the amount of the "free volume" [98].Although this model was repeatedly criticized in different directions [5,99], it nonetheless constitutes perhaps the most popular approach.Meanwhile, it is quite compatible with the IT.Indeed, one should recall that the volume change ∆V v upon vacancy formation is ∆V v /Ω = r v + 1, where r v is the corresponding relaxation volume and Ω is the volume per atom [56].Taking into account the volume change ∆V i occurring upon interstitial formation (see Section 5.2 above) and accepting, e.g., for Al, r i = 1.9 and r v = −0.38 [24], one can easily arrive to the ratio of the relative volume changes produced by interstitials and vacancies, (∆V/V) i /(∆V/V) v = (r i − 1)/(r v + 1) = 1.45.This means that the volume changes for interstitials and vacancies have the same sign and are quite comparable in the magnitude [100].Thus, a decrease of the defect concentration within both the free volume model and the IT should lead to the densification of the glass by about the same amount.Moreover, since the free volume in both models is proportional to the number of the defects, one can naturally expect a correlation of material properties with the amount of the free volume.However, there are major differences between these approaches.In the former one, it is the free volume that constitutes the principal source for the property changes.The IT considers the features of the interstitial-type defect (large shear susceptibility, high formation entropy, and specific strain fields) to be of major importance while defect-related volume changes are of secondary relevance.It is also very important that vacancy-like free volume appears to have the spherical symmetry, which is why it should not interact with the external shear stress.Conversely, an interstitial defect is strongly asymmetric (Figure 1b) and constitutes an elastic dipole displaying strong sensitivity to the external shear load.Moreover, if a metallic glass is formed by melting of a loosely packed crystalline structure (related to the generation of interstitial-type defects), the creation of the free volume is not necessary at all.In this case, the density of glass can be even bigger than that of the maternal crystal (see Section 5.2 above) that cannot be understood within the free volume approach.

Concluding Remarks
The Interstitialcy Theory (IT) of condensed matter states thermodynamically predicts that melting of simple metallic crystals is related to a rapid increase of the concentration of the interstitial defects in their most stable dumbbell form just near the melting temperature.Recently, rather convincing (although indirect) experimental arguments confirming this hypothesis were obtained.Computer simulations showed that these defects remain identifiable structural objects in the liquid state.Rapid melt quenching freezes them in the solid glass.In the liquid and glassy states, these defects do not have any clear topological pattern as in crystals (two atoms trying to occupy the same lattice site) but nonetheless display all the properties characteristic of dumbbell interstitials in crystals, as confirmed by computer modeling of mono-atomic metallic systems.These properties include strong susceptibility to the applied shear stress, specific strain fields and low-frequency modes in the terahertz vibration spectrum.It is found that rather numerous relaxation phenomena can be understood by assuming that these interstitial-type defects indeed exist in real multicomponent metallic glasses (MGs).
The mathematical formalism of the IT is quite simple and based on the two relations linking the unrelaxed shear modulus, which constitutes the basic thermodynamical parameter of the IT (being the second derivative of the Gibbs free energy with respect to the shear strain), with the formation enthalpy of interstitial-type defects and their concentration.A thermoactivated change of the defect concentration (due to structural relaxation below T g , in the supercooled liquid state or upon crystallization) leads to an alteration of the formation enthalpy and results in numerous heat effects, which are intrinsically related to the relaxation of the shear modulus, as verified by specially designed experiments.On the other hand, changes of the defect concentration lead to certain volume changes, which can also be monitored by the shear modulus relaxation.Accepting reasonable values of the material parameter for the interstitial-type volume relaxation (i.e., the relaxation volume), one can explain the kinetics and final volume changes occurring upon structural relaxation and crystallization.
The IT leads to the conclusion that the excess internal energy and enthalpy of the glassy structure with respect to the maternal crystal is mostly related to the elastic energy of the interstitial-type defect structure frozen-in from the melt upon glass production.The full dissipation of this elastic energy constitutes the heat of crystallization.On the other hand, since the total amount of the defects is determined by the melting, there arises a relationship between the latent heat absorbed upon melting and heat release occurring upon structural relaxation and crystallization, as verified experimentally.
Within the framework of the IT, the low temperature boson heat capacity peak originates from low-frequency resonance vibration modes of interstitial-type defects and, thus, its height should be proportional to the defect concentration.Specially performed experiments showed that this is indeed the case provided that the concentration is derived from measurements of the shear modulus, as assumed by the IT.The defect concentration determines the excess enthalpy of the glass and it is the origin, which defines the observed dependence of the boson peak height on the experimentally measured excess enthalpy.
In general, the obtained results convincingly demonstrate an intrinsic relationship of the shear modulus relaxation with the heat and volume effects occurring upon structural relaxation and crystallization of MGs.At that, the properties of the glass are tightly related with those of the maternal crystal, in line with the basic assumptions of the IT.It has also been shown that other models describing defects and properties of MGs are largely compatible with the IT.The ways for the development of the IT and its relation to other models on defects in MGs are discussed.

Figure 2 .
Figure 2. Histogram illustrating the distribution of the ratio of dilatation U bulk and shear U shear components of the elastic energy for dumbbell interstitials in 63 polycrystalline elemental metals.The same data for interstitial-type defects in 189 metallic glasses are also shown[27].With permission from JETP Letters, 2019.

Figure 3 .
Figure 3. Estimates of interstitial and vacancy concentrations in crystalline aluminum (a) and indium (b) derived from the diaelastic effect measurements.The melting temperatures are indicated[35,36].With permission from Pleiades Publishing, LTD, 2019.

Figure 4 .
Figure 4. Temperature dependences of the quantities ∆W and ∂∆G/∂T entering Equation (6) derived from calorimetric and shear modulus measurements[49].The data correspond to structural relaxation below the glass transition.With permission from Elsevier, 2019.

Figure 5 .
Figure 5. Experimental and calculated using Equation (7) temperature dependences of the shear modulus G of glassy Zr 65 Cu 15 Al 10 Ni 10 in the initial and relaxed states.Temperature dependence of the shear modulus µ after full crystallization is also shown.Calorimetric T g is indicated by the arrow [51].With permission from Elsevier, 2019.

Figure 6 .
Figure6.Dependence of the integral heat Q absorbed upon heating from 330 K to 610 K (supercooled liquid region) as a function of the shear modulus G 330 measured at 330 K just after heating onset.The points correspond to different preannealing treatment applied for shear modulus measurements and DSC tests as indicated.The solid line gives the lest square fit[52].With permission from Elsevier, 2019.

Figure 7 .
Figure 7. Dependence of the relative density change on the quantity ln(G/G 0 ), where G 0 is the initial room-temperature shear modulus and G is the room-temperature shear modulus after annealing of bulk glassy Pd 40 Cu 30 Ni 10 P 20 at T = 533 K [57].With permission from Elsevier, 2019.

11 Figure 8 .
Figure 8.The melting enthalpy vs. the absolute value of sum of the enthalpies of structural relaxation and crystallization.The numbers correspond to different Zr-, Pd-and La-based MGs [64].With permission from Elsevier, 2019.

Figure 9 .
Figure 9.The height of the boson peak as a function of the defect concentration c calculated using Equation(2).The line gives the least square fit[73].With permission from John Wiley and Sons, 2019.

Figure 10 .
Figure 10.Experimental and calculated height of the boson peak h B as a function of the excess enthalpy ∆H plotted according to Equation (15).The numbers near the data points indicate the corresponding preannealing temperatures in Kelvins [74].With permission from John Wiley and 2019.

Figure 11 .
Figure 11.Formation of a perfect icosahedron by the creation of dumbbell interstitials on the opposite faces of the FCC cell: (Left) elementary FCC cell where the arrows show how two atoms are inserted instead of one atom; and (Right) perfect icosahedron with < 0, 0, 12, 0 > Voronoi indexes formed by six dumbbell interstitials on the faces of the cell and one interstitial in the octahedral position (i.e., in center of the cell, see Figure 1a) [81].

Figure 12 .
Figure 12.Dumbbell [001]-oriented interstitial (two red circles) and the relative changes of the volume of the Voronoi polyhedra ∆V/V (indicated by the numbers) as compared with the ideal lattice for the atoms (shaded circles) in a copper crystal.∆V/V-changes are shown along the dumbbell axis, perpendicular to it and in the [111] direction.It is seen that the defect creates both positive and negative changes of the Voronoi polyhedra volume (i.e., changes of the local density) [27].With permission from JETP Letters, 2019.

Table 1 .
[49]rmination of the shear susceptiblity β on the basis of calorimetric and shear modulus data taken on MGs in the initial and relaxed state using Equation (6)[49].With permission from Elsevier, 2019.