Correlation between Viscoelastic Moduli and Atomic Rearrangements

School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287, United States. § I. Physikalisches Institut, Universität Göttingen, D...
0 downloads 24 Views 592KB Size
Subscriber access provided by CORNELL UNIVERSITY LIBRARY

Letter

Correlation between Viscoelastic Moduli and Atomic Rearrangements in Metallic Glasses Hai-Bin Yu, Ranko Richert, and Konrad Samwer J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.6b01738 • Publication Date (Web): 08 Sep 2016 Downloaded from http://pubs.acs.org on September 9, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Correlation between Viscoelastic Moduli and Atomic Rearrangements in Metallic Glasses Hai-Bin Yu1*, Ranko Richert2 and Konrad Samwer3* 1

Wuhan National High Magnetic Field Center, Huazhong University of Science and Technology, WuHan 430074, China 2

School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287, USA 3

I. Physikalisches Institut, Universität Göttingen, D-37077 Göttingen, Germany

AUTHOR INFORMATION Corresponding Author To whom correspondence should be addressed to H.-B.Y. Email: [email protected] or K.S. Email: [email protected]

ABSTRACT. Dynamical moduli, such as storage and loss moduli, characterize the viscoelasticity of materials (i.e., time-dependent elasticity) and convey important information about the relaxation processes of glasses and supercooled liquids. A fundamental question is what ultimately determines them in glassy materials. Here, for several model metallic glasses, we demonstrate that both the storage and loss moduli are uniquely determined by the most probable

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 16

atomic non-affine displacements, regardless of temperature or frequency. Moreover, the fastmoving atoms (which contribute to dynamical heterogeneity) does not contribute explicitly to the moduli. Our findings provide a physical basis for the origin of viscoelasticity in metallic glasses.

TOC GRAPHICS

The elastic properties of many glassy materials, like polymers, metallic glasses (MGs)1-4 and colloids5, are time dependent over a wide range of temperatures (or volume fraction for colloids), which are in contrast to the well-known Hook elasticity and are called viscoelasticity6. To describe viscoelasticity, it is necessary to study the time-dependent dynamical moduli, such as storage and loss moduli as probed by dynamical mechanical experiments. Specifically, the storage modulus measures the stored energy, representing the elastic portion of the mechanical deformation; while the loss modulus measures the energy dissipated as heat, representing the viscous or the damping portion7,8. These global dynamical moduli are a major concern in many applications of glassy materials9 whose mechanical properties change with time. They also

ACS Paragon Plus Environment

2

Page 3 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

convey important information about the relaxation processes of glasses and supercooled liquids which are the central topic in glassy physics7,8,10-13. The infinite frequency shear moduli are the essential parameters in theoretical understanding the glass transition14 and the mechanisms of plastic deformation of glasses15. For crystalline materials, elasticity can be understood from the reversible atomic displacements on specific crystallographic planes (e.g., bond stretching). However, similar arguments cannot be safely applied to glassy materials. For example, in MGs, a recent study shows that even in the apparent elastic deformation regime, about 25% atoms undergo irreversible rearrangements16. The moduli of MGs are also found to be spatially inhomogeneous17,18. Even perplexed by the disordered atomic structure19-22 and the complex (and sometimes materials specific23,24) relaxation processes7,8,25, the question what ultimately determines the viscoelasticity and dynamical moduli is still not very clear in glassy materials26,27. In this work, we study the origin of viscoelasticity by investigating the global dynamical moduli of a model MG over a wide range of temperature and time scales, via an approach of molecular dynamics simulation of dynamical mechanical spectroscopy (MD-DMS) together with structural analysis based on the statistics of atomic displacements28-30. We find that the viscoelasticity is determined predominantly by the most probable atomic displacements. The influence of temperature and frequency on dynamical moduli can be understood in a unified way as they contribute to the atomic displacements. Interestingly, dynamical heterogeneity as induced by atoms whose displacements larger than the most probable one does not contribute explicitly to the global dynamical moduli. The details of our model system and the protocols of MD-DMS are presented in Methods. Briefly, at a temperature T, we apply a sinusoidal strain ε (t ) = ε A sin( 2πt / tω ) , with a period tω

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry Letters

(related to frequency f = 1/ tω) and a fixed strain amplitude εA = 0.7 %, along the x direction of a model Ni80P20 MG and the resulting stress σ(t) is measured. To study the viscoelasticity and the dynamical moduli, we varied tω from 100 to 10,000 picoseconds (ps) and T from 20 to 780 K, covering both glass and supercooled liquid states. The storage- (E’) and loss (E”) moduli are calculated based on ε(t) and σ(t). a

c

10

10

4

d

0

3 2

10000 100

1

40 20

400

600 T (K)

e E" (GPa)

9

10-1

6 3 0 200

100 300 1000 3000 10000

0

100

12

0

tω (ps)

60

E' (GPa)

100 300 1000 3000 10000

up (Å)

40

p (u)

E' (GPa)

10

101

tω (ps)

20

b

10

tω (ps)

60

E" (GPa)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 16

400

600

800

10-2 -2 10

10-1

100

101

0 12 8 4 0

u (Å)

T (K)

10-1

100

101

up (Å)

Figure 1 Dynamical modulus and statistics on atomic displacements. (a) Storage E’ and (b) loss modulus E” as a function of temperature T for different periods tω as indicated. (c) Probability density function p(u) of atomic displacements u for different T at tω = 100 ps, T ranges from 50 to 650K in a step of 50K (from left to right). The inset of (c) the peak position up of p(u) as a function of T for two typical tω as indicated. The storage modulus (d) and loss modulus (e) as a function of up.

Figure 1 (a) and (b) show the E’ and E” respectively as a function of T for five different tω. For low temperatures (e.g., T < 400 K) both E’ and E” are insensitive to tω and changes slowly with T, indicating that the system is deep in the glassy state. However, for higher T (e.g., T > 500 K), E’ and E” depends strongly on tω and changes rapidly with T, indicating the viscoelastic nature of the deformation. Specifically, at higher T, the values of E’ drop in a

ACS Paragon Plus Environment

4

Page 5 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

sigmoidal manner while those of E” develop asymmetrical peaks. Both of them signal the dynamics of the transition from a glassy state to a supercooled liquid state, i.e., the α relaxation. We note that all these features are qualitatively consistent with experimental dynamical mechanical measurements7,8, which confirm the validity and reliability of our MD-DMS in studying the viscoelasticity of MGs. To investigate the structural origin of the viscoelasticity and what determines the dynamical moduli, we conducted structural analysis based on the atomic displacements u,

r r defined as u (t ) = r (t ) − r (0) for each atom during a time interval of ∆t = tω for all the combinations of T and tω as reported in Fig. 1(a) and (b). This choice of ∆t is meant to avoid atomic displacements due to the overall deformations applied by the MD-DMS. Figure 1(c) shows a typical probability density function p(u) based on the statistics of u for different T (ranging from 50 to 650K), at a fixed tω =100 ps. Obviously, the peak position up of p(u) represents the most probable atomic displacements of all the atoms for the given tω. The inset of Fig.1(c) quantitatively summarizes how the up varies with T for two typical tω = 100 and 10,000 ps, respectively. We find that up has the T and tω dependence that qualitatively similar to those of E’ and E” as presented in Fig. 1(a) and (b). For example, the up changes slowly with T and insensitive to tω at low temperatures but it increases dramatically with T and becomes sensitive to tω at higher T. This observation implies that there might be connections between up and E’ and E”. To reveal such connections, Fig. 1(d) and (e) plot E’ and E” as a function of up for all the studied T and tω. Remarkably, we find that all the data of E’ and E” collapses onto two unique curves as shown in Fig. 1(d) and (e) respectively. This result establishes correlations between up and dynamical moduli (E’ and E”) in our model MG. The collapse of the data reminds the so-

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry Letters

called “temperature-time superposition” where the spectra are shifted manually to construct a master curve31. Our result, however, does not involve any adjusting parameters. It suggests that the factors such as T and tω that affects E’ and E” can be understood in a unified manner as they contribute to up. This notion is validated by Fig.1 (d) and (e), which show that for any combinations of T and tω, if they have the same value of up, they would have uniquely the same dynamical moduli. Thus, the two major seemingly independent variables, T and tω, can be reduced into a single variable, up.

a

up = 0.38 Å

100

I (T = 560 K, tω = 100 ps) II (T = 520 K, tω = 1000 ps) III (T = 500 K, tω = 10000 ps)

p(u)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 16

10-1 10-2

I

0

1

II

2

3

III

4

5

u (Å)

b

c

6 u (Å) ≥4

2

0

Figure 2

Dynamical heterogeneity. (a) Probability density function p(u) for different pairs of

combinations of T and tω; from left to right (T = 560 K, tω = 100 ps), (T = 520 K, tω = 1,000 ps) and (T = 500 K, tω = 10,000 ps), respectively. They have the same value of up = 0.38±0.01 Å. (b) and (c) 2D slices of the atomic configurations correspond to (T = 560 K, tω = 100 ps) and (T = 500 K, tω = 10,000 ps), respectively.

ACS Paragon Plus Environment

6

Page 7 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Given the above analysis, one immediate question is that how does atoms with different mobility (i.e., dynamical heterogeneity) contribute to E’ and E”? This question is particularly meaningful because the atomic motions in glasses and supercooled liquid are intrinsically heterogeneous32-38. To address this issue, Fig.2 (a) selectively plots three curves of p (u). They are the different pairs of combinations between T and tω [from left to right in Fig.2(a), these pairs are (T = 560 K, tω = 100 ps), (T = 520 K, tω = 1000 ps) and (T = 500 K, tω = 10000 ps), respectively] but they have the same value of up = 0.38±0.01 Å and therefore the same global values of E’ = 33±1 GPa and E” = 5.6±0.5 GPa, according to Fig.1(d) and (e). We find that they show appreciable difference at the tails of p(u) when u > 1.5 Å. Especially, the sub-peaks of p(u) around about u = 2.5 Å grow higher in the order of decreasing T (or increasing tω), which suggest the increases in dynamical heterogeneity and the fraction of atoms with atomic displacements larger than up. Figure 2(b) and (c) compare two atomic configurations for (T = 560 K, tω = 100 ps) and (T = 500 K, tω = 10000 ps), respectively. The color code represents the magnitudes of u. Heterogeneous dynamics are observed in both configurations. Compared to Fig. 2(b), Fig. 2(c) is found to be more heterogeneous over a wider range of length scales, and it involves more particles with larger u than up. Additionally, the regions with heterogeneous dynamics are observed to be string-like structures in Fig. 2(b) but compact structures in Fig.2(c). Overall, the above results suggest that for some specified combinations of T and tω, the model MGs might have the same values of up and E’ and E”, however, their dynamical heterogeneity and the portion of fast-moving atoms (those move farther than up ) are not necessarily the same. In other words, dynamical heterogeneity as induced by the faster atoms does not contribute explicitly to the dynamical moduli.

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry Letters

E' (GPa)

a

75

50

tω (ps) 100 300 1000 3000 10000

25

0

b 12

E" (GPa)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 16

9 6 3 0 10-1

100

101

1/2 (Å)

Figure 3 Dynamical modulus and root mean square displacements. (a) Storage E’ and (b) loss E” moduli as a function of root mean square displacements 1/2 for different tω as indicated. The shadow region shows that E’ or E” do not collapse.

To provide more quantitative support for the above discussions, Fig. 3 reports another way to analyze the data. Instead of using up from statistical analysis, we calculated the root mean square atomic displacements 1/2 based on the u values, where < u 2 >= 1/ N ∑ ui2 is the mean square i

displacements and N = 32,000 is the number of atoms. Figure 3 (a) and (b) show E’ and E” as a function of 1/2, respectively. In contrast to Fig.1 (d) and (e), most of the data with different tω does not collapse onto unique curves [e.g., in the shadow region through Fig. 3(a) and (b)]. This suggests (i) that dynamical moduli could not be single variable functions of 1/2 and (ii) that only up fundamentally correlates with the viscoelasticity and the dynamical moduli. Since

ACS Paragon Plus Environment

8

Page 9 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

1/2 incorporates all the atomic displacements, including faster atoms with u > up, this result reinforces the above conclusion that faster atoms does not contribute to the viscoelasticity and dynamical moduli. By the way, a closer examine of the data in Fig.3 (a) and (b) reveals that the initial portion of the E’ and E” data indeed collapse when 1/2 is smaller than 0.5 Å. This is because at the conditions of small 1/2 (usually at low T or small tω), p(u) is narrowly distributed around up and therefore up ≈ 1/2. For other conditions up does not relate to 1/2 and lacks correlations between 1/2 and the dynamical moduli. We have also considered the mean atomic displacement < u >= 1 / N ∑ ui . It cannot collapse the data too, because of the same i

reason. In a colloidal glassy system39, Conrad et al identified that only the slow-moving particles contribute significantly to the dynamical shear moduli. This is consistent with the present findings, because in our model MG, atoms with u ~ up are indeed moving slowly, as can be seen in Fig.2. Therefore, the viscoelasticity of glassy materials and supercooled liquids could be considered as a consequence of the collective rearrangements of the atoms of the whole system. Possibly, the embedded fast-moving atoms are liquid like (with u >> up and E’ ~ E” ~ 0) and cannot restrain the stress and thus they do not contribute to the development of the viscoelastic modulus of a glass. We note that our findings are not limited to the Ni-P MG. In supporting online materials, Fig.S1-S3, we demonstrate similar conclusions hold for a Pd80Si20 MG (all the other conditions are same with the present Ni-P MG) as well. Additionally, our results are independent of the strain amplitudes (εA) conducted in the linear response regime (εA < 2%), Fig. S4. For the nonlinear response regime, plastic deformation is involved and the absolute values of E’ and E”

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 16

are difficult to obtain with accuracy (i.e., one has to consider multiple Fourier components). Since the dynamical moduli in experiments are measured in the linear response regime, we do not discuss the nonlinear issue here. Finally, the system size of the simulation affect the results very little, as shown in Fig. S5. One interesting extension of our work could be to study model MGs with secondary relaxations (e.g. the β relaxations), and to see if the relations between E’, E’’ ~ up remain valid. Such a study could provide information about the mechanism of secondary relaxations8,9,23,40,41. In summary, we have shown that both of the storage and loss moduli in a model metallic glass and its supercooled liquid are uniquely determined by the most probable atomic non-affine displacements. The temperature and time that affect dynamical moduli can be understood in a similar way as they contribute to the atomic displacements. Our findings bridge the gap between the dynamical mechanical properties of metallic glass and the atomic motions. Hopefully, they could offer a new benchmark for constructing viscoelasticity theory of glassy materials in the future.

SIMULATION METHODS An open source LAMMPS package42 was used for all the MD simulations. The model system contains N = 32000 atoms with the composition Ni80P20 and the constituting atoms are interacting via an embedded atom method (EAM) potentials20. For the sample preparations, the system was melted and equilibrated at T = 2000 K for 100ps, then cooled down to T = 1 K with a cooling rate of 1010 K/s, during which the cell-sizes were adjusted to give zero pressure with an NPT ensemble. Periodic boundary conditions were applied for all the calculations.

ACS Paragon Plus Environment

10

Page 11 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

As in the case of experimental DMS, we apply a sinusoidal strain ε (t ) = ε A sin( 2πt / tω ) along the x direction of the model MG, where tω is the period and is selected as 100, 300, 1,000, 3,000 and 10,000 ps, while εA is fixed at about 0.7 % to ensure all the deformations are in the linear response regime. For each MD-DMS, 10 full cycles were used, i.e., t in the range [0, 10tω]. We fitted the resulting stress as: σ (t ) = σ 0 + σ A sin(2π t / tω + δ ) where σ0 is a constant term and usually small (σ0 < 0.1σA in the glassy state), δ the phase difference between stress and strain. Storage (E’) and loss (E”) moduli are calculated according to E ' = σ A / ε A sin(δ ) and E " = σ A / ε A cos(δ ) , respectively. The MD-DMS were carried out during the cooling processes

of the sample-preparations, and NVT ensemble was applied during the cyclic deformations. The probability density function p(u) is defined as p(u ) = [ P(u + ∆u) − P(u )] / ∆u , where P(u) is the cumulative distribution that quantifies the probability of finding X ≤ u. We employed ∆u = 0.01 Å for all the calculations. The probability density is normalized according ∞

to ∫ p (u )du = P(∞) = 1 . 0

ACKNOWLEDGMENT We thank Prof. Yue Wu for discussions. H.-B.Y. thanks the support from Huazhong University of Science and Technology of China. K.S. acknowledges the supports from the German Science Foundation within the FOR 1394, P1. The computational work was carried out at Gesellschaft für wissenschaftliche Datenverarbeitung, Göttingen (GWDG), Germany and the TianHe-1(A) of National Supercomputer Center in Tianjin, China.

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 16

Supporting Information Available: We have considered other factors, such as the effects of different model MGs, the strain amplitude and the size of the system under simulation. They are shown in the supporting information.

REFERENCES (1)

Angell, C. A. Formation of glasses from liquids and biopolymers. Science 1995, 267,

1924-1935. (2)

Chen, M. W. A brief overview of bulk metallic glasses. NPG Asia Mater. 2011, 3, 82-90.

(3)

Wang, W. H.; Dong, C.; Shek, C. H. Bulk metallic glasses. Mater. Sci. Eng. R 2004, 44

45-89. (4)

Inoue, A. Stabilization of metallic supercooled liquid and bulk amorphous alloys. Acta

Mater. 2000, 48, 279-306. (5)

Hunter, G. L.; Weeks, E. R. The physics of the colloidal glass transition. Rep. Prog.

Phys. 2012, 75, 066501. (6)

Wang, W. H. The elastic properties, elastic models and elastic perspectives of metallic

glasses. Prog. Mater. Sci. 2012, 57, 487-656. (7)

Yu, H.-B.; Wang, W.-H.; Samwer, K. The beta relaxation in metallic glasses: an

overview. Mater. Today 2013, 16, 183-191. (8)

Yu, H. B.; Wang, W. H.; Bai, H. Y.; Samwer, K. The beta-relaxation in metallic glasses.

Natl. Sci. Rev. 2014, 1, 429-461.

ACS Paragon Plus Environment

12

Page 13 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(9)

Jancar, J.; Douglas, J. F.; Starr, F. W.; Kumar, S. K.; Cassagnau, P.; Lesser, A. J.;

Sternstein, S. S.; Buehler, M. J. Current issues in research on structure-property relationships in polymer nanocomposites. Polymer 2010, 51, 3321-3343. (10) Angell, C. A.; Ngai, K. L.; McKenna, G. B.; McMillan, P. F.; Martin, S. W. Relaxation in glassforming liquids and amorphous solids. J. Appl. Phys. 2000, 88, 3113-3157. (11) Maggi, C.; Jakobsen, B.; Christensen, T.; Olsen, N. B.; Dyre, J. C. Supercooled Liquid Dynamics Studied via Shear-Mechanical Spectroscopy. J. Phys. Chem. B 2008, 112, 1632016325. (12) Qiao, J. C.; Pelletier, J. M. Dynamic Mechanical Relaxation in Bulk Metallic Glasses: A Review. J. Mater. Sci. Tech. 2014, 30, 523-545. (13) Ediger, M. D.; Harrowell, P. Perspective: Supercooled liquids and glasses. J. Chem. Phys. 2012, 137, 4747326. (14) Dyre, J. C. Colloquium: The glass transition and elastic models of glass-forming liquids. Rev. Mod. Phys. 2006, 78, 953-972. (15) Johnson, W. L.; Samwer, K. A universal criterion for plastic yielding of metallic glasses with a (T/Tg)^(2/3) temperature dependence. Phys. Rev. Lett. 2005, 95, 195501. (16) Dmowski, W.; Iwashita, T.; Chuang, C. P.; Almer, J.; Egami, T. Elastic Heterogeneity in Metallic Glasses. Phys. Rev. Lett. 2010, 105, 205502. (17) Wagner, H.; Bedorf, D.; Kuechemann, S.; Schwabe, M.; Zhang, B.; Arnold, W.; Samwer, K. Local elastic properties of a metallic glass. Nat. Mater. 2011, 10, 439-442.

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 16

(18) Liu, Y. H.; Wang, D.; Nakajima, K.; Zhang, W.; Hirata, A.; Nishi, T.; Inoue, A.; Chen, M. W. Characterization of Nanoscale Mechanical Heterogeneity in a Metallic Glass by Dynamic Force Microscopy. Phys. Rev. Lett. 2011, 106, 125504. (19) Miracle, D. B.; Egami, T.; Flores, K. M.; Kelton, K. F. Structural aspects of metallic glasses. MRS Bull. 2007, 32, 629-634. (20) Sheng, H. W.; Luo, W. K.; Alamgir, F. M.; Bai, J. M.; Ma, E. Atomic packing and shortto-medium-range order in metallic glasses. Nature 2006, 439, 419-425. (21) Cheng, Y. Q.; Ma, E. Atomic-level structure and structure-property relationship in metallic glasses. Prog. Mater. Sci. 2011, 56, 379-473. (22) Ma, E. Tuning order in disorder. Nat. Mater. 2015, 14, 547-552. (23) Yu, H. B.; Samwer, K.; Wang, W. H.; Bai, H. Y. Chemical influence on beta-relaxations and the formation of molecule-like metallic glasses. Nat. Coummun. 2013, 4, 2204. (24) Ngai, K. L.; Paluch, M. Classification of secondary relaxation in glass-formers based on dynamic properties. J. Chem. Phys. 2004, 120, 857-873. (25) Lunkenheimer, P.; Schneider, U.; Brand, R.; Loidl, A. Glassy dynamics. Contemp. Phys. 2000, 41, 15-36. (26) Egami, T. Understanding the Properties and Structure of Metallic Glasses at the Atomic Level. JOM 2010, 62, 70-75. (27) Dubey, A. K.; Procaccia, I.; Shor, C. A. B. Z.; Singh, M. Elasticity in Amorphous Solids: Nonlinear or Piecewise Linear? Phys. Rev. Lett. 2016, 116, 085502.

ACS Paragon Plus Environment

14

Page 15 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(28) Yu, H. B.; Richert, R.; Maass, R.; Samwer, K. Unified Criterion for TemperatureInduced and Strain-Driven Glass Transitions in Metallic Glass. Phys. Rev. Lett. 2015, 115, 135701. (29) Yu, H.-B.; Richert, R.; Maass, R.; Samwer, K. Strain induced fragility transition in metallic glass. Nat. Coummun. 2015, 6, 7179. (30) Yu, H.-B.; Samwer, K. Atomic mechanism of internal friction in a model metallic glass. Phys. Rev. B 2014, 90, 144201. (31) Olsen, N. B.; Christensen, T.; Dyre, J. C. Time-Temperature Superposition in Viscous Liquids. Phys. Rev. Lett. 2001, 86, 1271-1274. (32) Kob, W.; Donati, C.; Plimpton, S. J.; Poole, P. H.; Glotzer, S. C. Dynamical heterogeneities in a supercooled Lennard-Jones liquid. Phys. Rev. Lett. 1997, 79, 2827-2830. (33) Kegel, W. K.; van Blaaderen, A. Direct observation of dynamical heterogeneities in colloidal hard-sphere suspensions. Science 2000, 287, 290-293. (34) Garrahan, J. P.; Chandler, D. Geometrical explanation and scaling of dynamical heterogeneities in glass forming systems. Phys. Rev. Lett. 2002, 89, 035704. (35) Dauchot, O.; Marty, G.; Biroli, G. Dynamical heterogeneity close to the jamming transition in a sheared granular material. Phys. Rev. Lett. 2005, 95, 265701. (36) Candelier, R.; Dauchot, O.; Biroli, G. Building Blocks of Dynamical Heterogeneities in Dense Granular Media. Phys. Rev. Lett. 2009, 102, 088001.

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 16

(37) Ediger, M. D. Spatially heterogeneous dynamics in supercooled liquids. Ann. Rev. Phys. Chem. 2000, 51, 99-128. (38) Wisitsorasak, A.; Wolynes, P. G. Dynamical Heterogeneity of the Glassy State. J. Phys. Chem. B 2014, 118, 7835-7847. (39) Conrad, J. C.; Dhillon, P. P.; Weeks, E. R.; Reichman, D. R.; Weitz, D. A. Contribution of Slow Clusters to the Bulk Elasticity Near the Colloidal Glass Transition. Phys. Rev. Lett. 2006, 97, 265701. (40) Zhu, Z. G.; Li, Y. Z.; Wang, Z.; Gao, X. Q.; Wen, P.; Bai, H. Y.; Ngai, K. L.; Wang, W. H. Compositional origin of unusual beta-relaxation properties in La-Ni-Al metallic glasses. J. Chem. Phys. 2014, 141, 084506. (41) Aji, D. P. B.; Johari, G. P. Kinetic-freezing and unfreezing of local-region fluctuations in a glass structure observed by heat capacity hysteresis. J. Chem. Phys. 2015, 142, 214501. (42) Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comp. Phys. 1995, 117, 1-19.

ACS Paragon Plus Environment

16