Abstract
The glass transition remains unclarified in condensed matter physics. Investigating the mechanical properties of glass is challenging because any global deformation that might result in shear rejuvenation would require a prohibitively long relaxation time. Moreover, glass is well known to be heterogeneous, and a global perturbation would prevent exploration of local mechanical/transport properties. However, investigation based on a local probe, i.e., microrheology, may overcome these problems. Here, we establish active microrheology of a bulk metallic glass, via a probe particle driven into host medium glass. This technique is amenable to experimental investigations via nanoindentation tests. We provide distinct evidence of a strong relationship between the microscopic dynamics of the probe particle and the macroscopic properties of the host medium glass. These findings establish active microrheology as a promising technique for investigating the local properties of bulk metallic glass.
INTRODUCTION
When a glass-forming liquid is rapidly cooled below its glass transition temperature, it falls out of equilibrium, and its viscosity, η, increases by more than 15 orders of magnitude (1, 2). Notably, the structural relaxation time of a glass, τα ∝ η, exceeds any other time scale in the realm of condensed matter physics and is comparable to cosmological time scales (3). The origin of this slowing down remains the primary challenge in the physics of glass transition, i.e., whether the long time scale originates from collective behavior (a thermodynamic description) (4), similar to critical phenomena, or the slowing down occurs at the microscopic level (a kinetic origin). A widely accepted picture of the slow progress of a supercooled liquid toward equilibrium is based on a trap-escape behavior: A particle is transiently trapped in a cage created by its neighbors; over a short time, τβ, the particle rattles inside the cage and escapes only after a long time, τα. This time scale diverges at the glass transition temperature, at which the cages become permanent. This interpretation is well supported by the characteristic two-step relaxation of a supercooled liquid, as observed in neutron scattering experiments and simulations (5). Several competing theories have aimed at explaining the origin of the slowing down of a glass. Among them, the mode-coupling theory (MCT) has successfully described power-law β-relaxation, stretched exponentials of α-relaxation, and the time-temperature superposition principle (6). An intriguing discovery in experiments (7, 8) and simulations (9, 10) has motivated a thermodynamic description of the glass transition: Localized particles, unable to diffuse away, form mesoscopic clusters whose typical size increases upon cooling (11). These mesoscopic structures, known as dynamic heterogeneities (DHs), whose origin is yet another controversy (static versus dynamic origin), suggest that a supercooled liquid can be treated as a complex fluid: An enlightening work by Furukawa and Tanaka (12) established the complex fluid view for a supercooled liquid. This interpretation contrasts with the original MCT, which was an extension of the simple liquid theory, although recent developments of the MCT have been reconciled with DHs (13).
Despite intense research in recent decades on the underlying origins of the glass transition, a consensus has not yet been established. Furthermore, understanding of the mechanical properties of glasses, which are investigated by shear deformation, is also elusive. The reason for this lack of consensus is fourfold: (i) Any global perturbation requires an extremely long time scale for relaxation; (ii) shear forces may rejuvenate the glass, thus relocating the glass to shallower energy minima (14, 15) and rendering the glass more prone to fluctuations; (iii) a realistic deformation rate for an atomistic glass is a technical impossibility in computer simulations (16, 17); and (iv), because of DHs, a glass is heterogeneous, and a global deformation allows for exploration of only an average mechanical response. Here, we overcome these problems by investigating the mechanical properties of a glass at the microscopic level. With this technique, known as active microrheology, a singled-out probe particle is pulled through a dense host liquid (18). Investigation of the probe particle dynamics reveals nonequilibrium properties of the host medium, such as viscosity.
RESULTS
Microrheology
We perform active microrheology of a probe particle with a constant velocity U. Density map snapshots of the system are given in fig. S1. The size of the probe particle mimics that of a typical length scale of DHs. We obtain the response of the host medium by measuring the force acting on the particle, i.e., the friction,
(A) Friction, F, versus the probe velocity, U, for various temperatures across the glass transition temperature. Each point is an ensemble average over 500 independent realizations. (B) Successive slope of the velocity-force curves
Scaling ansatz
To fully capture the nonlinear behavior of the velocity-force curves, we seek a scaling ansatz reminiscent of that for magnetization in an Ising model. In this approach, the friction, F, depends only on a scaling variable x = U/∣T − Tmc∣Δ via
Because the critical temperature Tmc is unknown, we fit Eq. 2 to the data in the critical region, 0.300 < T < 0.375. The resulting lines are depicted by solid curves in Fig. 1B, and we display the values of the fit parameters in Table 1.
The magnitude of the prefactor in Eq. 2 is crucial in this analysis because it determines the reliability of the fit. For high-quality data,
Equation 1 suggests that the inherent friction, Finh = F(U → 0, T), in the system scales as Finh ∝∣T − Tmc∣Γ. The numerical value of the exponent Γ within the measurement error bars is equal to the exponent β = 0.41 − 0.42 of the percolation theory in 3D. The latter describes the scaling of the order parameter P∞, the probability that a site belongs to the infinite cluster; it is zero below the critical point, and it scales as a function of the distance from the critical point P∞ ∝ (p − pc)β. This may indicate that the inherent friction is generated by percolation of immobile particles (DHs) near the glass transition point. To examine this possibility, a systematic percolation analysis would be required. Such an investigation is beyond the scope of this work, but we will address it in the near future.
Relaxations
We now focus on the properties of the time series of the force acting on the probe particle F = {F(t)}, and we investigate whether the structural relaxation process can be traced via a friction time series. In fig. S2 (A and B), we display the friction time series at T = 0.150 and 0.466, respectively, for a probe velocity of U = 0.0379269. As shown in fig. S2B, in the supercooled state, the friction changes rapidly with time; in contrast, in fig. S2A, corresponding to the glass phase, the friction exhibits a strong correlation. The time interval over which the friction is correlated is marked by an arrow as a guide to the eye (the shaded area). This apparent correlation may be related to the formation of cages. To investigate this possibility, in Fig. 2 (A and C), we plot the autocorrelation function for friction, CF(τ) = 〈F(τ)F(0)〉, for T = 0.200 and 0.466, respectively. 〈⋯〉 presents an ensemble average over 500 independent realizations (in fig. S3, we show a similar plot for T = 0.350). Each curve begins with an initial relaxation, which corresponds to the ballistic motion of particles at very small time scales. All the different curves with different U values are superimposed for a short time scale, thus indicating that the short-term dynamics is unaffected by the probe motion. This behavior is analogous to that of sheared liquids at very small time scales, which is similarly unaffected by shear (28). The curves level off with a dip followed by a damped oscillation and eventually a plateau. The dip and the damped oscillation correspond to the oscillatory nature of the initial system preparation, owing to shear annealing (29). The damped oscillation of the correlation function is known as quenched echoes (30). Kob and Barrat (31) have reported such oscillations for Lennard-Jones (LJ) glasses annealed with oscillatory protocols. However, Shiba et al. (32) have shown similar oscillations originating from the propagation of acoustic phonons. Kob and Barrat ruled out the possibility of this dip being related to the boson peak. After the damped oscillations, a plateau with a length of several orders of magnitude appears. This trend is akin to the β-relaxation process (5). After this rattling period, a long relaxation process terminates the curves. This final long relaxation is akin to an α-relaxation process, corresponding to the escape of a particle from a cage. The behavior of the system both below and above the glass transition temperature is consistent with that in shear-driven structural glasses (28, 33). Therefore, the relaxation of friction is described by a two-step relaxation process, thus revealing a strong relationship between dynamics of the probe particle and the two-step structural relaxation of the host medium glass. In the remainder of the paper, we provide more evidence of this intimate relationship.
The autocorrelation function CF(τ) = 〈F(0)F(τ)〉 of the time series of the force acting on the probe particle F versus the time window τ for T = 0.200 and 0.466 is depicted in (A) and (C), respectively. Each curve represents an average of 500 realizations. Each curve has an initial relaxation, which corresponds to ballistic transport at very small time scales. This fast relaxation is followed by a plateau that is akin to a β-process (rattling in cages). Each plateau starts with a damped oscillation, known as a quenched echo (30, 31), which corresponds to the oscillatory nature of the initial annealing of the system or phonon propagation (32). The plateau terminates with a long relaxation process akin to an α-process (escape). We plot CF(x) versus x = τ × Uν, in which ν = −(q − 1) = 0.95 ± 0.1 is the force thinning exponent in (B) and (D). The data collapse onto a compressed exponential function of the form Eq. 3 [dashed lines in (B) and (D)]. The collapse is akin to the time-shear superposition principle in spin glass and MCT (1, 33). The compressed (faster) exponential is a manifestation of the super-diffusive motion of the probe particle, corresponding to super-diffusive relaxation. In (E), the compressed exponent b is depicted versus the probe velocity for various values of T.
Time-force superposition
In Fig. 2 (B and D), CF(τ) is plotted versus x = τ × Uν, in which U is the velocity of the probe particle and ν = −(q − 1) = 0.95 ± 0.1 is the force thinning exponent (q is the critical exponent of the scaling ansatz in Eq. 1); here, an excellent collapse of the data into a master curve is obtained. The collapse of the correlation function is a time-force superposition, which is reminiscent of the “time-shear superposition principle” predicted by both the MCT (1) and spin glass theory (33). The force thinning exponent ν is consistent, within the error bars, with the shear thinning exponent reported by Furukawa et al. (12), ν = 0.8, and is very close to the values reported by Berthier and Barrat (16), ν = 2/3, and Varnik (28), ν = 1. Accordingly, our force thinning exponent, ν, and the critical exponent, q = 1 − ν, within the error bars are consistent with findings reported in the literature for macrorheology.
Compressed exponential
The master function (dashed lines in Fig. 2, B and D) is given by a compressed exponential function of the following form (dashed lines in Fig. 2, B and D)
Power spectral density
In Fig. 3, we display the power spectral density (PSD) of the friction time series F = {F(t)}, defined as
Absolute value of the Fourier transform of the autocorrelation function, PSD(ν) = ∣∫ CF (τ)e−2πiντdτ∣, the PSD, versus frequency ν for various T. The probe velocity is U = 10−3. At low frequency, a temperature-dependent scaling behavior is observed, owing to the temperature dependence of the compressed exponent b = b(T) of the autocorrelation function CF (τ) = exp (−(τ/τα)b). At very low T, the dependence converges to the well-known 1/ν colored-noise behavior because of the correlation in the data. As might be expected, the exponent of PSD decreases with increasing temperature. At very high T, the friction must be completely random, and therefore, the corresponding PSD must be a horizontal line of white noise, i.e., θ = 0.
Macrorheology
A conventional method of exploring the mechanical properties of supercooled liquids is macrorheology. We now compare the results obtained from active microrheology and macrorheology. We performed simulations of a simple shear applied to the host medium by using the Lees-Edwards boundary conditions. The resulting flow curves are shown in Fig. 4A. For most of the flow curves at low temperature, the shear stress saturates to a shear rate–independent regime at low
(A) Shear stress, σxy, versus the shear rate,
Bridge between micro- and macrorheologies
Creating a link between macro- and microrheology is crucial. To this end, we seek a direct comparison of microviscosity (ηm) and viscosity (η), defined as
A similar scaling ansatz can also be written for viscosity if ηm, Tmc, and U are replaced by η, Tc, and
To bridge the gap between micro- and macrorheology, we examine the scaling collapse of microviscosity (ηm bold symbols) and viscosity (η) [defined in Eq. 4]. We define an effective shear rate as
Distribution of microviscosity
We now demonstrate the superiority of microrheology to macrorheology, the latter of which yields only an average response. So far, we have demonstrated that the macrorheology can be attained by microrheology when the local properties, such as microviscosity, are spatially and temporally averaged. However, microrheology is superior to macrorheology in that it can probe the local properties of the system, such as local microviscosity and the local elastic modulus. To demonstrate this capability, in Fig. 6, using time series of microviscosity ηm, we plot the probability distribution function (PDF) of ηm. The horizontal axis is
Main: Probability distribution function of local microviscosity ηm, defined by Eq. 4, versus
CONCLUSION AND DISCUSSION
Here, we investigated the dynamics of a probe particle driven into a host glass-forming liquid. The friction acting on the probe exhibits a second-order phase transition. The microscopic origin of friction is one of the fundamental questions in tribology (37). Our results establish a nontrivial correspondence between the emergence of friction as a result of the approach to glass transition. We discovered a two-step relaxation process, a characteristic of the structural relaxation of a supercooled liquid, in the time autocorrelation function of the friction. Moreover, the autocorrelation function exhibits a time-force superposition reminiscent of the time-shear superposition principle in MCT (1) and spin glass theory (33). This finding strengthens the relationship between the dynamics of the probe particle and the structural relaxations of the host medium, on the basis of existing theoretical lines of argument. We further demonstrated that the autocorrelation function of friction is a compressed exponential, thus shedding light on a recent controversy regarding the origin of compressed exponential relaxation in active matter (34, 35). We compared the results obtained by active microrheology and macrorheology, revealing that spatiotemporally averaged microviscosity becomes equivalent to viscosity, thus providing clear evidence of the consistency of micro-macro rheologies. However, microrheology becomes superior to macrorheology because the former can provide the spatial distribution of microviscosity. We found that the skewness of PDF of microviscosity becomes larger as the temperature is decreased, a direct sign to the growth of heterogeneities at low temperatures.
The delocalization transition of the probe particle investigated by active microrheology has many common features with the so-called depinning transition. The depinning transition can be used to investigate collective behavior in many nonequilibrium systems. Specifically, it examines interacting particles driven over heterogeneous substrates, and it has abundant applications in vortex lattices in type II superconductors and driven colloids over optical periodic substrates. This topic has parallels with nonequilibrium phase transitions and directed percolation. More details can be found in an excellent recent review by Reichhardt and Reichhardt (38).
The following applications can be expected in microrheology studies of glasses or supercooled liquids. The spatial heterogeneities of the elastic modulus and viscosity can be observed in the systems. A previous study (39) has shown that elastic moduli are expected to be spatially heterogeneous in glass because the response on the strain fields for weak macroscopic flow becomes heterogeneous. In another study on supercooled liquids (12), a characteristic length scale has been observed from the wave number dependence of Fourier transformed viscosity. However, no studies have directly obtained the spatial distribution of the elastic modulus and viscosity. In the present study, we statistically confirmed the heterogeneous distribution of these mechanical properties. In the future, we expect to obtain the spatial distribution of these quantities. The heterogeneous distribution of these static mechanical properties suggests a relationship with the DH characterized by the particle-displacement field. Our analysis should be essential to understanding the origin of the DHs.
The low-temperature properties of glasses have recently been the subject of intensive investigations. Originally known as the Gardner transition in spin glass theory, the free energy landscape of a glass has been proposed to become fractal-like at low temperatures: A local minimum divides into many metabasins in the free energy landscape. As a result of this proliferation of the metabasins, local cages divide into several smaller cages upon cooling (or compression) (40). Therefore, the Gardner transition is associated with a change in the local morphology. As a result, a global deformation would not only hinder the exploration of such delicate local morphological changes but also might result in a loss of local information, thus potentially explaining why the existence of the Gardner transition at finite spatial dimensions has been debated. We believe that active microrheology provides a strong experimental and computational tool for investigating the low-temperature properties of glasses and, eventually, shedding light on the controversy regarding the existence of the Gardner transition at finite dimensions. Our work provides a basis for such explorations.
METHODS
Simulations
We performed large-scale 3D simulations of a 80 : 20 binary mixture (Kob-Andersen model) interacting via a smoothed LJ potential
Preparation of initial configurations
We prepared the initial configurations in three steps. (i) Particles are randomly inserted in the simulation box. (ii) A cosine soft pair potential of the form Uij(r) = A[1 + cos (πr/rc)] is applied to all pairs of interacting particles within the cutoff, rc, regardless of particle type. We then integrate the equations of motion of the particles by using the microcanonical ensemble (NVE) integration scheme followed by canonical ensemble (NVT) at T = 0.466ϵ/kB. This step is crucial to minimize the overlap between particles. After minimization, the soft cosine potential is replaced with the smoothed LJ potential given by Eq. 6. (iii) The system is annealed with cyclic-shearing deformation (15, 29) to obtain a stable glass less prone to fluctuations. The SLLOD algorithm was used to perform nonequilibrium molecular dynamics (MD) simulations with the Lees-Edwards boundary conditions.
Microrheology
We use a singled-out probe particle to perform microrheology. The probe particle is activated by a constant velocity U to mimic constant-velocity active microrheology (18). The probe particle interacts with the following expanded and shifted LJ pair potential
Unless otherwise stated, Δ = 2.0σ and rc = 21/6σ. Accordingly, the probe particle is fivefold larger than the A-type particle. This is a typical length scale for DHs in supercooled liquids (12). At t = 0, the probe particle is placed at the boundary, and after activation, the probe particle moves into the host at the prescribed velocity. We record the force acting on the probe particle in the direction of the velocity, F = −F ⋅ n where n = U/∣U∣. Here, F = ∑jfj, where fj is the force exerted by a neighboring particle, and the sum spans all neighbors of the probe particle. The force is averaged over time for 500 independent realizations. The total CPU time for the simulations is approximately 2 months for a cluster of 1656 cores. We gathered a dataset of 20 TB for all the simulations.
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/29/eaba8766/DC1
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).