Next Article in Journal
The Biological Effects of Complete Gasoline Engine Emissions Exposure in a 3D Human Airway Model (MucilAirTM) and in Human Bronchial Epithelial Cells (BEAS-2B)
Previous Article in Journal
Epidermolysis Bullosa-Associated Squamous Cell Carcinoma: From Pathogenesis to Therapeutic Perspectives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Fast Vibrational Modes and Slow Heterogeneous Dynamics in Polymers and Viscous Liquids

1
Dipartimento di Fisica “Enrico Fermi”, Università di Pisa, Largo B. Pontecorvo 3, I-56127 Pisa, Italy
2
Istituto per i Processi Chimico-Fisici-Consiglio Nazionale delle Ricerche (IPCF-CNR), via G. Moruzzi 1, I-56124 Pisa, Italy
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(22), 5708; https://doi.org/10.3390/ijms20225708
Submission received: 25 October 2019 / Revised: 7 November 2019 / Accepted: 8 November 2019 / Published: 14 November 2019
(This article belongs to the Section Macromolecules)

Abstract

:
Many systems, including polymers and molecular liquids, when adequately cooled and/or compressed, solidify into a disordered solid, i.e., a glass. The transition is not abrupt, featuring progressive decrease of the microscopic mobility and huge slowing down of the relaxation. A distinctive aspect of glass-forming materials is the microscopic dynamical heterogeneity (DH), i.e., the presence of regions with almost immobile particles coexisting with others where highly mobile ones are located. Following the first compelling evidence of a strong correlation between vibrational dynamics and ultraslow relaxation, we posed the question if the vibrational dynamics encodes predictive information on DH. Here, we review our results, drawn from molecular-dynamics numerical simulation of polymeric and molecular glass-formers, with a special focus on both the breakdown of the Stokes–Einstein relation between diffusion and viscosity, and the size of the regions with correlated displacements.

1. Introduction

When polymers, liquids, biomaterials, metals and molten salts are cooled or compressed, if the crystallization is avoided, they freeze into a microscopically disordered solid-like state, a glass [1,2,3]. On approaching the glass transition from states with high fluidity, the viscosity exhibits a huge increase of more than 10 orders of magnitude [1,2], along with the parallel decrease of the diffusivity [3,4]. Correspondingly, at microscopic level, solid-like behaviour becomes apparent, e.g., it is observed that a particle spends increasing time within the cage formed by the first neighbours where it rattles with amplitude u 2 1 / 2 on picosecond time scales [5]. This temporary trapping is rather persistent and the particle has average escape time, the structural relaxation time τ α , which increases from a few picoseconds in the low-viscosity liquid up to thousands of seconds close to the glass transition [6]. The quantity u 2 appears in the expression of the Debye–Waller (DW) factor, which, assuming harmonicity and isotropy of the thermal motion, takes the form exp q 2 u 2 / 3 , where q is the absolute value of the scattering vector [7]. Researchers investigating the cage motion in viscous liquids usually refer, as a metonym, to u 2 as the DW factor too, e.g., see the work in [8,9,10]. To keep maintain similarity with this literature, the same convention is adopted here.
The transition from a liquid to a glass is accompanied by the growth of transient domains which exhibit different mobility, e.g., see Figure 1. The phenomenon is usually dubbed “dynamical heterogeneity” (DH) and has been extensively studied, e.g., see the reviews in [4,6,11,12] and topical papers [13,14,15,16]. The size of the domains is relatively small involving approximately 10 molecule diameters [11], corresponding to a few nanometres [14]. On a more general ground, the size of DH domains is strictly related to the possible presence of characteristic length scales in glass-forming systems. Starting with the seminal paper by Adam and Gibbs, who invoked the presence of “cooperatively rearranging regions” in viscous liquids [17], increasing interest has been devoted to identifying possible growing length scales as mobility decreases [18,19]. A broad classification in terms of either static or dynamic length scales is usually used. Static (thermodynamic) length scales are determined by the free-energy landscape, whereas dynamic length scales are set by the rules governing the time evolution of the system and extracted from finite-time behaviour of time-dependent correlation functions and associated susceptibilities [6]. Even if growing static length scales have been reported by experiments [20] and simulations [21], there is still debate if they control the glass transition [22]. It is still not clear to what extent dynamic correlations are the consequence, or the primary origin of, slow dynamics occurring close to the glass transition [19].
Even if rooted at nanometric length scales, DH exerts clear influence at macroscopic level. One widely studied phenomenon is the breakdown of the Stokes–Einstein (SE) relation involving the diffusion coefficient D and the shear viscosity η (the more debated analogous phenomenon involving the rotational diffusion, where the breakdown is revealed [23,24] or missing [11], will not be considered here). For a single particle with radius R moving in a homogeneous fluid with viscosity η at temperature T, the SE relation states that
D = k B T ζ π η R
k B denotes the Boltzmann constant and ζ denotes a number depending on the boundary condition between the fluid and the particle [25]. Under a no-slip condition, ζ = 6 . Roughly, the SE law states that the quantity k B T / D η is a constant of the order of the size of the diffusing particle. Remarkably, despite its macroscopic derivation, SE also well accounts for the self-diffusion of many monoatomic and molecular liquids, provided the viscosity is low (≲10 Pa ·s) [26]. On the other hand, the finite diffusion coefficient of guest atoms in solid hosts, where viscous transport is missing strongly, suggest the SE failure close to the solidification occurring at the glass transition. In fact, a common feature of several fragile glass formers is the SE breakdown for increasing viscosity. The failure manifests itself as a partial decoupling between the diffusion and the viscosity, in the sense that D 1 increases less than η [4,11,27,28,29,30]. The decoupling is well accounted for by the fractional SE D η κ [31], where the non-universal exponent κ falls in the range 0.5 κ < 1 [15]. The usual interpretation of the SE breakdown relies on DH and the subsequent presence of a spatial distribution of characteristic relaxation times τ close to the glass transition [6,11,31]. The neat argument is that, although the viscosity is more sensitive to the longest relaxation times, the diffusivity is influenced by the shortest ones. As the shape of the distribution tends to widen on approaching the glass transition, the gap between D 1 and η increases as well, leading to the SE breakdown [4].
Diffusion, viscous transport, and structural relaxation involve time scales that are much longer that than the typical vibrational time t of the particle rattling in the cage of the first neighbours, typically a few picoseconds. The diffusion coefficient is expressed as D = 6 δ 2 / τ D , where τ D is the minimum time ensuring that the particle random displacements at a pace τ D are statistically independent with finite mean square value δ 2 [25]. On the other hand, the viscous flow requires the relaxation of the shear stress fluctuation, which occurs in a Maxwell time τ M = η / G , where G is the intermediate-time shear modulus [32]. On approaching the glass transition, t τ D , τ M , τ α .
Despite the huge difference in time scales, earlier [33] and later theoretical studies [5,34,35,36,37,38,39,40], and experimental ones [41], addressed the rattling process within the cage to understand the slow dynamics, rising a growing interest on the DW factor [8,29,30,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61]. Within this context, most interest has been devoted to the correlations between DW factor and the structural relaxation time τ α , which are found to be strong and encompassed by a universal master curve [47]:
τ α = F ( u 2 )
An analytical expression of the master curve is derived in Section 2. Alternative forms of the master curve are reported by Douglas and coworkers [8,9,10]. Correlations between DW factor and the structural relaxation time τ α are found in polymers in bulk [30,47,48,49] and thin films [61], binary atomic mixtures [48,55], colloidal gels [52], antiplasticised polymers [8,9], water [57] and water-like models [59,60]. The DW factor also provided an alternative interpretation of the so-called thermodynamic (or temperature/density) scaling [58]. The correlation between structural relaxation and DW factor has been inspected in the experimental data concerning several glass-formers in a wide range of fragility—the steepness index m defined by Angell [1] ( 20 m 191 ), including polymers, van der Waals and hydrogen-bonded liquids, metallic glasses, molten salts and the strongest inorganic glass-formers [47,50,51,55,56,57,58].
The structural relaxation time τ α is an average quantity which is certainly affected by DH but not in a straightforward way. Nonetheless, given the scaling expressed by Equation (2), it is legitimate to wonder if DH and fast vibrational dynamics exhibit correlations. Working in that direction, we have found positive answers and the present paper collects and reviews a selected part of our results, with a focus on the breakdown of the SE law [29,30]. Even if strictly related, we will not discuss here a study concerning ultrathin molecular films with strong mobility gradients analogous to DH, where the same scaling observed in bulk, Equation (2), has been revealed [61].
Our approach relies on the increasing evidence that the master curve, Equation (2), is a manifestation of a more fundamental correlation between the vibrational dynamics and the slow relaxation. It may be presented in the following terms. Let us consider a generic space- and time-dependent correlation function C ( x 1 , t 1 ; x 2 , t 2 ) , where x denotes a configuration of the liquid at a given time t, i.e., the set x of all the positions of the elementary microscopic particles (monomers, atoms, molecules, etc.). For steady states, C ( x 1 , t 1 ; x 2 , t 2 ) depends on the time difference t = t 2 t 1 . Let us set t 1 = 0 and define C ( x 0 ; x , t ) C ( x 0 , 0 ; x , t ) . If two states, labelled by a and b, have equal DW factor, the correlation function, when evaluated over the two states, has coincident time evolution at least between the typical vibrational time scale t and τ α [49]. Said otherwise, for t t τ α , it holds [32,47,49,50,51,53,54,55]:
u 2 ( a ) = u 2 ( b ) C ( x 0 ; x , t ) ( a ) = C ( x 0 ; x , t ) ( b )
In selected systems, Equation (3) holds beyond τ α and extends up to the diffusive regime, e.g., unentangled polymers and atomic binary mixtures [29,30,49,55].
Our studies were prompted by the finding by Widmer-Cooper and Harrowell that DH are predicted by particle displacements at short times [44]. However, it must be stressed that our DW factor is evaluated within the vibrational time scale t and not the time scale in [44], which is approximately one order of magnitude longer, a choice leading to differences for states with low viscosity.
The review outlines a model of the slow heterogeneous relaxation and transport in terms of vibrational dynamics in Section 2. The model is presented for completeness, but it is not essential to the understanding of the simulation results discussed in the rest of the paper. Later, a broad introduction to relaxation and transport in polymeric melts, and the correlation with the vibrational fast dynamics is given in Section 3 and Section 4, respectively. The signatures identifying the presence of heterogeneous dynamics are discussed in Section 5. The SE breakdown is presented in Section 6, with a final discussion on the length scale of the mutual influence between particle displacements in Section 7.

2. A model of the Slow Heterogeneous Relaxation and Transport in Terms of Vibrational Dynamics

An in-depth, microscopic understanding of the link between the fast and slow dynamics is still missing, even if the impact of anharmonicity has been noted [43,58,62]. Here, we present a model, extending first seminal ideas [34], where the key role is played by the DW factor u 2 , which is a single-particle quantity. Alternative pictures, in terms of the same quantity, are known [8,9,10]. Notice that, even if a single-particle quantity, the DW factor encodes information on collective dynamics and spatially extended cooperative phenomena [32,53,54,62,63,64].
At the present level of development, the model delivers expressions of the diffusion coefficient and the structural relaxation time in terms of the DW factor. It also accounts for the nonexponential character of the relaxation, an aspect which will be not presented here. However, even if it incorporates some consequence of DH, i.e., the presence of a wide distribution of relaxation times p ( τ ) , it does not cover any spatial aspect related to DH, which instead has been revealed by the simulations, as we will see in Section 5.1 and Section 7, and accounted for by Equation (3).

2.1. Relaxation Time

A first basis to connect fast and slow degrees of freedom was developed by Hall and Wolynes who, assuming that atomic motion is restricted to cells, pictured the glass transition as a freezing in an aperiodic crystal structure [34]. As a result, the viscous flow is described in terms of activated jumps over energy barriers Δ E k B T a 2 / u 2 , where a is the displacement to reach the transition state. The usual rate theory leads to the Hall–Wolynes equation:
τ α ( H W ) ( a 2 ) , η ( H W ) ( a 2 ) exp a 2 2 u 2
Equation (4) has the form of Equation (2). A very similar relation was derived by Buchenau and Zorn, in terms of soft vibrational modes [41]. Equation (4) is expected to fail when τ α becomes comparable to the typical rattling times of each atom in the cage, corresponding to picosecond timescales. This condition is quite mild, e.g., in selenium it occurs ~100 K above the melting temperature [41].
One basic assumption of Equation (4) is that the distance to reach the transition state has a characteristic value a. Actually, this length scale is dispersed. To constrain the related distribution, p ( a 2 ) , it is assumed that the latter does not depend on the state parameters such as the temperature, the density or the interacting potential. This complies with the spirit of the work in [34], where the a distance is said to be mostly controlled by the geometrical packings. It is also known that, irrespective of the relaxation time, τ α , the average distance moved by the relaxing unit within τ α is approximately the same, i.e., a fraction of the molecular diameter [1]. Averaging Equation (4) over the distribution p ( a 2 ) yields the structural relaxation time
τ α = τ α ( H W ) ( a 2 ) a 2
0 τ α ( H W ) ( a 2 ) p ( a 2 ) d a 2
Note that Equation (6) assumes that the distribution of the relaxation times is mainly due to the distribution of the displacement to reach the transition state in the different local environments, whereas the average DW factor u 2 is taken as homogeneous across the sample. This viewpoint relies on the picture that relaxation is related to long wavelength soft modes [41,46]. Support has been provided by the strong correlation observed in glass-formers between u 2 and the elastic modulus under quasi-static mechanical equilibrium [32].
As a suitable choice, the distribution of the squared distances p ( a 2 ) is taken as a truncated Gaussian form [47,48]
p ( a 2 ) = A exp ( a 2 a 2 ¯ ) 2 2 σ a 2 2 if   a > a m i n 0 otherwise
where A is the normalization ensuring 0 p ( a 2 ) d a 2 = 1 and a m i n 2 is the minimum displacement to reach the transition state. Given the weak influence, and to get rid of an adjustable parameter, one takes a m i n 2 = 0 [47,48]. The motivations behind the Gaussian form of p ( a 2 ) mainly rely on the Central Limit Theorem. In fact, a 2 ( r 0 2 in the notation in [34]) is the cumulative squared displacement of the N m o n particle that move [34].
Plugging Equation (7) into Equation (6) leads to the following generalized HW equation (GHW),
τ α = τ 0 exp a 2 ¯ 2 u 2 + σ a 2 2 8 u 2 2
τ 0 is a suitable constant. An analogous law is anticipated for the viscosity η , given the known near proportionality with τ α [3]. Equation (8) is the form of the master curve Equation (2) being adopted in our studies. Other variants useful in the comparison with numerical and experimental results are listed in Appendix A and Appendix B.
Obviously, if the distribution p ( a 2 ) is narrow and centred at a 0 2 , Equation (8) must reduce to the expression derived by Hall and Wolynes, Equation (4), τ α ( H W ) ( a 0 2 ) . For the specific choice of p ( a 2 ) , given by Equation (10), Equation (8) shows that this happens if σ a 2 2 / 8 u 2 2 a 2 ¯ / 2 u 2 , namely, the ratio R defined as
R σ a 2 2 / 4 a 2 ¯ u 2
is vanishingly small. Equation (9) depends on the magnitude of DW factor so that, being the parameters σ a 2 2 and a 2 ¯ independent of the physical state, the presence of the distribution p ( a 2 ) is negligible when the DW factor is large, thus leading to a very narrow distribution of relaxation times, a characteristic of homogeneous dynamics. This suggests to read the condition R = 1 as the crossover between homogeneous and heterogeneous dynamics, i.e.,
R 1 homogeneous dynamics R 1 heterogeneous dynamics
Finally, we notice that the distribution p ( a 2 ) in Equation (7) with a m i n 2 = 0 may be recast via Equation (4), as a log-normal distribution of relaxation times p ( ln τ )
p ( ln τ ) = B exp 2 u 2 2 σ a 2 2 ln τ τ ¯ 2 if   τ τ 0 0 otherwise
where B is the normalization ensuring p ( ln τ ) d ln τ = 1 , τ ¯ = τ 0 exp ( a 2 ¯ / 2 u 2 ) . An interesting feature of p ( ln τ ) is that its width σ a 2 / u 2 increases by decreasing the DW factor.

2.2. Diffusion Coefficient

The diffusion coefficient D may be expressed via the above model by the relation [30]
D = 1 6 a 2 τ α ( H W ) ( a 2 ) a 2
The above equation assumes that displacements as large as a occurring in a time τ α ( H W ) ( a 2 ) are statistically independent. Notice that, although Equation (6) signals that the structural relaxation time is affected by the larger a 2 values, i.e., the longest relaxation times of the distribution p ( ln τ ) , the diffusivity, according to Equation (12), is influenced by the shorter ones.
The explicit expression of the diffusion coefficient and an approximated version are given in Appendix A.2.

2.3. Stokes–Einstein Product

The Stokes–Einstein (SE) relation, Equation (1), states that the quantity D η / T is constant if the diffusing particle changes neither the size nor the boundary conditions with the liquid. As the numerical evaluation of the viscosity is a delicate point, proxies are often used [27,65]. As an example, as η T τ α in unentangled polymers [66], it is more suitable to study the breakdown of the SE law by considering the SE product
K S E = D M τ α
where M is the number of monomers. K S E is expected to be independent of the chain length, as D 1 / M in unentangled polymers [66] and the monomer relaxation at τ α poorly senses the chain connectivity. The above equation with M = 1 may be also used for liquids where the elementary units are atoms or small molecules, as the temperature factor in the ratio D η / T provides a change of approximately ∼20% in fragile glass-formers [3], much less than the observed increase of K S E on approaching the glass transition [4,11]. The explicit expression of the SE product K S E derived within the vibrational dynamics model and an approximated version K ˜ S E are given in Appendix A.3.

3. Transport and Relaxation in Polymeric Melts

The correlation between diffusivity, slow relaxation and fast vibrational dynamics has been studied by Molecular-Dynamics (MD) simulations of a coarse-grained model of a melt of linear unentangled polymer. Details about the model are given in Section 9. Even if rather crude, the model was proven to capture the universal aspects of the correlation and allowed an effective comparison with the experiment [47].
To provide a microscopic picture of the transport, the mean square displacement (MSD) of the monomer r 2 ( t ) is usually considered:
r 2 ( t ) = 1 N i x i ( t ) x i ( 0 ) 2
where x i ( t ) is the position of the i-th monomer at time t. In addition to MSD, with the purpose of characterizing the relaxation, the self part of the intermediate scattering function (ISF) is also considered [26]:
F s ( q , t ) = 1 N j N e i q · ( x j ( t ) x j ( 0 ) )
In an isotropic liquid, ISF depends only on the modulus of the wavevector q = | | q | | and features the rearrangements of the spatial structure of the fluid over the length scale ∼ 2 π / q , leading to a decaying profile in time starting from F s ( q , 0 ) = 1 . In our case, ISF was evaluated at q = q m a x , the maximum of the static structure factor ( 7.13 q m a x 7.55 ) corresponding to the length scale of the monomer size. F s ( q m a x , t ) vanishes when the monomer displacement in a time t largely exceeds the monomer diameter. The time needed to make F s ( q m a x , t ) small is a measure of the escape time of the monomer from the cage formed by the neighbours, also known as the structural relaxation time τ α , customarily defined by the relation F s ( q m a x , τ α ) = e 1 .
Figure 2 shows typical MSD and ISF curves of the polymeric monomers. At very short times (ballistic regime), MSD increases according to r 2 ( t ) ( 3 k B T / m ) t 2 and ISF starts to decay. The repeated collisions with the other monomers slow the displacement of the tagged one, as evinced by the knee of MSD at t 12 / Ω 0 0.17 , where Ω 0 is an effective collision frequency, i.e., it is the mean small oscillation frequency of the monomer in the potential well produced by the surrounding ones kept at their equilibrium positions [64,67]. At later times, a quasi-plateau region, also found in ISF, occurs when the temperature is lowered and/or the density increased. This signals the increased caging of the particle. Trapping is terminated after an average time τ α . For t τ α , MSD increases more steeply. The monomers of short chains ( M 3 ) undergo diffusive motion r 2 ( t ) t δ with δ = 1 . For longer chains, owing to the increased connectivity, the onset of the diffusion is preceded by a subdiffusive region ( δ < 1 , Rouse regime) [68]. At long time, the monomer displaces in a diffusive way with diffusion coefficient D = lim t r 2 ( t ) / 6 t .

4. Correlation between Vibrational Fast Dynamics and Slow Relaxation

4.1. Vibrational Caged Dynamics and Debye–Waller factor

In our model polymer, the term “vibrational dynamics” refers to the rattling of the trapped monomer within the cage formed by the closest monomers. It is crucial to provide a robust criterion to assess the presence of the cage, which is anticipated to lack in liquids with high molecular mobility and fast relaxation. Compelling evidence of the cage effect is provided by the time velocity correlation function, which, after a first large drop due to pair collisions, reverses the sign since the monomer rebounds from the cage wall [64]. As an alternative route to reveal the cage effect, we consider the slope of MSD in the log-log plot
Δ ( t ) log r 2 ( t ) log t
Representative plots of Δ ( t ) for the polymer system are given in the inset of Figure 2a. Δ ( t ) tends to 2 at short times, due to the ballistic motion, and reaches the plateau level 1 at long times, owing to the diffusive motion. In the absence of caging effect, Δ ( t ) decreases in a monotonous way on increasing time. Caging is indicated by the presence of a minimum of Δ ( t ) occurring, irrespective of the physical state in the present model polymer, at t = 1 . 0 ( 4 ) . In actual time units, t is ~1–10 ps [69].
The presence of the minimum paves the way to a robust definition of the DW factor u 2 , the mean square rattling amplitude of the monomer during the trapping period. In fact, the minimum, corresponding to the inflection point in the log-log plot of r 2 ( t ) ), separates two regimes. At short times, t < t , the inertial effects dominate, whereas for t > t , early escapes from the cage become apparent. Therefore, a convenient definition of the DW factor as a mean localization length is just MSD at t :
u 2 r 2 ( t = t )

4.2. Debye–Waller Scaling of the Slow Relaxation

The monomer dynamics depends in a complex way on the state parameters. Nonetheless, there is clear correlation between the DW factors and the long-time relaxation dynamics. First examples are shown in Figure 2 by considering MSD and ISF. Note that states with equal DW factor have coincident time evolution of both MSD and ISF at least between t and τ α [49]. In Section 5.1, it will be shown that these results are a manifestation of Equation (3).
It is seen that the coincidence of the MSD curves is lacking at times longer than τ α for states corresponding to polymer chains with different length. This effect is not a failure of the scaling at times exceeding τ α , but a mere consequence of the complex dependence of MSD on the chain length since it is affected by all the Rouse modes [66]. In fact, if the correlation function of the single Rouse mode with the slowest relaxation time is singled out, i.e., the one with characteristic relaxation time given by the average chain reorientation time τ e e [66], the scaling is still observed after proper account of the chain length dependence, see Figure 3. The finding proves that Equation (3) holds also at a time τ e e being much longer than τ α .
As a side product of the coincidence of the ISF curves in states with equal DW factor seen in Figure 2, one has that states with equal DW factor u 2 have equal structural relaxation time τ α too. This can be reformulated via the master curve Equation (2), which, according to the model detailed in Section 2.1, takes the form given by Equation (8), i.e., a simple parabolic law between log τ α and 1 / u 2 [47,48]. Figure 4 tests Equation (8), written in the form given by Equation (A1) for a wide variety of physical states of our model polymeric melt [48]. It is also shown that the scaling holds if one considers the end–end chain reorientation time τ e e , i.e., the time needed by the correlation function C e e ( t ) to drop to e 1 , see Figure 3; although, in this case, it is described by a different master curve.

5. Signatures of the Heterogeneous Dynamics

MSD and ISF well-expose the cage effect, whereas the possible DH influence on their shape is less apparent. Figure 1 shows that DH is characterized by the presence of clusters of monomers with rather different mobility [11,12]. We now present and discuss two quantities well tailored to provide quantitative insight into this aspect.

5.1. van Hove Function

One central quantity of the DH analysis is the self part of the van Hove function G s ( r , t ) [26]:
G s ( r , t ) = 1 N i = 1 N δ [ r + x i ( 0 ) x i ( t ) ]
where x i ( t ) is the position of the i-th monomer at time t, and δ [ · ] is the three-dimensional Dirac delta function. In isotropic liquids, the van Hove function depends on the modulus r of r . The interpretation of G s ( r , t ) is direct. The product G s ( r , t ) · 4 π r 2 is the probability that the monomer is at a distance between r and r + d r from the initial position after a time t. The moments of G s ( r , t ) are of interest:
r n ( t ) = 4 π 0 r n G s ( r , t ) r 2 d r
For n = 2 , one recovers the usual mean square displacement (MSD). If the monomer displacement is a Gaussian random variable, G s ( r , t ) reduces to the Gaussian form [26]:
G s g ( r , t ) = 3 2 π r 2 ( t ) 3 / 2 exp 3 r 2 2 r 2 ( t )
Equation (20) is the correct limit of G s ( r , t ) at very short (ballistic regime, r 2 ( t ) = 3 k B T / μ t 2 ) and very long times ( diffusion regime, r 2 ( t ) = 6 D t , where D is the monomer diffusion coefficient).
The spatial Fourier transform of the self part of the van Hove function yields the ISF function, Equation (15) [26].
Figure 5a presents the self-part of the van Hove function G s ( r , t ) , evaluated at τ α for the set of states with different mobility and relaxation shown in Figure 2. It is seen that if the relaxation and the mobility are fast, the shape of G s ( r , τ α ) decreases by increasing the displacement r from the initial position. On the other hand, the states belonging to the D and E set, the ones with slowest relaxation, exhibit a tendency toward a bimodal pattern, namely, in addiction to particles undergoing small displacements, a shoulder at r 1 (the monomer diameter) is observed. This signals the presence of particles exhibiting fast displacements by solid-state jump dynamics [27]. Said otherwise, the quasi-bimodal pattern of the van Hove function is clear signature of DH. Four other aspects are to be noted:
  • The self-part of the van Hove function is expressed by suitable correlation functions, see Appendix B. Then, the coincidence of G s ( r , τ α ) in states with equal DW factor observed in Figure 5a (the sets of states labelled as A, ⋯, E) is in harmony with Equation (3).
  • Equation (3) also holds if one inspects the spatial dependence of the correlation function, e.g., the van Hove function, at τ α . In particular, even in the presence of DH.
  • Given their relation with G s ( r , t ) , the coincidence of both MSD and ISF observed in Figure 2 for the sets of states labelled as A, ⋯, E is strictly linked to the one observed in Figure 5a.
  • The pattern of the D and E sets of states is not consistent with the Gaussian limit G s g ( r , τ α ) , Equation (20), predicting a progressive decay with r, i.e., the DH dynamics is not Gaussian;
To quantify the deviations of the self-part of the van Hove function G s ( r , τ α ) from the Gaussian limit, one defines the quantity [27,70]
N s ( r , τ α ) = G s ( r , τ α ) G s g ( r , τ α ) G s g ( r , τ α )
Figure 5b plots the ratio N s ( r , τ α ) . It exhibits increasing positive deviations at both short and large r values, evidencing the excess of nearly immobile and highly mobile monomers with respect to purely Gaussian behaviour, respectively. The analysis, in terms of the ratio N s ( r , τ α ) , reveals the wide distribution of mobilities pictured in Figure 1, right.

5.2. Non-Gaussian Parameter

An effective quantity to expose the time evolution of the non-Gaussian character of DH dynamics is the non-Gaussian parameter (NGP) [26]:
α 2 ( t ) = 3 5 r 4 ( t ) r 2 ( t ) 2 1
where r 2 ( t ) and r 4 ( t ) are the mean square and quartic displacements of the particle at time t, respectively. α 2 vanishes if the displacement is Gaussian, i.e., it follows from a series of independent elementary steps with finite mean and variance.
Figure 6 plots the NGP time evolution, Equation (22), for the set of states A, ⋯, E and additional states with very slow relaxation. It is seen that NGP vanishes at very short times, as the ballistic regime is Gaussian in nature. At intermediate times, a peak value α 2 m a x is observed increasing with the relaxation times [46,71,72]. The maximum occurs at a time slightly shorter than the structural relaxation time τ α , as in simpler molecular systems [27]. A snapshot of the microscopic mobilities in a lapse of time τ α , where DH is quite apparent, is plotted in Figure 1 (right). At later times, NGP decreases as the monomer dynamics enters the homogeneous diffusive regime, which is a Gaussian process [25].
It is seen that states belonging to the same set A, ⋯, E, i.e., with equal DW factor, have identical NGP in the time window [ t , τ α ] at least. This agrees with Equation (3), given the relation of NGP with the moments of the self part of the van Hove function G s ( r , t ) , Equation (19), and the expression of the latter in terms of suitable correlation functions, see Appendix B [73]. Note also the exponential increase of α 2 m a x with the ratio R defined in Equation (9) [47,48]. This is in harmony with the inequalities in Equation (10), stating that DH is characterized by R > 1 .

6. Breakdown of the Stokes–Einstein (SE) Law in the Presence of Dynamical Heterogeneity

6.1. SE Breakdown in Unentangled Polymers

The SE law is usually derived by considering the diffusivity of macroscopic bodies displacing in homogeneous viscous liquids [25]. The diffusion in the presence of strong DH does not comply with the SE law [15]. We have studied the SE breakdown in melts of unentangled polymers [29]. In these systems, helpful features are found [66]: (i) the diffusion coefficient D is inversely proportional to the chain length M, and (ii) the viscosity η is proportional to the end–end reorientation time which, in turn, is proportional to the structural relaxation time, e.g., see Figure 3, showing that states with equal structural relaxation time also have equal end–end reorientation time. Then, as discussed in Section 2.3, the study of the validity of the SE law is more efficiently carried out in terms of the product D M τ α , which is anticipated to be state-independent if the SE law holds.
Figure 7 shows that in states with homogeneous Gaussian dynamics, i.e., with small α 2 m a x values, the R values are comparable or less than the unit value and the product D M τ α is nearly constant, i.e., the SE law holds true. On the other hand, in the presence of significant DH, i.e., α 2 m a x > α 2 , c m a x = 0 . 40 ( 5 ) , one finds R > R c = 1 . 9 ( 1 ) and the product D M τ α tends to increase, i.e., the SE law fails [11,12,27]. The comparison between α 2 m a x and R substantiates the conclusion that the magnitude of the ratio R allows one to conclude whether DH is appreciable or not, as suggested by Equation (10). As the ratio R—apart from constants—depends only on DW, see Equation (9), the finding supports previous conclusions that the long-time DH is rooted in the fast dynamics [44].
It is seen that states with equal R (Figure 7a), i.e., states with equal DW factor according to Equation (9), exhibit nearly equal values of the product D M τ α . A similar result has been reported for atomic binary mixtures [55] and metallic alloys [30]. Recognising that the diffusivity D and the structural relaxation time τ α reflect the long time behaviour of MSD and ISF, respectively, Figure 7a reveals that Equation (3) is valid even in the diffusive regime which is entered in polymer melts at times fairly longer than τ e e , being τ e e τ α .

6.2. Quasi-Universal SE Breakdown of Fragile Glass-Formers

Having noted that the SE failure is tracked by the DW factor in unentangled polymers, we now pose the question if this finding exhibits universal features. To this aim, we consider the ratio K S E / K 0 with K S E defined in Equation (13) and K 0 a scaling factor to ensure the unit limit value at large DW factor.
In Figure 8, we plot the ratio K S E / K 0 as a function of u 2 / u g 2 . We complement the MD results on unentangled polymers already presented in Figure 7 with other MD data, considering the diffusion of Cu and Zr atoms in metallic alloy, A and B atoms in a Lennard–Jones binary mixture [30] together with experimental data concerning the popular fragile glass-former ortho-terphenyl (OTP) [74,75]. Figure 8 evidences the good collapse of the SE violation in terms of the reduced DW.
Figure 8 offers the opportunity to test the master curve predicted by the model of Section 2 with no adjustable parameters K ^ S E , Equation (A8), and its approximation K ˜ S E , Equation (A10). It is seen that K ^ S E predicts a stronger SE breakdown than actually observed. Larger deviations are exhibited by the approximant K ˜ S E . How to improve the agreement: The expression of the diffusion coefficient in Section 2.2, assuming that displacements as large as a are statistically independent, aims at a SE product K ^ S E , Equation (A8), with no additional adjustable parameters with respect to the ones of τ α , i.e., the ones of Equation (A2). This puts severe constraints on the shape of the distribution of the square displacements needed to overcome the relevant energy barriers p ( a 2 ) , Equation (7). The form of the distribution is adequate for large displacements to reach the transition state governing τ α , as proven by the effective fit of the MD data by the predicted master curve shown in Figure 4. However, the findings of Figure 8 suggest that it must be improved for small displacements affecting D. Alternatively, we may also state that the distribution p ( ln τ ) , Equation (11), should be refined as far as the short relaxation times are concerned.
Figure 8 shows that better agreement occurs by assuming the fractional SE form D τ α τ α 1 κ [15,31] with τ α as given from Equation (A2). The best fit is found with κ = 0 . 85 , which equals the universal value found by Mallamace et al. [76], deviating from the prediction of the “obstruction model” κ = 2 / 3 [15].

7. Displacement Correlation Length

Several results of the present paper suggest that the vibrational dynamics, as sensed by DW factor, provides insight into DH. A line of attack to understand how vibrational dynamics is related to slow relaxation is provided by the model of Section 2. The model is based on the distribution of the (squared) displacements needed by a particle to rearrange in the different local environments, Equation (7), leading in turn to the distribution of relaxation times, Equation (11). As noted in Section 6.2, the model needs further development. A further aspect to be improved is the absence of any detail on the localization of the particles with a given dynamics. This prevents any prediction concerning a peculiar aspect of DH, i.e., the existence of spatial domains with characteristic length scales where the particles undergo correlated motion, e.g., see Figure 1.
To make progress, it is worthwhile to preliminarily judge whether DW exhibits some correlation with possible dynamic length scales. To pursue this task, we studied the following monomer displacement–displacement correlation (DDC) functions [53,54]:
C u ( r , τ α ) = u ^ i ( t 0 , τ α ) · u ^ j ( t 0 , τ α ) ,
C δ u ( r , τ α ) = δ u i ( t 0 , τ α ) δ u j ( t 0 , τ α ) / [ δ u ( t 0 , τ α ) ] 2 .
An average over all the i-th and j-th monomers spaced by r is understood. u ^ k ( t 0 , t ) is the versor of the displacement vector of k-th monomer in a time interval from t 0 to t 0 + t , u k ( t 0 , t ) = r k ( t 0 + t ) r k ( t 0 ) and δ u k ( t 0 , t ) = | u k ( t 0 , t ) | | u ( t 0 , t ) | , where | u k ( t 0 , t ) | is the modulus of the displacement. Henceforth, C δ u ( r , τ α ) and C u ( r , τ α ) will be referred to as modulus (or mobility) and direction DDC functions, respectively. Local anisotropies and collective elastic solid-like response to the rattling of the monomers in the cage of their neighbours play a central role in the DDC build-up [64].
We consider DDCs of the states presented in part of the states in Figure 2. We remind that the states (i) exhibit different DH degree, e.g., see Figure 5 and Figure 6, and (ii) are grouped in sets labelled B through E, each set being characterized by a single value of the DW factor.
Figure 9a,b shows the spatial dependence of the direction and the modulus DDC functions, respectively, for the sets of states labelled B through E in Figure 2. Both correlation functions manifest damped oscillations in-phase with the pair correlation function g ( r ) , thus evidencing that the correlated motion of a tagged monomer and its surroundings is influenced by the structure of the latter. This agrees with previous work on DDCs in Lennard–Jones systems [72,77], hard-sphere [78] and experiments on colloids [79]. The highest correlations are reached at a distance corresponding to the bond length b = 0.97 which demonstrates the highly concerted dynamics of bonded monomers. The correlation peaks, located at the first-, second-,... neighbours shells, vanish approximately in an exponential way on increasing the distance from the tagged particle (see insets of Figure 9). In more detail, Figure 9a shows that the direction correlations do not show significant increase in their spatial extension on increasing the structural relaxation time. Figure 9b shows the modulus (mobility) correlations. Differently from the direction correlations, their spatial extension increase meaningfully with the structural relaxation time (see also the inset of Figure 9b).
Figure 9b clearly shows that physical states with equal DW factor, i.e., belonging to the same set of states (B, ⋯, E), exhibit the same spatial correlations. This provides further support that Equation (3) also holds if the spatial dependence of the correlation function is considered for a given time up to τ α at least. To provide additional insight, we evaluated the length scales of the exponential decays of the DDC maxima with the distance exp [ r / ξ X ( τ α ) ] with X = u , δ u , thus defining two distinct dynamic correlation lengths pertaining to direction and modulus DDCs, ξ u ( τ α ) and ξ δ u ( τ α ) , respectively. Figure 10 shows these quantities. It is seen that the spatial extension of the modulus DDC increases quite a lot with τ α and reaches distances beyond the next-nearest neighbours for the states with the slowest relaxation. Instead, the direction correlations are virtually independent of the structural relaxation. Irrespective of the correlation length under consideration, Figure 10 also shows that they are equal within the errors for states with equal DW factor, i.e., belonging to the same set of states (B, ⋯, E).
We are now in a position to compare our results with previous work on DDCs. Simulations of Lennard–Jones binary mixture (BM) observed that at time t α , corresponding to maximum dynamic heterogeneity, ξ δ u B M ( t α ) increases as the temperature decreases, whereas ξ u B M ( t α ) is almost constant [80]. This agrees with our findings in Figure 10 concerning unentangled polymers. As to the modulus DDC correlation length, one finds [53] that after suitable algebraic manipulation to allow comparison [79], our changes of ξ δ u ( τ α ) with τ α are in quantitative agreement with the results of Bennemann et al. reported in a study of the same polymer system investigated here [72].

8. Discussion

There is wide experimental, numerical and theoretical evidence that the fast vibrational dynamics, as sensed by the Debye–Waller factor u 2 , and the time scale τ α of the slow microscopic reorganisation of a liquid close to the transition to the glassy state are correlated in an universal way. Potential applicative implications concerning the quick characterisation of the stability of disordered structures with ultraslow relaxation are apparent. Less attention has been paid to a series of numerical MD simulation studies concluding in favour of strong correlations also between the vibrational dynamics and the dynamical heterogeneity, the spatial distribution of long-time mobility developing when approaching the disordered solid state. We reviewed these studies, mainly concerning melts of unentangled linear polymers, unifying all the results for the first time in terms of Equation (3). The latter has been tested both in space and time. In particular, we considered time-dependent quantities accounting for transport and relaxation like MSD, ISF and NGP and showed that they are related to the self-part of the van Hove function, which reduces to suitable correlation functions. In this respect, the correlation between the Debye–Waller factor and the breakdown of the SE law, a hallmark of DH presence, is seen as an ancillary consequence of the extension of Equation (3), which at times is much longer than τ α , where the motion is diffusive. We also inspected Equation (3) in space by considering both the van Hove function and DDC functions. Notably, DDC functions are collective in nature, differently from the self-part of the van Hove function, which is a single-particle quantity. This suggests that Equation (3) also holds for collective correlation functions. A further validation of this conclusion is offered by the collective stress–stress correlation function, which has been presented elsewhere [32] and not discussed in this review.
The understanding of the microscopic origin of the correlation between the vibrational dynamics and the heterogeneous dynamics close to the glass transition is still unsatisfactory in several respects. In particular, both the model discussed here, as well as other ones reported in the literature, even if successful in relating the Debye–Waller factor u 2 with the time scale τ α , are currently unable to account for the fact evidenced by the numerical simulations that the vibrational dynamics conveys also information on the spatial correlations between the mobility of different particles.

9. Methods

Most results discussed in this review concern a coarse-grained model of a linear polymer chain with M monomers is adopted. Bending and torsional potentials are neglected, i.e., the chain is fully flexible. While addressing the interested reader to the referenced papers for further details, we provide here some general aspect of the numerical model. We considered systems with total number of monomer N 2000 . Non-bonded monomers at a distance r interact via the truncated parametric potential:
U q , p ( r ) = ϵ p q p σ r q q σ r p + U c u t
where σ = 2 1 / 6 σ and the value of the constant U c u t are chosen to ensure U p , q ( r ) = 0 at r r c = 2.5 σ . The minimum of the potential U p , q ( r ) is at r = σ , with a constant depth U ( r = σ ) = ϵ . Note that U q , p ( r ) = U p , q ( r ) . Bonded monomers interact with a potential which is the sum of the Finitely Extendible Nonlinear Elastic (FENE) potential and the Lennard–Jones (LJ) potential [71]. The resulting bond length is b = 0.97 σ , within a few percent. We set σ = 1 and ϵ = 1 . The time unit is τ M D = ( m σ 2 / ϵ ) 1 / 2 , with m being the mass of the monomer. Temperature is in units of ϵ / k B , where k B is the Boltzmann constant. We set m = k B = 1 . All the data presented in this work are expressed in reduced MD units. It is interesting to map the reduced MD units to real physical units. The procedure involves the comparison of the experiment with simulations and provide the basic length ( σ ), temperature ( ϵ / k B ) and time ( τ M D ) units [69,71,81,82,83]. For polyethylene and polystyrene, it was found σ = 5.3 Å, ϵ / k B = 443 K, τ M D = 1.8 ps and σ = 9.7 Å, ϵ / k B = 490 K, τ M D = 9 ps, respectively [69]. For poly(vinyl alcohol) σ = 5.2 Å, ϵ / k B = 550 K and τ M D = 1.63 ps [83]. For polyisoprene σ = 6.7 Å, ϵ / k B = 307 K and τ M D = 10 ps [81]. The densities used in this and other studies are lower than the densities at atmospheric pressure, e.g., when mapping our model to polyethylene and polystyrene, one finds ∼0.5 and ∼0.7 g/cm 3 , to be compared to the actual values 0.78 and 0.92 g/cm 3 , respectively [69].

Author Contributions

F.P., A.T. and D.L. wrote the manuscript together. All authors have read and approved the final version.

Funding

This research was funded by the project PRA- 2018-34 (“ANISE”) from the University of Pisa.

Acknowledgments

A generous grant of computing time from IT Center, University of Pisa, and Dell EMC® Italia is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript.
DDCdisplacement–displacement correlation
DHdynamical heterogeneity
DWDebye–Waller
ISFIntermediate scattering function
MDMolecular-dynamics
MSDMean square displacement
NGPnon-Gaussian parameter
SEStokes–Einstein

Appendix A

Appendix A.1. Structural Relaxation

According to the model detailed in Section 2.1, the master curve relating the structural relaxation time τ α and the DW factor u 2 , Equation (2), takes the form given by Equation (8). Other variants of Equation (8) are of interest in the comparison with numerical and experimental results. As an example, the best fit of the master curve Equation (8) with the numerical data concerning the melt of unentangled polymer chains of interest here yields (in MD units) [47,48]:
log τ α = α + β 1 u 2 + γ 1 u 2 2
where α = 0.424 ( 1 ) , β = a 2 ¯ / ( 2 ln 10 ) = 2.7 ( 1 ) · 10 2 , γ = σ a 2 2 / ( 8 ln 10 ) = 3.41 ( 3 ) · 10 3 .
To recast Equation (A1) as a universal master curve removing system-dependent quantities, one considers the DW factor at the glass transition u g 2 (defined via τ α = 10 2 s) and introduces the reduced variable u 2 / u g 2 , so as to write Equation (A1),
ln τ α = α ^ + β ^ u g 2 u 2 + γ ^ u g 2 u 2 2
where α ^ = 2 ln 10 β ^ γ ^ . The ansatz is that Equation (A2) is system-independent and both β ^ and γ ^ are universal coefficients. To derive their numerical values, we use the value < u g 2 > 1 / 2 = 0.129 and the best-fit values of α , β , γ , drawn from the numerical simulations of the melt of unentangled polymer chains [47,48]. This yields β ^ = 3.7 ( 1 ) and γ ^ = 28.4 ( 2 ) [30].

Appendix A.2. Diffusion Coefficient

The diffusion coefficient is evaluated from Equation (12) via Equations (4) and (7). The result is
D ( u 2 ) = σ a 2 6 τ 0 exp 1 2 a 2 ¯ σ a 2 2 G ( u 2 ) 18 π 1 + erf ( a 2 ¯ / 2 σ a 2 )
where
G ( x ) = 1 π Λ ( x ) exp Λ 2 ( x ) erfc [ Λ ( x ) ] Λ ± ( x ) = σ a 2 2 ± 2 a 2 ¯ x 8 σ a 2 x .
As τ α ( H W ) ( a 2 ) depends on the parameter a in a much more marked way than a 2 , see Equation (4), an effective approximation of the diffusivity, as expressed by Equation (12), is D D ˜ with
D ˜ = a 0 2 6 1 / τ α ( H W ) ( a 2 ) a 2
where a 0 is a constant.

Appendix A.3. Stokes–Einstein Product

After suitable manipulation the Stokes–Einstein product K S E , Equation (13), takes the from ( M = 1 )
K S E ( u 2 ) = D ( u 2 ) × τ α ( u 2 )
where
K S E ( x ) = σ a 2 F ( x ) G ( x ) F ( x ) = 1 18 π exp Λ + 2 ( x ) exp a 2 ¯ σ a 2 2 1 + erf ( Λ + [ x ] ) 1 + erf ( a 2 ¯ / 2 σ a 2 ) 2 .
where the auxiliary functions G ( x ) and Λ ± ( x ) are defined in Appendix A.2.
When expressed in terms of the adimensional quantity z = u 2 / u g 2 , the product K S E takes the form K ^ S E ( z ) with
K ^ S E ( z ) = u g 2 F ^ ( z ) 1 π Λ ^ ( z ) exp Λ ^ 2 ( z ) erfc [ Λ ^ ( z ) ] Λ ^ ± ( z ) = 2 γ ^ ± β ^ z 2 γ ^ 1 / 2 z F ^ ( z ) = 2 3 γ ^ π exp Λ ^ + 2 ( z ) exp β ^ 2 2 γ ^ 1 + erf ( Λ ^ + [ z ] ) 1 + erf ( β ^ / 2 γ ^ 1 / 2 ) 2 .
The quantities β ^ and γ ^ are defined in Appendix A.1.
The Stokes–Einstein product K S E may be approximated as K S E K ˜ S E with K ˜ S E = D ˜ τ α , where D ˜ is defined in in Appendix A.2. This yields
K ˜ S E = a 0 2 6 τ α ( H W ) ( a 2 ) a 2 1 / τ α ( H W ) ( a 2 ) a 2
The explicit expression of the quantity K ˜ S E , in terms of z = u 2 / u g 2 , reads
K ˜ S E ( z ) = a 0 2 6 exp 2 γ ^ / z 2 1 + erf γ ^ 1 / 2 z + β ^ 2 γ ^ 1 / 2 erfc γ ^ 1 / 2 z β ^ 2 γ ^ 1 / 2 1 + erf β ^ 2 γ ^ 1 / 2 2

Appendix B

The van Hove function for a uniform fluid is defined as
G ( r , t ) = 1 N i = 1 N j = 1 N δ [ r + x i ( 0 ) x j ( t ) ]
Physically, G d r is the probability of finding a particle j in a region d r around a point r at time t if the particle i was at the origin at time 0. We may recast the van Hove function by resorting to the time-dependent, microscopic particle density [26]
ρ ( r , t ) = i = 1 N δ [ r x i ( t ) ]
We rewrite Equation (A11) as
G ( r , t ) = 1 N i = 1 N j = 1 N δ [ r + r x j ( t ) ] δ [ r x i ( 0 ) ] d r
= 1 N ρ ( r + r , t ) ρ ( r , 0 ) d r
= 1 ρ ρ ( r , t ) ρ ( 0 , 0 )
where ρ is the average number density. Equation (A15) shows that the van Hove function is proportional to the density correlation function. It is easily shown that the van Hove function may be written as [26]
G ( r , t ) = G s ( r , t ) + G d ( r , t )
where G s ( r , t ) and G d ( r , t ) are usually called the “self” and “distinct” parts. Equation (18) provides the explicit expression of G s ( r , t ) . The distinct part is written as
G d ( r , t ) = 1 N i = 1 N j i N δ [ r + x i ( 0 ) x j ( t ) ]
Finally, we show that both G s ( r , t ) and G d ( r , t ) may be expressed in terms of suitable correlation functions. To this aim, we define the auxiliary function B i ( r , t ) δ [ r x i ( t ) ] . By repeating the same passages leading from Equation (A11) to Equation (A13), one has
G s ( r , t ) = 1 N i = 1 N δ [ r + r x i ( t ) ] δ [ r x i ( 0 ) ] d r
= 1 N i = 1 N B i ( r + r , t ) B i ( r , 0 ) d r
= 1 ρ i = 1 N B i ( r , t ) B i ( 0 , 0 )
The last passage follows from the uniformity of the fluid. Analogously, one finds
G d ( r , t ) = 1 ρ i = 1 N j i N B i ( r , t ) B j ( 0 , 0 )

References

  1. Angell, C. Relaxation in liquids, polymers and plastic crystals-strong/fragile patterns and problems. J. Non-Cryst. Sol. 1991, 131–133, 13–31. [Google Scholar] [CrossRef]
  2. Debenedetti, P.G.; Stillinger, F.H. Supercooled liquids and the glass transition. Nature 2001, 410, 259–267. [Google Scholar] [CrossRef]
  3. Ngai, K.L. Relaxation and Diffusion in Complex Systems; Springer: Berlin, Germany, 2011. [Google Scholar]
  4. Sillescu, H. Heterogeneity at the glass transition: A review. J. Non-Cryst. Solids 1999, 243, 81–108. [Google Scholar] [CrossRef]
  5. Ngai, K.L. Why the fast relaxation in the picosecond to nanosecond time range can sense the glass transition. Philos. Mag. 2004, 84, 1341–1353. [Google Scholar] [CrossRef]
  6. Berthier, L.; Biroli, G. Theoretical perspective on the glass transition and amorphous materials. Rev. Mod. Phys. 2011, 83, 587–645. [Google Scholar] [CrossRef]
  7. Ashcroft, N.W.; Mermin, N.D. Solid State Physics; Holt, Rinehart and Wiston: New York, NY, USA, 1976. [Google Scholar]
  8. Simmons, D.S.; Cicerone, M.T.; Zhong, Q.; Tyagic, M.; Douglas, J.F. Generalized localization model of relaxation in glass-forming liquids. Soft Matter 2012, 8, 11455–11461. [Google Scholar] [CrossRef]
  9. Pazmiño Betancourt, B.A.; Hanakata, P.Z.; Starr, F.W.; Douglas, J.F. Quantitative relations between cooperative motion, emergent elasticity, and free volume in model glass-forming polymer materials. Proc. Natl. Acad. Sci. USA 2015, 112, 2966–2971. [Google Scholar] [CrossRef]
  10. Douglas, J.F.; Pazmiño Betancourt, B.A.; Tong, X.; Zhang, H. Localization model description of diffusion and structural relaxation in glass-forming Cu–Zr alloys. J. Stat. Mech. Theory Exp. 2016, 2016, 054048. [Google Scholar] [CrossRef]
  11. Ediger, M.D. Spatially heterogeneous dynamics in supercooled liquids. Annu. Rev. Phys. Chem. 2000, 51, 99–128. [Google Scholar] [CrossRef]
  12. Richert, R. Heterogeneous dynamics in liquids: Fluctuations in space and time. J. Phys. Condens. Matter 2002, 14, R703–R738. [Google Scholar] [CrossRef]
  13. Tarjus, G.; Kivelson, D. Breakdown of the Stokes–Einstein relation in supercooled liquids. J. Chem. Phys. 1995, 103, 3071–3073. [Google Scholar] [CrossRef]
  14. Tracht, U.; Wilhelm, M.; Heuer, A.; Feng, H.; Schmidt-Rohr, K.; Spiess, H.W. Length Scale of Dynamic Heterogeneities at the Glass Transition Determined by Multidimensional Nuclear Magnetic Resonance. Phys. Rev. Lett. 1998, 81, 2727–2730. [Google Scholar] [CrossRef]
  15. Douglas, J.; Leporini, D. Obstruction model of the fractional Stokes–Einstein relation in glass-forming liquids. J. Non-Cryst. Solids 1998, 235–237, 137–141. [Google Scholar] [CrossRef]
  16. Jung, Y.; Garrahan, J.P.; Chandler, D. Excitation lines and the breakdown of Stokes–Einstein relations in supercooled liquids. Phys. Rev. E 2004, 69, 061205. [Google Scholar] [CrossRef]
  17. Adam, G.; Gibbs, J.H. On the Temperature Dependence of Cooperative Relaxation Properties in Glass-Forming Liquids. J. Chem. Phys. 1965, 43, 139–146. [Google Scholar] [CrossRef]
  18. Starr, F.W.; Douglas, J.F.; Sastry, S. The relationship of dynamical heterogeneity to the Adam-Gibbs and random first-order transition theories of glass formation. J. Chem. Phys. 2013, 138, 12A541. [Google Scholar] [CrossRef]
  19. Karmakar, S.; Dasgupta, C.; Sastry, S. Length scales in glass-forming liquids and related systems: A review. Rep. Prog. Phys. 2015, 79, 016601. [Google Scholar] [CrossRef]
  20. Albert, S.; Bauer, T.; Michl, M.; Biroli, G.; Bouchaud, J.P.; Loidl, A.; Lunkenheimer, P.; Tourbot, R.; Wiertel-Gasquet, C.; Ladieu, F. Fifth-order susceptibility unveils growth of thermodynamic amorphous order in glass-formers. Science 2016, 352, 1308–1311. [Google Scholar] [CrossRef]
  21. Biroli, G.; Karmakar, S.; Procaccia, I. Comparison of Static Length Scales Characterizing the Glass Transition. Phys. Rev. Lett. 2013, 111, 165701. [Google Scholar] [CrossRef]
  22. Wyart, M.; Cates, M.E. Does a Growing Static Length Scale Control the Glass Transition? Phys. Rev. Lett. 2017, 119, 195501. [Google Scholar] [CrossRef]
  23. Andreozzi, L.; Schino, A.D.; Giordano, M.; Leporini, D. Evidence of a fractional DebyeStokesEinstein law in supercooled oterphenyl. Europhys. Lett. 1997, 38, 669–674. [Google Scholar] [CrossRef]
  24. De Michele, C.; Leporini, D. Viscous flow and jump dynamics in molecular supercooled liquids. II. Rotations. Phys. Rev. E 2001, 63, 036702. [Google Scholar] [CrossRef] [PubMed]
  25. Tyrrell, H.J.V.; Harris, K.R. Diffusion in Liquids; Butterworths: London, UK, 1984. [Google Scholar]
  26. Hansen, J.P.; McDonald, I.R. Theory of Simple Liquids, 3rd ed.; Academic Press: New York, NY, USA, 2006. [Google Scholar]
  27. De Michele, C.; Leporini, D. Viscous flow and jump dynamics in molecular supercooled liquids. I. Translations. Phys. Rev. E 2001, 63, 036701. [Google Scholar] [CrossRef] [PubMed]
  28. Lad, K.N.; Jakse, N.; Pasturel, A. Signatures of fragile-to-strong transition in a binary metallic glass-forming liquid. J. Chem. Phys. 2012, 136, 104509. [Google Scholar] [CrossRef] [PubMed]
  29. Puosi, F.; Leporini, D. Communication: Fast and local predictors of the violation of the Stokes- Einstein law in polymers and supercooled liquids. J. Chem. Phys. 2012, 136, 211101. [Google Scholar] [CrossRef] [PubMed]
  30. Puosi, F.; Leporini, D. Communication: Fast dynamics perspective on the breakdown of the Stokes–Einstein law in fragile glass-formers. J. Chem. Phys. 2018, 148, 131102. [Google Scholar] [CrossRef]
  31. Chang, I.; Fujara, F.; Geil, B.; Heuberger, G.; Mangel, T.; Sillescu, H. Translational and rotational molecular motion in supercooled liquids studied by NMR and forced Rayleigh scattering. J. Non-Cryst. Solids 1994, 172–175, 248–255. [Google Scholar] [CrossRef]
  32. Puosi, F.; Leporini, D. Communication: Correlation of the instantaneous and the intermediate-time elasticity with the structural relaxation in glass-forming systems. J. Chem. Phys. 2012, 136, 041104. [Google Scholar] [CrossRef]
  33. Tobolsky, A.; Powell, R.E.; Eyring, H. Elastic-viscous properties of matter. In Frontiers in Chemistry; Burk, R.E., Grummit, O., Eds.; Interscience: New York, NY, USA, 1943; Volume 1, pp. 125–190. [Google Scholar]
  34. Hall, R.W.; Wolynes, P.G. The aperiodic crystal picture and free energy barriers in glasses. J. Chem. Phys. 1987, 86, 2943–2948. [Google Scholar] [CrossRef]
  35. Angell, C.A. Formation of glasses from liquids and biopolymers. Science 1995, 267, 1924–1935. [Google Scholar] [CrossRef]
  36. Dyre, J.C.; Olsen, N.B.; Christensen, T. Local elastic expansion model for viscous-flow activation energies of glass-forming molecular liquids. Phys. Rev. B 1996, 53, 2171–2174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Martinez, L.M.; Angell, C.A. A thermodynamic connection to the fragility of glass-forming liquids. Nature 2001, 410, 663–667. [Google Scholar] [CrossRef] [PubMed]
  38. Sastry, S. The relationship between fragility, configurational entropy and the potential energy landscape of glass-forming liquids. Nature 2001, 409, 164–167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Ngai, K.L. Dynamic and thermodynamic properties of glass-forming substances. J. Non-Cryst. Solids 2000, 275, 7–51. [Google Scholar] [CrossRef]
  40. Xia, X.; Wolynes, P.G. Fragilities of liquids predicted from the random first order transition theory of glasses. Proc. Natl. Acad. Sci. USA 2000, 97, 2990–2994. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Buchenau, U.; Zorn, R. A relation between fast and slow motions in glassy and liquid selenium. Europhys. Lett. 1992, 18, 523–528. [Google Scholar] [CrossRef]
  42. Starr, F.; Sastry, S.; Douglas, J.F.; Glotzer, S. What do we learn from the local geometry of glass-forming liquids? Phys. Rev. Lett. 2002, 89, 125501. [Google Scholar] [CrossRef] [Green Version]
  43. Bordat, P.; Affouard, F.; Descamps, M.; Ngai, K.L. Does the interaction potential determine both the fragility of a liquid and the vibrational properties of its glassy state? Phys. Rev. Lett. 2004, 93, 105502. [Google Scholar] [CrossRef] [Green Version]
  44. Widmer-Cooper, A.; Harrowell, P. Predicting the Long-Time Dynamic Heterogeneity in a Supercooled Liquid on the Basis of Short-Time Heterogeneities. Phys. Rev. Lett. 2006, 96, 185701. [Google Scholar] [CrossRef] [Green Version]
  45. Zhang, H.; Srolovitz, D.J.; Douglas, J.F.; Warren, J.A. Grain boundaries exhibit the dynamics of glass-forming liquids. Proc. Natl. Acad. Sci. USA 2009, 106, 7735–7740. [Google Scholar] [CrossRef] [Green Version]
  46. Widmer-Cooper, A.; Perry, H.; Harrowell, P.; Reichman, D.R. Irreversible reorganization in a supercooled liquid originates from localized soft modes. Nat. Phys. 2008, 4, 711–715. [Google Scholar] [CrossRef] [Green Version]
  47. Larini, L.; Ottochian, A.; De Michele, C.; Leporini, D. Universal scaling between structural relaxation and vibrational dynamics in glass-forming liquids and polymers. Nat. Phys. 2008, 4, 42–45. [Google Scholar] [CrossRef]
  48. Ottochian, A.; De Michele, C.; Leporini, D. Universal divergenceless scaling between structural relaxation and caged dynamics in glass-forming systems. J. Chem. Phys. 2009, 131, 224517. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Puosi, F.; Leporini, D. Scaling between Relaxation, Transport, and Caged Dynamics in Polymers: From Cage Restructuring to Diffusion. J. Phys. Chem. B 2011, 115, 14046–14051. [Google Scholar] [CrossRef]
  50. Ottochian, A.; Leporini, D. Scaling between structural relaxation and caged dynamics in Ca0.4K0.6(NO3)1.4 and glycerol: Free volume, time scales and implications for the pressure-energy correlations. Philos. Mag. 2011, 91, 1786–1795. [Google Scholar] [CrossRef] [Green Version]
  51. Ottochian, A.; Leporini, D. Universal scaling between structural relaxation and caged dynamics in glass-forming systems: Free volume and time scales. J. Non-Cryst. Solids 2011, 357, 298–301. [Google Scholar] [CrossRef]
  52. De Michele, C.; Del Gado, E.; Leporini, D. Scaling between Structural Relaxation and Particle Caging in a Model Colloidal Gel. Soft Matter 2011, 7, 4025–4031. [Google Scholar] [CrossRef] [Green Version]
  53. Puosi, F.; Leporini, D. Spatial displacement correlations in polymeric systems. J. Chem. Phys. 2012, 136, 164901. [Google Scholar] [CrossRef]
  54. Puosi, F.; Leporini, D. Erratum: “Spatial displacement correlations in polymeric systems” [J. Chem. Phys. 136, 164901 (2012)]. J. Chem. Phys. 2013, 139, 029901. [Google Scholar] [CrossRef]
  55. Puosi, F.; Michele, C.D.; Leporini, D. Scaling between relaxation, transport and caged dynamics in a binary mixture on a per-component basis. J. Chem. Phys. 2013, 138, 12A532. [Google Scholar] [CrossRef]
  56. Ottochian, A.; Puosi, F.; Michele, C.D.; Leporini, D. Comment on “Generalized localization model of relaxation in glass-forming liquids”. Soft Matter 2013, 9, 7890–7891. [Google Scholar] [CrossRef]
  57. Novikov, V.N.; Sokolov, A.P. Role of Quantum Effects in the Glass Transition. Phys. Rev. Lett. 2013, 110, 065701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Puosi, F.; Chulkin, O.; Bernini, S.; Capaccioli, S.; Leporini, D. Thermodynamic scaling of vibrational dynamics and relaxation. J. Chem. Phys. 2016, 145, 234904. [Google Scholar] [CrossRef] [PubMed]
  59. Guillaud, E.; Joly, L.; de Ligny, D.; Merabia, S. Assessment of elastic models in supercooled water: A molecular dynamics study with the TIP4P/2005f force field. J. Chem. Phys. 2017, 147, 014504. [Google Scholar] [CrossRef] [Green Version]
  60. Horstmann, R.; Vogel, M. Common behaviours associated with the glass transitions of water-like models. J. Chem. Phys. 2017, 147, 034505. [Google Scholar] [CrossRef]
  61. Becchi, M.; Giuntoli, A.; Leporini, D. Molecular layers in thin supported films exhibit the same scaling as the bulk between slow relaxation and vibrational dynamics. Soft Matter 2018, 14, 8814–8820. [Google Scholar] [CrossRef] [Green Version]
  62. Bernini, S.; Puosi, F.; Leporini, D. Thermodynamic scaling of relaxation: Insights from anharmonic elasticity. J. Phys. Condens. Matter 2017, 29, 135101. [Google Scholar] [CrossRef]
  63. Puosi, F.; Leporini, D. The kinetic fragility of liquids as manifestation of the elastic softening. Eur. Phys. J. E 2015, 38, 87. [Google Scholar] [CrossRef] [Green Version]
  64. Bernini, S.; Leporini, D. Cage effect in supercooled molecular liquids: Local anisotropies and collective solid-like response. J. Chem. Phys. 2016, 144, 144505. [Google Scholar] [CrossRef] [Green Version]
  65. Allen, M.P.; Tildesley, D.J. Computer Simulations of Liquids; Oxford University Press: New York, NY, USA, 1987. [Google Scholar]
  66. Doi, M.; Edwards, S.F. The Theory of Polymer Dynamics; Clarendon Press: Oxford, UK, 1988. [Google Scholar]
  67. Boon, J.P.; Yip, S. Molecular Hydrodynamics; Dover Publications: New York, NY, USA, 1980. [Google Scholar]
  68. Gedde, U.W. Polymer Physics; Chapman and Hall: London, UK, 1995. [Google Scholar]
  69. Kröger, M. Simple models for complex nonequilibrium fluids. Phys. Rep. 2004, 390, 453–551. [Google Scholar] [CrossRef]
  70. 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. [Google Scholar] [CrossRef] [Green Version]
  71. Baschnagel, J.; Varnik, F. Computer simulations of supercooled polymer melts in the bulk and in confined geometry. J. Phys. Condens. Matter 2005, 17, R851–R953. [Google Scholar] [CrossRef]
  72. Bennemann, C.; Donati, C.; Baschnagel, J.; Glotzer, S.C. Growing range of correlated motion in a polymer melt on cooling towards the glass transition. Nature 1999, 399, 246–249. [Google Scholar] [CrossRef]
  73. Ottochian, A.; De Michele, C.; Leporini, D. Non-Gaussian effects in the cage dynamics of polymers. Philos. Mag. 2008, 88, 4057–4062. [Google Scholar] [CrossRef] [Green Version]
  74. Tölle, A. Neutron scattering studies of the model glass former ortho -terphenyl. Rep. Prog. Phys. 2001, 64, 1473. [Google Scholar] [CrossRef] [Green Version]
  75. Chang, I.; Sillescu, H. Heterogeneity at the Glass Transition: Translational and Rotational Self-Diffusion. J. Phys. Chem. B 1997, 101, 8794–8801. [Google Scholar] [CrossRef]
  76. Mallamace, F.; Branca, C.; Corsaro, C.; Leone, N.; Spooren, J.; Chen, S.H.; Stanley, H.E. Transport properties of glass-forming liquids suggest that dynamic crossover temperature is as important as the glass transition temperature. Proc. Natl. Acad. Sci. USA 2010, 107, 22457–22462. [Google Scholar] [CrossRef] [Green Version]
  77. Donati, C.; Glotzer, S.C.; Poole, P.H. Growing Spatial Correlations of Particle Displacements in a Simulated Liquid on Cooling toward the Glass Transition. Phys. Rev. Lett. 1999, 82, 5064–5067. [Google Scholar] [CrossRef] [Green Version]
  78. Doliwa, B.; Heuer, A. Cooperativity and spatial correlations near the glass transition: Computer simulation results for hard spheres and disks. Phys. Rev. E 2000, 61, 6898–6908. [Google Scholar] [CrossRef]
  79. Weeks, E.R.; Crocker, J.C.; Weitz, D.A. Short and long-range correlated motion observed in colloidal glasses and liquids. J. Phys. Condens. Matter 2007, 19, 205131. [Google Scholar] [CrossRef]
  80. Narumi, T.; Tokuyama, M. Simulation study of spatial‚Äìtemporal correlation functions in supercooled liquids. Philos. Mag. 2008, 88, 4169–4175. [Google Scholar] [CrossRef]
  81. Kremer, K.; Grest, G.S. Dynamics of entangled linear polymer melts: A molecular-dynamics simulation. J. Chem. Phys. 1990, 92, 5057–5086. [Google Scholar] [CrossRef]
  82. Paul, W.; Smith, G.D. Structure and dynamics of amorphous polymers: Computer simulations compared to experiment and theory. Rep. Prog. Phys. 2004, 67, 1117–1185. [Google Scholar] [CrossRef]
  83. Luo, C.; Sommer, J.U. Coding coarse grained polymer model for LAMMPS and its application to polymer crystallization. Comput. Phys. Commun. 2009, 180, 1382–1391. [Google Scholar] [CrossRef]
Figure 1. Monomer arrangements at a time t 0 of two states of a polymer melt with (a) homogeneous and (b) heterogeneous dynamics. Bonds are removed for clarity reasons. Particles are coloured according to their squared displacements in the time interval [ t 0 τ α , t 0 ] . Bright yellow particles have squared displacements no less than 1. Notice that the two states have comparable mean square displacement (∼0.21, homogeneous state; ∼0.28, heterogeneous state) but rather different relaxation times τ α (∼9, homogeneous state; ∼1550, heterogeneous state). Homogeneous, i.e., position-independent, dynamics of the monomers is an aspect of systems with fast relaxation. Conversely, in the presence of heterogeneous dynamics, clusters of particles with extremely high mobility coexist with nearly immobile ones, slowing down the relaxation.
Figure 1. Monomer arrangements at a time t 0 of two states of a polymer melt with (a) homogeneous and (b) heterogeneous dynamics. Bonds are removed for clarity reasons. Particles are coloured according to their squared displacements in the time interval [ t 0 τ α , t 0 ] . Bright yellow particles have squared displacements no less than 1. Notice that the two states have comparable mean square displacement (∼0.21, homogeneous state; ∼0.28, heterogeneous state) but rather different relaxation times τ α (∼9, homogeneous state; ∼1550, heterogeneous state). Homogeneous, i.e., position-independent, dynamics of the monomers is an aspect of systems with fast relaxation. Conversely, in the presence of heterogeneous dynamics, clusters of particles with extremely high mobility coexist with nearly immobile ones, slowing down the relaxation.
Ijms 20 05708 g001
Figure 2. Monomer dynamics in the polymer melt. (a) Mean square displacement (MSD) for polymers in selected states (see below for details). For clarity reasons, MSDs are multiplied by indicated factors. Inset: corresponding MSD slope Δ ( t ) , Equation (16); the position of the minimum at t = 1.0 ( 4 ) is signalled by the arrow in the inset and the dashed line in the main panel. (b) corresponding ISF curves. The figure shows a manifestation of Equation (3), see Section 5.1 for details, i.e., if states have equal DW factor u 2 , both the MSD and ISF curves coincide at least in the time window [ t , τ α ] ( τ α is marked with dots on each curve). The physical states are labelled by the string (M, ρ , T, q, p) where M is the number of monomers per chain, ρ the number density, T the temperature and the pair ( q , p ) refers to the characteristic parameters of the non-bonding potential, Equation (25). The six sets of states are as follows. Set A: (2,1.086,0.7,7,6), (3,1.086,0.7,7,6), (10,1.086,0.7,7,6), (10,1.033,0.7,8,6). Set B: (2,1.033,0.7,10,6), (3,1.039,0.7,11,6), (3,1.041,0.7,11,6). Set C: (2,1.033,0.5,10,6), (3,1.056,0.7,12,6), (5,1.033,0.6,12,6), (10,1.056,0.7,12,6). Set D: (3,1.086,0.7,12,6), (5,1.086,0.7,12,6), (10,1.086,0.7,12,6). Set E: (2,1.0,0.7,12,11), (3,1.1,1.1,15,7). Data from [48].
Figure 2. Monomer dynamics in the polymer melt. (a) Mean square displacement (MSD) for polymers in selected states (see below for details). For clarity reasons, MSDs are multiplied by indicated factors. Inset: corresponding MSD slope Δ ( t ) , Equation (16); the position of the minimum at t = 1.0 ( 4 ) is signalled by the arrow in the inset and the dashed line in the main panel. (b) corresponding ISF curves. The figure shows a manifestation of Equation (3), see Section 5.1 for details, i.e., if states have equal DW factor u 2 , both the MSD and ISF curves coincide at least in the time window [ t , τ α ] ( τ α is marked with dots on each curve). The physical states are labelled by the string (M, ρ , T, q, p) where M is the number of monomers per chain, ρ the number density, T the temperature and the pair ( q , p ) refers to the characteristic parameters of the non-bonding potential, Equation (25). The six sets of states are as follows. Set A: (2,1.086,0.7,7,6), (3,1.086,0.7,7,6), (10,1.086,0.7,7,6), (10,1.033,0.7,8,6). Set B: (2,1.033,0.7,10,6), (3,1.039,0.7,11,6), (3,1.041,0.7,11,6). Set C: (2,1.033,0.5,10,6), (3,1.056,0.7,12,6), (5,1.033,0.6,12,6), (10,1.056,0.7,12,6). Set D: (3,1.086,0.7,12,6), (5,1.086,0.7,12,6), (10,1.086,0.7,12,6). Set E: (2,1.0,0.7,12,11), (3,1.1,1.1,15,7). Data from [48].
Ijms 20 05708 g002
Figure 3. Correlation function of the end-to-end vector joining the two ends of a polymer chain. Each group of curves corresponds the physical states A, ..., E with identical DW factor detailed in Figure 2. Polymer states contributing to one cluster of scaled curves have not necessarily equal chain length. However, the scaled time removes the chain length dependence. Dots mark the time 4 τ e e / M 2 . The results prove that Equation (3) holds also at times τ e e much longer than τ α . Data from [48].
Figure 3. Correlation function of the end-to-end vector joining the two ends of a polymer chain. Each group of curves corresponds the physical states A, ..., E with identical DW factor detailed in Figure 2. Polymer states contributing to one cluster of scaled curves have not necessarily equal chain length. However, the scaled time removes the chain length dependence. Dots mark the time 4 τ e e / M 2 . The results prove that Equation (3) holds also at times τ e e much longer than τ α . Data from [48].
Ijms 20 05708 g003
Figure 4. The structural relaxation time τ α and the scaled reorientation time τ e e of the polymer chains vs. the DW factor u 2 . Empty circles highlight the cases plotted in Figure 2. The other states are detailed in Ref. [48]. The dashed line across the τ α curve is Equation (A1). The dashed curve across the chain reorientation time curve is a guide for the eyes. Data from [48].
Figure 4. The structural relaxation time τ α and the scaled reorientation time τ e e of the polymer chains vs. the DW factor u 2 . Empty circles highlight the cases plotted in Figure 2. The other states are detailed in Ref. [48]. The dashed line across the τ α curve is Equation (A1). The dashed curve across the chain reorientation time curve is a guide for the eyes. Data from [48].
Ijms 20 05708 g004
Figure 5. (a) Self part of the van Hove function G s ( r , t ) of the states of Figure 2 at the structural relaxation time t = τ α . The curves are multiplied by indicated factors. The sets of clustered curves A–E show that, if states have equal DW factor, they have coincident van Hove functions too. As G s ( r , t ) may be expressed in terms of correlation functions, the coincidence reflects Equation (3). Data from [49]. (b) The ratio N s ( r , τ α ) , Equation (21), of the states of Figure 2. On increasing the structural relaxation time from A states to E states, the system tends to increase the fractions of monomers with either much lower or much higher mobility with respect to the fraction predicted by the Gaussian approximation. Data from [53].
Figure 5. (a) Self part of the van Hove function G s ( r , t ) of the states of Figure 2 at the structural relaxation time t = τ α . The curves are multiplied by indicated factors. The sets of clustered curves A–E show that, if states have equal DW factor, they have coincident van Hove functions too. As G s ( r , t ) may be expressed in terms of correlation functions, the coincidence reflects Equation (3). Data from [49]. (b) The ratio N s ( r , τ α ) , Equation (21), of the states of Figure 2. On increasing the structural relaxation time from A states to E states, the system tends to increase the fractions of monomers with either much lower or much higher mobility with respect to the fraction predicted by the Gaussian approximation. Data from [53].
Ijms 20 05708 g005
Figure 6. Non-Gaussian parameters (NGPs) of states with different relaxation times τ α (marked with grey dots). The physical states A, ⋯, E are the states with identical DW factor detailed in Figure 2. Note that they have coinciding NGPs in the time window [ t , τ α ] at least, in agreement with Equation (3). The curve labelled as F is the state (M, ρ , T, q, p) = (3, 1.1, 0.65, 12, 6) with τ α 2 · 10 3 , see Figure 4. Inset: the NGP maximum α 2 m a x vs. the ratio R, Equation (9). The dot with the largest α 2 m a x value corresponds to the state with the longest structural relaxation time τ α in Figure 4 with parameters (M, ρ , T, q, p) = (3, 1.2, 0.95, 6,12). Data from [29].
Figure 6. Non-Gaussian parameters (NGPs) of states with different relaxation times τ α (marked with grey dots). The physical states A, ⋯, E are the states with identical DW factor detailed in Figure 2. Note that they have coinciding NGPs in the time window [ t , τ α ] at least, in agreement with Equation (3). The curve labelled as F is the state (M, ρ , T, q, p) = (3, 1.1, 0.65, 12, 6) with τ α 2 · 10 3 , see Figure 4. Inset: the NGP maximum α 2 m a x vs. the ratio R, Equation (9). The dot with the largest α 2 m a x value corresponds to the state with the longest structural relaxation time τ α in Figure 4 with parameters (M, ρ , T, q, p) = (3, 1.2, 0.95, 6,12). Data from [29].
Ijms 20 05708 g006
Figure 7. (a) The product D M τ α vs. the ratio R, Equation (9) (b) the same product vs. α 2 m a x , the maximum of the non-Gaussian parameter, Equation (22). The onset of the Stokes–Einstein (SE) violation for α 2 m a x > α 2 , c m a x and R > R c , respectively, is indicated with the full vertical lines (uncertainty marked by dashed lines). The thick line in the panel (a) is the master curve between log τ α and the DW factor, Equation (A1), recast in terms of R and the thin line, is the corresponding linear approximation for small R values. Note that the SE violation is apparent where the linear approximation is poor. Data from [29].
Figure 7. (a) The product D M τ α vs. the ratio R, Equation (9) (b) the same product vs. α 2 m a x , the maximum of the non-Gaussian parameter, Equation (22). The onset of the Stokes–Einstein (SE) violation for α 2 m a x > α 2 , c m a x and R > R c , respectively, is indicated with the full vertical lines (uncertainty marked by dashed lines). The thick line in the panel (a) is the master curve between log τ α and the DW factor, Equation (A1), recast in terms of R and the thin line, is the corresponding linear approximation for small R values. Note that the SE violation is apparent where the linear approximation is poor. Data from [29].
Ijms 20 05708 g007
Figure 8. Stokes–Einstein product K S E , normalised by its high temperature value K 0 ( τ α 1 ps), as a function of the reduced DW factor u 2 / u g 2 , u g 2 being the DW factor at the glass transition. In addition to unentangled polymers, the plot also considers MD data concerning atomic binary mixtures (atoms labelled as A and B) and metallic alloys made by Cu and Zr atoms, as well as experimental data for ortho-terphenyl (OTP) [74,75]. Two predictions of the master curve are presented in terms of the quantity K ^ S E and K ˜ S E , Equations (A8) and (A10), respectively. Both quantities have no adjustable parameters. K ^ 0 and K ˜ 0 are suitable constants to ensure the unit limit value at large u 2 / u g 2 . A third master curve, drawn from the fractional SE law τ α 1 κ with κ = 0.85 (orange curve), is superimposed to the other curves. For numerical data, u g 2 is obtained according to the procedure outlined in [47]. Data from [30].
Figure 8. Stokes–Einstein product K S E , normalised by its high temperature value K 0 ( τ α 1 ps), as a function of the reduced DW factor u 2 / u g 2 , u g 2 being the DW factor at the glass transition. In addition to unentangled polymers, the plot also considers MD data concerning atomic binary mixtures (atoms labelled as A and B) and metallic alloys made by Cu and Zr atoms, as well as experimental data for ortho-terphenyl (OTP) [74,75]. Two predictions of the master curve are presented in terms of the quantity K ^ S E and K ˜ S E , Equations (A8) and (A10), respectively. Both quantities have no adjustable parameters. K ^ 0 and K ˜ 0 are suitable constants to ensure the unit limit value at large u 2 / u g 2 . A third master curve, drawn from the fractional SE law τ α 1 κ with κ = 0.85 (orange curve), is superimposed to the other curves. For numerical data, u g 2 is obtained according to the procedure outlined in [47]. Data from [30].
Ijms 20 05708 g008
Figure 9. Radial dependence of the correlation of the direction (a) and the mobility (b) displacements occurring in a time range as wide as the structural relaxation time τ α . For comparison, the radial distribution function g ( r ) (dashed line) of the state with { M = 2 , ρ = 1.086 , T = 0.7 , q = 7 , p = 6 } is plotted. Note that g ( r ) is virtually state-independent. The insets are semi-log plots of the corresponding main panels. Note the approximate exponential decay of the peak amplitudes with slopes ξ u ( τ α ) and ξ δ u ( τ α ) , respectively. Data from [53].
Figure 9. Radial dependence of the correlation of the direction (a) and the mobility (b) displacements occurring in a time range as wide as the structural relaxation time τ α . For comparison, the radial distribution function g ( r ) (dashed line) of the state with { M = 2 , ρ = 1.086 , T = 0.7 , q = 7 , p = 6 } is plotted. Note that g ( r ) is virtually state-independent. The insets are semi-log plots of the corresponding main panels. Note the approximate exponential decay of the peak amplitudes with slopes ξ u ( τ α ) and ξ δ u ( τ α ) , respectively. Data from [53].
Ijms 20 05708 g009
Figure 10. The direction ξ u ( τ α ) (full symbols) and modulus ξ δ u ( τ α ) (open symbols) correlation lengths vs. the structural relaxation time τ α of selected set of states of Figure 2. Dashed lines are guides for the eyes. States with equal DW factor, i.e., belonging to the same set B, ⋯, E exhibit equal directional and mobility correlation lengths. Data from [53].
Figure 10. The direction ξ u ( τ α ) (full symbols) and modulus ξ δ u ( τ α ) (open symbols) correlation lengths vs. the structural relaxation time τ α of selected set of states of Figure 2. Dashed lines are guides for the eyes. States with equal DW factor, i.e., belonging to the same set B, ⋯, E exhibit equal directional and mobility correlation lengths. Data from [53].
Ijms 20 05708 g010

Share and Cite

MDPI and ACS Style

Puosi, F.; Tripodo, A.; Leporini, D. Fast Vibrational Modes and Slow Heterogeneous Dynamics in Polymers and Viscous Liquids. Int. J. Mol. Sci. 2019, 20, 5708. https://doi.org/10.3390/ijms20225708

AMA Style

Puosi F, Tripodo A, Leporini D. Fast Vibrational Modes and Slow Heterogeneous Dynamics in Polymers and Viscous Liquids. International Journal of Molecular Sciences. 2019; 20(22):5708. https://doi.org/10.3390/ijms20225708

Chicago/Turabian Style

Puosi, Francesco, Antonio Tripodo, and Dino Leporini. 2019. "Fast Vibrational Modes and Slow Heterogeneous Dynamics in Polymers and Viscous Liquids" International Journal of Molecular Sciences 20, no. 22: 5708. https://doi.org/10.3390/ijms20225708

APA Style

Puosi, F., Tripodo, A., & Leporini, D. (2019). Fast Vibrational Modes and Slow Heterogeneous Dynamics in Polymers and Viscous Liquids. International Journal of Molecular Sciences, 20(22), 5708. https://doi.org/10.3390/ijms20225708

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop