Merging of relaxations and step-like increase of accompanying supercooled liquid region in metallic glasses via ultrafast nanocalorimetry

Glassy materials under external stimuli usually display multiple and complex relaxations. The relaxations and the evolution paths of glassy materials significantly affect their properties and are closely related to many key issues in glass physics, such as glass transition and thermoplastic forming. However, until now, the relaxation dynamics in the presence of external stimuli and the microscopic atomic motion of glassy materials have been unclear due to the lack of structural information. By combining Flash and conventional differential scanning calorimetry (DSC), we applied a very large range of heating rate of six orders of magnitude and investigated the relaxation dynamics of three typical metallic glasses. We discovered the merging of distinct relaxation events with increasing rate of heating. Most interestingly, the experiments revealed new behaviors with step-like increases in the supercooled liquid region and excess heat capacity during the merging of multiple relaxations. A comprehensive scheme was proposed for the evolution of the thermal relaxation spectrum, the heterogeneity of the corresponding atomic motion and the potential energy landscape with rate of heating. These experimental results shed light on the mechanism of atomic rearrangement during heating and provided a new approach to regulate the physical properties of amorphous materials by controlling their intrinsic relaxation dynamics. Glasses are non equilibrium materials that develop by cooling of a supercooled liquid, where the rapidly increasing viscosity results in a kinetic arrest of long range atom rearrangements. During heating from glass state, different relaxations are thermally activated and display different relaxation spectrum. Here, we used one ultrafast nanocalorimetry to examine the evolution of multiple relaxations and discover the merging of the relaxation modes with increasing heating rates, resulting in step-like increases in both the supercooled liquid region and excess heat capacity. Our findings provide new insights on the evolution of the relaxation spectrum and the associated heterogeneous atomic motion.


Introduction
The relaxation phenomena of metallic glasses (MGs) have attracted a great deal of attention from researchers who seek to unveil the nature of the glass transition and other key issues in glass physics, such as physical aging, the memory effect and thermoplastic forming [1][2][3][4][5] . Compared to their crystalline counterparts, MGs display diverse relaxation modes over wide ranges of temperature and timescale, such as the primary α relaxation, the secondary β (Johari-Goldstein) relaxation and the faster relaxation modes called β' or γ relaxations (shown in the upper part of Fig. 1a) [6][7][8][9][10][11] . Different relaxation modes gradually decouple during vitrification. At a sufficiently high temperature, there is only the α relaxation; when the liquid quickly quenches into the supercooled liquid region, the relaxation mode splits into α and β relaxations. In the glass state below the glass transition temperature, the α relaxation disappears, and the relaxation mode further decouples into β and β' relaxations 5,12 . However, few studies have focused on the evolution of the relaxation mode during the glass-to-liquid transition (GLT) by fast heating to bypass crystallization [13][14][15][16] . From the view of practical applications, the GLT process is one of the critical issues in the net-shape thermoplastic forming processes of MGs, such as injection molding, blow molding, and microreplication 14,17,18 . Meanwhile, researchers have recently reported that the macroscopic GLT is induced by the activation and percolation of numerous flow units that are closely correlated with local β and β' relaxations 9,15 . Thus, understanding how atomic rearrangements govern these relaxation processes and the detailed evolution paths of different relaxations is of great importance in fundamental research and practical applications for MGs.
For thermal activation processes in MGs, such as crystallization and different relaxations, including the primary α relaxation, the relationship between the typical temperature (crystallization temperature or glass transition temperature) and the heating rate roughly follows the Kissinger equation or the Arrhenius equation 9,10,[19][20][21][22] . Thus, different thermal activation processes can be separated under high heating rates. For example, the hidden glass transition (α relaxation) in marginal Al-based MGs has been separated from primary crystallization, and the weak β relaxation in fragile MGs is separated from the α relaxation when the heating rate reaches a critical value 19,20 . Furthermore, the activation processes of different relaxations during the GLT with increasing temperature in MGs are usually accompanied by obvious thermal signals, such as endothermic peaks before the glass transition, which usually correspond to the appearance of β and β' relaxations and the unfreezing of the corresponding local regions 9,15,16,20 . Considering that the activation energies of β and β' relaxations are much smaller than that of the α relaxation, a relaxation mode coupling should appear for different relaxations with increasing heating rates. Moreover, recent studies have established that the secondary β relaxations in MGs are governed by string-like clusters of particles with cooperative motion 23,24 and that fast β' or γ relaxations result from caging-breaking motion driven by internal stress in nonequilibrium states (shown in the lower part of Fig. 1a) 9,10,25,26 . Therefore, the above possible relaxation mode coupling with increasing rate of heating may imply the coupling of various atomic rearrangement motions, which provides critical information for understanding hidden physical process during the GLT. However, due to the limited experimental capacity of calorimetry instruments, until now, there has been no report about this kind of mode coupling for different relaxation events with heating. Recently, one advanced commercial chip-based instrument for fast differential scanning calorimetry (DSC) (Mettler Toledo, Flash DSC) enabled the use of heating rates for thermoanalytical measurements over a range of six orders of magnitude, which provides a good opportunity to systematically study relaxation events during heating and related relaxation dynamics in amorphous materials [27][28][29][30] .
In this paper, we systematically investigate the relaxation dynamics of three typical MG systems with different relaxation patterns in a large range of heating rates, six orders of magnitude, by combining advanced Flash and conventional DSC instruments. First, we find the step-like emergence of three distinct relaxation events as the temperature and heating rate increase in a LaCe-based MG system. Moreover, we observe and discover a large increase in the supercooled liquid region during the merging of multiple relaxations. Additionally, the corresponding excess heat capacity exhibits a similar step-like increase with the merging of multiple relaxations. In addition, similar relaxation merging behaviors with heating rate take place in the other two MG systems. These results reveal the novel activation path of different relaxation events during the GLT with the heating rate and provide a new approach to tune the relaxation dynamics by ultrafast heating to increase the supercooled liquid region.

Materials
Ingots of La 55 Ce 15 Ni 20 Al 10 , La 60 Ni 25 Al 15 and Al 90 Ca 10 were prepared by arc melting the elemental components several times in a Ti-gettered argon atmosphere to ensure homogeneity. Then, ribbon samples were prepared by single-roller melt spinning on a copper wheel with a tangential speed of 55 m/s. The ribbons had a cross section of approximately 2 mm × 20 μm and were approximately several meters in length. The cooling rates for preparing the ribbon-like samples in this study were estimated as approximately 2.5 × 10 6 K/s. The glassy nature for all ribbon-like samples was ascertained by X-ray diffraction (D8 Discover Diffraction with Cu K α radiation, Bruker, United States) and DSC (Diamond DSC, Perki-nElmer Inc., United States).

Calorimetry measurements
The samples for all the DSC tests, including conventional DSC and Flash DSC tests, were cut from the above as-cast ribbons in sizes that were suitable for different DSC instruments. A high-rate differential scanning calorimeter with chip sensors (Flash DSC 2, with maximum heating and cooling rates of 4 × 10 4 and 1 × 10 4 K/s, respectively, Mettler Toledo, Switzerland) was used to investigate the thermal behaviors of various MG systems under different heating rates from 5 K/s to 10,000 K/s. The measurement temperature range for the Flash DSC 2 was between −90°C and 1000°C. The as-cast ribbon-like samples were cut into tiny rectangular pieces of approximately 150 μm × 150 μm × 20 μm (length × width × thickness) and then loaded onto the Flash DSC chip. As a reference, the heat flow curve at a low heating rate of 0.33 K/s was also measured by Diamond DSC in this study. For conventional DSC tests, the as-cast ribbon-like samples were cut into small pieces by scissors; each piece was approximately 5 mm × 2 mm × 20 μm (length × width × thickness), which was suitable to be loaded into an aluminum pan for conventional DSC measurements. The masses of the samples in conventional DSC measurements were weighed. For the mass of the tiny Flash DSC samples, the mass was estimated from the corresponding volume (measured by optical microscope) and density (estimated by the Archimedean technique).
Heat capacity measurement based on the Flash DSC platform and confirmation of excess heat capacity at the glass transition temperature Based on the Flash DSC platform, the heat capacity c p was measured 31 . The detailed temperature-time program for the step response analysis included a temperature jump of 2 K (far less than the thermal fluctuations) at one heating rate φ h and isotherms of duration Δt. The value of Δt was set to guarantee a frequency of 10 Hz. Thus, for different applied heating rates, the value of Δt was adjusted. The complex heat capacity, c p *(ω) = c p ' (ω) − ic p "(ω), was calculated by a Fourier transform of the heat flow rate HF(t) and the instantaneous heating rate φ h (t), , where c p * is the complex heat capacity, c p ' and c p " are the real part and imaginary part of the complex heat capacity, ω is the frequency, HF is the heat flow and φ h is the heating rate. From the obtained complex heat capacity, the reversible heat capacity c p was calculated byc p ðwÞ ¼ ðc 0 p ðwÞ 2 þ c 00 p ðwÞ 2 Þ 0:5 . For each applied heating rate φ h , the temperature of the onset of the glass transition was obtained from the corresponding heat flow curve. Then, the excess heat capacity Δc p at the glass transition temperature was calculated by subtracting the linear fit of the heat capacities of the supercooled liquid state and glass state.

Observation of multiple relaxation events via the combined DSC platform in LaCe-based MG
Two typical heat flow curves by conventional DSC with a heating rate of 0.33 K/s and by Flash DSC with a heating rate of 100 K/s for the LaCe-based MG are exhibited in Fig. 1b, c, respectively. The conventional heat flow curve displays only the glass transition signal corresponding to the primary α relaxation. In contrast, the Flash DSC heat flow curve exhibits two distinct endothermic peaks before the glass transition. For LaCe-based MG, the primary α relaxation, the secondary β relaxation and the fast β' relaxation have been reported based on results of dynamic mechanical analysis (DMA) 9 . However, the conventional calorimetry method with limited heating rates of several K/s cannot observe the activation of the secondary β relaxation and the fast β' relaxation except when the thermal signal of secondary β relaxation is induced by a thermal annealing treatment 32,33 . By comparison, in this study, based on ultrafast Flash DSC, the secondary β relaxation and the fast β' relaxation are observed for the first time by the calorimetry method, as shown in Fig. 1c. Previous studies have reported that the secondary β relaxation and fast β' relaxation commonly observed in rare earth-based MGs are separately correlated to stringlike motion and more local cage-breaking motion, which can be attributed to the presence of highly 'mobile' atomic pairs, such as the rare earth element-Ni pair 9,10,16 . Therefore, from this perspective, ultrafast nanocalorimetry provides a new and powerful tool to investigate dynamic behaviors in glasses, such as the GLT and various relaxation behaviors during heating, and the accompanying local atomic rearrangement motions during these thermally activated processes in glassy materials.
To investigate the evolution of various relaxation events with heating rates in the LaCe-based MG, a series of heating rates ranging over six orders of magnitude were applied to measure the corresponding heat flow curves. All heat flow curves were normalized by sample mass and are shown in Fig. 2a, b. From Fig. 2a, b, when the heating rate φ h increased to approximately 10 K/s and 50 K/s, the thermal signals of the endothermic peaks for the secondary β relaxation and fast β' relaxation appeared, respectively. Considering that the thermal signals for the secondary β relaxation and fast β' relaxation correspond to the activation of local string-like and cage-breaking motions, as shown in Fig. 1a, the thermal signal intensity was relatively weak, and it was difficult to detect these signals at low heating rates by conventional DSC instruments 20 . In contrast, when the heating rates were larger than the corresponding critical heating rates for the appearance of the β and β' relaxations, the thermal signals corresponding to the activation of the β and β' relaxations were strong enough to be observed by Flash DSC; they appear as two obvious endothermic peaks before the glass transition in Fig. 1c. From the perspective of recovery enthalpy during heating, when the heating rate is very small, the effect of the decrease in the recovery enthalpy induced by structural relaxation during heating is very large, and the thermal signals from the β and β' relaxations are relatively weak. When the applied heating rates become large enough, the effect of the decrease in the recovery enthalpy can be minimized, and thus, the thermal signals from the β and β' relaxations become strong enough to be detected 16,20 . Moreover, it is apparent that the thermal signals for the primary α relaxation (glass transition), the secondary β relaxation and the fast β' relaxation evolve with increasing rate of heating. In particular, it is intriguing to observe in Fig. 2b that the secondary β relaxation started to merge with the primary α relaxation when the heating rate increased into one critical value φ β→α of approximately 2400 K/s. The merging of the initial α relaxation (marked as α I ) and the secondary β relaxation led to the second stage of the α relaxation (marked as α II ). Furthermore, when the heating rate continued to increase to another critical value φ β'→α of approximately 15,000 K/s, the fast β' relaxation merged with the α II relaxation and entered the third stage of α relaxation (marked as α III ). Therefore, the above results directly show that the thermal relaxation spectra for MGs at different heating rate ranges are not the same and that there appear to be different activated relaxation modes for one MG system within different heating rate ranges.

Activation energies for various relaxation events
To further illustrate the evolution paths of the different relaxation events for LaCe-based MG, the detailed evolutions of the onset temperature T onset and the end temperature T end for the different relaxation events with heating rates were obtained based on Fig. 2a, b and plotted in Fig. 2c. The α relaxation process within the whole experimental heating rate range clearly displayed three stages: stage I (heating rate below 2400 K/s), stage II (heating rate between 2400 K/s and 15,000 K/s) and stage III (heating rate above 15,000 K/s). When the heating rate increased to 2400 K/s, the end temperature of the secondary β relaxation reached the onset temperature of the α relaxation, which implied the merging of the secondary β relaxation and α relaxation (shown in the upper-right plot of A in Fig. 2c). When the heating rate increased to 15,000 K/s, the end temperature of the fast β' relaxation reached the onset temperature of the α relaxation, which implied the merging of the fast β' relaxation and the α relaxation (shown in the lower-right plot of B in Fig. 2c).
Notably, the shift of all three relaxations to higher temperatures with increasing heating rate in Fig. 2a and b directly implied that these relaxations were thermally activated processes. Although the representation of the temperature and heating rate dependence has been presented differently in the literatures 9,10,20-22 , such as the Kissinger equation or Arrhenius equation, the values of the effective activation energy for the relationships obtained by fitting the Kissinger equation and Arrhenius equation within a small temperature range were very close 20,34 . Thus, in this study, the effective activation energy for relaxations was obtained by simply fitting the onset values of the relaxation events in Fig. 2c by the Kissinger equation, as shown in Fig. 2d. The values of the effective activation energies for the α, β and β' relaxations were 129.3 ± 2.6, 75.6 ± 1.5 and 41.6 ± 1.2 kJ/mol (1.34 ± 0.03, 0.78 ± 0.02 and 0.43 ± 0.01 eV), respectively, which were consistent with the previous results obtained by the DMA method for LaCe-based MG 9 . Moreover, two step-like decreases in the onset temperature for the α relaxation were observed, 10 K and 19 K, when the β and α relaxations merged and the β' and α relaxations merged (shown in the right plots of Fig. 2c). Meanwhile, the corresponding activation energies for the α relaxation showed two step-like decreases, from 129.3 ± 2.6 kJ/mol (1.34 ± 0.03 eV) to 106.8 ± 2 kJ/mol (1.11 ± 0.02 eV) and from 106.8 ± 2 kJ/mol (1.11 ± 0.02 eV) to 84.2 ± 3.2 kJ/mol (0.87 ± 0.03 eV). Thus, for the three stages of the α relaxation, the corresponding activation energies differed. These results indicated that the α relaxation was more easily activated at a large heating rate than a slow heating rate, considering the decrease in the α relaxation activation energy with heating rate.

Supercooled liquid region evolution with heating rates
In general, for glassy materials, the glass transition behavior during heating (the GLT process) can be characterized by a single endothermic reaction, where the specific heat, c p , increases abruptly to a maximum value and then remains constant or slightly decreases to the crystallization onset temperature, T x . The temperature range between the glass transition and crystallization is called the supercooled liquid region (ΔT = T x − T g ), which is closely related to glass formation ability and thermoplastic deformation ability. Moreover, considering the different evolution paths for the glass transition and crystallization with temperature during heating, the supercooled liquid region should continuously change with heating rate. Here, from Fig. 2c, the extent of the supercooled liquid region ΔT for LaCe-based MG is plotted against the heating rate in Fig. 3a. The supercooled liquid region displays three stages: I, II and III, which are consistent with the three stages of the α relaxation in Fig. 2c. When the heating rate was below 2400 K/s, ΔT slowly increased with increasing heating rate; and when the heating rate increased to 2400 K/s, ΔT exhibited an obvious step-like increase from 20 K to 70 K. Similarly, when the heating rate was between 2400 K/s and 15,000 K/s, ΔT slowly increased with the heating rate; when the heating rate increased to 15,000 K/s, ΔT developed a second step-like increase from 70 K to 110 K. The inserted graphs in Fig. 3a give the thermal relaxation spectrum corresponding to heating rates of 800 and 3000 K/s. Most interestingly, these two step-like increases in ΔT were accompanied by the merging of the secondary β relaxation and α relaxation and the merging of the fast β' relaxation and α relaxation, respectively. For MGs, the larger the supercooled liquid region is, the greater the resistance to crystallization, which is one of the prerequisites for thermoplastic forming applications 17 .

Excess heat capacity evolution with heating rates
From the perspective of statistical physics, the heat capacity of matter in an equilibrium state reflects the possible motions of the basic atoms in the system 35 .
Compared to the absolute heat capacity for one state, the excess heat capacity Δc p of supercooled liquids relative to their glassy state contains contributions from the influence of changes in configuration states on the thermodynamic contributions of atomic rearrangement processes, which also give rise to secondary relaxations below the glass transition [36][37][38][39] . Thus, the excess heat capacity between the glass and liquid can reveal the kinds of atomic rearrangement motions that are involved in dynamic processes. Here, the heat capacities corresponding to different heating rates for LaCe-based MG were measured (discussed in detail in the Methods section), and one heat capacity curve for 20,000 K/s is shown in the inset of Fig. 3b. The excess heat capacity between the glassy state and the supercooled liquid state with different heating rates is plotted in Fig. 3b. Similar to the step-like increase in the supercooled liquid region, the excess heat capacity displayed two step-like increases with the merging of multiple relaxations. For stage I, before the β relaxation merged into the α relaxation, the value of Δc p was almost constant, 12.4 ± 0.3 Jmol −1 K −1 . For MGs with only the thermal α relaxation, previous studies reported that the values of excess heat capacity were almost constant, 13.7 ± 2.1 Jmol −1 K −1 38,39 , which was consistent with the above results of stage I in this study. After the merging of the β, β' relaxations and α relaxation in stages II and III, the values of Δc p jumped to 15 ± 0.1 Jmol −1 K −1 and 16.5 ± 0.1 Jmol −1 K −1 , respectively. This behavior can be used to understand the possible atomic rearrangement motions during the GLT. When the heating rate was in stage I, only the atomic rearrangement motions related to the α relaxation contributed to the GLT. When the heating rate increased into stage II, both rearrangement motions related to the α and β relaxations contributed to the GLT. When the heating rates further increased into stage III, all the rearrangement motions related to the α, β and β' relaxations participated in the GLT.

Two additional typical MG systems
For MGs, the relaxation patterns are closely related to their chemical compositions and structural heterogeneities. Various MG samples with different chemical compositions and microscopic structures display different relaxation behaviors [7][8][9][10][11]21,32,33 . Specifically, for the LaCe-based MG in this study, three different relaxation events developed during heating 9 . Thus, to verify whether the enlargements of the supercooled liquid region and the excess heat capacity with multiple relaxation merging via ultrafast  Fig. 4a but not for Al 90 Ca 10 MG (Fig. 5a). In contrast, the thermal signal corresponding to the α relaxation for Al 90 Ca 10 MG can be observed by Flash DSC (Fig. 5b, c). For marginal MGs, such as Al-based MGs, the weak glass transition signal is prone to overlap due to  strong primary crystallization 40 . Thus, conventional DSC with a heating rate of only several K/s is not a suitable calorimetric measurement method to detect the glass transition of Al-based MGs. However, considering that the glass transition and crystallization are significantly different kinetic processes, they follow different evolution paths with heating rates. For higher heating rates, the thermal signals corresponding to the glass transition and crystallization evolve along two different paths, and the glass transition and crystallization separate when the heating rate increases to one critical heating rate 19 . Thus, the glass transition for marginal Al 90 Ca 10 MG can be detected by Flash DSC with a large heating rate. Moreover, compared to the LaCe-based MG with three relaxation events in Fig. 1c, only two relaxation events, the secondary β relaxation and the α relaxation, appear in Figs. 4b, c and 5b, c. For Al-based and La-based MGs, with increasing rates of heating, the evolution rate of the β relaxation is much faster than that of the α relaxation, and finally, the β relaxation and the α relaxation merge when the heating rate increases to a critical value. In addition, based on the heat flow curves in Fig. 4b, c and Fig. 5b, c, the corresponding activation energies of the secondary β relaxation and the α relaxation for the Al-and La-based MGs were calculated by fitting the Kissinger equation; the detailed values are shown in Fig. 6a and d. Similar to the LaCe-based MG, the activation energy of the α relaxation decreases after the merging of the β relaxation and the α relaxation. Moreover, the supercooled liquid regions for the Al-and La-based MGs underwent a large increase from 34 K to 73 K and from 104 K and 221 K, respectively. These results indicate that the supercooled liquid region can be enlarged by multiple relaxation merging via ultrafast heating.
In addition, the excess heat capacity exhibited a similar step-like increase with the merging of the secondary β relaxation and α relaxation for the La-and Al-based MGs in Fig. 6c, f. The excess heat capacity Δc p increased from 12.6 ± 0.2 Jmol -1 K −1 to 15.2 ± 0.1 Jmol −1 K −1 and from 12.7 ± 0.2 Jmol −1 K −1 to 15.1 ± 0.1 Jmol −1 K −1 for the Labased and Al-based MGs, respectively. Compared to the LaCe-based MG in Fig. 3b, it is interesting to note that the increases in Δc p induced by the merging of the β relaxation and the α relaxation were very similar for all three MG systems, approximately 2.5 ± 0.1 Jmol −1 K −1 . This result indicated that the similar value of the increase in Δc p resulted from the contribution of the atomic rearrangement motions corresponding to the β relaxation. To illustrate the contribution of the fast relaxation to the increase in the excess heat capacity during the α relaxation, a series of schemes for the heat capacity corresponding to different heating rate ranges are displayed in Fig. 3c. These schemes clearly show a common physical mechanism for the merging of the β relaxation and the α relaxation during heating in these MG systems, which implies that different atomic rearrangement motions became more uniform with the merging of different relaxations.

Discussion
One of the major questions in glass physics is whether vitrification from a viscous liquid into a nonequilibrium glass and devitrification from a glass to a liquid by heating are related exclusively to the primary structural relaxation process, the α relaxation process, which is attributed to the cooperative motion of several structural units, or whether other atomic rearrangement motions play a role. The conventional view, based on experimental evidence for various glassy materials, indicates that the cooling rate dependence of the glass transition temperature exhibits the same behavior as the temperature dependence of viscosity or the structural α relaxation time [1][2][3] . However, recent experiments on polymer glasses under geometrical confinement and MGs directly indicate that the vitrification kinetics and the α relaxation are decoupled, and the role of other atomic motion mechanisms of fast relaxation play a key role 31 . In this work, we used Flash DSC with a heating rate range of six orders of magnitude and discovered the enlargement of the step-like supercooled liquid region enlargement and the excess heat capacity induced by multiple relaxations merging in three MGs. In particular, the increase in the excess heat capacity between the supercooled liquid and the glass during the GLT comes from the contribution of the atomic rearrangement motions corresponding to fast relaxations with the merging of the relaxation events. Therefore, these results indicate that the atomic rearrangement motions related to fast relaxation modes in addition to α relaxation may be engaged in the GLT during heating at a fast heating rate.
Moreover, for glassy materials, the endothermic signals before the glass transition are actually the calorimetric features of the unfreezing of atomic rearrangement in local regions 16 . Here, in this work, the observed endothermic thermal signals for different relaxation events in Fig. 2a and b imply that the atomic rearrangement motions corresponding to the activation modes of multiple local structural heterogeneities in MGs are heterogeneous 9,10,16,25,26 . Therefore, from the viewpoint of atomic rearrangement, the merging of multiple relaxations implies the merging of different modes of atomic motion, and the atomic rearrangement during the GLT evolves from spatially heterogeneous modes to more homogeneous modes with increasing rates of heating. In addition, for the merging of relaxations with the heating rate, it is interesting to see whether these relaxations interact. If the relaxations undergo noninteraction mixing, then these relaxations may separate when the heating rate continues to increase to a critical value. Unfortunately, according to the evolution of the onset and end points of different relaxations in Figs. 2d, 4d and 5d, there are no new merging or separation behaviors after the merging of the β', β and α relaxations for these three MG systems within the experimental heating rate limit of the Flash DSC instrument. New calorimetric technology with a faster heating rate is needed to investigate this question. However, considering that the transition of atomic rearrangement motion with the merging of different relaxations is usually irreversible, the merging of the relaxations may be due to interaction mixing. This point will be the next focus of our research.
However, from the perspective of the potential energy landscape 5,15 , when the volume of the glass sample is constant, the potential energy landscape should be fixed. The energy landscape contains multiple metabasins separated by high barriers. For the relaxation events in glasses, the β relaxation is considered to involve stochastic and reversibly activated hopping events across "subbasins" confined within the inherent "metabasin" (intrabasin hopping), and the α relaxation is considered to involve irreversible hopping events extending across different metabasins in the landscape (as shown in Fig. 1a) 5 . Considering the fractal nature of the potential energy landscape in glasses, the hopping barriers corresponding to the α relaxation should be multiple rather than single, as indicated for a trimodal scheme at the top of Fig. 7. However, it is still unclear how a glass sample selects a corresponding energy landscape under different external stimuli. In this study, based on the above experimental results for the activation energy and the excess heat capacity for the α relaxation within different heating rate ranges, three different potential energy landscapes corresponding to different heating rate ranges are displayed, as shown in Fig. 7. In the low heating rate range (stage I), the potential energy landscape displays a rough pattern with a three-level metabasin pattern as shown in Fig. 1a. In the medium heating rate range (stage II), the potential energy landscape displays a less rough pattern with a two-level metabasin pattern including the β' relaxation and the α relaxation. In the high heating rate range (stage III), the potential energy landscape displays a smooth pattern with only a onelevel pattern including the α relaxation. Thus, the above evolution of the relaxation pattern with heating rate reflects the change in sampling of the potential energy landscape for MGs.
The above results for new relaxation behavior by ultrafast heating provide a new possibility to tune the relaxation pattern and enlarge the supercooled liquid region by selectively designing applied heating rates. The activated relaxation pattern and the corresponding activated energy landscape are selected according to the range of heating rate. To clearly display the evolution of the thermal relaxation spectrum, the atomic mobility heterogeneity and the activated energy landscape within different heating rate ranges, a simple scheme was proposed; it is shown in Fig. 7. Here, the atomic mobility refers to the dispersion of the atomic mobility compared to the average atomic mobility at a given heating rate, which is indicated by different colors in Fig. 7. Strong atomic mobility means a larger dispersion of the atomic mobility compared to the average atomic mobility, and weak atomic mobility corresponds to a smaller dispersion of the atomic mobility compared to the average atomic mobility. Thus, with the experimental strategy of probing the evolution of the relaxation dynamics over a larger hearing rate range, new details about the potential energy landscape for a glass sample can be obtained. These new details can be extended to allow the study of the chemical composition and the effects of external stimuli on the evolution of the energy landscape. Fig. 7 Scheme of the evolution of the thermal relaxation spectrum, atomic dynamic heterogeneity and activated energy landscape at different heating rates.