arxiv:cond-mat/ v1 [cond-mat.soft] 5 Nov 2002

Transcription

1 Simulation tudy of Non-ergodicity Tranition: Gelation in Colloidal Sytem with Short Range Attraction Antonio M. Puerta Department of Phyic, Univerity of Almeria, 42 Almeria, Spain arxiv:cond-mat/287v [cond-mat.oft] 5 Nov 22 Matthia Fuch Intitut Charle Sadron, 6783 Strabourg Cedex, France Michael E. Cate School of Phyic, The Univerity of Edinburgh, JCMB King Building, Edinburgh EH9 3JZ, UK (Dated: October 3, 28) Computer imulation were ued to tudy the gel tranition occurring in colloidal ytem with hort range attraction. A colloid-polymer mixture wa modelled and the reult were compared with mode coupling theory expectation and with the reult for other ytem (hard phere and Lennard Jone). The elf-intermediate cattering function and the mean quared diplacement were ued a the main dynamical quantitie. Two different colloid packing fraction have been tudied. For the lower packing fraction, α-caling hold and the wave-vector analyi of the correlation function how that gelation i a regular non-ergodicity tranition within MCT. The leading mechanim for thi novel non-ergodicity tranition i identified a bond formation caued by the hort range attraction. The time cale and diffuion coefficient alo how qualitatively the expected behaviour, although different exponent are found for the power-law divergence of thee two quantitie. The non-gauian parameter wa alo tudied and very large correction to Gauian behaviour found. The ytem with higher colloid packing fraction how indication of a nearby high-order ingularity, cauing α-caling to fail, but the general expectation for non-ergodicity tranition till hold. I. INTRODUCTION Colloidal upenion are often referred to a model ytem for tudying fundamental problem in condened matter phyic[]. Mot of the propertie of colloidal ytem are imilar to thoe of imple liquid, except for the difference in the time cale involved in the procee in liquid or colloid, making the latter more ueful in the tudy of ome baic quetion. Moreover, the interaction force between particle in a colloidal ytem are eaily tailored (e.g. by adding alt or polymer). However, there are ome feature found only in colloid, uch a aggregation or gelation, which make the tudy of thee ytem even more facinating. Gel formation, or gelation, i found in ytem with trong hort-range attraction, and i a univeral phenomenon oberved experimentally in many different ytem, ranging from colloid-polymer mixture [2, 3] to charged ytem [4], or to globular protein ytem [5]. Gelation i the formation of a percolating network (typically fractal) of dene and more dilute region of particle with void which coaren up to a certain ize and freeze when the gel i formed. Thi proce i oberved in the tructure factor a a low-q cattering peak which move to lower q, increaing it height, and then arret [6, 7, 8]. Decription of thi phenomenon ha been attempted with percolation theorie, theorie of phae eparation for tate inide the liquid-ga binodal (which i meta-table with repect to fluid-olid coexitence for hort interaction range) or in term of a gla tranition of cluter of particle [8, 9]. Recently, acknowledging it non-equilibrium character, gelation ha been interpreted uing the formalim of mode coupling theory, MCT, for non-ergodicity tranition [,, 2]. Thi approach view the gel a particle trapped by a network of bond which hinder the particle motion, reulting in a non-ergodic tate. Thu, gelation i caued by formation of long lived bond, whoe collective arret i decribed a a a normal non-ergodicity tranition. (Thi i ditinct from many earlier approache whereby the bond were aumed to form irreveribly from the outet.) In the preent imulation tudy, we want to tet thi uggetion critically, thereby etablihing the exitence or otherwie of a non-ergodicity tranition correponding to bonding network formation. Alo preent in colloidal ytem i the equivalent of the uual gla tranition in imple liquid, which occur at high denitie, and i driven by teric imprionment. Thi tranition ha been tudied experimentally and compared to MCT thoroughly [3, 4, 5, 6]. When two different non-ergodicity tranition are oberved in a ytem, MCT predict a high order ingularity in the region where the driving mechanim for both tranition are preent [7, 8, 9]. Therefore, a higher order tranition i expected at high attraction trength and high denity in colloidal ytem with attractive interaction [, 2, 2]. Computer imulation have been ued to tet the expectation from MCT in many different ytem, uch a a Lennard-Jone liquid [22, 23, 24, 25], water [26, 27], trong gla former [32, 33, 34] and polymer [28, 29, 3, 3]. The tet have hown that the prediction from MCT are correct, not only qualitatively but alo, in part, quantitatively [33, 35]. However, they have alo pointed out ome difference, epecially in the patial correlation

2 2 of particle mobility [36, 37, 38, 39]. In none of thee imulated ytem, however, did gelation occur, preumably becaue the attraction were not hort-ranged enough. In thi work, we have ued molecular dynamic imulation to tudy the propertie of the gel tranition, and compared them with the prediction from MCT. (Thi wa initiated in [4] where ome further reult may be found.) We take the numerou univeral prediction of the theory to tet the cenario qualitatively. Comparing with quantitative prediction available for ytem of hard phere[4, 42], phere with hort range attraction [,, 2, 2], and the mentioned imulation tudie, we identify the novel mechanim driving the non-ergodicity tranition which i the caue of gelation for moderately dene upenion. Molecular dynamic were ued intead of Brownian dynamic becaue the choice of microcopic dynamic doe not affect the relaxational dynamic of a ytem cloe to a non-ergodicity tranition [24]. By mean of the Aakura-Ooawa interaction potential [43], we imulate the behaviour of a colloid-polymer mixture, which i a well-undertood ytem [44, 45, 46, 47]. For hort interaction range, thi ytem exhibit a fluidcrytal tranition, at intermediate denitie and increaing attraction trength, with a liquid-ga tranition metatable to the fluid-crytal one. In our imulation, the ytem wa modified to prevent both of thee phae tranition from occurring, in order to be able to tudy the tranition from the fluid to the non-equilibrium tate. The paper i organized a follow: Section II decribe ome reult from MCT which will be ued in the ubequent analyi of the imulation reult. In ection III the imulation method i preented and the detail are given. Section IV, deal with the reult and i divided into four ubection tudying i) the correlation function, ii) the time cale and the diffuion coefficient, iii) the mean quared diplacement and iv) a higher colloid concentration. Finally, in ection V, we preent the concluion of thi work. II. MODE COUPLING THEORY In thi ection we will preent the mot important MCT reult on non-ergodicity tranition. MCT attempt a decription of the denity correlator and it elf part, in term of a fluctuating-force correlator [48, 49]. In thi paper, only the elf part of the denity correlator will be tudied, defined a: Φ q(t) = exp{iq(r j (t) r j ())} () where the bracket denote average over particle j and time origin, and q i the wave-vector. The equation of motion of Φ q in Brownian (coare grained) dynamic, i given by: τ q t Φ q(t)+φ q(t)+ t m q (t t ) t Φ q(t )dt = (2) where τ q i a ingle particle diffuive time cale and m q (t) i a mode coupling kernel which decribe the cage effect [5]. Within MCT, gla tate are given by nonzero olution of thi equation for the long time limit of Φ q(t ) = fq, the o-callednon-ergodicityparameter. It decribe the gla tructure and may alo be called Lamb-Möbauer factor. The gla tranition i marked by a (generally) dicontinuou tranition from the unique trivialolution in the liquid, fq =, to multiple olution in the gla, fq >, where only the highet olution i phyical. Gla tranition can be claified according to the number, l, of non-trivial olution merging with the highet one, and the type of tranition i noted a A l. For liquid tate cloe to the gla, a two tep decay i oberved for the correlator; the plateau i at fq and ignal the proximity of the gla tranition. Around thi plateau, Φ q how ome univeral propertie, depending on the type of tranition. For the mot common type of tranition, A 2, the decay to the plateau, and that from the plateau, can both be expreed a power law expanion. In particular, the decay from the plateau i given by: Φ q(t) = f q h () q (t/τ) b +h (2) q (t/τ) 2b +O((t/τ) 3b ) (3) with h () q and h (2) q amplitude and τ the final or α- relaxation time cale. b i known a the von Schweidler exponent, and depend on the detail of the interaction potential. Expreion (3) implie time caling for the decay from the plateau, called α-decay, for different tate cloe to the gla tranition. The time cale, τ, diverge a the gla tranition i approached according to a power law, with an exponent γ, which can be related to the von Schweidler exponent: τ σ γ, with σ the ditance to the tranition [48, 49]. On the other hand, the wave-vector dependence of the non-ergodicity parameter and amplitude give ome non-univeral propertie of the tranition, providing information about the mechanim cauing the non-ergodicity tranition. For high order ingularitie, the fluid tate cloe by how again a two tep decay in the correlation function, but the decay to and from the plateau are no longer power law expanion. Intead, logarithmic law are obtained [8, 9]. A alient feature i that a logarithmic decay around the plateau i predicted: Φ q(t) = f A q C q log(t/t ) (4) where fq A i the non-ergodicity parameter of the high order ingularity, C q i an amplitude and t i a time cale (the time when the correlator lie on the plateau). The mean quared diplacement (MSD) can be tudied intead of the correlation function, obtaining a imilar two tep behaviour. Similar aymptotic law to decribe the decay to and from the plateau can be derived, and the parameter and exponent can be related to thoe of the correlation function [42]. The value of the plateau in

3 3 the MSD define the localization length and i a meaure ofthe ize ofthe cage. However,it hould be noticed that the cage, a formed by other particle, i contantly retructuring cooperatively. Only when the particle have broken free of their cage, diffuive motion i oberved, with a elf-diffuion coefficient, D, that tend to zero a the gla tranition i approached a D σ γ for the uual A 2 tranition. Two different non-ergodicity tranition have been found in colloidal ytem with a hort range attraction [, 2]: a teric hindrance driven gla tranition and an attraction driven gel tranition [,, 2]. While the firt i found at high denitie and i qualitatively imilar to the glatranitionin the hardphereytem(hss) or Lennard-Jone ytem (LJS), the gel tranition occur at high attraction trength for all volume fraction. Different propertie for thee two tranition are predicted, the main difference ariing from the driving mechanim: the localization length i horter in the gel than in the gla, reulting in higher non-ergodicity parameter. Alo, a maller von Schweidler exponent for the gel than for the gla i expected, implying a higher value of γ, i.e. the tranition a oberved by τ q or D i more abrupt. The actual hape of the non-ergodicity tranition line depend on the detail of the interaction potential, although ome general feature can be found. From lower to higher interaction trength, the gla line, i lanted to higher concentration, howing that a weak attraction fluidize the gla. However, at even higher interaction trength, the gel tranition occur at lower colloid denity the higher the attraction trength. A a reult, a re-entrance tranition i obtained at high colloid volume fraction. The line may be wedge-haped or curved in thi region, depending on the range of the interaction. If the line i wedge-haped a high-order tranition (generically A 3 ) i preent near the corner, wherea none exit if the line i moothly continuou. An A 4 ingularity appear right at the vanihing of the A 3 point when the line firt become mooth [2]. III. SIMULATION DETAILS Equilibrium molecular dynamic imulation mimicking a colloid-polymer mixture were performed for a ytem compoed of oft-core polydipere colloidal particle. The core-core interaction between particle wa modeled by: ( ) 36 r V c (r) = k B T (5) a 2 where a 2 = a +a 2, with a and a 2 are the radii of the interacting particle. A flat ditribution of radii with a width of δ =.a, where a i the mean radiu, wa ued. The exponent in V c wa elected high enough to avoid problem related to the oftne of the potential[5]. The polymer induce an attractive depletion interaction between the colloidal particle, which wa modeled by the Aakura-Ooawa interaction potential [43, 47]. The extenion of thi potential to take polydiperity into account read [52]: {[ V AO (r) = k B Tφ p ( η +) 3 3r ] 4ξ ( η +)2 + r3 6ξ ξ [ 4r (η η 2 ) 2 ( η +) r ] } 2 2ξ (6) for r 2(a 2 +ξ) and for larger ditance. Here, η i = a i /ξ; η = (η + η 2 )/2, and φ p i the volume fraction of the polymer. Note that the range of the potential i given by the polymer ize, ξ, and it trength by φ p. Thi potentialwamodifiedaroundr = a 2, toenurethatthe minimum of the total potential (V c +V A ) occur at thi point: forr 2a 2 +ξ/5aparabolicform,whichconnect analytically to V AO at 2a 2 + ξ/5 and ha a minimum in 2a 2, wa ued. In our imulation, the range of the interaction, 2ξ, wa et to.2a, which would correpond to polymer with R g /a =. where R g i the radiu of gyration. A long-range repulive barrier wa added to the interaction potential in order to prevent liquid-ga eparation (a hown below). The barrier had a maximal height of k B T, according to a fourth-order polynomial: V bar (r) = k B T { ( r r r r ) 4 ( ) 2 r r 2 +} r r (7) for r r r and zero otherwie. The limit of the barrier were et to r = 2(a 2 +ξ), and r = 4a, which wa enough to prevent phae eparation. The maximum height of the barrier equal the depth of the depletion interactionatcontactforφ p =.625,muchlowerthan the value where the gel tranition take place. The reulting total interaction potential, V tot = V c +V AO +V bar, i analytical everywhere. It i hown in figure, where in order to indicate the pread induced by polydiperity, the potential among three different pair with differing radii are plotted. In our imulation, length were meaured in unit of themeanradiu, a, andtimein unitof 4a 2 /3v 2, where thethermalvelocityv waetto 4/3. Equationofmotion were integrated uing the velocity-verlet algorithm, in the canonical enemble(contant NTV), to mimick the colloidal dynamic. Every n t time tep, the velocity of the particle wa re-caled to aure contant temperature. No effect of n t wa oberved for well equilibrated ample. The time tep wa et to.25. Equilibration of the ytem wa teted by monitoring the total energy, and other order parameter (ee below), and by meauring Φ q(t) and the MSD at different initial time. When the order parameter were contant and the Φ q (t) and

4 4 V tot V tot r / a Ψ r / a FIG. : Total pair interaction potential V tot a function of the radial ditance r = r r 2 for three different particle pair; a pair of particle with minimal radii a = a 2 = a δ, one with average radii a = a 2 = a, and one with maximal a = a 2 = a + δ (from left to right). The inet how the enlarged region of the attractive minimum. Croe mark where the parabolic minimum moothly matche to eq. (6). MSD curve howed no dependency on the initial time (ageing), the ytem wa conidered to be equilibrated. The( volume fraction of the colloidal particle, φ c = 4 3 πa3 + ( ) δ 2 ) a n c, with n c the colloid number denity, and the polymer volume fraction, φ p, were the control parameter ued to identify the tate in the phae diagram. In order to explore the whole φ p φ c plane in earch of the gel tranition, phae tranition which forbid acce (in equilibrium) to important part of the plane mut be prevented. Several order parameter were ued to identify different kind of ordering in our ytem and to monitor whether unwanted liquid-ga or fluid-crytal tranition were taking place. Firt, the onet of phae eparation involving tate of different denity can be detected by dividing the ytem into n 3 boxe and meauring the denity in every box. The demixing order parameter i defined a the tandard deviation of the ditribution of denitie: Ψ n = (ρ k ρ) 2 (8) n 3 k= where ρ k i the denity of particle in box k, and ρ i the mean denity. Thi parameter i cloe to zero for an homogeneou ytem, and increae if it demixe into phae of different denity. In our cae, n ha been et to 4, implying 64 boxe, and a box edge of about 5a (depending on φ c ). On the other hand, the orientational order parameter, Q 6, a defined by Steinhardt et al. [53, 54], ignal the preence of an ordered phae, and i ued to detect crytallization. The phae diagram wa probed uing thee parameter. In figure 2 the reult are preented for a bare Q φ p FIG. 2: Demixing (Ψ 4) and orientational (Q 6) parameter for φ c =.4 and increaing polymer fraction, φ p, for different ytem: monodipere without long-range barrier (quare), polydipere without barrier (circle), and polydipere with long-range barrier (croe). ytem (monodipere and without the long-range barrier), a polydipere ytem without the long-range barrier, and the final ytem with both polydiperity and barrier. In thi figure, the colloid volume fraction i contant, φ c =.4, and the polymer concentration varie; an iochore i tudied. The udden increae in both Ψ 4 and Q 6 occurring at φ p =.2 for the bare ytem, ignal the crytallization boundary, in accordance with Dijktraetal. [47]. Becaueofthe hortrangeofthe potential, thi ytem ha no liquid phae; i.e. the liquid-ga coexitence i meta-table with repect to the crytal-ga tranition. When polydiperity i introduced in the ytem, crytallization i prevented, a indicated by the contant trend of both parameter cloe to φ p =.2. However, a the ytem now doe not crytallize, the liquid-ga tranition can be reached upon increaing the trength of the interaction. Thi demixing i ignalled by an increae in Ψ 4, not involving local ordering. In order to avoid thi eparation, the long-range barrier ha been introduced in the interaction potential. The energy of a dene phae i raied, and demixing i thu energetically unfavourable. Figure 2 how that liquid-ga eparation i indeed inhibited by the repulive barrier. Intead, individual void of finite ize are created in the ytem, cauing a low-q peak in the tructure factor, S(q), preented in Figure

5 5 S(q) g(r) r / a q a FIG. 3: Structure factor for different polymer fraction at φ c =.4: φ p = gray line, φ p =.2 dahed line, φ p =.35 dotted line and φ p =.425 olid black line. Note the low-q peakraie aφ p iincreaed. Inet: Pairditributionfunction, g(r) for the ame tate. 3. There, S(q) i hown for different polymer fraction, ranging from no attraction (φ p = ) to the cloet tate to the gel we have acceed (φ p =.425). In the inet to Figure 3, the pair ditribution function, g(r), i preented for the ame tate a the tructure factor. The value at contact, r = 2, increae continuouly a the atraction trength grow, ignalling increaed local contact probabilitie. Thi proce will be hown to be reponible for the gel tranition. In S q it become evident a an increae in the ocillation for large q. The low-q peak in the tructure factor reemble the low-angle peak oberved in light cattering experiment with colloidal gel [6, 8]. However, wherea the peak in our ytem i an equilibrium property, induced by the pecific hape of the interaction potential, the experimental peak ha a non-equilibrium origin. We alo checked for the poibility of microphae eparation, which in ome cae can be induced by a repulive barrier [55]. In our cae, the mall angle peak continuouly increae with φ p, but tay finite and maller than the neighboring peak a we approach the gel tranition. We interpret thi to indicate that we do not have microphae eparation, andwealoobervednootherignofuchordering. Furthermore, ince the relevant wave vector in the MCT calculation of the gel tranition are the high one(around 2π/ξ), the change in low-q region in the tructure factor i expected to have little effect on the gel tranition. IV. RESULTS AND DISCUSSION Thi ytem ha been previouly hown to undergo both the gla and gel tranition a tated by MCT. It alo exhibit a logarithmic decay in the correlation function at high colloid concentration, indicating a highorder ingularity in that region [4]. In thi ection we will dicu the propertie of the gel tranition, and compare them with MCT and with thoe of the HSS and other ytem, which are imilar to the gla tranition at high colloid concentration. We tet for difference by comparing quantitatively the non-univeral feature of the tranition, which will aid in the identification of the driving mechanim. The gel line i predicted to extend to low packing fraction with the ame qualitative propertie. In order to tet thee propertie we have performed imulation at two different colloid concentration, φ c =.4 and φ c =.5, where the gel line i far away from the percolation one. At high concentration, the higher order ingularity i expected to affect the equilibrium tate, diturbing ome feature of the gel tranition. A. Self-Intermediate Scattering Function The caling prediction for the α-decay of tate cloe to a non-ergodicity tranition i teted in Figure 4, for contant colloid packing fraction, φ c =.4. Two different repreentative wave-vector are preented in thi Figure, q = 6.9 and q = 5. A oberved at the gla tranition in the HSS and many other different ytem [22, 23, 27, 3, 34, 4], the α-caling property hold. In comparing thee correlation function with thoe typical for the HSS or LJS, it i noticed that in Figure 4 the α-decay of the correlator i more tretched, implying a maller von Schweidler exponent at the gel tranition than at the gla tranition. Becaue of thi tretching in the α-decay, a clear plateau i not oberved, although a lowing down of more than four decade i tudied. Neverthele, extrapolating the relaxation curve to extract plateau value, much higher value are found than at the gla tranition in the HSS or the LJS. We have analyzed the tate φ c =.4 and φ p =.42 in more detail, which how four decade of lowing down compared to the purely repulive ituation upon turning on the attraction. Becaue caling i oberved in Figure 4, tudying only one tate i enough to analye the α-decay of the correlation function. The lowet tate, φ p =.425 wa not choen becaue it trongly deviate from the expected behaviour of τ q v. φ p (ee Figure 8 and dicuion thereafter). The correlation function at different wave-vector for tate φ c =.4 and φ p =.42 are preented in Figure 5. The range of wave-vector tudied, where the plateau height change i much wider thanthe rangeforaimilarchangeinf q attheglatranition of hard phere or Lennard-Jone particle. Thi feature indicate that the relevant ditance for the gel tranition are much horter than for the uual gla tranition. The correlation function were meaured until the average particle diplacement wa 5a, which i one fourth of the box ize(2.95 a). Thu, extending thi meaurement to longer time in order to oberve the whole α-decay at low q i troubleome. If the diffuion coefficient diverge at the ame rate a the α-time cale (a predicted by

6 6.8.8 Φ q.6.4 q = 6.9 Φ q (t) t Φ q.4.2 q = 5 FIG. 5: Intermediate cattering function (elf part), Φ q, for the tate φ c =.4 and φ p =.42, for different wave-vector; from top to botton: q = 3.9, 6.9, 9.9, 5, 2, 25, 3, 35, 4, 5. The dahed line are fitting from eq. (3) up to econd order, with the ame von Schweidler exponent for all q t / τ FIG. 4: Intermediate cattering function (elf part), Φ q, v. re-caled time with the α-time cale, τ q, for different tate: φ c =.4 and φ p =.375,.39,.4,.4,.45,.42 and.425 from right to left. Two different wave-vector are tudied: q = 6.9 (upper panel) and q = 5 (lower panel), with the KWW fit (dahed line) included. f q, hq (), hq (2) MCT), thi problem would not appear. Thu, we are alo oberving a dicrepancy between both time cale divergence, that will be further dicued below. The impoibility of oberving a clear plateau, a mentioned above, make it more difficult to analye the correlator, ince fq cannot be fixed a priori. Therefore, the α-decay of the correlation function ha been fitted uing expreion (3), with fq, h () q and h (2) q a fitting parameter. The von Schweidler exponent wa alo fitted but wa kept identical for different wave-vector. It wa found a b =.37, and the other reult for the fitting parameter, are hown in Figure 6. The trend of thee parameter are imilar to that of the gla tranition in both HSS and LJS, but over a wider q-range in the gel cae. Thi indicate that the localization length i quite different in the preent ytem. The non-ergodicity parameter exhibit a bell haped curve, wherea the firt order amplitude decribe a maximum. The latter i determined from the fit up to a prefactor which depend on the choice of τ in eq. (3). A an etimate we have ued τ q for q = 9.9 (τ q i defined by Φ q (τ q) = f q /e), which yield value that are imilar (in magnitude) to the HSS. The econd amplitude how a monotonouly increaing behaviour with q, in accordance with the HSS, but it i alway poitive, unlike the HSS where it goe through q FIG. 6: Non-ergodicity parameter, f q (open circle), and firt (filled circle) and econd order (diamond) amplitude from the fitting in Figure 6 with b =.37 for all wave-vector. The line give the Gauian approximation from eq. (9) fitted to all wave-vector (olid line) and the three lowet wave-vector (dahed line). zero at the peak of h () q. The non-ergodicity parameter, f q, can be approximated uing the Gauian expreion: f q exp { q 2 r 2 l/6 } (9) where r l i the localization length. Thi approximation i known to be valid for low wave-vector, and important deviation from the Gauian behaviour are expected cloe to the gla tranition. However, the value for the localization length obtained from fitting thi curve (olid line in Fig. 6), can be ued a an etimate of the one in the MSD.

7 7 A q and f q, βq and b τ q K and τq q q FIG. 7: Parameter ued in the KWW fitting and comparion with the von Schweidler fitting. Main figure: A q (cloed circle) and f q (open circle), β q (croe) and b (horizontal dahed line). Inet: τ K q (cloed circle) and τ q (open circle). The localization length o obtained i r 2 l =.26a 2, much maller than for the HSS or the LJS, where r l i of the order of the Lindemann ditance. Thi feature how that the proce cauing the non-ergodic tranition in our cae ha a typical ditance much maller than in the cae of gla tranition in the HSS or the LJS. Thi agree with the obervation and dicuion about the height of the plateau, and of the different q-range covered by f q in Figure 6. Wherea the gla tranition in the HSS i driven by core-core repulion, the gel tranition i caued by the hort range attraction, therefore by bond between particle, (ee inet to Figure 3) whoe ize i of the order of the interaction range. An intereting analogy ha been etablihed between the mechanim driving the formation of gel and glae, and the freezing tranition [56]. The α-decay of near-non-ergodic tate can be alo tudied uing the Kohlrauch-William-Watt (KWW) tretched exponential. The KWW expreion i given by: { ( ) } βq t Φ K q (t) = A q exp τ K q () where β q i known a Kohlrauch exponent, which ha been hown to coincide with the von Schweidler exponent at high wave-vector [57]. Thi expreion ha been fitted to very different ytem, and decribe the α-decay down to zero. We have fitted thi expreion to the α- decay in our ytem. However, ince the correlator in Figure 5 do not how the complete α-decay, we have fitted expreion() to the mater curve, obtained from the α-recaling. Two of thee fitting are preented in Figure 4 by the dahed line, howing that the KWW tretched exponential decribe well the α-decay in thi ytem. The fitting parameter, A q, β q and τ q, are preented in Figure 7, and compared with the correponding parameter in the von Schweidler formalim. In uch a way, A q i compared with the non-ergodicity parameter, β q with the von Schweidler exponent, and the τq K with τ q. A expected, the height of the plateau can be determined equally well both by the KWW or von Schweidler analyi. The ame hold for the time cale, τq KWW and τ q. The Kohlrauh exponent i expected to tend to at low wave-vector, and to approach the value of the von Schweidler exponent at high q. The low-q limit i explained becaue diffuion i the dominant proce over long ditance, wherea at hort ditance(comparable to the cage ize) the dynamic i dominated by the cooperative local rearrangement. Thi behaviour i predicted from MCT [57], and ha been oberved in different ytem, uch a molecular gla former [58], and in imulation of polymer melt [3], and of water [59]. In our cae, the low-q limit i not oberved, but β rie a the wave-vector decreae, indicating that the expected behavior may appear at lower q below the mall angle peak in S(q). At high wave-vector, the Kohlrauh exponent croe the von Schweidler value, but tay cloe to it. Although an exact agreement i not oberved, we may conclude that the correct general trend i obtained. B. Time Scale and Diffuion Coefficient An important univeral prediction of MCT i the exitence of power law divergence for both the time cale, τ, and the invere of the elf diffuion coefficient D, with the ame exponent in both cae, γ: τ q ( φ G p φ p) γ and D ( φ G p φ p) γ () where φ G p i the polymer volume fraction where the gel tranition occur. The relation between exponent γ and the von Schweidler exponent, b, i alo univerally etablihed by MCT [49]. Teting of the power law divergence (and meauring of γ) i uually carried out plotting τ q a a function of φ G p φ p for different value of φ G p, looking for a traight line. Thi method i cumberome, even more a deviation from it are expected for tate cloe to the tranition, and precie value for γ and φ G p cannot be given. To avoid thi difficulty, we have calculated γ from b, a given by MCT, and with thi particular value of the exponent looked for the power law divergence. In uch a way, we are teting the compatibility of MCT prediction with our data. Figure 8 how the wave-vector dependent time cale τ q v. φ G p φ p for different wave-vector. For everywavevector, τq /γ wa extrapolated to zero, yielding a value for the polymer fraction at the gel tranition, φ G,q p. The final value of φ G p, ued in Figure 8, wa calculated a the average value for all wave-vector tudied. The linear trend in Figure 8 for φ G p φ p > 5 3, how the powerlaw behaviour predicted by MCT, with exponent γ = 3.

8 8 4 φ p G = γ = D τ q τ q 2 γ = 3.25 D φ p γ =.23-2 γ = φ p G - φp φ p G - φp FIG. 8: Wave-vector dependent time cale, τ q, v. φ G p φ p for the iochore φ c =.4 for different wavevector; ymbol from top to botton correpond to q = 3.9, 6.9, 9.9, 5, 2, 25, 3. The line are power law fitting to q = 3.9 and q = 3. In all cae γ keep cloe to thi two value. FIG. 9: Self diffuion coefficient, D, v. φ G p φ p for two value of φ G p : φ G p =.459 (cloed circle) and φ G p =.4265 (open circle). The dahed line are the power law fitting to the data, with the exponent hown in the figure. Inet: D τ q v. polymer fraction for different wave-vector. Symbol a in Figure 8. and φ G p = The cloet tate to the gel tranition, φ p =.425 deviate from the power-law behaviour oberved for lower polymer fraction. Similar deviation have been oberved in the HSS and LJS and can tentatively be attributed to thermally activated procee (or hopping event) [6]. A hown in eq. (), MCT predict a power-law for the elf diffuion coefficient, D, with the ame exponent a the divergence of the time cale. Simulation on HSS and LJS have hown that a power law divergence i indeed obtained, but with a different exponent than in the cae of τ q. Uing the ame procedure a decribed above (calculating γ from b and extrapolating D /γ to obtain φ G p ) yield a value for φg p =.459, with the ame γ a for the time cale. Thi value of φ G p i too far from that obtained uing τ q. Therefore, we cannot have imilar φ G p and γ to explain the behaviour of both τ q and D, implying that the MCT prediction, eq. () i violated. In Figure 9, we preent D v. φ G p φ p uing for φ G p both the value etimated from τ q and that from D. We conider more deirable to have imilar φ G p to explain the behaviour of τ q and D, even though thi implie two different γ: γ = 3. for τ q and γ =.23 for D. A obtained in other non-ergodicity tranition [22, 3, 34], the γ exponent i lower in the diffuion coefficient than in the time cale, although the difference between both value of γ i bigger in our cae. In order to tre the different γ exponent in the divergence of the time cale and /D, we have plotted D τ q a a function of φ p for different wave-vector in the inet to figure 9. Thi product, that hould be contant according to MCT, diverge a the polymer fraction approache φ G p. The divergence follow a power law with the exponent equal the difference between both value of γ. The maximum in the elf diffuion coefficient (upperrightcorneroffigure9)iaconequenceofthere-entrant gla tranition at high packing fraction[6, 62]. A weak hort range attraction at firt detabilize the cage and thu the gla tranition move to higher particle concentration initially a the polymer fraction i increaed. At contant colloid concentration the diffuion thu firt peed up with increaing φ p, until for intermediate attraction trength the gel line i approached, where the oppoite trend then dominate. At φ c =.4, the gla tranition i rather far removed and thu ha little effect, but the increae of D i till meaurable and the diffuion coefficient can be ued a a meaure of the ditance to the cloet tranition. The maximum thu indicate the re-entrant hape of the non-ergodicity line. The wave-vector dependence of the time cale τ q can alo be compared with theoretical prediction. At low q, the time cale i expected to behave a q 2, correponding to a diffuive proce over large ditance. Yet, becaue the imulated cattering function exhibit non exponential relaxation even for the mallet wavevector, thi imple theoretical cenario i not expected to appear in our cae. At intermediate wave-vector, where the Kohlrauch exponent become comparable to the von Schweidler one the theory predict a decreae a q /b, wherea at even higher q, the ditance involved are dominated by the microcopic dynamic, and correction to thi behaviour are expected [63]. The inet to Figure how τ q for different tate cloe to the gel tranition. In order to make clear common propertie the curve have been caled vertically to collape (main figure). It can be een in thi figure that the behaviour of τ q at low wave-vector (below q = ), indeed how a q 2 behavior, which however i not the one explained by MCT.

9 9 2 q -2 q -3.3 Gla tranition Recaled τ q q τ q δr Gauian approximation q t FIG. : Inet: Time cale, τ q a a function of the wavevector, q, for different tate cloe to the gel: φ p =.42 (circle), φ p =.45 (quare), φ p =.4 (diamond), φ p =.4 (upward triangle), φ p =.39 (left-ward triangle) and φ p =.375 (downward triangle). Main body: Same data, re-caled to collape in the low-q power-law behavior. The dahed line how power-law behaviour with exponent 2 (gray line) and 3.3 (black line). FIG. : Mean quared diplacement of the particle a the gel tranition i approached. Black curve from left to right: φ p =.3,.325,.35,.375,.39,.4,.4,.45,.42,.425. Grey curve: φ p =. The dahed horizontal line indicate the localization length correponding to the gla tranition (hort dah) in the HSS and the etimate from Figure 6 (long dah). Dotted curve: mean quared diplacement for a ingle particle in a frozen environment at φ p =.425. At higher q, another power-law trend i oberved, with a higher exponent: q 3.3. The crooverfrom the low-q behaviour to the high-q one, compare nicely with the wavevector where the Kohlrauh exponent become equal to the von Schweidler one (Figure 7). The exponent of the high-q region yield b =.3, lower than the value obtained from the analyi in Figure 5. However, thi value i quite cloe to the meaured von Schweidler value and much maller than the HSS one. Deviation from thi power-law behaviour are oberved at high q for the lowet φ p preented in the figure. Thee deviation are caued by the microcopic dynamic, a they occur when τ q i lower than a certain value, regardle the polymer fraction. Thi value, preented in the inet a an horizontal line, i t.6, which agree with the time one would etimate from the correlator in Figure 5. C. Mean Squared Diplacement We turn now our attention to the MSD curve, that were partially analyed to obtain the diffuion coefficient preented in Figure9. We areonly intereted in the lowing down cloe to the gel tranition and thu we do not how the MSD for low polymer fraction, where the attraction peed up the dynamic and increae the diffuivity (ee figure 9 and [6]). The MSD, after a hort initial regime of free flight, δr 2 t 2, low down becaue of the particle interaction and take longer and longer to reach the long-time regime diffuive, where δr 2 = 6D t. AnimportantfeaturethatcanbeobtainedfromtheMSD of the particle in the ytem, i the localization length, where the particle interaction hinder particle motion mot trongly, and in the idealized gla tate, arret it. It can be compared with the etimate uing the Gauian approximation (ee Fig. 6). In Figure we preent the MSD for increaing polymer volume fraction. A the gel tranition i approached, the localization length how up a an indication of a plateau, ignaling the bond formation. A already dicued above, r l i much horter than in the HSS gla tranition (upper dahed line in Figure ), becaue of the driving mechanim. The lower dahed line in thi figure i the localization length, a etimated from the non-ergodicity parameter uing the Gauian approximation (r 2 l =.26). Although a clear plateau ha not fully developed in our curve, it height eem to be above that etimate, by a factor.5 2. Since the Gauian approximation work very well in the cae of the HSS, thi ugget big non- Gauian correction at the gel tranition. Before teting the Gauian approximation, we tre that the localization length give a typical ize of the meh of bond formed between neighbouring particle, and that the low tructral unit are continuouly and cooperatively rearranging. In order to tet thi idea about a correlated region which cooperatively rearrange with and around each particle, a ingle mobile particle i conidered in a fixed environment. A well equilibrated ytem with φ p =.425 i frozen, and only one particle i allowed to move. Thi mobile particle now explore a frozen environment, providing the tructural ize of the region it i confined to. The mean quared diplacement o obtained igivenin Figure(dotted line). Someparticle(.6%) were able to break their bond and diffue freely in the frozen environment. For the particle that tay localized, it can be oberved that the length of the frozen bond i much maller than the localization length. Thi fact demontrate that the tructure of bond, like the repul-

10 ive cage at the gla tranition in the HSS or LJS, i dynamic, and contantly rearrange cooperatively. Thi collective retructuring of the ytem fluidize it, and retore ergodicity, which cannot occur in the frozen ytem, where the particle are not able to diffue even at very long time. We turn now back to the Gauian approximation, and it accuracy. Uually, thi i teted by meauring the non- Gauian parameter, defined a: α 2 = 3 r4 (t) 5 r 2 (t) 2 (2) where the average imply enemble averaging. Thi parameter meaure the deviation of the probability denity function for the ingle particle motion from Gauian behaviour, and vanihe for diffuive motion. Special care mut be taken when performing the enemble average in polydipere ytem, a pointed out in [64]. The non- Gauian parameter mut be calculated for every particle (the average in the definition above thu implying timeorigin averaging only), and particle averaging i taken on the value of α 2 (o long a long enough time interval are tudied, each particle will ample the ditribution relevant to it own ize in an ergodic fahion). The non- Gauian parameter for tate with increaing φ p are preented a a function of time in Figure 2. At hort time α 2 tend to zero, ince the ytem how Gauian behaviour during it unhindered ballitic regime. At long time, when the particle break free from their bond andhydrodynamicdiffuionhold, α 2 againgoebacktoward zero. At intermediate time, correponding to the plateau in both the correlation function and the mean quared diplacement, α 2 grow, ince the ingle particle motion hindered by bonding i not Brownian. A a reult, α 2 how a maximum, whoe height and poition grow in time, becaue the particle take longer and longer to break free and tart diffuing. The behaviour of the imulated α 2 obey the general expectation [36, 37, 39], but important difference are oberved in the comparion with the reult for the HSS or LJS. Wherea in thoe cae the height of the maximum for imilar (or even higher) α-relaxation time i around 2, at the gel tranition much higher value are meaured. Another intereting difference i the failure of the hort time caling, oberved both in the HSS and LJS. Both effect can be rationalized conidering that the cage i indeed a network of bond in the cae of a gel, rather than a cavity. The trength of thee bond i given by the intenity of the interaction, and thu, it i modified for different tate, diabling the hort time collape. Becaue the bond are hort ranged, they affect the particle motion from very hort time onward, o that the particle feel the hindrance much longer in the gel cae. It can be concluded that the non-gauian correction are very important in the gel tranition. Therefore, the localization length etimated from the non-ergodicity pa- α t FIG. 2: Non-Gauian parameter α 2 a a function of time for tate approaching the gel tranition at the ame tate a Figure. The maximum increae with increaing φ p. rameter may be inaccurate, a dicued above. However, it till provide an indication of how mall the localization lengthi. Abetterindicationofr l canbeobtainedwithin the Gauian approximation if only low wave-vector are ued in fitting expreion (9). The fitted curve i preented in Figure 6 by the dahed line, where only the three lowet q are fitted. The etimated f q deviate from the data at higher wave-vector, howing high non- Gauian correction. The localization length i higher than the previou value: r 2 l =.62. Thu, thi fitting provide data more conitent with the MSD curve and the non-gauian parameter. D. Higher Colloid Volume Fraction We move now to a higher colloid volume fraction: φ c =.5. Thee reult are preented to upplement the finding at the lower packing fraction and tet for the prediction of tronger tretching cloer to the higher order ingularity. A indicated in the theoretical ection, MCT predict a higher order ingularity in the vicinity of the junction of the gel and gla line; i.e. at high polymer and colloid denitie. In thi particular ytem we found clear indication of thi ingularity in imulation at φ c =.55 and φ p =.375 [4]. The iochore under tudy now, φ c =.5, could be cloe enough to the higher order ingularity to how ome effect. In Figure 3 we preent the correlation function for increaing polymer fraction at the ame wave-vector a Figure 4, re-caled to collape in the long-time decay. It i intereting to note that the polymer concentration tudied in thi cae are lower than thoe tudied at the lower colloid volume fraction. In accordance with experiment and theory, thi indicate that the gel tranition take place at lower polymer fraction the higher the colloid concentration. In Figure 3, it can be oberved that the correlator

11 .8.8 Φ q.6.4 q = 6.9 Φ q (t) t Φ q.4.2 q = 5 FIG. 4: Intermediate cattering function (elf part), Φ q, for the ame wave-vector a Figure 5 wave-vector for φ c =.5 and φ p =.39. The dahed line are fitting from eq. (3) up to econd order, with the ame von Schweidler exponent a φ c = t / τ FIG. 3: Correlation function for φ c =.5 and different polymer concentration. From left to right: φ p =.35,.36,.375,.38,.385,.39. Two wave-vector are tudied, a labeled in the figure. The dahed line repreent the KWW fitting to the φ c =.4 correlation function. do not collape over the whole α-decay, but only in the end. Thee deviation are expected becaue of the higher order ingularity, which i at higher denitie. However, we tre that although thi ingularity ha clear effect on the correlation function, they do not how o clear ignature a that of the φ c =.55 iochore [4]. At thi high concentration, a logarithmic decay wa oberved, with a wave-vector dependent extenion. Becaue thee correction affect the early α-decay, analyzing the correlation function i difficult. Furthermore, the plateau i not oberved, and the von Schweidler analyiithuextremelydifficult. Inordertoanalyetheelf intermediate cattering function, we compare the tretching of the curve at φ c =.4 and φ c =.5; in Figure 3 the α-decay mater function of the φ c =.4 tate, a parametrized by the KWW fitting i included. It can be een that thi curve can be recaled to collape onto the α-decayofthe correlatorat φ c =.5for both wavevector at long time. Thi indicate that the von Schweidlerexponentiveryimilarinbothcae,butalopoint out the effect of the high order ingularity. According to MCT, b hould decreae a the ingularity i approached, but thi behaviour i not oberved in our cae. Comparion of the α-decay by fitting the KWW tretched exponential to the mater function i troubleome, ince only the late decay i obtained unambiguoly. The imilarity of both α-decay wa ued in the von Schweidler analyi of the correlation function, and only the non-ergodicity parameter and amplitude were fitted. Since the upper part of the decay i known to be affected by the higher order ingularity cloe-by, that part mut be dicarded in the fitting. The correlation function and fitting are preented in Figure 4 for the tate φ c =.5 and φ p =.39, for the ame wave-vector a Figure 5. The main concluion i that the late α-decay at all wavevector can be correctly decribed by the von Schweidler decay, with the ame exponent a the tate at φ c =.4. The non-ergodicity parameter obtained from the fitting are lightly lower than thoe of φ c =.4, but imilar within the error bar. According to MCT, fq decreae when approaching the gla part of the non-ergodicity line(ignaling an increae in the localization length). Our reult i thu conitent with thi prediction. With thee value of the non-ergodicity parameter one can define alo the wave-vector dependent time cale, τ q, a dicued above. In order to tet the value of the von Schweidler exponent, uing eq. () we have performed a three parameter fitting to obtain γ and φ G p. In Figure 5, τ q i preented a a function ofφ G p φ p for different wavevector. The power-law fitting for two wave-vector are alo plotted and the critical polymer fraction φ G p i given. The value of γ obtained from thi analyi for different wave-vector range from γ = 3.37 to γ = 3.82, the mean value being γ = 3.7. Thi value of γ implie a maller von Schweidler exponent, b =.33, in diagreement with our previou etimate, but backing the MCT prediction. Uing the ame value of φ G p, the vanihing of the elf-diffuion coefficient, D, can be analyed, and i preented in the inet to thi figure. A power-law i oberved in thi cae, with an exponent, γ =.92, which, again in contradiction to MCT, leave u with a big difference between the two value of γ. The diffuion coefficient in the inet of figure 5 again

12 2 τ q γ = 3.8 φ p G =.3975 D φ q G - φp -2 - φ q G - φp γ = 3.4 γ=.92 FIG.5: Wave-vectordependenttimecale, τ q, v. φ G p φ p for different wave-vector; ymbol from top to botton correpond to q = 3.9, 6.9, 9.9, 5, 2, 25, 3. The line are power law fitting to q = 3.9 and q = 3. φ G p and γ were fitted. indicate the re-entrat gla tranition. They decribe a maximum, more pronounced than that oberved in Figure 9 becaue the gla line i cloer to the φ c =.5 iochore. The minimum in τ q, which i oberved only for q = 3.9, in an equivalent way indicate the hape of the non-ergodicity tranition line. At higher wave-vector, theglatranitioncaueverylowf q andthetimecale merge with the microcopic tranient and thu thi feature i uppreed. The wave-vector dependence of τ q can alo be tudied, a done for the lower concentration, yielding another etimate of b. In thi cae, a imilar plot a Figure i obtained, wherethe lowq regionicompatible with aq 2 behaviour, and a higher exponent at higher q, yielding a value of b =.38. Thi value i in agreement with the nice comparion between the φ c =.4 and φ c =.5 iochore, but not with γ or the MCT prediction. We may then conclude that analyi of thi tate i extremely difficult, but our indication tate that the von Schweidler exponent i imilar for both packing fraction, but probably lightly lower in the higher concentration. Finally, we would like to point out that the non- Gauian parameter at thi packing fraction how a behaviour imilar to that hown in Figure 2, i.e. the peak i a high, and no hort-time caling i oberved. V. CONCLUSIONS In thi paper, by mean of imulation, we have teted the univeral prediction of MCT for gelation in colloidal ytem, viewed a an attraction driven gla tranition. The elf part of the intermediate cattering function for tate cloe to thi tranition have been analyed and the reult were compared with the theoretical prediction. For the φ c =.4 iochore, which i far enough from the high order ingularity, the correlation function can be α-caled. The time cale of the α-decay wa hown to obey a power law divergence, with an exponent, γ, related to the von Schweidler exponent, obtained from the early α-decay. Both feature are predicted by MCT for all non-ergodicity tranition. Alo, the wave-vector analyi of the time cale follow the behaviour predicted by MCT, with a mall difference in the value of the von Schweidler exponent. The wave-vector analyi of the correlation function depend on detail of the interaction potential, and thu provide information about the mechanim, leading to the tranition. In our cae, it etablihe that the gel tranition i driven by a hort-range mechanim, namely, bond formation, a oberved in the pair ditribution function. Additionally, it ha been hown that the KWW tretched exponential can account for the α-decay of the correlation function, a in other non-ergodicity tranition. We have alo teted the Gauian approximation, which work very well for the HSS. The non-gauian parameter, α 2, etablihe that thi approximation i much wore in the cae of the gel tranition than for the gla tranition. It wa alo teted when comparing the etimated localization length from the non-ergodicity parameter with the MSD of the particle. The diffuion coefficient ha been alo tudied. It tend to zero a the tranition i approached following a power-law, with an exponent much lower than γ, in accordance with imulation of gla tranition in other ytem, but in diagreement with MCT, where both exponent are equal. Finally, when the colloid concentration i increaed, the ytem how ignature of the high order ingularity nearby and little can be dicued about the exponent b or γ. However, only light change in the number are expected, ince the qualitative behaviour i reproduced, except for the α-caling. Alo, the diffuion coefficient follow a power law with a different exponent and the non-gauian parameter reache value imilar to the φ c =.4 cae. Therefore, our main concluion i that MCT account for mot feature of the imulated ytem on approach to the gel tranition, but the dicrepancie already found in other non-ergodicity tranition(uch a the repuliondriven gla tranition in hard phere ytem) are alo obtained here. Acknowledgement The author thank W. Kob and R. Sear for ueful dicuion. A.M.P. acknowledge the financial upport by the CICYT (project MAT2-55-CO3-2). M.F. wa upported by the DFG under grant Fu 39/3.