Structural origins of Johari-Goldstein relaxation in a metallic glass

Abstract

Johari-Goldstein or β relaxation, persisting down to glassy state from a supercooled liquid, is a universal phenomenon of glassy dynamics. Nevertheless, the underlying micromechanisms leading to the relaxation are still in debate despite great efforts devoted to this problem for decades. Here we report experimental evidence on the structural origins of Johari-Goldstein relaxation in an ultra-quenched metallic glass. The measured activation energy of the relaxation (~26 times of the product of gas constant and glass transition temperature) is consistent with the dynamic characteristics of Johari-Goldstein relaxation. Synchrotron X-ray investigations demonstrate that the relaxation originates from short-range collective rearrangements of large solvent atoms, which can be realized by local cooperative bonding switch. Our observations provide experimental insights into the atomic mechanisms of Johari-Goldstein relaxation and will be helpful in understanding the low-temperature dynamics and properties of metallic glasses.

ARTICLE

Received 14 Dec 2013 | Accepted 10 Jan 2014 | Published 3 Feb 2014

Johari-Goldstein or b relaxation, persisting down to glassy state from a supercooled liquid, is a universal phenomenon of glassy dynamics. Nevertheless, the underlying micromechanisms leading to the relaxation are still in debate despite great efforts devoted to this problem for decades. Here we report experimental evidence on the structural origins of Johari-Goldstein relaxation in an ultra-quenched metallic glass. The measured activation energy of the relaxation (B26 times of the product of gas constant and glass transition temperature) is consistent with the dynamic characteristics of Johari-Goldstein relaxation. Synchrotron X-ray investigations demonstrate that the relaxation originates from short-range collective rearrangements of large solvent atoms, which can be realized by local cooperative bonding switch. Our observations provide experimental insights into the atomic mechanisms of Johari-Goldstein relaxation and will be helpful in understanding the low-temperature dynamics and properties of metallic glasses.

DOI: 10.1038/ncomms4238

Structural origins of Johari-Goldstein relaxation in a metallic glass

Y.H. Liu1,*, T. Fujita1,*, D.P.B. Aji1,*, M. Matsuura2 & M.W. Chen1,3

1 WPI Advanced Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan. 2 Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan. 3 State Key Laboratory of Metal Matrix Composites, School of Materials Science and Engineering, Shanghai Jiao Tong University, Shanghai 200030, China. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to M.W.C.(email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 5:3238 | DOI: 10.1038/ncomms4238 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4238

Johari-Goldstein (JG) relaxation has intrinsic correlation with glass transition and many fundamental properties of glassy materials15. It is generally explained in term of hindered

intra- and/or inter-molecular motion that remains active even when the whole system is frozen into the glassy state1. However, the underlying micromechanisms leading to JG relaxation have not been fully understood213, mainly because of the lack of well-dened experimental characterization on the structural evolution of the disordered glasses. It remains indistinct that the relaxation is essentially involved all atoms/molecules that happens spatially uniformly or originates from local atomic rearrangement associated with structural and chemical heterogeneity2,3,1116. Different from organic and net-work oxide glasses, metallic glasses are basically atomic glasses with fundamental components of individual atoms and can be treated as randomly packed hard spheres without complex inter- and intra-molecular interactions. Besides their technical importance for a wide range of applications, metallic glasses are also an ideal system for understanding the basic physical process of structural relaxation in glassy solids on an atomic scale.

Since the relaxation of conductive metallic glasses cannot be studied using dielectric spectroscopy, JG relaxation in metallic glasses is mainly characterized by differential scanning calorimetry (DSC)7,8,16 and dynamical mechanical spectroscopy9,10,17,18. Unlike spectroscopic method, the DSC measurements are an out-of-equilibrium and nonlinear approach. JG relaxation is superimposed on the broad sub-Tg exothermic relaxation patterns that are the direct outcome of combined secondary and primary relaxations, while JG relaxation appears at the low-temperature side as the tails of the primary relaxation bands7,8,16. It has been suggested that the relaxation of metallic glasses depends on fragility, cooling rate, chemical bonding, atomic size disparity and constitute atom diffusivity3,8,9,1619. Since these kinetic and chemical variables

are essentially associated with structure of metallic glasses, directly investigating the correlation between atomic structure and JG relaxation may shed a light on the atomic origins of the low-temperature relaxation. However, metallic glasses usually have a weak JG relaxation in proportion to trivial structure evolution because the high atomic mobility gives rise to fast structure relaxation during cooling from liquid to glass state. Although it has been noticed that the intensity of JG relaxation can be enhanced by high cooling rates16,20, the intensity enhancement from the currently feasible rapid quenching techniques, such as melting spinning, is still insufcient for detectable structure changes because of the maximum cooling rate is only B106 K s 1.

In this study we developed an ultra-quenched metallic glass with an ultrahigh cooling rate of B3.6 109 K s 1 which is

B36 orders of magnitude higher than conventional metallic glasses. The ultrahigh cooling rate effectively freezes the glass at a high-energy state with a high ctive temperature and, hence, amplies JG relaxation for structure study. The measured activation energy for the low-temperature relaxation of the ultra-quenched metallic glass, B1.39 eV (B26RTg, where R is the gas constant and Tg is the glass transition temperature), is consistent with the characteristics of JG relaxation. Systematic structural characterization using synchrotron X-ray scattering demonstrates that the JG relaxation mainly arises from short-range conguration transformation of solvent atoms by local cooperative bonding switch.

ResultsThermal measurements. Figure 1a presents the apparent heat capacity, Cp, during the rst up-scan of the ultra-quenched lm, ribbon and bulk glass samples, as well as that of the second up-scan of the slowly cooled specimen within DSC at the rate of

0.8

0.20 Ultra-quenched

Bulk

0.7

0.15

1 K1 )

0.6

Ribbon

1 K1 )

C p(Jg

Standard

Bulk Ribbon

0.5

C p (J g

0.10

0.4

0.05

0.3

Ultra-quenched

0.00 450 500 550 600 650 700

Tg

0.2 450 500 550 600 650 700

T (K)

T (K)

0.20

14

0.15

10 K min1

1 K1 )

20 K min1

13

40 K min1

0.10

C p (J g

80 K min1

ln (T2 / Q ) (Ks)

160 K min1

Ebeta=1.39 eV

12

0.05

0.00 450 500 550 600 650 700

11 1.7 1.8 1.9 2.0 2.1

T (K)

1,000/T (K1)

Figure 1 | Thermal characterizations with differential scanning calorimetry. (a) Apparent heat capacity Cp of the ultra-quenched, ribbon, bulk and standard samples. (b) Activation spectra of the three glass samples that experienced high cooling rates. (c) Activation spectra of the ultra-quenched glass measured at different heating rates. (d) Kissinger plot of the onset temperature versus heating rate for the secondary relaxation.

2 NATURE COMMUNICATIONS | 5:3238 | DOI: 10.1038/ncomms4238 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4238 ARTICLE

a

0.33 K s 1 (dened as the standard sample as the reference). The DSC curves indicate that there is a weak cooling rate dependence of the DSC glass transition temperatures (Tg) of the four samples upon heating. The mean value of Tg is B690 K. Different from the standard sample, broad exothermic peaks can be seen from the lm, ribbon and bulk glass samples that experienced faster cooling. The release of the trapped enthalpy by the JG relaxation starts at an onset temperature of B500 K, which is independent of cooling rates. However, the amplitudes of the JG relaxation show the strong dependence on the cooling history of the samples. Compared with the ribbon and bulk samples, the lm exhibits a much more pronounced DCp peak (Fig. 1b), where DCp characterizes the heat capacity difference between the as-prepared and the standard samples. As illustrated in Fig. 1c, the relaxation peak of the lm shifts to higher temperatures with heating rates increasing, indicating the strong rate dependence of the JG relaxation. Figure 1d is the Kissinger plot of the onset temperatures (Tc) of the relaxation at different heating rates. The measured activation energy, Ebeta, is determined to be 1.39 eV, which is well consistent with the reported values of Zr-based glasses obtained by the enthalpy relaxation approach8. It is worth noting that the Ebeta value is in accordance to the universal activation energy (B26RTg) of JG relaxation, the dynamic characteristics widely observed in glassy materials8,9,17,18,21,22.

Although the broad sub-Tg exothermic peak in the DSC traces of the ultra-quenched metallic glasses may be the outcome of combined JG and primary relaxations together with other possible sources, such as irreversible structural relaxation, the dynamic identity of the activation energy provides the solid evidence that the low-temperature side of relaxation spectra and the onset of Cp mainly originate from JG relaxation, while the contribution of the structural relaxation is minor.

We evaluated the energy states of the glassy samples by calculating their ctive temperatures, Tf, at which the structure of the equilibrium liquid is frozen into a glass state. According to the energy matching method23, Tf can be derived by equalling the area hatched by the heat capacity curves of the rst (Cp1) and

second (Cp2) up-scans to that enclosed by Cpl and Cpg from Tg to Tfthat is, R

Teq

TC Cp2 Cp1

12 Y = 37 37X

10 Film

Q =3.6109

Ribbon

6 Q =1.8106

Rod

Q =5.8103

8

Log (Q/Q s)

4

2

s =1.34

s =1.21

s =1.11

0

T f/T f

T f/T f

T f/T f

2 0.5 0.6 0.7 0.8 0.9 1.0 1.1

Tfs/Tf

b

m~40

5

4

lnQ(K min1 )

3

2

dT, where Tc is

the temperature at which energy release starts; Teq is the

temperature where the equilibrium liquid is reached; Cpl and

Cpg are the heat capacities of the supercooled liquid and glass, respectively. On the basis of the equation, the ctive temperatures of the ribbon and bulk glasses are calculated to be 1.21Tg and1.11Tg, respectively, in line with the typical values of conventional metallic glasses. For the ultra-quenched lm glass, the Tf is as large as 1.34Tg, approaching the melting temperature Tm of the alloy (Tf 0.83Tm), verifying that the lm glass is trapped at a

high-energy state far away from equilibrium.

Figure 2a shows the Arrhenius plot of the scaled cooling rate Q/Qs and ctive temperature Tf/Tsf by the standard sample.

Since the standard sample experienced a constant slow rate of Qs 0.33 K s 1 on cooling and heating, its ctive temperature Tsf

is approximately equal to Tg ( 690 K). The nominal cooling

rates of the glasses can be determined by the scaled Arrhenius plot on the basis of their estimated ctive temperatures in a calibration line24. In addition to the standard sample with the cooling rate of 0.33 K s 1 and ctive temperature of 690 K, the calibration line is constructed by the designed cooling rates of8.33 and 3.33 K s 1 using DSC. According to the derived Tf values, the cooling rate of the lm is as high as B3.6 109 K s 1,

about 3 and 6 orders of magnitude higher than that of the ribbon (1.8 106 K s 1) and bulk (5.8 103 K s 1) samples.

The slope of the plot in Fig. 2a yields the fragility (m) of B37, which is related to the activation energy Ea for glass transition by m Ea/2.303RTg (ref. 24). The fragility (m 37) is in the good

agreement with that (m 40) obtained by measuring the heating

rate dependence of Tg (Fig. 2b)25, validating that the Tf and cooling rates determined from the calibration line are dependable.

Synchrotron X-ray diffraction characterization. We utilized synchrotron radiation X-ray diffraction (XRD) and the extended X-ray absorption ne structure (EXAFS) spectroscopy to experimentally investigate the structural evolution of the glassy samples caused by JG relaxation. Figure 3a displays the pair distribution functions (PDFs), g(r), of ultra-quenched, bulk and standard glasses, which were derived from XRD data. The PDFs reveal short to medium-range ordering up to B15 , beyond which conspicuous correlation peaks cannot be observed. Close inspection indicates that the rst peaks of the Zr55Cu30Ni5Al10

glasses split into two sub-peaks locating at B2.80 and 3.15 (Fig. 3b) due to the discrete bonding length distribution of the nearest neighbour atoms. According to the Goldschmidt radii and weight factors (Table 1), the rst sub-peak can be assigned mainly to the Zr-Cu pairs, while the Zr-Zr pairs dominate the second

dT R

TfTg Cpl Cpg

1 680 690 700 710 720 730

Tg (K)

Figure 2 | Determination of ctive temperatures. (a) Scaled Arrhenius plot of ctive temperature and cooling rate. (b) Dependence of glass transition temperature on heating rate. The solid line is the VogelFulcher Tamman equation tting to the data points.

NATURE COMMUNICATIONS | 5:3238 | DOI: 10.1038/ncomms4238 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4238

2.5 Ultra-quenched

Bulk Standard

2.5 Ultra-quenched

Bulk Standard

2.0

2.0

1.5

1.5

g(r)

g(r)

1.0

1.0

0.5

0.5

0.0 2 4 6 8 10 12 14 16 18 20

0.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8

r ()

r ()

1.4

1.2

1.2

1.1

1.0

g(r)

g(r)

1.0

0.8

0.6

Ultra-quenched BulkStandard

0.9

Ultra-quenched BulkStandard

0.4 4.0 4.5 5.0 5.5 6.0 6.5

0.8 6.8 7.2 7.6 8.0 8.4 8.8

r () r ()

Figure 3 | Pair distribution functions derived from synchrotron XRD. (a) The overall tendency of the pair distribution functions. (bd) Close inspections of the rst, second and third peaks of PDFs.

Table 1 | Inter-atomic bond lengths and weight factors of atomic pairs in the Zr55Cu30Ni5Al10 metallic glass.

i-j rij () Weight factor Heat of mixing (kJ mol 1)

Zr-Zr 3.16 0.377

Zr-Cu 2.85 0.345 23

Zr-Ni 2.86 0.080 49

Zr-Al 3.01 0.050 44

Cu-Cu 2.54 0.079

Cu-Ni 2.55 0.036 1

Ni-Ni 2.56 0.004

Cu-Al 2.70 0.023 1

Ni-Al 2.71 0.005 22

Al-Al 2.86 0.002

sub-peak. The large weight factor of the Zr-Cu pairs is consistent with the fact that the Cu-centred clusters are the dominant short-range order as demonstrated by molecular dynamics simulation26. Compared with the standard and bulk samples, the ultra-quenched lm shows obviously higher g(r)Zr-Zr/g(r)Zr-Cu ratio,

indicating more Zr-Zr pairs and lesser Zr-Cu pairs in the sample. By integrating the rst peak of the PDFs, the average coordination numbers are found to be B13.6 for all the three glasses, indicating very small variation in the global average coordination numbers. The cooling rate-induced difference in the coordination number is probably too small to be identied. However, the expanded width of the rst peak of the lm suggests that the atomic packing in the nearest neighbours is looser compared with the standard and bulk samples (Fig. 3b), which coincides with the high fraction of the longer Zr-Zr pairs. For the second and third peaks of the lm glass, substantial peak broadening cannot be seen. Instead, their intensities are noticeably lower than those of

the standard sample (Fig. 3c,d). This indicates the loss of medium-range order in the ultra-quenched lm glass, which is presumably due to the expanded nearest neighbours that cause the structure frustration beyond the rst shells27.

EXAFS results. The local atomic structure of the glass samples is further characterized by EXAFS, which is highly sensitive to local elemental specicity. Figure 4a,b shows the EXAFS k2w(k) spectra of Cu-K and Zr-K edges from ultra-quenched, bulk and standard glass samples, respectively. The spectra of the three samples at the Cu-K edge are nearly identical, indicating that the huge cooling rate difference among the samples does not cause detectable changes in the nearest neighbours of Cu atoms. However, at the Zr-K edge the spectrum of the ultra-quenched glass shows detectable difference from the bulk and standard samples, particularly at the low-k region. The Fourier transform proles of k2w(k) spectra of Cu and Zr-K edges are shown in Fig. 4c,d. Two sub-peaks, corresponding to Cu-Cu and Cu-Zr pairs, can be seen in the Cu-K-edge spectra (Fig. 4a), signifying that the rst shell of the Cu-centred clusters includes both Cu and Zr solvent atoms. While the standard specimen is almost fully relaxed in comparison with the ultra-quenched glass, its peak intensity ratio between Cu-Cu and Cu-Zr pairs remains the nearly same as those of the samples with higher cooling rates. This indicates that the pronounced JG relaxation of the lm glass does not entail visible alteration of the local chemical environments of the prevalent Cu-centred polyhedra in the metallic glass. In contrast, the Fourier transform of the Zr-K edge k2w(k) spectra shows the noticeable difference of the ultra-quenched glass from the bulk and standard ones in the magnitude ratios of the peaks of Zr-Zr and Zr-Cu pairs (Fig. 4d). As the peak magnitude ratio corresponds to the number fraction of neighbouring atoms, considering the inconspicuous difference in the Fourier-

4 NATURE COMMUNICATIONS | 5:3238 | DOI: 10.1038/ncomms4238 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4238 ARTICLE

0.8

Zr-K edge

0.8 Ultra-quenched BulkStandard

Cu-K edge

0.6

0.6

Ultra-quenched BulkStandard

0.4

0.4

EXAFS(k)k2

0.2

0.2

0.0

EXAFS (k)k2

k (1) k (1)

0.0

0.2

0.2

0.4

0.4

0.6

0.6

0.8 0 2 4 6 8 10 12

0.8 0 2 4 6 8 10 12 14

0.8

0.7

Cu-K edge

Zr-K edge

0.7

Ultra-quenched BulkStandard

0.6

Ultra-quenched BulkStandard

0.6

Cu-Zr

0.5

Zr-Cu

|(r)| (3 )

0.5

0.4

|(r)| (3 )

r () r ()

0.4

Cu-Cu

0.3

Zr-Zr

0.3

0.2

0.2

0.1

0.1

0.0 0 1 2 3 4 5

0.0 0 1 2 3 4 5

Figure 4 | EXAFS spectra of the ultra-quenched, bulk and standard samples. (a) k2w(k) spectra of Cu-K-edge; and (b) Zr-K-edge. (c) Fourier transform proles of the experimental EXAFS spectra of Cu-K-edge, and (d) Zr-K-edge.

transformed Cu-K edges of the three samples (Fig. 4c), it can be concluded that the fraction of Zr-Cu pairs (2.85 ) evidently increases, while that of longer Zr-Zr pairs (3.16 ) decreases in the relaxed standard and bulk sample. This local chemical measurements are fully consistent with the statistic geometry analysis by XRD (Fig. 3b) in which the Zr-Zr pairs have a higher fraction than that of Zr-Cu in the ultra-quenched lm.

To further verify the local structural changes mainly originating from JG relaxation, we investigate the atomic conguration evolution of the isothermally annealed ultra-quenched glass at 623 K, close to the peak temperature of JG relaxation, for 30 min and 12 h. Since this temperature is far below Tg (690 K), the primary relaxation (glass transition) and irreversible structural relaxation should not take place. The EXAFS k2w(k) spectra of

Cu-K and Zr-K edges from the annealed ultra-quenched glass lms are presented in Fig. 5a,b, while their Fourier transform proles are displayed in Fig. 5c,d. As seen from the Cu-K edge, the peak positions as well as the magnitude ratios of the Cu-Cu and Cu-Zr pairs remain constant with the annealing time increasing from 30 min to 12 h (Fig. 5c). Nonetheless, the isothermal annealing-induced structure changes are visible from the Zr-K edge. The intensity magnitude of the Zr-Zr peak detectably decreases with annealing time (Fig. 5d), further conrming that JG relaxation mainly involves the reduction of solvent Zr-Zr pairs.

DiscussionBoth XRD and EXAFS characterizations demonstrate that very local chemical ordering around the solvent Zr atoms takes place during JG relaxation. The reduced number fraction of Zr-Zr pairs during the JG relaxation is expected to lead to the formation of more Cu-centred polyhedra. Although we cannot detect the Al- and Ni-related atomic pairs because of their low weight

factors (Table 1), the EXAFS analysis demonstrates that the local structure changes during JG relaxation mainly arise from the local rearrangement of the solvent Zr atoms with a larger size, whereas the chemical constitution of the dominant Cu-centred clusters remain nearly unaffected. On the basis of the XRD and EXAFS measurements, the key structure changes of the lm glass upon JG relaxation can be summarized as the following: (1) the packing density of the nearest neighbours increases as evidenced by the difference in the rst peak widths of the PDF proles (Fig. 3b); and (2) the number fraction of Zr-Zr atomic pairs in the nearest neighbours of Zr atoms decreases by substitution of Zr-Cu pairs, while the prevailing Cu-centred clusters do not show too much change in the number of solvent Zr atoms. Therefore, the physical process of JG relaxation of the ultra-quenched glass is apparently associated with a local atomic rearrangement beyond the dominant Cu-centred clusters and only solvent atoms, mostly Zr atoms, are rearranged by a short-range collective process. There are many possible dynamic paths to realize the local structure changes during JG relaxation. The most feasible one is probably by a simple bonding switch as depicted in Fig. 6. This easy cooperative atomic rearrangement can give rise to the key structure changes revealed by XRD and EXAFSthat is, the decreased Zr-Zr inter-atomic pairs (correspondingly, increased Zr-Cu pairs) and denser atomic packing. Since Zr atoms have a large atomic size and typically locate at the shells of densely packing clusters or adjunct sites between these clusters, some of them may be dynamically too slow to travel to a low-energy basin during fast quenching and hence are trapped in a high-energy state by bonding with other Zr atoms (a of Fig. 6). The JG relaxation process can be accomplished by a simple bonding switch from a long Zr-Zr pair to a tight Zr-Cu one through a cooperative short-range atomic motion as shown in b of Fig. 6. This process leads to the formation of a densely packed Cu-centred cluster and the loss of a long Zr-Zr pair. It is

NATURE COMMUNICATIONS | 5:3238 | DOI: 10.1038/ncomms4238 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4238

0.8 30 min at 623 K

12 h at 623 K

Cu-K edge

0.8 30 min at 623 K

12 h at 623 K

k (1)

Zr-K edge

0.6

0.6

0.4

0.4

0.2

0.2

EXAFS (k)k2

0.0

EXAFS (k)k2

0.0

0.2

0.2

0.4

0.4

0.6

0.6

0.8 0 2 4 6 8 10 12

0.8 0 2 4 6 8 10 12 14

k (1)

0.8 30 min at 623 K

12 h at 623 K

Cu-K edge

0.7 30 min at 623 K

12 h at 623 K

Zr-K edge

0.7

0.6

0.6

Zr-Cu

Zr-Zr

Cu-Zr

0.5

0.5

0.4

|(r)|(3 )

0.4

Cu-Cu

|(r)|(3 )

0.3

0.3

0.2

0.2

0.1

0.1

0.0 0 1 2 3 4 5

0.0 0 1 2 3 4 5

r ()

r ()

Figure 5 | EXAFS spectra of the ultra-quenched glass annealed at 623 K. (a) k2w(k) spectra of Cu-K-edge; and (b) Zr-K-edge. (c) Fourier transform proles of the experimental EXAFS spectra of Cu-K-edge, and (d) Zr-K-edge.

Thermal relaxation

Mechanical rejuvenation

Zr Cu

Figure 6 | Schematic diagram of secondary relaxation by local bonding switch. The break of the long Zr-Zr pair distance and formation of the short Zr-Cu bond lead to the denser local atomic packing.

analogous to shear-induced bonding switch28 but, however, the difference is also obvious since the bonding switch during JG relaxation is solely assisted by thermal uctuation at high temperatures and can be revisable by being quenched from the same ctive temperature of the glass. The short-range collective rearrangement of solvent atoms leads to the increased packing density, while the chemical environments of individual Cu-centred clusters are not affected except the increase in their number density, which will not lead to any changes in the EXAFS Cu-K edge proles.

As illustrated by this study, JG relaxation mainly originates from the local structure transformation associated with large solvent atoms. This appears to be in contradiction of the recent report on the coupling of small atoms diffusion with JG relaxation in metallic glasses9. In fact, the two observations tell the same scenario but from different aspects that are emphasized. Our structure characterization reveals the atomic conguration changes upon JG relaxation, while the activation energy measurement by Yu et al.9 discloses the dynamic process of

high mobility of small atom resulted from JG relaxation. As shown in the bonding switch model (Fig. 6), the resultant structure changes, corresponding to JG relaxation, are accomplished by the collective rearrangements of large solvent atoms. However, the dynamic process of the collective atomic transformation is mainly realized by the motion of smaller atoms over a short rangefor example, Cu in this studybecause of their high mobility and coordination capacity of lling the excess open space among the large solvent atoms.

It has been noticed that JG relaxation at a temperature below Tg gives rise to signicant change in the volume of shear transformation zones (STZs)10, the basic ow units in metallic glasses. Although it has been long known that STZs are associated with heterogeneous atomic packing in metallic glasses29, the physical process of STZ formation is very vague and lacks distinct physical description. As the barrier energy of STZs is also equal to the activation energy of JG relaxation, it is rational to deduce that STZs share the same physical process with JG relaxation12. STZs may initiate from the regions with loosely packed solvent atoms, consistent with previous molecular dynamics simulations30, and are dynamically manipulated by stress-induced atom motion. Accordingly, the STZ formation may be considered as the result of dynamic coupling between applied stress and inherent structure relaxation.

In summary, we developed an ultra-quenched metallic glass that undergoes a cooling rate on the order of B1010 K s 1. The glass shows pronounced sub-Tg relaxation compared with conventional ribbon and bulk glasses. The measured activation energy demonstrates that the relaxation has the dynamic characteristics of JG relaxation. Atomic structural characterization suggests that the relaxation originates from short-range collective rearrangements of solvent atoms beyond densely packed atomic clusters. This study may shed light on the atomic

6 NATURE COMMUNICATIONS | 5:3238 | DOI: 10.1038/ncomms4238 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4238 ARTICLE

mechanisms of JG relaxation and bear important implications on the structural process of plastic deformation in glassy metals.

Methods

Sample preparation. The ultra-quenched Zr55Cu30Ni5Al10 metallic glass was fabricated by means of rf-magnetron sputtering deposition31. To acquire sufcient DSC and XRD signals, B30-mm thick glassy lms were prepared with the well-controlled sample temperature below 50 C. Ribbon and bulk glasses with the same chemical composition were made by conventional melt-spinning and copper-mold-casting methods, respectively. Slowly cooled samples were prepared by reheating the glassy sample to the supercooled liquid region and then cooling to ambient temperature at the rate of 20 K min 1.

DSC measurements. Thermal analysis was conducted by using a Perkin-Elmer 8500 DSC with which controllable cooling rate can be achieved up to 500 K min 1.

The samples were rst heated to their supercooled liquid region at the rate of20 K min 1 to obtain the heat capacity traces (Cp1) of the as-prepared samples, and then slowly cooled at the same rate to room temperature. The slowly cooled samples were then subjected to the second up- and down-scan following the same procedure to establish the standard heat capacity curves (Cp2).

Structural characterization. The structural characterization was performed at SPring-8 synchrotron radiation facility, Japan. High-energy XRD with the beam energy of 113.3 keV was carried out in transmission geometry. The diffraction data were analysed following the procedures in Kohara et al.32 The local atomic structure of the samples was characterized by the EXAFS spectroscopy33. EXAFS measurements of the Cu and Zr-K edges were conducted in a transmission mode at the bending magnet beam line BL01B1 of the synchrotron radiation facility (SPring-8) in Hyogo, Japan. Incident X-ray energy was selected by a xed-exit double-crystal monochromator using Si(111) planes for the Cu-K-edge and the Zr-K-edge. EXAFS oscillation functions w(k) and their magnitude of Fourier transforms w(r) are deduced from the EXAFS spectra using a coded program of Athena34. The energy edge (E0) was set to a set fraction (50%) of the edge step. The k-range of forward Fourier transform is from 2 to 11 1 for both Cu-edges and

Zr-edges. It is common practice in EXAFS to multiply w(k) by different k-weights. In this study, we used k2 to calculate w(r) ( 3).

References

1. Johari, G. P. & Goldstein, M. Viscous liquids and the glass transition. II. secondary relaxations in glasses of rigid molecules. J. Chem. Phys. 53, 23722388 (1970).

2. Johari, G. P. Localized molecular motions of beta-relaxation and its energy landscape. J. Non-crystal. Solids 307, 317325 (2002).

3. Ngai, K. L. Relaxation and Diffusion in Complex Systems (Springer, 2011).4. Thayyil, M. S., Capaccioli, S., Prevosto, D. & Ngai, K. L. Is the Johari-Goldstein beta-relaxation universal? Philos. Mag. 88, 40074013 (2008).

5. Yu, H. B., Wang, W. H. & Samwer, K. The b relaxation in metallic glasses: an overview. Mater. Today 16, 183191 (2013).

6. Stevenson, J. D. & Wolynes, P. G. A universal origin for secondary relaxations in supercooled liquids and structural glasses. Nat. Phys. 6, 6368 (2010).

7. Hornboll, L. & Yue, Y. Z. Enthalpy relaxation of hyperquenched glasses and its possible link to a and b-relaxations. J. Non-Crystal. Solid 354, 350354 (2009).

8. Hu, L. N. & Yue, Y. Z. Secondary relaxation in metallic glass formers: its correlation with the genuine Johari-Goldstein relaxation. J. Phys. Chem. C 113, 1500115006 (2009).

9. Yu, H. B., Samwer, K., Wu, Y. & Wang, W. H. Correlation between beta relaxation and self-diffusion of the smallest constituting atoms in metallic glasses. Phys. Rev. Lett. 109, 095508 (2012).

10. Pan, D. et al. Correlation between structural relaxation and shear transformation zone volume of a bulk metallic glass. Appl. Phys. Lett. 95, 141909 (2009).

11. Yu, H. B. et al. Tensile plasticity in metallic glasses with pronounced beta relaxations. Phys. Rev. Lett. 108, 015504 (2012).

12. Yu, H. B., Wang, W. H., Bai, H. Y., Wu, Y. & Chen, M. W. Relating activation of shear transformation zones to beta relaxations in metallic glasses. Phys. Rev. B 81, 220201 (2010).

13. Tanaka, H. Origin of the excess wing and slow beta relaxation of glass formers: a unied picture of local orientational uctuations. Phys. Rev. E 69, 021502 (2004).

14. Fujita, T. et al. Atomic-scale heterogeneity of a multicomponent bulk metallic glass with excellent glass forming ability. Phys. Rev. Lett. 103, 075502 (2009).

15. Fujita, T. et al. Coupling between chemical and dynamic heterogeneities in a multicomponent bulk metallic glass. Phys. Rev. B 81, 140204 (R) (2010).16. Chen, H. S. & Coleman, E. Structure relaxation spectrum of metallic glasses. Appl. Phys. Lett. 28, 245247 (1976).

17. Hachenberg, J. & Samwer, K. Indications for a slow beta-relaxation in a fragile metallic glass. J. Non-crystal. Solids 352, 51105113 (2006).

18. Hachenberg, J. et al. Merging of the alpha and beta relaxations and aging via the Johari-Goldstein modes in rapidly quenched metallic glasses. Appl. Phys. Lett. 92, 131911 (2008).

19. Zhao, Z. F., Wen, P., Shek, C. H. & Wang, W. H. Measurements of slow beta- relaxations in metallic glasses and supercooled liquids. Phys. Rev. B 75, 174201 (2007).

20. Hu, L. N. & Yue, Y. Z. Secondary relaxation behavior in a strong glass. J. Phys. Chem. B 112, 90539057 (2008).

21. Liu, X. F., Zhang, B., Wen, P. & Wang, W. H. The slow beta-relaxation observed in Ce-based bulk metallic glass-forming supercooled liquid.J. Non-Cryst. Solids 352, 40134016 (2006).22. Ngai, K. L. Correlation between the secondary beta-relaxation time at Tg with the Kohlrausch exponent of the primary alpha relaxation or the fragility of glass-forming materials. Phys. Rev. E 57, 73467349 (1998).

23. Yue, Y. Z., Christiansen, J.deC. & Jensen, S. L. Determination of the ctive temperature for a hyperquenched glass. Chem. Phys. Lett. 357, 2024 (2002).

24. Velikov, V., Borick, S. & Angell, C. A. The glass transition of water, based on hyperquenching experiments. Science 294, 23352338 (2001).

25. Busch, R., Bakke, E. & Johnson, W. L. Viscosity of the supercooled liquid and relaxation at the glass transition of the Zr46.75Ti8.25Cu7.5Ni10Be27.5 bulk metallic glass forming alloy. Acta Mater. 46, 47254732 (1998).

26. Ma, D. et al. Nearest-neighbor coordination and chemical ordering in multicomponent bulk metallic glasses. Appl. Phys. Lett. 90, 211908 (2007).

27. Hirata, A. et al. Geometric frustration of icosahedron in metallic glasses. Science

341, 376379 (2013).

28. Iwashita, T. & Egami, T. Atomic mechanism of ow in simple liquids under shear. Phys. Rev. Lett. 108, 196001 (2012).

29. Argon, A. S. Plastic deformation in metallic glasses. Acta Metall. 27, 4758 (1979).

30. Falk, M. L. & Langer, J. S. Dynamics of viscoplastic deformation in amorphous solids. Phys. Rev. E 57, 71927205 (1998).

31. Liu, Y. H. et al. Deposition of multicomponent metallic glass lms by single-target magnetron sputtering. Intermetallics 21, 105114 (2012).

32. Kohara, S. et al. Structural studies of disordered materials using high-energy X-ray diffraction from ambient to extreme conditions. J. Phys. Condens. Matter 19, 506101 (2007).

33. Stern, E. A., Newville, M., Ravel, B., Yacoby, Y. & Haskel, D. The UWXAFS analysis package-Philosophy and details. Physica B 208, 117120 (1995).

34. Ravel, B. & Newville, M. ATHENA, ARTEMIS, HEPHAESTUS: data analysis for X-ray absorption spectroscopy using IFEFFIT. J. Synchrotron Rad. 12, 537541 (2005).

Acknowledgements

We thank Professor Takeshi Egami and Professor Weihua Wang for stimulating discussions and the support of JASRI/SPring-8 under Proposals of 2010A1078 and 2010B1364. This work is sponsored by the Grant-in-Aid for Young Scientists and Grant-in-Aid-S, JSPS; National Natural Science Foundation of China (grant 51271113); and the State Key Laboratory for Advanced Metals and Materials (2013-ZD03).

Author contributions

Y.H.L. and M.W.C. conceived the research. Y.H.L. prepared the samples. Y.H.L. and D.P.B.A. performed DSC measurements. Y.H.L., D.P.B.A., T.F. and M.M. conducted synchrotron experiments. Y.H.L. and M.W.C. wrote the manuscript. M.W.C. supervised the project. All author contributed to data analysis and discussion.

Additional information

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsand

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsand

Web End =reprintsandpermissions/

How to cite this article: Liu, Y. H. et al. Structural origins of Johari-Goldstein relaxation in a metallic glass. Nat. Commun. 5:3238 doi: 10.1038/ncomms4238 (2014).

NATURE COMMUNICATIONS | 5:3238 | DOI: 10.1038/ncomms4238 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2014 Macmillan Publishers Limited. All rights reserved.