Poer la rheology of generalised Maxell and Kelvin-Voigt models arxiv:56.927v3 [physics.bio-ph] 2 May 28, 2, P. Patrício ISEL - Instituto Superior de Engenharia de Lisboa, Instituto Politécnico de Lisboa, 959-7 Lisboa, Portugal. 2 CFTC - Centro de Física Teórica e Computacional, Universidade de Lisboa, Campo Grande, Edifício C8, P-749-6 Lisboa (Dated: January 5, 28) We analyse the rheological behaviour of the generalised Maxell and Kelvin-Voigt models ith a poer la spectrum of relaxation or retardation times, H /τ α or L τ α, respectively. When < α <, both models have the same eak poer la behaviour, hich occur in a variety of soft materials, including living cells. Beyond this region, for other values of α, the to models have different and distinctive rheological behaviours, converging continuously to the -element Maxell and Kelvin-Voigt models, respectively, for α < or α > 2. I. INTRODUCTION In a large variety of soft materials, the elastic and viscous moduli exhibit approximately the same eak poer la behaviour, ith the angular frequency of the oscillatory deformations, G () G () α ( < α < ). In diluted unentangled polymers, the poer la exponent is associated ith the relaxation modes of each molecular chain [ 4]. In pre-gel and post-gel polymers, it is associated ith percolation and netork self-similarity [5 2] (an extensive revie is found in [3]). The universality of the eak poer la behaviour, even for systems devoid of any obvious fractal structure, like foams, emulsions, pastes and slurries, led Sollich and coorkers to propose the phenomenological Soft Glassy Materials (SGM) model [4, 5]. Based on the idea of structural disorder and metastability, this model relates the poer la exponent to a mean-field noise temperature, closely above the glass transition. In the context of cell mechanics, the eak poer la behaviour of the cytoskeleton has been recently associated ith the Glassy Worm Like Chain model [6, 7], hose relaxation modes follo essentially the Worm Like Chain model for semiflexible polymers, of hich the cytoskeleton is composed, but retain nevertheless some analogy ith the SGM model (see [8] for a revie). To understand this behaviour, the different models use a microscopic description of the materials to propose a specific spectrum of relaxation or retardation modes. Thus, they may be directly associated ith the generalised Maxell or Kelvin-Voigt models (see Fig. ) from the linear theory of viscoelasticity [9, 2]. In fact, it is knon that e may mathematically translate the relaxation spectrum of the generalised Maxell model to the retardation spectrum of the generalised Kelvin-Voigt model. The general relation is hoever non trivial and, in some situations, a spring and/or a dashpot should be added, in parallel for the generalised Maxell model and Electronic address: pedro.patricio@adf.isel.pt in series for the generalised Kelvin-Voigt model [2]. In this article, e ill firstly highlight the fundamental aspects of these interesting and universal eak poer la behaviours, in the frameork of the linear theory of viscoelasticity. In the folloing sections, e ill make a parallel excursion to the rheological behaviour of the generalised Maxell and Kelvin-Voigt models, ith a poer la spectrum of relaxation or retardation times, H /τ α or L τ α, respectively. When < α <, e ill see that these models lead to approximately the same eak poer la behaviour, so often encountered in nature. For other values of α, the to models have different and distinctive rheological behaviours, converging continuously to the -element Maxell and Kelvin-Voigt models, respectively, for α < or α > 2. II. WEAK POWER LAW RHEOLOGY In the linear theory of viscoelasticity [2], the relaxation function G(t) is defined through the force response σ(t) = ε G(t) () to an applied step deformation ε = ε Θ(t) (Θ(t) is the step function). Analogously, the creep compliance J(t) is defined through the deformation response ε(t) = σ J(t) (2) to an applied step force σ = σ Θ(t). More generally, if e impose an arbitrary deformation ε(t), hich may be vieed as an infinite sum of small successive step deformations dε(τ), ith < τ < t, then Boltzmann superposition principle states that the total force exerted on the system may be determined by: σ(t) = t G(t τ)dε(τ) = t G(t τ) dε(τ). (3) Similarly, for an imposed general force σ(t), the total deformation of the system is given by: ε(t) = t J(t τ) dσ(τ). (4)
2 η η 2 η N G G 2 G N function is necessarily a decreasing function of time, implying α >. For α <, its Laplace transform is given by G(s) s α, the complex modulus becomes G () (i) α = α e iπα/2, yielding G () G () α, G G = tan πα 2 (9) It should be noted that, in these conditions ( < α < ), J(s) /s α+ and the creep compliance, hich is necessarily an increasing function of time, is given by G η FIG. : Generalised Maxell (top) and Kelvin-Voigt (bottom) models. The generalised Maxell model is composed of a collection of Maxell elements in parallel, each made out ( < α < ) of a dashpot of drag coefficient η i in series ith a spring of stiffness G i. Each Maxell element corresponds to a different relaxation mode, ith relaxation time τ i = η i/g i. The generalised Kelvin-Voigt model is composed of a collection of Kelvin-Voigt elements in series, each made out of a dashpot in parallel ith a spring. Each Kelvin-Voigt element corresponds to a different retardation mode, ith retardation time τ i = η i/g i. G N η N J(t) t α. () These are the fundamental features of the eak poer la rheology. A system ith this type of behaviour is considered to be solid like henever < α <.5 (because G /G < ) and fluid like henever.5 < α < (as G /G > ). In the folloing sections, e ill see that these eak poer la behaviours, ith ( < α < ), are found in systems described by the generalised Maxell or Kelvin- Voigt models ith a poer la spectrum of relaxation or retardation times, H /τ α or L τ α, respectively. These spectra are hoever ell defined for all values of α. For α <, or α >, e ill see that these to to models have different and distinctive rheological behaviours, converging continuously to the -element Maxell and Kelvin-Voigt models, respectively, for α < or α > 2. If the system is subjected to an oscillatory deformation ε(t) = ε e it or force σ(t) = σ e it, to important rheological assays, its response is given, respectively, by σ(t) = ε G e it, (5) here G = G + ig is the complex modulus, or, ε(t) = σ J e it, (6) here J = J ij is the complex compliance. G (J ) and G (J ) are defined as the elastic and viscous moduli (compliances). The Boltzmann superposition principle, expressed by Eq. (3) and (4), allos us to relate the different rheological functions through: G () = (i) G(i) = (i) J(i) = J () (7) here the Laplace transform is given by f(s) = f(t)e st dt. From these relations it may be made clear that any relaxation function ith a eak poer la behaviour: G(t) t α ( < α < ) (8) ill necessarily lead to the eak poer la behaviours of the elastic and viscous moduli. Indeed, the relaxation III. GENERALISED MAXWELL MODEL We consider first a generalised Maxell model (see Fig., top) composed of a parallel array of N Maxell viscoelastic elements. Each Maxell element is composed of a dashpot and a spring in series. These are the basic elements of the linear theory of viscoelasticity. A spring of stiffness G follos Hooke s la, σ = Gε, and a dashpot of drag coefficient η follos Neton s la, σ = η ε = ηdε/dt. As a result, the Maxell element behaves elastically on small times scales and as a fluid on long ones. Its relaxation function is G(t) = Ge t/τ, here τ = η/g is the element relaxation time. Because the Maxell elements are in parallel, the total relaxation function is just the sum: G(t) = G i e t/τi () This sum may be approximated by an integral by multiplying it by di =. If di is considered to be small, e have: G(t) τn τ G τ e t/τ di (2) The derivative di/ corresponds to the density of modes (per unit relaxation time). Choosing the change of variables ζ = /τ, and for t beteen the minimum and
3 maximum relaxation times, and τ max, this integral approaches a Laplace transform. To establish the eak poer la behaviour G(t) /t α (ith < α < ), the relaxation spectrum must scale as H(τ) G τ di τ τ α (3) As ith other rheological functions, the relaxation spectrum of a viscoelastic fluid characterises its rheological behaviour. Because H(τ) is defined in terms of an arbitrary choice of stiffnesses, G τ, and density of modes, di/, the relaxation spectrum is not subjected to any constraint, and α may assume any value. To understand the rheological behaviour of a material ith a general poer la spectrum of relaxation times, e ill determine the corresponding elastic and viscous moduli. After performing a Laplace transform and using relations (7), e have G () G () τ α (τ) 2 τ α + (τ) 2 τ (4) τ + (τ) 2 τ. (5) For /τ max, the integrals (4) and (5) are much simplified, and e obtain G 2 and G. For /, e have simply G and G. These scalings correspond to the scalings of the -element Maxell model, both for small and large. For /τ max /, e may extend the limits of the integrals (4) and (5) to and τ max, hich allo us to obtain the results: G () π 2 α csc πα 2 (if < α < 2) (6) one Maxell element, corresponding respectively to the maximum or minimum relaxation times. To plot representative elastic and viscous moduli vs, for different values of α, e used the discrete generalised Maxell model, ith ell defined N relaxation modes. The poer la relaxation spectrum of Eq. () may be implemented from to independent functions, G τ and di/. So, one of them may be chosen arbitrarily. α =.5 3 ' 2 α =.25 3 α =.5 2 ' ' α =.5 3 ' α =.25 ' α =.5 3 ' G () π 2 α sec πα 2 (if < α < ). (7) When < α <, both integrals are ell defined and e recover the eak poer la behaviours of Eq. (9), as expected. If α <. the integral (4) for the elastic modulus diverges as. The elastic modulus is then dominated by the smallest relaxation time, and e have G. On the contrary, if α > 2, then this integral is dominated by the largest relaxation time, and e have G 2. By the same line of reasoning, e may determine from Eq. (5) the poer la behaviours G for α < and G for α >. Applying these results to each interval of α, e obtain the interesting poer la behaviours: 4 α =.75 2 ' α =.5 ' 5 α =.75 ' 6 α =.5 ' 6 8 G G (α < ) (8) G G α ( < α < ) (9) G α G α ( < α < ) (2) G α G ( < α < 2) (2) G 2 G (2 < α) (22) We note that for α < or α > 2, the rheological response of the material is entirely dominated by only FIG. 2: Elastic and viscous moduli, G () (solid line) and G () (dashed line) (in units of G) obtained from the generalised Maxell relaxation spectrum of Eq. (23) (left) and the generalised Kelvin-Voigt retardation spectrum of Eq. (39) (right), for α =.5,.25,.5,.75,.5 (top to bottom), = 3, τ max = 3 (in arbitrary units), and N = 3.
4 An interesting choice consists of assuming the poer la behaviour G i = Gi β for the stiffnesses. In this case, the poer la spectrum of relaxation times is recovered through: τ i = τ i (+β)/α (G i = Gi β ), (23) independently of our choice of β (the density of modes di/ is easily obtained from Eq. (23)). Mathematically, this choice allos us to define not only the time relaxation limits, τ max,, but also the exact number of elements N e ant to use, for any α. We just have to take τ = τ max and β = (α ln(τ max / )/ ln N). If e use the Laplace transform of Eq. (), e obtain, ithout approximation, the ell-knon elastic and viscous moduli of the generalised Maxell model, for any choice of discrete relaxation spectra: G () = G () = G i (τ i ) 2 + (τ i ) 2 (24) G i τ i + (τ i ) 2. (25) Figure 2 (left) shos the elastic and viscous moduli, G () and G () for α =.5,.25,.5,.75,.5, = 3, τ max = 3 (in arbitrary units), and N = 3. To obtain these plots, e have used the sums of Eq. (24) and (25), together ith the discrete relaxation spectrum defined by Eq. (23). As it may be seen from this figure, e recover the global behaviour previously described in the analysis of the continuous relaxation spectrum. Hoever, the referred poer la exponents (including those of Eq. (9)) correspond only to an approximation (they ould be valid for and τ max ), e sho in Fig. 3 the numerically converged (e have increased N to verify it) computed values, established in the interior range < <. It is important to mention that the exponent of G is alays greater than the exponent of G, even in the region < α <, here they assume almost identical values. This is a distinctive feature of the generalised Maxell model. It should be noted that e have performed the same calculations for different numbers of elements. The possibility of changing β allos us to have a good approximation of the poer la behaviours even for a reduced number of elements (N > ). Nevertheless, as the number of elements diminishes, some oscillations of the elastic and viscous moduli start to appear. IV. GENERALISED KELVIN-VOIGT MODEL We no turn to the generalised Kelvin-Voigt model (see Fig., bottom) composed of a serial array of N 2.5.5.5 2 2 3.5 2 2 3 FIG. 3: Exponents x and y vs α of the poer las G x and G y, obtained from the generalised Maxell relaxation spectrum of Eq. 23 (top) and the generalised Kelvin- Voigt retardation spectrum of Eq. 39 (bottom), ith = 3, τ max = 3 (in arbitrary units), and N = 3. These exponents ere computed numerically, for every α, in the interior range < <. Kelvin-Voigt viscoelastic elements. Each Kelvin-Voigt element is composed of a dashpot of drag coefficient η i in parallel ith a spring of stiffness G i. It behaves elastically on long times scales and its dynamics comes from the viscous element. Its creep compliance is J i (t) = J i ( e t/τi ), here J i = /G i and τ i = η i /G i is the element retardation time. Because the Kelvin-Voigt elements are in series, the total creep compliance is just the sum: J(t) = x y y x α α J i ( e t/τi ) (26) This sum may be approximated by the integral: J τn τ J τ ( e t/τi ) di (27) If e ant to establish the eak poer la dependence J GK t α, the retardation spectrum must then scale as L(τ) J τ di τ τ α (28) As before, the exponent α of the retardation spectrum may assume any value. To understand the rheological behaviour of a material ith a general poer la spectrum of retardation times, e ill determine the corresponding elastic and viscous moduli. After performing a Laplace transform and using
5 relations (7), e first obtain J () J () τ α + (τ) 2 τ τ α τ + (τ) 2 τ (29) (3) The elastic and viscous moduli may then be determined from the relations: G = J J 2 + J 2, J G = J 2 + J 2 (3) For /τ max, the integrals of Eq. 29 and 3 are much simplified, and e obtain the scalings J and J, ith J J. The elastic and viscous moduli scale as G /J and G J /J 2. For /, the same integrals are again simplified, yielding the scalings J 2 and J, ith J J. We obtain then the same scalings for the elastic and viscous moduli: G J /J 2 and G /J. The coefficients of G and G are dependent of α but their scalings are not. In fact, the exponents of coincide ith the scalings of a simple Kelvin-Voigt model. For /τ max /, e may in some cases extend the limits of the integrals of Eq. (29) and (3) to and τ max, hich allos us to obtain the results: J () π 2 α csc πα 2 J () π 2 α sec πα 2 (if < α < 2) (32) (if < α < ) (33) When < α <, both integrals are ell defined and e recover the eak poer la behaviours of Eq. (9). If α <, the integral for the elastic compliance (Eq. (29)) diverges as. The elastic compliance is then dominated by the smallest retardation time. In this case, e have τ n, and J. On the contrary, if α > 2, then this integral is dominated by the largest retardation time, τ max. We have τ, and the poer la behaviour J 2. By the same line of reasoning, e may determine from Eq. (3) the poer la behaviours J for α < and J for α >. Applying these results to each interval of α, and using the approximations J J for α < and J J for α > (hich e may infer from Eq. (32) and (33)), e obtain the poer la behaviours: G G (α < ) (34) G G α ( < α < ) (35) G α G α ( < α < ) (36) G α G ( < α < 2) (37) G G (2 < α) (38) We note that for α < or α > 2, e have a single poer la behaviour for G and G, for all values of the angular frequency. Indeed, for these ranges of α, the hole structure is entirely dominated by only one Kelvin-Voigt element, corresponding respectively to the minimum (α < ) or the maximum (α > 2) retardation times. To plot representative elastic and viscous moduli vs, for different values of α, e used the discrete generalised Kelvin-Voigt model, ith N ell defined retardation modes. Again, the poer la retardation spectrum of Eq. (28) may be implemented by to independent functions, J τ = /G τ and di/. So, one of them may be chosen arbitrarily. If e choose the poer la behaviour J i = Ji β the condition given by Eq. 28 implies the folloing discrete spectrum of retardation times: τ i = τ i (+β)/α (J i = Ji β ) (39) For a given set of τ max,, α, e are able to construct the eak poer la behaviours of Eq. (9) from a generalised Kelvin-Voigt model ith N modes by making τ = and β = (α ln(τ max / )/ ln N). The exact elastic and viscous compliances are given through the Laplace transform of Eq. (26): J () = J () = J i + (τ i ) 2 (4) J i τ i + (τ i ) 2 (4) In order to determine the elastic and viscous moduli e use Eq. (3). Figure 2 (right) shos the elastic and viscous moduli, G GK () and G GK () for α =.5,.25,.5,.75,.5, = 3, τ max = 3 (in arbitrary units), and N = 3. To obtain them, e have used the sums of Eq. (4) and (4), together ith the discrete retardation spectrum defined by Eq. (39). As it may be seen, e recover the global behaviour previously described in the analysis of the continuous retardation spectrum. As these poer la exponents correspond only to an approximation, only valid as and τ max, e sho in Fig. 3 the numerically converged computed values, established in the interior range < <. We note that the exponent of G is alays greater than the exponent of G. This is a distinctive feature of the generalised Kelvin-Voigt model. This inversion in the values of the slopes is reflected also in the intersection points of G and G (see Fig. 2, right and left). If the intersection takes place at small values of, for the generalised Maxell model (hen α <.5), it occurs at large values of for the generalised Kelvin-Voigt model (for the same value of α) and vice versa. V. CONCLUSION In this article, e have focused on the interesting eak poer la behaviours G () G () α ( < α < )
6 that occur in a large variety of soft materials. They may appear in diluted, fluid systems, encompassing both the Bronian motion of elastic chains in a solvent, modelling unentangled polymers, or the important SGM model, that describes fluid behaviour above the glass transition. The eak poer la behaviours may also occur in gelled, elastic systems, made out of polymeric netorks, hich retain to important characteristics of the cell s cytoskeleton, namely the idea of an existing prestress on its fibres, and the idea of a possible fractal structure, characterised by a poer la distribution of elastic element lengths and retardation times [22]. The cytoskeleton is hoever a very complex active structure, and exhibits other important rheological properties that are not considered here, such as stress stiffening [23] or a ne poer la exponent (of 3/4) for higher frequencies, usually associated ith the rheological behaviour of diluted semi-flexible filaments [24]. Particularly interesting, e find that beyond the region corresponding to the eak poer la behaviours (α < or α > ), the generalised Maxell and Kelvin- Voigt models give very different and distinctive rheological behaviours, converging continuously to the - element Maxell and Kelvin-Voigt models, respectively, for α < or α > 2. Acknoledgements We thank C. R. Leal, J. Duarte, C. Januário, J. M. Tavares and P. I. C.Teixeira for stimulating discussions. We thank P. I. C. Teixeira for linguistic assistance. [] P. E. Rouse Jr, The Journal of Chemical Physics 2, 272 (953). [2] B. H. Zimm, The journal of chemical physics 24, 269 (956). [3] M. Doi and S. F. Edards, The theory of polymer dynamics, vol. 73 (oxford university press, 988). [4] M. Rubinstein and R. Colby, Polymers Physics (Oxford, 23). [5] P.-G. De Gennes, Journal de Physique Lettres 37, (976). [6] M. Cates, Journal de Physique 46, 59 (985). [7] M. Muthukumar, The Journal of chemical physics 83, 36 (985). [8] H. H. Winter and F. Chambon, Journal of Rheology (978-present) 3, 367 (986). [9] F. Chambon and H. H. Winter, Journal of Rheology (978-present) 3, 683 (987). [] M. Muthukumar, Macromolecules 22, 4656 (989). [] N. Goldenfeld and P. Goldbart, Physical Revie A 45, R5343 (992). [2] F. Tanaka, Polymer Physics: Applications to Molecular Association and Thermoreversible Gelation (Cambridge University Press, 2). [3] H. H. Winter and M. Mours, in Neutron spin echo spectroscopy viscoelasticity rheology (Springer, 997), pp. 65 234. [4] P. Sollich, F. Lequeux, P. Hébraud, and M. E. Cates, Physical revie letters 78, 22 (997). [5] P. Sollich, Physical Revie E 58, 738 (998). [6] K. Kroy and J. Glaser, Ne Journal of Physics 9, 46 (27). [7] L. Wolff, P. Fernandez, and K. Kroy, Ne Journal of Physics 2, 5324 (2). [8] R. H. Pritchard, Y. Y. S. Huang, and E. M. Terentjev, Soft matter, 864 (24). [9] R. Christensen, Theory of viscoelasticity: an introduction (Elsevier, 22). [2] N. Phan-Thien, Understanding Viscoelasticity: An Introduction to Rheology (Springer Science & Business Media, 22). [2] R. S. Lakes, Viscoelastic solids, vol. 9 (CRC press, 998). [22] P. Patrício, C. Leal, J. Duarte, and C. Januário, Physical Revie E 92, 472 (25). [23] P. Fernández, P. A. Pullarkat, and A. Ott, Biophysical journal 9, 3796 (26). [24] L. Deng, X. Trepat, J. P. Butler, E. Millet, K. G. Morgan, D. A. Weitz, and J. J. Fredberg, Nature materials 5, 636 (26).