arxiv: v2 [cond-mat.soft] 28 Sep 2014

Transcription

1 Disordered Solids Without Well-Defined Transverse Phonons: The Nature of Hard Shere Glasses Xieng Wang, Wen Zheng, Lijin Wang, and Ning Xu CAS Key Laboratory of Soft Matter Chemistry, Hefei National Laboratory for Physical Sciences at the Microscale, and Deartment of Physics, University of Science and Technology of China, Hefei 2326, Peole s Reublic of China. (Dated: Setember 3, 24) arxiv:49.596v2 [cond-mat.soft] 28 Se 24 We robe the Ioffe-Regel limits of glasses with reulsions near the zero-temerature jamming transition by measuring the dynamical structure factors. At zero temerature, the transverse Ioffe- Regel frequency vanishes at the jamming transition with a diverging length, but the longitudinal one does not, which excludes the existence of a diverging length associated with the longitudinal excitations. At low temeratures, the transverse and longitudinal Ioffe-Regel frequencies aroach zero at the jamming-like transition and glass transition, resectively. As a consequence, glasses between the glass transition and jamming-like transition, which are hard shere glasses in the low temerature limit, can only carry well-defined longitudinal honons and have an oosite ressure deendence of the ratio of the shear modulus to the bulk modulus from glasses beyond the jamminglike transition. PACS numbers: 63.5.Lm,64.7.v,6.43.Bn Uon comression, colloidal systems undergo the glass transition when the relaxation time (or the viscosity) exceeds the measurable value [ 3]. In the absence of the thermal energy, the comression leads to the jamming transition of ackings of reulsive articles with the sudden formation of rigidity and static force networks [4, 5]. Although the initial jamming hase diagram [4] rooses that the glass transition of urely reulsive systems collases with the jamming transition in the zero temerature (T = ) limit, or equivalently in the hard shere limit [6], recent studies have evidenced that the T = or hard shere glass transition haens at a acking fraction φ g lower than the critical value φ j of the T = jamming transition at the so-called oint J [7 2]. The dearture of φ g from φ j leaves φ (φ g,φ j ) a secial region. At T =, the systems in this region are unjammed and unable to sustain the shear or comression. However, they are by definition glasses when being thermally excited, because articles are unable to diffuse freely. While glasses are believed to be mechanically rigid solids, it seems erlexing that the systems in (φ g,φ j ) are rigid in the glass ersective but not in the T = jamming ersective, which makes T = singular. For soft and reulsive sheres, recent simulations have indicated that the glass transition temerature T g is roortional to the ressure and vanishes at φ g [6 9]. Both colloidal exeriments and simulations have also shown that inside the glass regime, at fixed temerature, there exists a crossover ressure j at which the first eak of the air distribution function reaches the maximum height g max [7, 8, 3 5], reminiscing the structural signature of the T = jamming transition [6]. At such a crossover, the ressure deendence of the temerature T j () is different from T g (). T j () is lower than T g () and vanishes at φ j [7, 8]. In other words, when comressed at a fixed low temerature, the system undergoes the glass transition first and then the emergence of g max. The whole icture is reroduced here as well in Fig. 3. Arecentstudyhasdemonstratedfurtherthatat > j and at a fixed low temerature some tyical scaling relations observed in the T = jammed solids are recovered, while they break down at < j, so it has been suggested that the emergence of g max at j signifies the jamming-like transition [8]. Solid roerties of glasses on both sides of the jamming-like transition must be distinct. For instance, the thermal energy softens jammed solids, but it hardens the systems that are unjammed at T = by increasingthe yield stresswith the temerature [2]. We will show that the underlying mechanism of such distinctions challenges our understanding of glasses with reulsions, esecially the nature of hard shere glasses. By measuring the dynamical structure factors of soft glasses with reulsions, we obtain both the transverse and longitudinal Ioffe-Regel frequencies, i.e., the critical frequencies at which the honon wavelength is comarable to the mean free ath [7]. At T =, the transverseioffe-regelfrequencyir T vanishesatointj, while the longitudinal one IR L remains nonzero. Our analysis of the low-temerature glasses suggests that IR T and IR L vanish with a diverging length at the jamming-like transition and glass transition, resectively. Therefore, glasses below the jamming-like transition do not have well-defined transverse honons, which should be the intrinsic cause of their significantly different material roerties from glasses above the jamming-like transition. Our systems are L L boxes in two dimensions with eriodic boundary conditions. To avoid crystallization, we ut in the box N bidiserse disks with an equal mass m and a diameter ratio.4. If not secified, the results shown are for N = 496 systems. Particles i and j interact via urely reulsive otential U ij = ǫ α ( r ) α ij Θ( r ) ij, () σ ij σ ij

2 2 where r ij is their searation, σ ij is the sum of their radii, and Θ(x) is the Heaviside function. We show here the results of Hertzian reulsion (α = 5/2), which has been verified as a good aroximation to the interaction of PNIPAM colloids [7]. In our simulations, the length, energy, and mass are in units of small article diameter σ, characteristic energy of the interaction ǫ, and article mass m. The time is in units of mσ 2 /ǫ. The temerature is in units of ǫ/k B with k B the Boltzmann constant. We generate jammed solids at T = by finding the local otential energy minima using the fast inertial relaxation engine minimization algorithm [8]. The normal modes of vibration are obtained by diagonalizing the Hessian matrix using ARPACK [9]. The dynamical structure factors at T = are measured from the modes [2]: S (k,) = k2 m 2 F n, (k)δ( n ), (2) n wherethe sum is overallmodes, denotest (transverse) or L (longitudinal), and F n,l = j F n,t = j ( e n,j ˆk)ex(i k r j ) 2, (3) ( e n,j ˆk)ex(i k r j ) 2, (4) where the sums are over all articles, e n,j is the olarization vector of article j in mode n, r j is the location of article j, and ˆk = k/k with k satisfying the eriodic boundary conditions. For thermal systems, we erform molecular dynamics simulations at constant temerature and ressure. The structure and dynamics are evaluated from the air j i δ(r distribution function g(r) = (4L 2 /N 2 ) i r ij ) and intermediate scattering function F s (q,t) = (2/N) j ex(i q [ r j(t) r j ()]) for large articles, where the sums are over all large articles, q is chosen in the x direction with the amlitude aroximately equal to the value at the first eak of the static structure factor, and. denotes the ensemble average. The relaxation time τ satisfies F s (q,τ) = e F s (q,). By fitting the relaxation time measured at fixed temerature with the Vogel-Fulcher function, τ = τ ex[m/( g )], where τ and M are fitting arameters, we estimate the glass transition ressure g. The dynamical structure factors are measured from the time correlation of current [2, 2]: S (k,) = k2 2πN 2 where J L (k,t) = j J T (k,t) = j dt J (k,t) J ( k,) ex(it), (5) [( v j (t) ˆk)ˆk]ex[i k r j (t)], (6) [ v j (t) ( v j (t) ˆk)ˆk]ex[i k r j (t)], (7) S T (k,) Ω T -/2, πγ T -/2 IR, IR T L k=.95 k=.69 k= (c).. k -/6 (e) S L (k,) Ω L -/2, πγ L -/2 R 6 (b) 4 2 k=.89 k=.2673 k= (d).. k -/3 (f) k -2/9 FIG. : (color online). Analysis of the T = jammed solids. and (b) Examles of the dynamical structure factors S T(k,) and S L(k,) at = 5 4 with the lines fittings with Eq. (8). (c) and (d) Pressure scaled disersion relation Ω /2 (emty symbols) and sound attenuation coefficient πγ /2 (solid symbols) against k /6 for the transverse ( = T) and k /3 for the longitudinal ( = L) excitations at = 4 (circles), 5 4 (squares), 3 (diamonds), and 5 3 (triangles). The solid and dashed lines have a sloe of and 2, resectively. (e) Pressure deendence of the Ioffe-Regel frequencies T IR (circles) and L IR (squares). The solid line is the fitting with T IR /2. (f) Ratio of the Rayleigh comonent to the Brillouin comonent, R, against k 2/9. The symbols are defined as in (c) and (d). The solid line has a sloe of 3/2. where the sums are over all articles, and v j is the velocity of article j. For simlifications, here we use the same symbol S (k,) as that for T = [22]. Figures and (b) show examles of the dynamical structure factors at T = calculated from Eqs. (2)-(4), which can be well fitted with [2] S (k,) = S,B (k,)+s,r (k,) (8) A (k) = [ 2 Ω 2 (k)]2 + 2 Γ 2 (k) + B (k) 2 +(D k 2 ) 2, where S,B (k,) and S,R (k,) are the Brillouin (solid) and Rayleigh (liquid) comonents, A (k) and B (k) are fitting arameters, Ω (k) gives the disersion relation, Γ (k) is the sound attenuation coefficient, and D is the thermal diffusivity. For S T (k,), the Brillouin comonent itselfcanfit the datawell [2, 23], sos T,R (k,). However, for S L (k,), the Rayleigh comonent gives rise to the low-frequency increase with the decrease of the frequency, esecially at large k and low.

3 3 The Rayleigh comonent of the longitudinal dynamical structure factor has been overlooked in most of revious studies of disordered solids. Its resence imlies that disordered solids still exhibit liquid characters at certain length scales. Figure (f) shows the ratio of the Rayleigh contribution to the Brillouin contribution, R(k) = ds L,R (k,)/ ds L,B (k,), for the T = jammed solids at different ressures. At each ressure, R(k) k 3/2. Interestingly, all curves collase when we lot R(k) against k 2/9, which imlies a diverging length at oint J ( = ): l 2/9. This is consistent with the icture of jamming transition: The Rayleigh comonent dominates over infinitely large length scale, so that the jammed solids are on the verge of instability. In Figs. (c) and (d), we illustrate how the Ioffe-Regel limit is determined. As reorted for some model glasses, Ω k, and Γ k 2 [2]. The Ioffe-Regel limit is achieved when Ω (kir ) = πγ (kir ), from which we obtain the Ioffe-Regel frequency IR = Ω (kir ). When > IR, a honon undergoes multile collisions before moving to a distance comarable to its wavelength and is thus ill-defined. At fixed ressure, IR L > T IR. For the T = jammed solids, various characteristic frequencies have been defined, e.g., the onset frequency of the lateau in the density of vibrational states [24 27] or energy diffusivity [28]. It has been shown that these characteristic frequencies are scaled well with the ressure: (φ φ j ) (α )/2 /2 [24 26, 28], so at oint J decays to zero. It has also been argued that relates to the boson eak or Ioffe-Regel frequencies [27, 28]. However, no direct measure of the Ioffe-Regel frequencies has ever been erformed for jammed solids. In Fig. (e), we exlicitly show the ressure deendenceoftheioffe-regelfrequenciesforthe T = jammed solids. For the transverse one, IR T /2, which vanishes at oint J. Because IR T = c TkIR T, where the transverse seed of sound c T = G/ρ (α 3/2)/(2α 2) [5] with G the shear modulus and ρ the mass density, for systems with Hertzian reulsion, IR T /3 kir T. Therefore, there exists a diverging length aroaching oint J: lir T (kt IR ) /6. As verified in Fig. (c), Ω T (k) and πγ T (k) measured at different ressures collase onto the same master curves with the same intersection when we lot Ω T /2 and πγ T /2 against k /6. The longitudinal Ioffe-Regel frequency, however, shows the tendency to aroach a nonzero value at oint J, as shown in Fig. (e). It has been argued that the vanishing of the characteristic frequency at oint J imlies two ossible diverging length scales, associated with the transverse and longitudinal excitations, resectively [24, 25, 29]. Here we show that IR L > at oint J, so the length lir L (kl IR ) associated with the longitudinal excitations remains finite at oint J. The longitudinal seed of sound c L = (B +4G/3)/ρ B (α 2)/(2α 2), because the bulk modulus B G near oint J. If IR L /2, we would observe scaling collase S T (k,) T IR IR k=.844 k=.94 k=.688 S L (k,)..2.3 (c) (e) c 2-5 (b) L IR Ω T, πγ Τ (d) (f) k=.888 k=.267 k= k FIG. 2: (color online). Analysis of glasses at T >. and (b) Examles of the dynamical structure factors S T(k,) and S L(k,) at T = 6 and = 3 with the lines fittings with Eq. (8). (c) and (d) Pressure deendence of the Ioffe- Regel frequencies T IR and L IR at T = 7 (circles), 5 7 (squares), 6 (diamonds), and 5 6 (triangles) with the lines fittings with IR ( c) /2. (e) Scaling collase of IR when lotted against c with the lines having a sloe of /2. The uer and lower branches are L IR and T IR, resectively. (f) Disersion relation Ω T(k) (circles) andsound attenuation coefficient πγ T(k) (squares) for a glass comosed of N = 6384 articles at T = 6 and = 4 in ( g, j) (located by the star in Fig. 3). The lines are to guide the eye. of Ω L (k) and Γ L (k) by lotting Ω L /2 and πγ L /2 againstk /3, whichishowevernotthecaseinfig.(d). Then, where will IR L vanish? At T =, mode analysis is inaccessible for unjammed systems below oint J. We can alternately measure excitations in low-temerature systems using Eqs. (5)-(7) and redict behaviors in the T = or hard shere limit. Figures 2 and (b) show examles of the dynamical structure factors measured at T = 6 and = 3, which can also be well fitted with Eq. (8). We again estimate the Ioffe-Regel frequencies by locating the intersection between Ω (k) and πγ (k). As shown in Figs. 2(c) and (d), at fixed ressure, both IR T and L IR decrease when increasing the temerature. The decrease is faster at lower ressures, leading to the tendency that IR = at a nonzero ressure c (T). We find that all the curves can be well fitted with IR ( c) /2, as illustrated by the scaling collase of all the T > data in Fig. 2(e). Does the vanishing of IR at c relate to any transitions? As mentioned above, a thermal system with reulsion

4 4 T Liquid Glass L Glass TL G/B- -2 T= -6 T=5-7 T= -7 T=5-8 T= (G/B) T /5 5.5 (b). T 3/ FIG. 3: (color online). Phase diagram including the glass transition (circles with the line T g ) and jamming-like transition (squares with the line T j 5/3 ). The diamonds and triangles locate the crossover temeratures T L ( L IR = ) and T T ( T IR = ), resectively. The star marks the location of the state shown in Fig. 2(f). undergoes the glass transition and jamming-like transition in sequence under comression. Figure 3 is the T hase diagram with both transitions. As reorted before, the glass transition temerature T g and jamminglike transition temerature T j are scaled well with the ressure: T g and T j α/(α ) [6 8]. For Hertzian reulsion, T j 5/3. In Fig. 3, we lot as well the temeratures T T and T L at which IR T and L IR tend to vanish, i.e., (T ) = c. Surrisingly, T T() collases with T j (), while T L () agrees well with T g (). This finding reveals the most significant underlying hysics of the jamming-like transition: At fixed temerature, it is the lower limit for transverse honons to survive. In the region g ( L c ) < < j( T c ), only longitudinal honons can be well-defined. In the low-temerature limit this region is just the territory of hard-shere glasses. To further evidence that transverse honons are ill defined in ( g, j ), we show in Fig. 2(f) Ω T (k) and πγ T (k) of a state in the middle of this regime. Γ T (k) is no longer scaledwithk 2 atsmallk andπγ T (k) > Ω T (k)alltheway down to k =. The jamming-like transition thus divides the whole glass regime into two arts. In Fig. 3, we use Glass TL and Glass L to denote glasses with both transverse and longitudinal honons at > j and with only longitudinal honons at g < < j. In the solid ersective, the absence of well-defined transverse honons should be the intrinsic cause of the distinct material roerties of Glass L from Glass TL [8]. In Fig. 4 we show the ratio of the shear modulus to the bulk modulus, G/B = (c 2 L /c2 T 4/3), at various temeratures and ressures. In the T = jammed solids, it is known that G/B (φ φ j ) /2 /(2α 2) [5, 3 32]. For Hertzian reulsion, G/B /3. Figure 4 showsthat G/B /3 when > j, while such ascaling relation fails at < j. Interestingly, all the curves in Fig. 4 collase well onto the same master curve with two distinct regimes when the scaling function G/B = T /5 H(/T 3/5 ) is alied, as shown in Fig. 4(b). The FIG. 4: (color online). Pressure deendence of the ratio of the shear modulus to the bulk modulus, G/B, at various temeratures. (b) Scaling collase of by alying G/B = T /5 H(/T 3/5 ). The solid lines have a sloe of /3, while the dashed line has a sloe of /2. scaling function is constructed from the known scaling j T 3/5 [7, 8, 3]. In contrast to G/B /3 at > j, G/B (/T) /2 at < j and low temeratures with a short transition region near j. The oosite ressure deendence of G/B on both sides of j directly reflects the unusual vibrational roerties of Glass L. In the hard shere limit, G/B of hard shere glasses decreases with increasing the acking fraction or ressure and decays to zero at oint J, where /T, consistent with the divergence of the mean squared dislacement aroaching oint J [3]. Our major findings that Glass L can only suort longitudinal honons and its distinct ressure deendence of G/B from Glass TL aeal for the reconsideration of the nature of hard shere glasses. They also validate the jamming-like transition, which converges to oint J in the T = limit, to be a hysically meaningful transition and a ossible substitute of the glass transition as the boundary of the jamming hase diagram in the T φlane[4]: Itisthetransitionto normal disordered solids with both transverse and longitudinal honons. The conclusions drawn here are valid for lowtemerature systems with reulsions near oint J. They may not be simly generalized to systems with long-range attractions, because the jamming transition is inaccessible to these systems. Furthermore, recent studies have demonstrated that reulsive and attractive systems differ in their glassy dynamics and vibrational roerties at low densities [33 35], which magnify when lowering the density and are thus the direct manifestations of their distinct awareness of the existence of oint J. We are grateful to Ludovic Berthier for the critical reading. This work is suorted by National Natural Science Foundation of China Grant No , National Basic Research Program of China (973 Program) Grant No. 22CB825, CAS -Talent Program Grant No. 2324, and Fundamental Research Funds for the Central Universities Grant No

5 5 [] P. G. Debenedetti and F. H. Stillinger, Nature (London) 4, 259 (2). [2] L. Berthier and G. Biroli, Rev. Mod. Phys., 83, 587 (2). [3] M. D. Ediger and P. Harrowell, J. Chem. Phys. 37, 89 (22). [4] A. J. Liu and S. R. Nagel, Nature 396, 2 (998). [5] C. S. O Hern, L. E. Silbert, A. J. Liu, and S. R. Nagel, Phys. Rev. E 68, 36 (23). [6] N. Xu, T. K. Haxton, A. J. Liu, and S. R. Nagel, Phys. Rev. Lett. 3, 2457 (29). [7] Z. Zhang, N. Xu, D. T. N. Chen, P. Yunker, A. M. Alsayed, K. B. Atowicz, P. Habdas, A. J. Liu, S. R. Nagel, and A. G. Yodh, Nature 459, 23 (29). [8] L. Wang and N. Xu, Soft Matter 9, 2475 (23). [9] L. Berthier and T. A. Witten, Phys. Rev. E 8, 252 (29); Eurohys. Lett. 86, (29). [] F. Krzakala and J. Kurchan, Phys. Rev. E 76, 222 (27). [] G. Parisi and F. Zamoni, Rev. Mod. Phys. 82, 789 (2). [2] A. Ikeda, L. Berthier, and P. Sollich, Phys. Rev. Lett. 9, 83 (22). [3] H. Jacquin, L. Berthier, and F. Zamoni, Phys. Rev. Lett. 6, 3572 (2); L. Berthier, H. Jacquin, and F. Zamoni, Phys. Rev. E 84, 53 (2); A. Ikeda, L. Berthier, and G. Biroli, J. Chem. Phys. 38, 2A57 (23). [4] X. Cheng, Phys. Rev. E 8, 33 (2). [5] P. J. Yunker, K. Chen, M. D. Gratale, M. A. Lohr, T. Still, and A. G. Yodh, Re. Prog. Phys. 77, 566 (24). [6] L. E. Silbert, A. J. Liu, and S. R. Nagel, Phys. Rev. E 73, 434 (26). [7] A. F. Ioffe and A. R. Regel, Prog. Semicond. 4, 237 (96). [8] E. Bitzek, P. Koskinen, F. Gähler, M. Moseler, and P. Gumbsch, Phys. Rev. Lett. 97, 72 (26). [9] htt:// [2] H. Shintani and H. Tanaka, Nat. Mater. 7, 87 (28). [2] J. P. Hansen and I. R. McDonald, Theory of Simle Liquids (Elsevier, Amsterdam, 986). [22] Note that S (k,) in Eq. (5) includes the contribution of thethermal energyandis roortional tok BT, while itin Eq. (2) is scaled by k BT. This difference only affects the magnitude of S (k,). It does not affect the measure of the disersion relation and sound attenuation coefficient. [23] H. Mizuno and R. Yamamoto, Phys. Rev. Lett., 959 (23). [24] L. E. Silbert, A. J. Liu, and S. R. Nagel, Phys. Rev. Lett. 95, 983 (25). [25] M. Wyart, L. E. Silbert, S. R. Nagel, and T. A. Witten, Phys. Rev. E 72, 536 (25); M. Wyart, S. R. Nagel, and T. A. Witten, Eurohys. Lett. 72, 486 (25). [26] N. Xu, M. Wyart, A. J. Liu, and S. R. Nagel, Phys. Rev. Lett. 98, 7552 (27). [27] E. DeGiuli, A. Laversanne-Finot, G. Düring, E. Lerner, and M. Wyart, Soft Matter, 5628 (24). [28] N. Xu, V. Vitelli, M. Wyart, A. J. Liu, and S. R. Nagel, Phys. Rev. Lett. 2, 38 (29); V. Vitelli, N. Xu, M. Wyart, A. J. Liu, and S. R. Nagel, Phys. Rev. E 8, 23 (2). [29] S. S. Schoenholz, C. P. Goodrich, O. Kogan, A. J. Liu, and S. R. Nagel, Soft Matter 9, (23). [3] A. Basu, Y. Xu, T. Still, P. E. Arratia, Z. Zhang, K. N. Nordstrom, J. M. Rieser, J. P. Gollub, D. J. Durian, and A. G. Yodh, Soft Matter, 327 (24). [3] B. P. Tighe, Phys. Rev. Lett. 7, 5833 (2). [32] C. P. Goodrich, A. J. Liu, and S. R. Nagel, Nat. Phys., 578 (24). [33] L. Berthier and G. Tarjus, Phys. Rev. Lett. 3, 76 (29); J. Chem. Phys. 34, 2453 (2); Phys. Rev. E 82, 352 (2). [34] Z. Zhang, P. J. Yunker, P. Habdas, and A. G. Yodh, Phys. Rev. Lett. 7, 2833 (2). [35] L. WangandN.Xu, Phys.Rev.Lett.2, 557 (24).