Structure-dynamics relation in metallic glass revealed by 5-dimensional scanning transmission electron microscopy

Introduction
Glass transition is a universal phenomenon observed in a wide range of materials, including metals, polymers, and ceramics, which were cooled from the melt quickly enough to avoid crystallization and freeze the constituent atoms or molecules in disordered positions1,2. Since the glass phase has different physical properties compared with a liquid or a crystal, and the glass transition is strongly related to crystallization and melting, this phenomenon is important for controlling physical and technological properties such as moldability3. However, the process of freezing the atomic dynamics is still not understood. In the glass transition phenomenon, the viscosity changes rapidly near the glass transition temperature (Tg), even though the average structure remains almost unchanged4. Hence, to explain the viscosity variations, researchers focused on the microscopic structure of glasses and supercooled liquids. Dynamical heterogeneity was discovered in simulations and experiments on colloidal systems5,6, in which the motion of particles is heterogeneous and results in clustering. The dynamical heterogeneity shows a divergent behavior near Tg and is expected to explain the drastic change in viscosity during the glass transition. Furthermore, it has been revealed by electron tomography and other methods that there is a spatial heterogeneity of atomic arrangement in glasses7,8,9,10,11. This heterogeneity is called structural heterogeneity. The relation between dynamical and structural heterogeneities is important. If a common atomic structure of the slow-dynamics region in the dynamical heterogeneity is detected, it may become possible to identify the structures that cause the slow dynamics (high viscosity) and are necessary for the glass transition. Since atomic motion is affected by the surrounding atomic structure, there must be a relation between them. In fact, the relation between the structure and the temperature dependence of viscosity (fragility) has already been reported12,13, and a relation between heterogeneities has been suggested,14,15,16,17,18 but has not been confirmed experimentally so far, possibly because of the scarcity lack of observation methods that can simultaneously probe both the structural and dynamical heterogeneities. Those methods require a spatial resolution of a few nanometers, a time resolution of a few seconds, and the ability to measure the local structure of glass.
Experimental methods for characterizing dynamics at the atomic scale in glass began with measuring average values using X-ray photon correlation microscopy (XPCS)19,20,21, progressed to probing local dynamics with electron correlation microscopy (ECM)22,23, and advanced to visualizing dynamical heterogeneities using dark-field electron correlation microscopy (DF-ECM)24. These methods did contribute to the understanding of glass dynamics, but they suffered from low spatial resolution (XPCS), inabilities to observe the spatial distribution (ECM) or local structures (DF-ECM). However, a novel technique of 5-dimensional scanning transmission electron microscopy (5D-STEM)25,26,27 that involves convergent electron beam diffraction (CBED) does meet the requirements mentioned above. 5D-STEM is an extension of 4D-STEM28. It consists of 2D real-space scanning and 2D diffraction patterns, with the addition of a time dimension, and enables the measurement of the spatiotemporal distribution of CBED patterns. These patterns can yield local information that is hard to retrieve in the conventional parallel-beam diffraction mode. When a parallel electron beam irradiates a macroscale volume of a glass, a halo pattern is observed. However, when the probed volume shrinks to the nanoscale, intensity fluctuations, called speckles, appear instead29,30. The speckle pattern reflects the local order of the atomic structure31,32,33,34, and its temporal change reflects the motion of the local atomic structure22. Thus, as illustrated in Fig. 1, the local structure of a glass can be analyzed from the spatial distribution of diffraction patterns, and the local atomic rearrangement can be measured from the corresponding time sequences.

(Vleft(Qright)) is calculated from each diffraction pattern and averaged over time to make fluctuation map. (gleft(Q,tright)) is calculated from sequential CBED images at each position according to equations (1) and (2). By fitting the KWW function to (gleft(Q,tright)) at each position, relaxation time map can be obtained. Examples of sequential CBED images whose relaxation times are slow and fast are shown in bottom left and their (gleft(Q,tright)) are shown in bottom center panel.
In this study, we used 5D-STEM to simultaneously observe dynamical and structural heterogeneities and analyze the relation between them in a Zr50Cu40Al10 metallic glass. We also conducted an in situ heating experiment to measure the temperature dependences of the dynamical and structural heterogeneities and the structure-dynamics relation below Tg.
Materials and methods
Zr50Cu40Al10 metallic glass was fabricated by the tilt-melting method. Its glass transition and crystallization temperatures were 673 K and 750 K at the heating rate of 1.4 min/K35. TEM samples were fabricated by focused ion beam milling. STEM observation was conducted at 200 kV in an aberration-corrected JEM-ARM200F microscope (JEOL. Ltd.) equipped with a cold field emission gun and a 4DCanvas camera (JEOL. Ltd.). In the normal operation condition, the convergence semi-angle of the electron probe was in the range of 10 to 30 mrad (which corresponds to 6.0 to 18.0 nm–1). This value is much larger than the typical diffraction angle of speckles, and hence it was not possible to observe them. Therefore, the condenser lens and one of the transfer lenses of the corrector were adjusted to achieve an appropriate convergence semi-angle of 1.6 mrad (0.95 nm–1) that results in a probe diameter of 0.78 nm (full width at half maximum). The probe current was 8.9 pA. We used a beam stop to save the detector from the transmitted electron beam. The shadow of the beam stop was excluded from the analysis by excluding this area from the calculation of eqs. (1) and (5). To investigate temperature dependences of dynamical and structural heterogeneities and structure-dynamics relations, we heated the sample from 633 to 673 K in 10 K increments. Each diffraction pattern was acquired within 1 ms at 633 and 643 K and within 0.5 ms at 653, 663, and 673 K. The number of scan points and the step were 60 × 120 and 0.16 nm, respectively. The total observation area was 9.7 × 19.4 nm. Since 5D-STEM requires a sample edge for data calibration, we set the observation area as a horizontal rectangle. A total of 106 diffraction patterns were acquired for each scan point with a 9.2 s time resolution at 633 and 643 K and a 4.8 s time resolution at 653, 663, and 673 K. The total observation time was 1048 s at 633 and 643 K and 524 s at 653, 663, and 673 K. The pixel size and number of pixels for diffraction patterns were 0.058 nm–1 and 264 × 264. The diffraction patterns were 4 × 4 binned to increase the signal-to-noise ratio. The resultant pixel size and the number of pixels were 0.23 nm–1 and 66 × 66. The diffraction patterns near the sample edge (the left part of the sample) were too noisy for reliable measurements because of their thin thickness. Therefore, we used only the right half of the observed area whose thickness range is 40–50 nm in this article (see Fig. S1) to reduce thickness dependence and eliminate spurious signals.
Results and discussion
The local dynamics were evaluated from the temporal change of the local diffraction patterns using the following equations:16,20,21,36,37
Here ({boldsymbol{Q}}) is the position vector in reciprocal space, its absolute value (Q) is the spatial frequency and (t), ({t}_{1}) and ({t}_{2}) are acquisition times. (Gleft(Q,{t}_{1},{t}_{2}right)) measures the correlation between diffractions acquired at times ({t}_{1}) and ({t}_{2}). (Ileft({boldsymbol{Q}},{t}_{1}right)) is the signal intensity at ({boldsymbol{Q}}) and ({t}_{1}), and 〈⋯〉ϕ the subscript (phi) indicates that the averaging is performed on the ensemble of pixels over a range of spatial frequency. In this study, we limited the frequency to the range of 4.2–4.6 nm–1 that covered the first sharp diffraction peaks. Strong speckles are observed when the electron beam irradiates an ordered atomic structure along its symmetry axis. Speckles are formed when the electron beam irradiates an ordered atomic structure along its symmetry axis. Thus, this method probes the dynamics of the ordered atomic structures whose crystal symmetry axis is aligned to the beam direction. (gleft(Q,tright)) takes the average of (Gleft(Q,{t}_{1},{t}_{2}right)) over t1, at a fixed delay time, (t={t}_{2}-{t}_{1}). The Kohlrausch-Williams-Watts (KWW) function was fitted to (gleft(Q,tright)) to measure the dynamics.
Here (beta), (tau) and (gamma) are a scaling constant, the relaxation time, and the stretch factor, respectively. The stretch factor represents the nonlinearity of relaxation and indicates whether relaxation is compressed ((gamma > 1)) or stretched ((gamma < 1)) as compared to the simple exponent. Relaxation time is an indicator of dynamics; long relaxation time indicates slow dynamics and vice versa. We measured the relaxation time from all scan points to obtain a relaxation time map that directly reflects the dynamical heterogeneity. The relaxation time map for 633 K is presented in Fig. 2a. As seen from the map, the relaxation time spatially varies from 200 s to 5000 s revealing a heterogeneity of local dynamics. As shown in Fig. 2b–e, although the relaxation times vary, the dynamical heterogeneity has been observed at different temperatures. The average relaxation time is plotted in Fig. 2f; it shortens as the temperature approaches Tg, indicating an increase in atomic mobility. These values show good agreement with the previous research of Zr–Cu–Al systems, and such temperature dependence was observed in a series of metallic glasses using various measuring methods38,39,40. In samples with similar composition (Cu59Zr41), the relaxation process below Tg is thought of as (beta)-relaxation. From the energy landscape perspective, (beta)-relaxation is recognized as hopping across sub-basins inside an identical mega-basin. The activation energy to cross the sub-basins can be estimated from the temperature dependence of the relaxation time below Tg. The relation between relaxation time and activation energy is expressed by the Arrhenius law.
Where τ∞ is the relaxation time at infinite temperature, ({E}_{beta }) is the activation energy and (R) is the gas constant. The activation energy was calculated as 97(pm ,) kJ/mol. This value falls within the range of the previously reported values, namely, 46 kJ/mol by internal friction measurements in Cu59Zr4141, 96–159 kJ/mol by MD simulations of the Zr-Cu system42,43,44, and 62–174 kJ/mol by dynamic mechanical analysis of Zr-Cu-Al alloys45,46,47.

Maps of relaxation time at a 633 K, b 643 K, c 653 K, d 663 K, and e 673 K. f Temperature dependence of the spatially averaged relaxation time. The error bar shows the distribution of relaxation time.
To evaluate the structural heterogeneity, we measured the intensity fluctuation along the azimuthal direction around the halo pattern. This intensity fluctuation mainly originates from the speckle pattern. The fluctuation is evaluated by the following equation:
Since the speckle pattern reflects the local atomic structure, so does the value of (Vleft(Q,{boldsymbol{r}}right))—a higher (Vleft(Q,{boldsymbol{r}}right)) corresponds to a more ordered structure32,48. It should be noted that the value of (Vleft(Q,{boldsymbol{r}}right)) is also influenced by the direction of the symmetry axis. However, since glass is isotropic, such effects can be greatly reduced by statistical sampling, making (Vleft(Q,{boldsymbol{r}}right)) a reliable indicator for measuring overall trends and correlations. In this experiment, (Vleft(Q,{boldsymbol{r}}right)) was measured from diffractions whose spatial frequency ranged from 3.0 to 5.8 nm–1. We measured (Vleft(Q,{boldsymbol{r}}right)) for all scanning points and averaged it over time. Figure 3a shows the calculated fluctuation map at 633 K. The values range from 0.45 to 0.75 and show heterogeneous distribution, indicating the coexistence of different local atomic structures. This image directly visualizes the structural heterogeneity of metallic glass. Thus, it was confirmed that the local atomic structure was heterogeneous. As shown in Fig. 3b–e, although the (Vleft(Q,{boldsymbol{r}}right)) values vary, the structural heterogeneity has been observed at different temperatures. The temperature dependence of average (Vleft(Q,{boldsymbol{r}}right)) is shown in Fig. 3f. The average (Vleft(Q,{boldsymbol{r}}right)) slightly decreased with temperature from 0.56 to 0.52, reflecting the decrease in atomic order. A similar slight decrease in the average local order was theoretically predicted for Zr47.5Cu47.5Al5 glass49.

Maps of (Vleft(Q,{boldsymbol{r}}right)) at a 633 K, b 643 K, c 653 K, d 663 K, and e 673 K. f The average (Vleft(Qright)) at all temperatures. The errorbar shows the distribution of (Vleft(Qright)). (Q) is the spatial frequency and ranged from 3.0 to 5.8 nm−1 in this calculation.
As mentioned above, the dynamical and structural heterogeneities can be measured simultaneously by 5D-STEM. The structure-dynamics relation was evaluated from these heterogeneities by correlating the maps of relaxation time and (Vleft(Q,{boldsymbol{r}}right)). Since the local order was probed by the speckle pattern and the local relaxation time was measured from the temporal change of the same pattern, the local atomic order and atomic rearrangements were measured from the same nanostructure, and the correlation between dynamical and structural heterogeneities directly reflects the structure-dynamic relation. The degree of relaxation depends both on the relaxation time and the value of γ50. To simplify the analysis, considering that (gamma) values are close to 1 (see Figs. S1 and S2), we fixed (gamma) at 1.0. The correlation was calculated by Spearman’s rank correlation, where the correlation value can range from −1 to 1. Positive and negative values mean positive and negative correlations between the maps of relaxation time and (Vleft(Q,{boldsymbol{r}}right)), while 0 means lack of correlation. A scatter plot is shown in Fig. 4 with (Vleft(Q,{boldsymbol{r}}right)) on the vertical axis and relaxation time on the horizontal axis. At 633 K, the distribution of the points is stretched from the center to the upper right, and the correlation value is 0.39, statistically indicating that regions with a long relaxation time and regions with a partial order tend to coincide. Thus, ordered regions tend to have slow dynamics. This tendency is also observed in simulations of polydisperse particles14. It is noteworthy that not all structural heterogeneities exhibit correlations with dynamical heterogeneities. Structural heterogeneities calculated by density or potential energy fields do not exhibit a positive correlation with dynamical heterogeneities14,51. Simulation studies have demonstrated that the multi-body effects such as flexibility17, deviations from sterically favored structures14, or softness18 are significant in explaining dynamical heterogeneity, as the relaxation process is affected by surrounding atoms or particles. The diffraction or speckle pattern formed by the interference of electrons scattered by atoms naturally includes information about surrounding atoms31. Therefore, a positive correlation has been observed between structural heterogeneity measured by (Vleft(Qright)) and dynamical heterogeneity. The same tendency is observed at higher temperatures, as shown in Fig. 4b–e, and especially in Fig. 4f, which superimposes the scatter plots for all temperatures. Figure 4f reveals that the relaxation time shortens at higher temperatures. Meanwhile, the shift in the distribution toward the lower left corner indicates that the number of less regular structures that exhibit higher atomic mobility increases at high temperatures.

Scatter plots of relaxation time and (Vleft(Qright)) at a 633 K, b 643 K, c 653 K, d 663 K, e 673 K. f Superimposed scatter plots at all temperatures.
In summary, our results can be interpreted as follows: at low temperatures, the atomic rearrangement was slow. As the temperature increased, the atomic motion sped up, and the degree of local atomic order decreased. The atomic motion was slow for ordered regions regardless of temperature. The atomic motion accelerated near Tg for two reasons: (1) an increase in the number of structures that have a low degree of order and hence a high atomic mobility, and (2) an increase in the atomic mobility in structures with any degree of order.
Structural and dynamical heterogeneities are important for elucidating the mechanism of glass transition. However, the relationship between these heterogeneities is difficult to analyze because of the lack of experimental methods that can simultaneously measure these heterogeneities. In this study, we simultaneously visualized the dynamical and structural heterogeneities in a Zr50Cu40Al10 metallic glass by recording the spatiotemporal distribution of diffraction patterns via 5D-STEM. We also measured the temperature dependencies of relaxation time and structural order. The results directly reveal a positive correlation between dynamical and structural heterogeneities, indicating that atomic rearrangement is slower in more ordered atomic structures. This study demonstrates the high potential of 5D-STEM in monitoring the dynamical and structural heterogeneities in glasses, yet the reported calculation of (Vleft(Q,{boldsymbol{r}}right)) can only hint at the presence of the ordered structure. Therefore, a better correlation of structure and dynamics would require the development of new analytical methods for diffraction patterns that can identify translational symmetry, rotational symmetry, and whether the structural order is energy-driven or entropy-driven, as proposed for atomic structures in real space52.
Responses