1. Introduction
The non-crystallinity of glasses necessarily leads to the heat release or heat absorption upon heat treatment. This applies in full to metallic glasses (MGs) [
1]. Exothermal effects in MGs take place upon structural relaxation below the glass transition temperature
as well as upon crystallization. Strong endothermal reaction is observed in the supercooled liquid state, i.e., at temperatures
, where
is the crystallization onset temperature [
1]. Besides that, any physical impact (thermal cycling, plastic deformation, irradiation, etc.) leads to certain additional heat effects in MGs [
2,
3,
4].
Current literature most often relates the heat effects in MGs to the changes of the free volume frozen-in upon glass production, as was originally suggested in the 1980s [
5]. Later, the enthalpy relaxation in MGs was linked to the amount of the free volume in the simplest linear form [
6] and subsequent experiments seemed to confirm this relationship [
7,
8]. As a matter of fact, heat effects in MGs are often considered as an indirect measure of the free volume since then [
8].
However, it was repeatedly mentioned that the free volume concept has definite theoretical shortcomings [
9] and the enthalpy release observed during relaxation cannot be simply interpreted as a measure of the free volume change [
10]. On the other hand, exothermal reaction occurring upon crystallization is currently considered solely from a general viewpoint relating it to the appearance of crystalline phases, without discussing any specific details that might lead to the release of heat. We are unaware of any physical model, which could explain all heat effects taking place in the glassy state as well as upon crystallization within a unified physical concept, with the only exception discussed below.
Heat effects can be naturally explained if the defect system of glass is considered within the framework of the Interstitialcy theory (IT) suggested by Granato [
11,
12]. The IT argues that melting of metals is associated with the rapid increase of the concentration of interstitial defects in their most stable dumbbell form, in line with recent experimental observations [
13,
14]. These “defects” remain identifiable structural units in the liquid state [
15] while rapid melt quenching freezes them in the solid glass. They retain all basic properties of dumbbell interstitials in crystals—high sensitivity to the applied shear stress, specific shear strain fields as well as characteristic low- and high-frequency vibration modes—although do not have any characteristic geometrical image like in crystals (two atoms trying to occupy the same lattice site) [
16]. These entities can be considered as interstitial-type “defects” (quotation marks are omitted hereafter). Heat effects can be then interpreted in terms of the changes of the defect concentration. This approach provides quantitative explanations for quite a few other relaxation phenomena in MGs [
17].
A description of the heat effects within the framework of the IT is based on an expression for the formation enthalpy of interstitial-type defects [
11,
12],
, where the dimensionless
is related to the defect strain field [
18],
is the volume per atom and
G is the unrelaxed shear modulus. The latter is related to the defect concentration
c as
, where
is the shear modulus of the maternal crystal (i.e., the one, which was melted and then used for glass production by melt quenching) measured at the same temperature and the dimensionless shear susceptibility
characterises the anharmonicity of the interatomic potential and by the definition equals the ratio of the fourth-order non-linear shear modulus to the second-order (i.e., “usual”) shear modulus [
17]. Typically,
depending on MGs’ chemical composition [
17].
Thus, any change of the defect concentration leads to heat release/absorption depending on the
c-change sign. For the glassy state, this mechanism was first suggested in Ref. [
19]. It was later found that the whole excess enthalpy (heat content) of glass with respect to the maternal crystal (i.e., the difference between the heat contents of the glassy and crystalline states) within the IT framework is related to the elastic energy of interstitial-type defect system and this energy is fully released as heat when the defect concentration drops down to zero as a result of crystallization [
20]. Thus, all exo- and endothermal heat effects occurring upon structural relaxation and crystallization can be considered as different sides of the same process—a change of interstitial-type defect concentration. In this case, the heat flow upon structural relaxation and crystallization of glass should be described by the same kinetic law, which relates the heat effects with the defect concentration. At that, the latter can be monitored by shear modulus measurements.
The two main equations of the IT given above lead to the difference in the heat flow of the glassy and crystalline states (per unit mass) conditioned by this relaxation mechanism [
18]:
where
is the enthalpy of glass,
is the heating rate,
is the density and
m is the sample’s mass. It is seen that the heat flow
is fully controlled by the shear moduli of glass and maternal crystal,
and
. All other quantities are constants while any fitting parameters are absent. Equation (
1) was successfully tested for the heat effects occurring in the glassy state [
18] (A similar heat flow law was first derived in Ref. [
19]. It differs from Equation (
1) by slightly another expression in the square brackets:
, where
and
are the shear moduli of glass and maternal crystal at room temperature, respectively. This heat flow law provides the results very close to those given by Equation (
1). The latter, however, is more general [
18]). However, as mentioned above, the same mechanism should be valid for the crystallization as well, as first argued in Ref. [
18] and recently confirmed experimentally on a particular metallic glass [
21]. In this work, we show that Equation (
1) is valid both for the relaxation in the supercooled liquid state and the whole crystallization kinetics of four Zr-based MGs displaying diverse physical properties.
2. Experimental
X-ray non-crystalline
,
,
and
(at.%, labelled as
,
,
and
hereafter) produced as
mm
bars were used for the investigation. The choice of these glasses was conditioned by the following reasons. First, all above Zr-based glasses (either in the initial or preannealed states) display a level of the mechanical damping in the supercooled liquid state (i.e., above
), which is small enough to ensure automatic measurements of the shear modulus (see below). As noted earlier, a high enough damping level above
results in the loss of automatic resonant frequency tracking [
21], which strongly limits the range of possible glass compositions. Second, these glasses have diverse and, in some sense, particular properties. Specifically, the
glass is very resistant to oxidation [
22] that strongly favours precise calorimetric measurements on relatively small samples and their subsequent comparison with calculation results. The glasses
and
display, respectively, big and small compression plasticity [
23] that could potentially affect their high temperature behavior. The
glass is being produced industrially for different applications and demonstrates enhanced corrosion resistance [
24]. Besides that, the choice of the above MGs is conditioned by the fact that they can be fully crystallized (without any additional phase transformations) below the maximal temperature (≈900 K) achievable upon standard calorimetric measurements. Thus, it is important to check the IT approach sketched above on MGs displaying diverse physical properties.
Heat effects were measured by a Hitachi DSC 7020 instrument in flowing
(99.999%). The mass of the samples was 50–70 mg. Every DSC run on a glassy sample was taken up to the temperature of the full crystallization. It was next followed by the 2nd run on the same sample and the difference between the two runs was then calculated. It is this difference, which is shown below in
Figure 2 and compared with
-calculations performed using Equation (
1).
The electromagnetic acoustic transformation (EMAT) method (see Ref. [
25] for the method’s details) was applied to measure the transverse resonant frequencies
f (500–700 kHz) of
mm
samples at temperatures of up to 810 K in a vacuum of ≈0.01 Pa. For this purpose, frequency scanning was automatically performed every 10–15 s and the resonant frequency was determined as a maximal signal response received by the pick-up coil upon scanning. The half-width of the resonant curve was used to calculate the mechanical quality factor (equal to the inverse damping) and the latter was used to estimate the precision of resonant frequency determination. The shear modulus was then calculated as
, where
and
are the vibration frequency and shear modulus at room temperature, respectively. This way of
G-calculation ignores possible density changes that can occur upon heating (usually less than 1%). The errors for the absolute
-values were accepted to be 1–2%. Then, the errors in the absolute
data are about the same while the error in the measurements of
-changes was estimated to be 5 ppm near room temperature and about 100 ppm near
. A heating/cooling rate of 3 K/min was accepted in all shear modulus and DSC measurements.
Three of the MGs under investigation (, and ) were tested in the initial state. This was impossible for the fourth glass () because of the large damping near , which results in the loss of EMAT automatic signal tracking. To avoid this effect, this MG was first annealed by heating into the supercooled liquid region and cooling back to room temperature. This MG is labelled as “relaxed” below. Both initial and relaxed MGs were first tested from room temperature up to temperatures of 780–810 K, which in all cases lead to the full crystallization. The second run for each glass was performed on the same sample in order to measure shear modulus in the crystalline state.
3. Results and Discussion
Figure 1a–d shows temperature dependences of the shear modulus
for all MGs under investigation. The room-temperature shear moduli
listed in this Figure are taken from Ref. [
26]. It is seen that
-patterns for all MGs are quite similar. Heating of glassy samples leads to a monotonous decrease of
G while the slope
decreases by several times near
, which constitutes a typical behavior upon high-frequency
G-measurements [
17]. Heating by 50–70 K above
results in the beginning of crystallization, which leads to an increase of the shear modulus by 22% (
) to 40% (
). In the crystalline state (2nd run), the shear modulus
demonstrates a featureless decrease with temperature.
Figure 2 shows DSC traces of the MGs under investigation at the same heating rate of 3 K/min. Initial glasses (
,
and
) demonstrate (
i) exothermal reaction below
(not seen in the scale of
Figure 2a–c), (
ii) heat absorption in the supercooled liquid state (i.e., at temperatures
), which is bigger compared with that of the glassy and crystalline phases, and (
iii) large heat release due to the crystallization. In the relaxed MG (
), the feature (
i) is absent, as one would expect.
Figure 2 also gives the heat flow
-curves calculated with Equation (
1) at
K/min using the corresponding temperature dependences of the shear moduli
and
given in
Figure 1. The way of the determination of the shear susceptibility
from the experimental data and its values together with the densities for the MGS under investigation are given in Ref [
26]. The insets in
Figure 2 show experimental and calculated
-curves in the supercooled liquid and crystallization temperature regions on enlarged scales.
Figure 2 in general demonstrates that the calculation provides a rather good quantitative reproduction of the experimental
-curves. First of all, this applies to the heat absorption for all glasses in the supercooled liquid state. Temperature position of the crystallization exothermal peak is reproduced within ≈10 K for
and ≤6 K for other glasses. The height of this peak is reproduced within 10%–15% accuracy with the exception of
for which the calculation gives about
of the experimental height (not shown in
Figure 2). It should be emphasized that the density of data points provided by our EMAT system in the range of fast crystallization (i.e., upon a rapid change of the shear modulus) is significantly smaller than that for temperatures
(see
Figure 1) and this constitutes a significant source for the calculation errors when using Equation (
1), which contains both derivatives of the shear moduli and their difference. Nonetheless, one can conclude that this equation provides a good description of the heat effects on the basis of shear modulus data. It is interesting to note that the calculation reproduces the fine details of the crystallization kinetics as exemplified by the data around
K shown in the inset of
Figure 2d.
The obtained results confirm the basic idea of the IT sketched above: the origin of the heat absorption and heat release occurring in the glass transition range as well as upon crystallization can be understood as a result of the change of the concentration of interstitial-type defects frozen-in from the melt upon glass production. In particular, structural relaxation below
provides relatively small decrease of the defect concentration and constitutes the reason of moderate exothermal heat effect (feature (
i) mentioned above). The defect concentration above
increases with temperature leading to the heat absorption (feature (
ii)). Finally, fast crystallization at
results in a rapid drop of the defect concentration down to zero providing quick dissipation of the defect elastic energy into heat, which is released as a strong crystallization peak (feature (
iii)) [
20].
As mentioned above, the interstitial defects in crystalline structures are unambiguously geometrically defined as two atoms trying to occupy the same lattice site. However, the situation in metallic non-crystalline structures is more complex. Molecular static simulation performed on a monoatomic glassy structure revealed localized nano-regions, which display the properties similar to those of dumbbell interstitials in simple metals [
16]. These findings imply that interstitial-type defects can indeed exist in a monatomic non-crystalline structure. The issue whether such structural entities can be found in polyatomic glassy structures constitutes a major challenge for further research in this direction.