Fractional rheology of colloidal hydrogels with cellulose nanofibers

Colloidal gels are soft solids composed of particles dispersed in a fluid phase. Their rheological behavior highly depends on the particle concentration, but establishing a relationship can be challenging. This article showcases the potential of fractional rheology to model and predict linear viscoelastic responses of colloidal hydrogels containing TEMPO-oxidized cellulose nanofibers. Cellulose nanofiber hydrogels are soft solids whose rheology is directly related to the particle concentration. Therefore, this work defined the rheological behavior of the hydrogels using a fractional order derivative analytically solved to determine rheological responses in frequency, stress relaxation, and creep. Using two parameters, it evaluated the rheology of cellulose nanofiber hydrogels and established tests that predict rheological behaviors for given particle concentrations. The findings suggested that the fractional approach could become a standard method for characterizing cellulose nanofiber hydrogels in the reported concentration regime. The two parameters of the fractional model build a comparison framework to assess the rheology of different viscoelastic materials.


Introduction
Colloidal gels are viscoelastic materials with a soft solid texture commonly used in the food, cosmetics, and pharmaceutical industries.Colloidal gels are crucial for dispersing solid components that may otherwise settle quickly.As metastable solids, colloidal gels have a low particle density, forming a rigid structure with an interconnected fractal-like network (Jiang and Seto 2023;Angelini et al. 2023;Rouwhorst et al. 2020).Describing the rheology of colloidal gels is challenging because their stress-strain response is greatly influenced by their complex microstructures (Bantawa et al. 2023;Faber et al. 2017a, b).In real-world situations, not fully grasping the rheology of colloidal gels creates obstacles in material processing and product design.Simply put, there is no arbitrary way to connect rheological behavior with material composition.
Typically, the rheology of viscoelastic materials can be described phenomenologically.Existing rheological models consider the viscoelastic behavior as a mixture of elasticity and viscosity using two rheological elements: springs and dashpots.A spring represents an elastic solid with a stress-strain relationship given by Hooke's law.In contrast, a dashpot stands for a fluid with a Newtonian stress-strain relationship.Although the spring-dashpot approach yields rational representations of rheological phenomena, for complex materials such as colloidal gels, it requires many parameters (Bonfanti et al. 2020).Assigning physical meaning to the resulting parameters constrains the potential universality of the spring-dashpot approach.
Cellulose nanofiber (CNF) hydrogels are an example of colloidal gels with complex rheological behavior.These gels are made from cellulose nanofibers, nanocelluloses with a slender shape that can extend to several microns (Isogai et al. 2011).Despite their small size and semicrystalline structure, CNF particles exhibit remarkable elasticity, making them helpful in, for example, designing hydrogels for biomedical applications, where they serve as a dispersion medium for drugs and other nanoparticles (Fneich et al. 2019;Hubbe et al. 2017;De France et al. 2017;Kontturi et al. 2018;Zeng et al. 2021;Reichler et al. 2021).Characteristically, CNFs are obtained as aqueous colloidal systems that tend to percolate into hydrogels even at low concentrations (Hill 2008;Pääkkö et al. 2007).Their rheology is known for displaying powerlaw signatures over a wide range of conditions and concentrations ( ) , which can be challenging to describe mathematically without using recursive constitutive 1 3 Vol.: (0123456789) equations and numerous spring-dashpot arrangements (Pääkkö et al. 2007;Hill 2008;Koponen 2020;Naderi and Lindström 2014;Šebenik et al. 2019).
Particularly, describing the rheology of CNF hydrogels is constrained to empirical equations (Pääkkö et al. 2007;Naderi and Lindström 2014).Therefore, characterizing the rheology of CNF hydrogels requires a different mathematical framework than the common spring-dashpot approach.Scott Blair, a pioneer in establishing the foundations of rheology, suggested that certain viscoelastic materials show a rheological response intermediate between elasticity and viscosity (Scott Blair et al. 1947;Scott Blair 1969).For those materials, Scott Blair proposed using fractional calculus and introduced the concept of "principle of intermediacy".Using fractional calculus has allowed the creation of models that better align with experimental rheological data.This, in turn, has improved the understanding of viscoelastic phenomena on a physical level (Ciancio et al. 2023;Tarasov and Tarasova 2020).
The principle of intermediacy defines the viscoelastic behavior with a fractional order derivative that generalizes Hooke's and Newton's laws.Nowadays, the principle of intermediacy translates to fractional rheology, the application of fractional calculus to the viscoelasticity of materials (West et al. 2003).Fractional rheology has great potential for describing viscoelastic systems with complex microstructures (Bantawa et al. 2023;Bonfanti et al. 2020;Faber et al. 2017b;Owens et al. 2022).Yet, to our knowledge, no work has used fractional rheology to comprehend the viscoelasticity of CNF hydrogels.Describing the rheology of CNF hydrogels with constitutive equations rather than with empirical relationships would shed light on understanding the hierarchical nature of colloidal systems.
This article proposes using fractional rheology to characterize the rheological properties of colloidal hydrogels with TEMPO-oxidized CNFs.For that, the principle of intermediacy introduces a rheological element called springpot.The springpot constitutive equation describes the degree of viscoelasticity ( ) of a material by varying the fractional order of the derivative between 0 and 1 (Koeller 1984).In particular, the springpot can accurately describe timescales for materials exhibiting power-law signatures.For instance, casein food gels were characterized by Faber et al. (2017b) and Keshavarz et al. (2017); the authors demonstrated that the power-law rheological behavior of such protein gels could be described using the principle of intermediacy.Similar to the former example, the literature showcases multiple materials whose rheology was successfully characterized by including a springpot in the mathematical framework (Alcoutlabi and Martinez-Vega 1998;Bantawa et al. 2023;Bonfanti et al. 2020;Faber et al. 2017a, b;Reyes-Melo et al. 2008).
This article hypothesizes that one springpot can describe and predict the viscoelastic behavior of CNF hydrogels within the linear viscoelastic regime (LVR).In its simplest form, the springpot constitutive equation is called the Scott-Blair (SB) model.This paper uses the SB model to fit the material functions storage [ G � ( ) ] and loss [ G �� ( ) ] modulus measured from dynamic mechanical analysis (DMA) with only two parameters: and . is a quasi-property with Pa s units related to the gel firmness, while defines the degree of viscoelasticity.To show the potential of the principle of intermediacy, we computed the stress relaxation and creep response of the CNF hydrogels from the frequency spectrum, indicating that the data obtained with the SB model approximated the experimental results.We demonstrated that the SB model allowed for building a framework to predict rheological responses given a CNF concentration.Our findings showed us at a glance that the SB model could be used as a standard model to characterize the rheology of CNF hydrogels in the reported concentration regime and therefore ease material design and processing operations.

Experimental section
The CNF hydrogels are made from TEMPO-oxidized cellulose nanofibers obtained from Kraft softwood pulp according to the methods described in the literature (Leppiniemi et al. 2017;Pöhler et al. 2020;Skogberg et al. 2017).Figure S1 shows the structure of a CNF hydrogel under polarized optical microscopy; the image highlights the broad particle size distribution of the fibers.VTT Technical Research Centre of Finland kindly provided the CNFs as a hydrogel with 1.16 wt% fiber concentration (Koochi et al. 2023).To assess the effect of CNF concentration ( ) in the hydrogel, a dilution series of 0.90, 0.70, 0.50, and 0.30 wt% was prepared by adding deionized water to the 1.16 wt% CNF hydrogel.All the hydrogels Vol:.( 1234567890) were homogenized mechanically and then rested at 4 • C for three days to coarsen any bubble formed during the homogenization.After this time, we evaluated the rheology of the samples.

Small strain shear rheology
Rheological characterization of the CNF hydrogels consisted of three time-dependent experiments: dynamic mechanical analysis (DMA), stress relaxation, and creep.These three experiments were conducted on hydrogels with 0.30, 0.50, 0.70, 0.90, and 1.16 wt% CNF concentration.The experiments took place on an Anton Paar (Austria) modular compact rheometer (MCR) 302 equipped with a Couette geometry.The geometry comprises a lower cup and upper bob tool with serrated walls (Viitanen et al. 2023).For the experiments, the tested CNF hydrogel was slowly poured into the cup, avoiding bubble formation.Then, the bob tool was lowered to the measuring position.Before starting an experiment, the sample was allowed to rest for 15-20 min at 10 • C to ensure its thermal equilibrium and structure recovery.All the tests were carried out at a constant temperature since modifying the temperature showed a neglectable influence on the viscoelastic properties of the CNF hydrogels, see Fig. 1a.Accordingly, to avoid drying out the CNF hydrogels, all experiments were at a temperature of 10 • C.

Dynamic Mechanical Analysis (DMA):
The rheometer measured the material functions storage and loss modulus ( G ′ and G ′′ ) in response to angu- lar frequency ( ) in a range from 0.1 to 100 rad s −1 .The tests were performed at a constant temperature T = 10 • C and constant strain amplitude 0 = 0.2 % .The parameter 0 was within the linear viscoelastic regime (LVR), and its magnitude was chosen from strain amplitude tests on the CNF hydrogels with the lowest and highest concentrations, see Fig. 1b.During DMA, the angular frequency was increased logarithmically from = 0.1 rad s −1 to = 100 rad s −1 , and the rheometer set the data acquisition rate.We discarded the experimental points that the rheometer tagged as unsteady.For some authors, a common practice is to measure from high to low frequencies rather than from low to high. Figure 1c shows a DMA result from high to low frequencies.We use Fig. 1c to demonstrate that the material functions were independent of the measuring protocol.
Stress relaxation: On the rheometer, the sample was submitted to a step strain of 0 = 0.2 % for a time of t = 100 s.The device measured the relaxa- tion modulus [G(t)], a property that represents the ability of a material to relieve stress.G(t) is defined as (t)∕ 0 , where (t) is the stress signal as a func- tion of time.The tests were carried out at a constant temperature T = 10 • C. Data acquisition was loga- rithmic from t = 0.1 s to t = 100 s.The lower time Fig. 1 Pre-characterization of the cellulose nanofiber (CNF) hydrogels for determining measuring conditions.a Dynamic mechanical thermal analysis on a CNF hydrogel with concentration ( ) of 0.90 wt%.The results show that the material functions G � (T) , G �� (T) show a neglectable dependence on the temperature (T).No clear effects of hysteresis are observed from heating and cooling cycles at 1 • C min −1 at a constant angular frequency ( ) of 10 rad s −1 and strain amplitude ( 0 ) of 0.2 %. b Strain tests on samples with = 0.30 wt% (empty and solid circles of lower magnitudes) and = 1.16 wt% (empty and solid circles of higher magnitudes), = 10 rad s −1 .The plot also shows the chosen 0 (0.2 %) value for dynamic mechanical analysis (DMA) and stress relaxation.c DMA measured from high to low frequencies on a CNF hydrogel with = 1.16 wt% ( 0 = 0.2 %) Vol.: (0123456789) limit corresponds to the timescale that the rheometer can accurately record.
Creep: Similar to the stress relaxation protocol, creep tests were carried out on the rheometer to measure the creep compliance [J(t)] of the CNF hydrogels.A step stress ( 0 ) of 5 Pa or 20 Pa (depending on the sample concentration: 5 Pa for the lowest concentration and 20 Pa for the highest) was generated for t = 100 s, followed by a recovery phase of t = 100 s, where the stress was relieved.Since the magnitude of 0 highly affects the creep response, we chose a small 0 to deform the sample elastically (Owens et al. 2022).J(t) is defined as (t)∕ 0 , and the test was carried out at a constant temperature T = 10 • C. Data acquisition was loga- rithmic from t = 0.1 s to t = 100 s, following the same protocol used during stress relaxation.

Principle of intermediacy: Scott-Blair model
To describe the rheology of the CNF hydrogels, we used the Scott-Blair (SB) model based on the principle of intermediacy.As Fig. 2a illustrates, the SB model defines the shear stress ( ) signal of a viscoelastic material as an intermediate behavior between a spring (Hooke's law for elastic solids) and dashpot (Newton's law for viscous liquids).Such intermediate behavior is represented by a rheological element called springpot, and Eq. 1 presents its constitutive equation as a fractional order derivative (Koeller 1984), Equation 1 generalizes Hooke's and Newton's laws.The fractional order ( ), which ranges between 0 and 1, specifies whether the rheology of a (1) = d dt .sample approximates a solid-like ( ≈ 0 ) or liquid- like ( ≈ 1 ) behavior.Furthermore, indicates the time-dependence response of the viscoelastic sample (Faber et al. 2017a).In Eq. 1, t is the time, and is the shear strain.The factor is a quasi-property, namely quasi-modulus, with units of Pa s .When = 0 , has units of elasticity; has units of viscos- ity when = 1 .is correlated with the material firm- ness.It is important to remark that others have defined the springpot as a three-parameter element, introducing the concept of relaxation time (Friedrich et al. 1999;Schiessel et al. 1995).However, it is not possible to design an experiment to measure the relaxation time, which constrains the practical use of the threeparameter springpot (Bonfanti et al. 2020;Jaishankar and Mckinley 2013).Accordingly, the two-parameter version of the springpot is preferred as it eases fitting the springpot parameters to experimental data.Next, we present a succinct description of how Eq. 1 can be solved to describe the material properties of viscoelastic materials.The reader may refer to textbooks by Mainardi (2010) and West et al. (2003) to obtain more details about the springpot fractional derivative.
Equation 1 has analytical solutions for describing G � ( ) and G �� ( ) in DMA, G(t) in stress relaxation, and J(t) in creep (Mainardi 2010).To find the analytical solutions, Eq. 1 can be redefined using the Caputo derivative as: where Γ is the Gamma function (Surguladze 2002).Equation 2 represents a response given by a Boltzmann memory integral with a power-law relaxation (Faber et al. 2017a;Koeller 1984;Torvik and Bagley 1984).The Caputo derivative states the assumption that the system is at rest before time 0, which reflects the reality of the test.
In DMA, G � ( ) and G �� ( ) result from shear oscil- lations; the material functions are frequency-dependent.The condition can be defined by applying the Fourier transform to the Caputo derivative of (t) , as Eq. 3 (Schiessel et al. 1995), The Fourier transform results in the complex modulus function G * ( ) given by Eq. 4, (2) Then, the same equation can be separated, providing a mathematical framework to characterize G � ( ) with Eq. 5, and G �� ( ) with Eq. 6 (Friedrich et al. 1999), On the other hand, the Laplace transform in Eq. 7, allows for establishing the G(t) and J(t) functions as given in Eq. 8, and Eq. 9, respectively (Bird et al. 1987;Faber et al. 2017a, b).The Laplace transform is a mathematical tool that eases obtaining time-dependent analytical solutions.

Modeling and prediction of materials functions
The DMA results of the CNF hydrogels, except the one with = 0.90 wt%, were fitted using Eq. 5 and Eq. 6.It is possible to fit only the material function G * ( ) ; however, fitting G � ( ) and G �� ( ) reduces the bias (Faber et al. 2017b).Fitting the DMA results yields two parameters, and , for each of CNF hydrogel.Then to show the potential of the Scott-Blair model, both parameters were replaced in Eq. 8 and 9 to compute the materials functions G(t) and J(t).These computations were compared to the experimental results measured during stress relaxation and creep.Finally, we use the classic scaling law approach to describe the dependence between the Scott-Blair model parameters and .From describing how and vary as a function of , we attempt to predict all the rheological responses for a given CNF concentration of = 0.90 wt%.

Rheological models: a comparison
We first present in Fig. 2b-d  In comparison, the SB model depicts in frequency a power-law behavior whose viscoelastic nature can be assessed using the parameter .In the case of Fig. 2d, it shows a more elastic response given by = 0.30.The response depicted in Fig. 2d is common in soft solids, such as cheese, bitumen, and some natural gums (Bonfanti et al. 2020;Faber et al. 2017b;Fagioli et al. 2018).Mainly, cellulose-based materials, which tend to percolate, also exhibit DMA responses that can be characterized using a springpot (Miranda-Valdez et al. 2023).Accordingly, we present in Fig. 2e one of the DMA results measured for a CNF hydrogel with = 0.50 wt%.One can see that the rheological response of the CNF hydrogel resembles the SB model curves in Fig. 2d.Recently, Bantawa et al. (2023) correlated the rheological behavior characterized using a springpot with the fractal and hierarchical microstructure in soft particulate gels.CNF hydrogels meet both features; they are soft particulate gels and display microstructures with fractal dimensions (Serra-Parareda et al. 2021).Intuitively, such structural characteristics of CNF hydrogels may reveal why the SB model approaches the experimental rheological response (Fig. 2d,e).
All in all, a qualitative comparison between the theoretical and experimental material functions, G � ( ) and G �� ( ) in Fig. 2d,e, showed at a glance that the linear rheology of CNF hydrogel can be described using the SB model.The efficacy of the SB model is based on the fact that the springpot parameterizes the relaxation spectrum, which in turn determines the distribution of characteristic timescales, unlike the spring-dashpot models that define specific timescales (Bonfanti et al. 2020;Suzuki et al. 2021).In other words, the springpot can be seen as an infinite fractal representation of spring-dashpot elements which can be correlated to the fractal scaling properties of gel clusters (Schiessel and Blumen 1993;Suzuki et al. 2021;Larson 1999).Furthermore, the concept aligns with the fractal structure of CNF hydrogels, as it suggests a broad range of lengths and timescales.In the following sections, we confirmed the relevance of using the principle of intermediacy to characterize other rheological responses, focusing on the relaxation modulus [G(t)] and creep compliance [J(t)].

Scott-Blair model: fitting and computation of material functions
All the results from DMA, stress relaxation, and creep are presented in Fig. S2  Vol:.( 1234567890) complement the discussion.The material functions of the CNF hydrogels agreed with the physical dispersion stability precondition since G � ( ) > G �� ( ) .In other words, the CNF systems met some of the characteristics of a viscoelastic solid (Winter and Mours 1996).Furthermore, we observed that G � ( ) > G �� ( ) by one decade, which is a characteristic behavior in CNF hydrogels (Dimic-Misic et al. 2013;Nechyporchuk et al. 2014).This precondition can be seen in the DMA results exemplified in Fig. 2e and Fig. 3a,b.What complicates understanding the rheology of CNF hydrogels is the slight dependency of G ′ and G ′′ on , meaning that the CNF hydrogels devi- ate from the ideal behavior of true gels (Winter and Mours 1996).Generally, when the material functions are frequency-dependent, they indicate an imperfect gel network (Rubinstein and Colby 2003).This is expected since CNF hydrogels, as colloidal gels, form a metastable solid structure rather than a chemically crosslinked one.
Regarding the results obtained using the SB model, we present a subset of DMA experiments on CNF hydrogels.Figure 3a,b show the SB model fittings of the DMA results on hydrogels with = 0.70 wt% and = 1.16 wt%.In addition to these DMA results, the SB model also fitted the experimental data for systems with = 0.50 wt% and = 0.30 wt%.All the experimental data were well described using the SB model since the mean absolute error (MAE) was small; see Table 1.In previous work, we showcased how the MAE using a springpot can be four times smaller than the analysis with a generalized Maxwell model with 200 parameters (Miranda-Valdez et al. 2023).Furthermore, Fig. 3c shows how the fitting parameters and changed as a function of .One must remark that at high frequencies, G �� ( ) deviated from the SB model; G �� ( ) slightly increased, indicat- ing a gradual change of the dangling gel structure as a function of frequency.Despite this deviation, the SB model offered an analytical description of G � ( ) and G �� ( ) for the CNF hydrogels.This allowed for com- puting other material functions in the CNF hydrogels and built a systematic approach to capturing their rheology.For example, the SB model could replace empirical equations that have been suggested in the literature to describe G � ( ) in CNF hydrogels and systems with similar rheological responses (Jampala et al. 2005;Kim and Yoo 2006;Naderi and Lindström 2014).What is remarkable from the results in Fig. 3a,b is that, to the best of our knowledge, this is the first time that the SB model has been used successfully to fit the material functions G � ( ) and G �� ( ) of CNF hydrogels.This demonstrates that the SB model can capture the behavior of colloidal systems containing high aspect ratio particles, as CNFs are.
Knowing that the hydrogel concentration directly affects the rheological properties, we analyzed how the parameters changed as a function of .Table 1 summarizes the fitted values for the SB model parameters.We observed two characteristic behaviors in Fig. 3c.On the one hand, scaled as a function of following a power-law behavior.Alternatively, decreased exponentially as a function of , then reached a plateau.We should mention that this work only evaluates the effect of on the rheological properties of CNF hydrogels.However, there are more experimental conditions that may influence the rheology.For example, the pH of the liquid phase, the -potential, and particle size distribution are other parameters that can be controlled and change the rheological behavior of a colloidal system.The effect of these different experimental conditions on the rheology could be interconnected with the SB model parameters by using experimental design approaches that allow predictions in high-dimensional spaces (response surface methodology or Bayesian optimization) (Miranda-Valdez et al. 2022).
As the parameter is closely related to the hydrogel stiffness, the observed power-law increment in means that the CNF hydrogels became stiffer when increased.Such power-law behavior of is analogous The parameters were used for establishing the scaling law of the cellulose nanofiber hydrogels.The MAE column shows the mean absolute error of the SB model.*Scaled values from the results in Fig. 3c.For the hydrogel with = 0.90 wt%, the fitted parameters using the Scott-Blair model were = 407 Pa s and = 0.056
The drawback of the latter approach is its bias raised from neglecting the effect of frequency on the material functions G � ( ) and G �� ( ) .A more systematic solution can be obtained from determining the scaling law of the quasi-modulus ∝ n .From fitting the data in Fig. 3c to a power-law function, we found that ≈ 578 2.45 ( R 2 = 0.999).
Regarding the behavior of as a function of , from Fig. 3c, we observed that the CNF hydrogels took a viscoelastic solid conformation since readily approached a minimum value close = 0.053 , a number in the order of magnitude reported for other polymer gels with similar G � ( ) and G �� ( ) (Kes- havarz et al. 2017).To characterize ior of as a function of , we fitted the results in Fig. 3c using an exponential function.This is, ≈ 0.053 + 0.322 ⋅ 0.0045 .However, one should remark that the function for describing as a function of depends on case-by-case.Here we chose the function that best fitted the -relationship.Moreover, and were determined only as a function of , but they are sensitive to varying conditions (e.g., pH) that may lead the colloidal cluster to form a different network.The next section shows the relevance of the defined functions of and as a function of .Given a value, one can predict the rheological response of the CNF hydrogel in frequency, stress relaxation, and creep.Further, as shown in Fig. S4, predicting flow curves is possible for a given .
To show how the scaling law ∝ n can be used, we estimated the fractal dimension d f of the CNF hydrogels following the approach reported by Shih et al. (1990).The authors reported that d f of a strong-gel network, as the CNF hydrogels studied here, can be extracted from the storage modulus scaling law, G � ∝ n (Shih et al. 1990).We have demonstrated that G � ∝ n and ∝ n are analogous; furthermore, they share the same exponent n.Accordingly, one can estimate d f from n = (d + x)∕(d − d f ) .There, d = 3 represents a three-dimensional network, and x is the backbone fractal dimension of flocs for strong-linked CNF hydrogels (Fneich et al. 2019), usually ranging from 1 to 1.4.Recalling that n = 2.45 , the fractal dimen- sion can be approximated to a constant value ranging from d f = 1.37 to d f = 1.20 , depending on the value of x.This value is slightly lower than the one reported in the literature for CNF hydrogels (Fneich et al. 2019), implying weaker aggregation.In the weak link regime, where n = (d − 2)∕(d − d f ) , the Shih model would give a much higher fractal exponent ( d f = 2.59 ), but one can explore intermediate link strengths by employing the Wu-Morbidelli model (Wu and Morbidelli 2001).This model is given by n = ∕(d − d f ) where the exponent = d − 2 + (2 + x)(1 − z) includes a parameter z which interpolates between the strong ( z = 0 ) and weak link ( z = 1 ) regimes.Then, with a quite small deviation from the strong link regime ( z = 0.2 ) one can achieve d f ≈ 1.5 , which was previously found for CNF hydrogels (Fneich et al. 2019) in the dilute regime.As this is only an estimation, we suggest that future research attempts to correlate the fractional parameters with the structure of CNF hydrogels.The latter is out of the scope of this article, but it could be achieved by combining rheological characterization with confocal microscopy.
The SB model can be analytically solved for stress relaxation and creep experiments.Similar to the analytical equations of G � ( ) and G �� ( ) , the model could describe the material functions of the relaxation modulus [G(t)] and creep compliance [J(t)] using the same two parameters: and .Stress relaxation and creep tests are critical to assessing the durability and shape stability of CNF hydrogels.For example, creep behavior is desired when a material needs to flow through a pipe.In particular, for CNF hydrogels, creep is favorable for injectability or pipe flow.Therefore, characterizing G(t) and creep compliance J(t) allows a thorough comprehension of the rheological behavior.
For all the concentrations, we replaced the previously obtained values of and in Eq. 8 and Eq. 9.Then, we compared the computed responses of G(t) and J(t) with the experimental results obtained for each CNF hydrogel.Figure 4a,b demonstrate that the SB computation approximates the experimental results, confirming that the SB model could reduce 1 3 Vol:. ( 1234567890) the experimental burden during the rheological characterization of CNF hydrogels.One must highlight that the SB model can only describe the powerlaw regime during stress relaxation and creep.The latter was less noticeable for stress relaxation since the rheometer cannot record the G(t) response at short timescales.On the other hand, during creep, one can see that only the second stage of creep can be accurately described.Since the CNF hydrogels were highly elastic, it took longer to induce creep than relaxation.In general, with the SB model, as shown in Fig. 4c, and its analytical equations, one can obtain a holistic framework that characterizes the rheology of CNF hydrogels in the reported concentration regime and could be expanded to characterize a broader range of concentrations.

Scaling law of quasi-modulus and degree of viscoelasticity: prediction of material functions
In practice, the relevance of rheological models relies on their potential to predict rheological responses given a material composition (Bonfanti et al. 2020).In that context, for CNF hydrogels, this translates to tailoring a hydrogel with the desired rheological properties based on the CNF concentration.Hence, to validate the relevance of the SB model, we used the functions defined for and (i.e., the dashed lines in Fig. 3c) to predict the magnitude of both parameters given = 0.90 wt%.Then, we replaced the predicted values of and at = 0.90 wt% in Eqs. 5 and 6. Figure 5a shows that the predicted material functions, G � ( ) and G �� ( ) , corresponded to the experimen- tal results.Additionally, knowing that the prediction in Fig. 5a was in agreement with the experimental results, one could further predict the material functions G(t) and J(t) by replacing and in Eqs. 8 and 9.The predicted responses, G(t) and J(t), for a hydrogel with = 0.90 wt% are shown in Fig. 5b,c.The experimental results agreed again with the predicted responses.In addition, the constitutive equation of the springpot can be extended to describe other experiments.For example, the springpot can be analytically solved to predict the engineering stress-strain curve of the CNF hydrogels (Meng et al. 2019;Meng 2021).Furthermore, when the Cox-Merz rule is valid, the SB parameters can be used to predict the nonlinear response * ( ) , i.e., complex viscosity (Jais- hankar and Mckinley 2014).Figure S4 compares the experimental data of * ( ) with the predicted values from scaling the SB parameters to = 0.90 wt%.The results in Fig. S4 confirmed that it is possible to compute non-linear responses from the linear behavior using the SB model.Also, it demonstrated that the SB model could describe the dynamic viscosity of materials with yield point such as CNF hydrogels.

Conclusions
Fractional rheology successfully described the viscoelastic phenomena in colloidal hydrogels containing TEMPO-oxidized cellulose nanofibers, proving our hypothesis that the Scott-Blair (SB) model could characterize the power-law behavior depicted by the material functions.The model accurately describes the stress relaxation, creep, and even complex viscosity response from the information obtained from the frequency spectrum of CNF hydrogels.The Scott-Blair model uses only two parameters: the quasimodulus ( ) and the degree of viscoelasticity ( ).Assessing the effect of particle concentration ( ) on the magnitude of and showed that the quasimodulus followed the same scaling law that has been suggested empirically for the storage modulus ( G ′ ) of CNF hydrogels.Regarding the magnitude of , the degree of viscoelasticity dropped exponentially as a function of .The latter was evidence of the highly elastic character of the hydrogels.The SB model provided constitutive equations for material functions previously described empirically for CNF hydrogels.
The results demonstrated that fractional rheology can be used to describe materials with complex viscoelastic responses.Accordingly, our findings suggest that the SB model could be used as a standard method to characterize CNF hydrogels in the studied regime.For example, the and parameters characterized for a CNF hydrogel could be compared with those characterized for another material.Offering a comparison framework is crucial for easing material selection and design.In particular, since CNF hydrogels can be produced with high particle concentrations (e.g., = 7 wt%), such colloidal systems could be easily diluted to achieve a rheological behavior suiting the target application.We showcased an example of how after scaling the and as a function of , the material functions [ G � ( ) , G �� ( ) , G(t), J(t) and * ( ) ] can be predicted for a given .We found that the predictions matched the experimental results.The results highlight the SB model robustness since all the experimental observations were consistently accurate with the predictions.
Vol:. ( 1234567890) original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Fig. 2 a
Fig. 2 a Scheme of the principle of intermediacy depicted by a spring, springpot, and dashpot (from top to bottom order).The springpot represents the stress ( ) response of viscoelastic materials as an intermediate behavior between Hooke's (spring) and Newton's (dashpot) laws.The constitutive equation of the springpot generalizes both laws with a fractional order derivative.is the fractional order of the derivative (between 0 and 1) of strain ( ) respective to time (t).represents the degree of viscoelasticity, and is a quasi-property with Pa s units.Signature rheological response of b Maxwell, the material functions G � ( ) and G �� ( ) for the Maxwell, Voigt-Kelvin, and Scott-Blair (SB) models.In brief, Maxwell and Voigt-Kelvin are two of the simplest linear viscoelastic traditional models; the first assembles a spring and dashpot in series, and the second in parallel.While the Maxwell model can characterize rheological behaviors in frequency for cross-linking materials, Voigt-Kelvin fails to represent typical rheological responses in viscoelastic materials (Zelenkova et al. 2017; Webber and Tibbitt 2022).
Fig. 3 Rheological responses of cellulose nanofiber (CNF) hydrogels with a concentration ( ) of 0.70 wt% and b of 1.16 wt% in dynamic mechanical analysis (DMA).c Parame-

Fig. 4
Fig. 4 Rheological responses of cellulose nanofiber (CNF) hydrogels with a concentration ( ) of 1.16 wt% and of 0.70 wt% in a stress relaxation and b creep experimental.The dashed lines are the results of replacing the Scott-Blair (SB)