i-Rheo: determining the linear viscoelastic moduli of colloidal dispersions from step-stress measurements

We report on the application of a Fourier transform-based method, ‘i-Rheo’, to evaluate the linear viscoelastic moduli of hard-sphere colloidal dispersions, both in the fluid and glass states, from a direct analysis of raw step-stress (creep) experimental data. We corroborate the eﬃcacy of i-Rheo by comparing the outputs of creep tests performed on homogenous complex fluids to conventional dynamic frequency sweeps. A similar approach is adopted for a number of colloidal suspensions over a broad range of volume fractions. For these systems, we test the limits of the method by varying the applied stress across the materials’ linear and non-linear viscoelastic regimes, and we show that the best results are achieved for stress values close to the upper limit of the materials’ linear viscoelastic regime, where the signal-to-noise ratio is at its highest and the non-linear phenomena have not appeared yet. We record that, the range of accessible frequencies is controlled at the higher end by the relative weight between the inertia of the instrument and the elasticity of the complex material under investigation; whereas, the lowest accessible frequency is dictated by the extent of the materials’ linear viscoelastic regime. Nonetheless, despite these constrains, we confirm the eﬀectiveness of i-Rheo for gaining valuable information on the materials’ linear viscoelastic properties even from ‘creep ringing’ data, confirming its potency and general validity as an accurate method for determining the material’s rheological behaviour for a variety of complex systems.


Introduction
In rheology, the determination of the linear viscoelastic (LVE) properties of materials is of great importance for the determination of their states; i.e., whether a material could be classified as a viscoelastic solid or fluid.This is because, over the years, the LVE properties of materials have been successfully correlated both theoretically and experimentally to their microscopic structure and dynamics occurring at different time and length scales. 1 These are fully described by the frequency-dependent viscoelastic moduli, G 0 (o) (storage modulus) and G 00 (o) (loss modulus), which are classically measured over a finite set of frequencies via oscillatory measurements called dynamic frequency sweeps (DFSs). 2 However, due to the mechanical inertia carried by the electric motor of rheometers, the frequency of oscillation cannot exceed typically 100 Hz, whereas the lowest accessible frequency is limited by the measurement duration and the stability of the sample.Typically, frequencies lower than 0.1 Hz require exceedingly extensive measuring times.In a recent work 3 Tassieri et al. proposed a novel rheological 'tool' to derive the frequency-dependent viscoelastic moduli via a discrete Fourier transform 4 of step-strain and relaxation experiments.This tool has been made available to the scientific community in the form of an open access executable called 'i-Rheo' and it has been recently adapted to analyse atomic force microscopy measurements 5,6 and molecular dynamic simulations. 7An advantage of using i-Rheo is that there is no need of assuming any rheological model for data interpretation nor idealisation of the real measurements.Moreover, thanks to the fast sampling of modern rheometers, this tool is particularly useful for accessing the materials' high frequency responses, especially for strain-controlled rheometers, for which the motor inertia is less disruptive than in the case of stress-controlled rheometers; but, would still limit anyway the range of accessible frequencies to 100 Hz for oscillatory measurements.In this work, we exploit the efficacy of i-Rheo for the analysis of step-stress (creep) measurements that suit better the stress-controlled rheometers.In particular, we first validate the use of i-Rheo for the analysis of creep measurements performed on homogeneous complex fluids and then we exploit it for the investigation of the LVE properties of hardsphere colloidal suspensions, over a broad range of colloidal volume fractions f including non-equilibrium glass states.Despite their complex rheological behaviour, colloidal suspensions are a class of soft-materials widely used in daily applications including paints, food products, cosmetics and cleaning products. 1,8Moreover, due to their inherent complexity, model suspensions, such as one-component hard-sphere colloids (for which the volume fraction f is the only control parameter), have been extensively studied for a better understanding of the rheological behaviour of more complex suspensions. 9,10It is known that for such systems the equilibrium state diagram shows fluid states for f o 0.494, fluid-crystal coexistence for 0.494 o f o 0.545 and a crystalline solid for f 4 0.545. 9,10oreover, when the width of the size distribution of the colloidal particles exceeds 7% circa, the increase of f above f g E 0.58 leads to the formation of a non-equilibrium amorphous solid (glass) 11,12 instead of a crystal. 13Glass formation is attributed to the trapping of a particle by the nearest neighbors (the so-called cage formation) 9 and is detected due to the appearance of a solid-like response of the suspension. 14,15][16] The experimental determination of these moduli by means of conventional DFS measurements is time consuming especially at low frequencies, where the relaxation processes of these samples eventually occur.2][23][24][25][26][27][28] The transition to flow can be affected by shear banding, [29][30][31][32] ageing and history effects. 33Here, we demonstrate that the evaluation of the materials' LVE moduli by means of i-Rheo applied to creep measurements is a reliable and faster alternative to conventional DFS measurements, also in the case of colloidal suspensions.The linear viscoelastic moduli obtained via i-Rheo are in good agreement with conventional DFS measurements of hard-sphere dispersions in both fluid and glass states.We discuss the results obtained for a range of different applied stresses across the materials' linear and non-linear response regimes, identifying the best choice of the applied stress.We also discuss the limits of the accessible range of frequencies for such systems, in particular in the non-equilibrium glass state, where non-homogeneous flows are observed for the majority of the experimental conditions.

Linear rheology in the frequency-domain: oscillatory measurements
The LVE properties of a material are fully described by its frequency-dependent complex shear modulus G*(o), which is a complex number whose real and imaginary parts provide information on the elastic and the viscous nature of the material, respectively.The complex shear modulus is defined as the ratio between the Fourier transforms of the stress s(t) and the strain g(t), whatever their temporal forms and regardless of whichever has been imposed or measured: 34 where o is the angular frequency, i is the imaginary unit (i.e., i 2 = À1), and G 0 (o) and G 00 (o) are the material storage (elastic) and loss (viscous) moduli, respectively.The standard method of measuring G*(o) over a finite set of frequencies is based on the imposition of an oscillatory stress s(o,t) = s 0 sin(ot) (where s 0 is the amplitude of the stress function) and the measurement of the resulting oscillatory strain, which would have a form of g(o,t) = g 0 (o)sin(ot + j(o)), where g 0 (o) is the strain amplitude and j(o) is the phase-shift between the stress and the strain.From eqn (1) it follows that: from which the expressions of the viscoelastic moduli can be inferred.Nonetheless, conventional rotational rheometers are limited in the range of accessible frequencies, mainly at the top end, because of the inertia carried by the electric motor that generates the torque (proportional to s(t)) necessary to induce the material deformation.At low frequencies the motor inertia is negligible and the measurements are mainly limited by the 'patience' of the operator (when of course material related issues such as ageing, evaporation and mutational effects 35 have been already pondered).For instance, to perform an oscillatory measurement at o = 10 À3 rad s À1 (and considering 'only' 1.5 of a cycle) a rheometer would take (1.5 Â 2p)/10 À3 C 10 4 s (i.e., circa three hours) to read a single value of G*(o).

Fourier series of a square wave
In this section, we wish to remind the reader of an important mathematical equivalence that underpins the effectiveness of creep-measurements for the evaluation of the materials' LVE properties over a wide frequency spectrum.This is related to the equipollence of a square wave to an infinite sum of sinusoidal functions having different frequencies and weighted amplitudes; i.e., a Fourier series.In particular, in the case of real experiments, where a step-measurement cannot last for infinite time (for obvious reasons), the temporal behaviour of the step-function assumes the form of a half square wave u(t) with a relatively long period T (as drawn in Fig. 1).It follows that u(t) could be expressed as: where o = 2p/T is the fundamental angular frequency of the square wave and k identifies the order of the higher harmonics.
The first four terms of this series are shown in Fig. 1.Therefore, in real experiments it follows that, by applying a stress having a temporal form of a half square wave u(t) with a sufficiently long period T, it is like testing the sample with an infinite number of sinusoidal stress functions with angular frequencies ranging from o min = 2p/T to o max = 2p/t sam , where t sam is the acquisition time of the instrument.

Linear rheology in the time-domain; step-stress measurements
In the case of creep measurements in simple shear, where a constant shear stress of magnitude s 0 (i) is applied to the sample within a brief period of time e (Fig. 2), (ii) is achieved with a constant rate of stress _ s ¼ s 0 =e and (iii) holds for an unlimited period of time, the following constitutive equation simplifies: where J(t) is the material's shear creep compliance and t 0 is the time at which the stress reaches the constant value of s 0 .By appealing to the mean value theorem, eqn (4) solves as: which is indistinguishable from when t 0 is set equal to 0 and for t c e. Eqn (6) has been accepted by the rheology community because the data acquisition rates of rheometers were (until few years ago) slower than e.However, nowadays, thanks to the development of sensitive and fast transducers, it is possible to acquire a significant number of experimental data within the time window defined by e.Nevertheless, these data are commonly discarded because of the condition on which eqn (6) relies on, i.e. t c e.

Evaluating the Fourier transform of raw data
An analytical procedure for the evaluation of the Fourier transform of any generic function sampled over a finite time window was proposed by Evans et al., 36,37 to convert creep compliance J(t) into G*(o) directly, without the use of Laplace transforms nor fitting functions.This method is based on the interpolation of the finite data set by means of a piecewise-linear function.In particular, the general validity of the proposed procedure makes it equally applicable to find the Fourier transform g ˆ(o) of any time-dependent function g(t) that vanishes for negative t, sampled at a finite set of data points (t k ,g k ), where k = 1. ..N, which extend over a finite range, and need not be equally spaced: 37 where : g N is the gradient of g(t) extrapolated to infinite time and g(0) is the value of g(t) extrapolated to t = 0 from above.
This method has been improved by Tassieri et al., 38 while analysing microrheology measurements performed with optical tweezers.The authors found that a substantial reduction in the size of the high-frequency artefacts, from which some highfrequency noise tends to spill over into the top of the experimental frequency range, can be achieved by an over-sampling technique.The technique involves first numerically interpolating between data points using a standard non-overshooting cubic spline, and then generating a new, over-sampled data set, by sampling the interpolating function not only at the exact data points but also at a number of equally-spaced points in between.Notice that, over-sampling is a common procedure in  signal processing and it consists of sampling a signal with a sampling frequency f s much higher than the Nyquist rate 2B, where B is the highest frequency contained in the original signal.A signal is said to be oversampled by a factor of b f s /(2B). 39he advantages of an analytical tool, such as i-Rheo, which allows the evaluation of the Fourier transform of raw experimental data, become apparent when both eqn (4) and the relationship between the material's complex shear modulus and the creep compliance are considered together in the frequency domain: Therefore, the materials' LVE properties can be deduced directly from the ratio of the Fourier transform of the experimental strain and stress functions, without any assumption on the time scales (i.e.t c e) nor on the functional form of s(t) within e; as long as the experiments are performed within the materials' LVE regime.
3 Materials and methods

Samples
3.1.1Homogeneous complex fluids.In order to validate the efficacy of i-Rheo for analysing creep measurements, we prepared two homogeneous complex fluids with significantly different LVE properties: (I) a polyacrylamide (PAM, from Polysciences Inc.) with a nominal molecular weight of M w = 18 MDa was used to prepare an aqueous solution at a mass concentration of 1% w/w.The solutions were mixed with a magnetic stirrer at 200 rpm for 48 h at room temperature.(II) A standard calibration fluid was provided with the purchase of the rotational rheometer MCR 302 (Anton Paar).This fluid is a polydimethylsiloxane (PDMS) sample that is expected to show a low frequency crossover of the moduli at o = 100 rad s À1 and T = 25 1C.
3.1.2Colloidal suspensions.We prepared dispersions of polymethylmethacrylate (PMMA) spheres of radius R = 150 nm (polydispersity 12%), sterically stabilized with a layer of polyhydroxystearic acid (PHSA), in cis-trans-decalin. 40In this solvent the particles behave like hard spheres. 41Dispersions with different colloidal volume fractions f were prepared by diluting a sediment obtained by centrifugation of a dilute suspension.][44] We prepared dispersions with f = 0.40, 0.48, 0.50, 0.52, 0.54, 0.56, 0.58, and 0.60, going from the fluid to the glass state.After dilution, the samples were mixed in a rotating wheel for a least 12 h before the experiments were started.

Rheology
Rheological measurements on homogeneous complex fluids were performed with a stress-controlled rotational rheometer MCR 302 (Anton Paar) having the lowest torque value of B0.1 mN m.The instrument was equipped with a cone and plate geometry of D = 50 mm in diameter and 11 angle.For the PDMS sample, a stress of s 0 = 500 Pa was applied to the sample for a duration of 10 3 s.For the 1% weight polyacrilamide solution a stress of s 0 = 0.5 Pa was applied to the sample for a duration of 10 4 s.Conventional dynamic frequency sweep measurements were performed with strain amplitudes of g = 5% and g = 20%, respectively, in the frequency range of 10 À2 o o o 10 2 rad s À1 .All the tests performed on these two samples were carried out at a temperature of T = 5 1C.This temperature was chosen to shift the crossing point of the viscoelastic moduli well within the explored range of frequencies.
Measurements on the colloidal suspensions were performed with a DHR3 stress-controlled rotational rheometer (TA instruments), using a plate-plate geometry of diameter D = 30 mm and a gap d = 0.5 mm.We used a solvent trap to minimize solvent evaporation.The temperature was set to T = 20 1C and controlled within 0.1 1C via a Peltier system.For all samples, a dynamic strain sweep (DSS) at frequency o = 1 rad s À1 was performed directly after loading to determine the linear response regime and, for glassy samples, the strain amplitude at which fluidization occurs.For these last samples, corresponding to dispersions with f = 0.58 and 0.60, the effects of sample loading and ageing were reduced by performing the following rejuvenation procedure before each test.First, a dynamic time sweep with a large strain amplitude, typically g = 500%, well above the yield strain, with a duration of 200 s and a frequency o = 1 rad s À1 was performed in order to fluidize the sample.After that, another dynamic time sweep with the same frequency but strain amplitude within the linear viscoelastic regime, typically g = 1% was applied until a steady-state response, i.e. a timeindependent storage G 0 (o) and loss modulus G 00 (o), was observed, which typically took about 150 s.This indicated that no further structural changes occurred and hence a reproducible state of the sample was obtained.The linear viscoelastic moduli where measured using a conventional DFS with a strain amplitude of g = 0.1-1%, depending on the sample, in the frequency range 10 À1 o o o 10 2 rad s À1 .Creep experiments for different applied stress s values were measured for all the samples for a duration of 10 3 s using the fast sampling option, which provides a sampling time t sam E 1.5 Â 10 À3 s.

Results and discussion
Fig. 3 shows the front panel of the executable i-Rheo, which reads the experimental raw data (i.e., [t k ,s k ,g k ]) in the form of a tab-separated .txtfile and generates a new oversampled set of data (with a sufficiently high value of b f s t sam /2 c ot sam , usually f s D 10 MHz); then it applies eqn (7) to this new data set and returns the viscoelastic moduli of the material under investigation.
It is important to highlight that the values of the two parameters g(0) and : g N required by eqn (7) have been set equal to zero for g(0) = 0, and _ s 1 ¼ 0; whereas, s(0) and _ g N are obtained from the experimental data as follows: (i) s(0) is the first experimental value of the stress (e.g., in Fig. 3 s(0) = 426 Pa); (ii) _ g N is obtained from a linear regression of the long-time behaviour of g(t).Notably, the assumptions made for g(0) and _ s 1 are actually exact, because the sample is initially undeformed and when the stress is constant its time derivative is null.In the cases of s(0) and _ g N , they are sensible approximations.Indeed, for viscoelastic fluids, whose linear response is characterised by the existence of the terminal region at low frequencies (where G 0 (o) p o 2 and G 00 (o) p o for o -0), the long-time behaviour of the strain is expected to be: g(t) p t/Z for t -N and therefore _ g N p 1/Z for t -N, where Z is the fluid's steady state shear viscosity.A similar conclusion is achieved in the case of viscoelastic solids (e.g., gels and rubbers), for which J(t) tends to a finite value J 0 at long times and for which _ g N D 0. The effectiveness of i-Rheo to evaluate the LVE properties of complex materials from conventional bulk rheology step-stress measurements has been initially corroborated by investigating the rheological properties of two homogeneous complex fluids (Section 4.1), then applied to the rheological study of concentrated suspensions of colloidal particles, including glasses (Section 4.2).

Creep experiments of homogeneous complex fluids
In order to corroborate the applicability of i-Rheo to creep measurements, preliminary tests were performed on two homogeneous complex fluids; i.e., a polydimethylsiloxane (PDMS) sample and a solution of polyacrylamide at a concentration of 1 wt%.The results are reported in Fig. 4 and 5, respectively.In particular, in the main plot of both the figures, the material strain g is reported as a function of time t.The inset shows the temporal behaviour of the ratio t/g(t), which in the case of viscoelastic fluids should approach a plateau value over long time frames.This is because at sufficiently longer time the fluids' compliance should assume the temporal form of: J(t) D t/Z; as the elastic character of the complex fluid should vanish towards zero (i.e., G 0 (1/t) p (1/t) 2 -0 for t -N).Therefore, it is a valid approximation to extract the value of _ g N from the terminal (plateau) region of the measurements.Once all the parameters required for the application of eqn (7) have been assigned, then the raw experimental data of both stress and strain functions can be fed to i-Rheo for their conversion into the materials' LVE properties via eqn (8).The efficacy of i-Rheo is validated by the results reported in both Fig. 4 and 5 (right), where the transformed data are directly compared with conventional oscillatory measurements performed on the same samples prior to performing the creep test.In particular, in the case of PDMS (Fig. 4, right), the two experimental methods agree over a range of frequencies spanning four decades, with a net gain of half decade at frequencies higher than 100 rad s À1 in favour of i-Rheo.In the case of the polyacrylamide solution (Fig. 5, right), it is interesting to highlight the effects of instrument inertia at relatively short time scales (i.e., for t o 1 s in this case), where one could identify a characteristic time t I (dotted lines in Fig. 5, left) before which the instrument is not able to account for the inertia of the electric motor and the strain rapidly increases from zero up to a region where the 'creep ringing' 45,46 phenomenon is observed (which is caused by the coupling of the instrument inertia with the elasticity of the sample).As we show also later for colloidal suspensions, the creep ringing data are retained and actually contain useful information on the fluid's LVE properties.Instead the data affected by the motor inertia at time scales shorter than t I (i.e., for o \ 10 rad s À1 ) cannot be used to determine the viscoelastic moduli.This also implies that, This journal is © the Owner Societies 2020 even if the sampling rate t sam is shorter than t I , 1/t sam will not define the largest accessible frequency, which is not the case for step-strain experiments. 3If these initial data affected by inertia are discarded, the two methods show a noteworthy agreement over three decades of frequencies, with a net gain of two decades at low frequencies in favour of i-Rheo, which would be extremely time consuming to achieve with oscillatory measurements.The initial range of creep data to be discarded can be identified by those values of t/g(t) showing a monotonic decrease at short time scales (i.e., all those data on the left side of the dotted line drawn in the inset of Fig. 5, left).Notice that, as a means of comparison, in Fig. 5 (right) we also report the viscoelastic moduli derived via i-Rheo applied to the entire set of data, including those affected by the instrument inertia; the results endorse the arguments discussed above.
Moreover, with regards to the lowest accessible frequencies, typically the limiting factor is the signal-to-noise (SNR) ratio of the data measured after a long time.In a creep experiment though, large strain g values are measured at long times, meaning that the SNR value should be large.This is different from the use of i-Rheo to determine viscoelastic moduli from stress relaxation after a step strain.Indeed, in that case the stress tends to zero at long times and therefore SNR is small.Thus creep tests are in principle a better choice to determine viscoelastic moduli at small frequencies.However, as shown hereafter, the sample behaviour at long times might impose additional constraints.

Creep experiments on colloidal suspensions
After corroborating the efficacy of i-Rheo for determining the LVE moduli of homogeneous complex fluids, we apply the same approach to concentrated suspensions of hard-sphere like colloids with volume fraction 0.40 r f r 0.60.In these concentrated colloidal suspensions additional complexity is present due to hydrodynamic effects and in particular nonequilibrium effects associated with the glass transition.
Hereafter we report the results for a selected sample with f = 0.54 to illustrate the procedure used to determine the linear viscoelastic moduli by means of i-Rheo and investigate the dependence on the applied stress.Then, we report the linear viscoelastic moduli obtained via i-Rheo for all the investigated samples and their comparison with those of conventional DFS.
The results of a dynamic strain sweep (DSS) test performed on a sample with f = 0.54 at o = 1 rad s À1 are reported in Fig. 6.At this frequency, the sample behaves like a viscoelastic solid, with the storage modulus larger than the loss modulus within the material's linear response regime, which extends up to a strain amplitude of g 0 E 1%, with a corresponding oscillatory stress amplitude of circa 0.3 Pa.In the non-linear regime of the material response the moduli cross at g 0 E 12%, which corresponds to a stress amplitude of approximately 1.5 Pa, after which a fluid-like response is observed.This is consistent with recent studies, which indicate the crossing point of the frequency dependent viscoelastic moduli as the yielding point of soft glassy materials. 47We have used the DSS tests to identify for each sample the stress values that fall within the linear and non-linear regimes of the materials.We then selected a series of stresses across the two regimes to investigate the effects of the applied stresses on the frequency-dependent moduli evaluated by means of i-Rheo.Fig. 7 shows the creep measurements performed at different applied stresses (i.e., s = 0.1, 0.2, 0.3, 0.4, 0.5 and 1 Pa) on a sample having a volume fraction of f = 0.54.
As for the polyacrylamide solution, also the response of this system shows inertia effects in the time interval 10 À3 o t o 10 À1 s, followed by a 'creep ringing' regime for 10 À1 o t o 1 s.At longer times, all curves show a sub-linear increase of the strain g(t) with time, which is characteristic of solid-like materials for applied stresses lower than their yield stress s y , which for this system can be estimated from Fig. 6 to be s y D 1.5 Pa.Moreover, from Fig. 7 it can be seen that the measured strain increases with the increase of the applied stress, as expected.

Results of i-Rheo for different applied stresses
In Fig. 8 we report the comparison between the frequencydependent linear viscoelastic moduli obtained by means of i-Rheo applied to the creep measurements shown in Fig. 7 (lines) and those obtained from a DFS (symbols).In particular, it can be seen that for the smallest applied stress of s = 0.1 Pa the moduli obtained from i-Rheo overestimate the solid-like behaviour of the material, with a higher G 0 (o) and a lower G 00 (o) than those measured with a DFS, whereas, for s = 0.2 Pa, while the storage modulus shows a good agreement with that obtained from the DFS measurement, G 00 (o) still shows a significant deviation.Interestingly, for s = 0.3 Pa the agreement between the two methods is very good for both moduli, especially in the frequency range of 0.1 o o o 10 rad s À1 .Notably, this value of the applied stress is very close to the limit of the material's linear response, see Fig. 6.For larger values of the applied stress we observed a detrimental effect on the agreement between the two methods.This is expected as i-Rheo assumes a linear relation between stress and strain, which of course does not hold in the materials' non-linear regime.Therefore, the onset of deviations between the moduli obtained by i-Rheo and those obtained by conventional oscillatory measurements could be used to estimate the onset of the materials' non-linear response.In particular, the moduli derived by means of i-Rheo increasingly deviate from the expected values at small frequencies (i.e.long relaxation times), with a crossover of the moduli that shifts towards increasingly larger frequency values.This is consistent with the progressive shortening of the structural relaxation time of the sample associated with the structural rearrangements induced by the shear deformation, which leads to yielding and fluidisation of the sample when the applied stress falls in the materials' non-linear response regime.It is interesting to notice that, of the two moduli G 0 (o) is less sensitive than G 00 (o) to the variations of the applied stress.A similar behaviour has been observed for all samples investigated in this work.Specifically, we found that the best agreement between the two approaches was obtained for stress values that were closely approaching from the bottom the upper limit of the materials' linear response regime.This might be due to the fact that noise contributions are smaller for larger stress, i.e. for stresses close to the limit of the materials' linear response regime the signal to noise ratio is close to its maximum value, but yet within the materials' LVE regime, where the effects of shear-induced structural changes are not observed.

Comparison between i-Rheo and DFS for different / values
In Fig. 9 we report the comparison between the frequencydependent linear viscoelastic moduli obtained by means of i-Rheo (lines) by converting creep measurements (see Fig. 10) and those obtained from DFS (symbols) for different colloidal volume fractions, ranging from the fluid state at f = 0.40 to the glass state at f = 0.60.Notably, for all the investigated samples the agreement between the two experimental procedures is very good, confirming the effectiveness of i-Rheo.It must be noted that, while all the creep measurements were performed over the same experimental time window, the related viscoelastic moduli were computed over different frequency ranges, depending on the values of f, because of the instrument inertia 45,46 as already discussed earlier.In particular, for f = 0.40 the instrument inertia was affecting the measurement up to a time of t 1 E 1.5 s (Fig. 10), which implies a maximum frequency of o E 5 Â 10 À1 rad s À1 .For f = 0.48, 0.50 and 0.52 the influence of the inertia is less detrimental than the previous case (Fig. 10), and data were truncated at a time of t 2 E 1.5 Â 10 À1 s, which translates into a maximum accessible frequency of the order of 7 rad s À1 .For higher volume fractions (or higher material stiffness), i.e. f = 0.54, 0.56 and 0.58 the creep data were discarded (because of inertia) up to a time of   t 3 E 5 Â 10 À2 s (Fig. 10) and therefore achieving frequencies up to o E 20 rad s À1 .Finally, for the largest volume fraction f = 0.60 (with the highest stiffness) the maximum frequency was o E 10 2 rad s À1 .As a preliminary conclusion, we could infer that by increasing the volume fraction of the dispersion we increase the overall stiffness of the material, which mitigates the disruptive effects of instrument inertia on the measurements of the materials' LVE properties via creep tests, as corroborated by the results reported in Fig. 4 and 10.9][50] It is interesting to mention that the regions of the creep measurements affected by the ringing effects were not discarded during the analysis with i-Rheo; they actually provide output data that are in good agreement with oscillatory measurements, as shown in Fig. 5 and 10, confirming that the ringing regime contains useful information on the LVE properties of materials.
In terms of the lowest accessible frequency, we can argue that it is affected by the physics of the suspension rather than by the technical limitations of the rheometer.32]51 This was not a problem for the lowest volume fractions explored in this work, i.e. f = 0.40, 0.48 and 0.50, which did not manifest flow inhomogeneities within the measurement time-window (Fig. 10) and therefore the viscoelastic moduli were evaluated down to o E 4 Â 10 À3 rad s À1 .In the cases of f = 0.40 and 0.48 some oscillations observed in G 00 (o) could be attributed to a low signal and noise of the measurements.Moreover, for all f Z 0.52 the samples showed an anomalous behaviour of the strain g(t) for t \ 10 2 s (Fig. 10), which was limiting the range of the lowest accessible frequencies by means of i-Rheo.

Conclusions
In this work we report experimental evidence of the applicability of i-Rheo to the analysis of time-dependent step-stress (creep) measurements for the evaluation of the linear viscoelastic properties of materials.I-Rheo was recently introduced for the analysis of time-dependent step-strain measurements and successfully exploited for the analysis of both molecular dynamics simulations and atomic force microscopy measurements.
Here we have validated the applicability of i-Rheo to analyse creep measurements performed on homogeneous complex fluids, and then exploited its effectiveness for the evaluation of the viscoelastic moduli of suspensions of hard sphere-like colloids at different volume fractions, spanning the state diagram from fluids to glasses.For all the explored samples, we found a very good agreement between the outputs obtained by means of i-Rheo and those obtained via conventional oscillatory measurements.We found that the transformation of creep measurements by means of i-Rheo works best for values of the applied stress that are close to the limit of the suspension's linear response regime.By entering the materials' non-linear response regime, we found an increasingly larger deviation from oscillatory measurements, which can be associated with the progressive structural rearrangement of the samples induced by the shear flow.The range of accessible frequencies for the systems investigated in this work was limited by two factors: at high frequencies, by the inertia of the instrument, which becomes less important with increasing the sample stiffness (or increasing f); and at low frequencies, by physical phenomena occurring due to the nature of the samples, such as shear banding.For values of the volume fraction approaching and crossing the glass transition, shear induced inhomogeneities appearing at long times (revealed by an irregular behaviour of the time-dependent strain g(t)) limited the lowest accessible frequency.Nonetheless, we confirm the effectiveness of i-Rheo for gathering valuable information on the materials' linear viscoelastic properties even from 'creep ringing' data, confirming its potency and general validity as an accurate method for determining the material's rheological behaviour for a variety of complex systems.Indeed, we anticipate the use of i-Rheo for the rheological investigation of novel and more complex colloidal systems, with particular attention to the system-specific limitations for the evaluation of the low frequency linear viscoelastic properties.Moreover, besides obvious differences in the experimental protocols, i-Rheo could be successfully used to determine the viscoelastic properties of harder materials, such as for example metallic glasses. 52

Fig. 1 A
Fig. 1 A schematic representation of half square wave u(t) and its approximation obtained by considering only the first four terms of the Fourier series shown in eqn (3).

Fig. 2
Fig. 2 Schematic representation of a real step-stress measurement.(top) A semi-idealised temporal behaviour of the stress function s(t); (bottom) description of the material's strain function g(t) in response to the applied constant stress.

Fig. 5
Fig.5(left) The temporal behaviour of the measured strain g when a step-stress of s 0 = 0.5 Pa is applied to a solution of polyacrylamide at a concentration of 1 wt%.The inset shows the temporal behaviour of t/g.The dotted lines in both the graphs indicate the minimum time (t I ) after which the data are transformed.(right) Comparison between the solution's viscoelastic moduli derived by means of i-Rheo applied to the step-stress data shown on the left (with and without the initial data affected by the instrument inertia, dashed and continuous lines, respectively) and those obtained via a conventional oscillatory measurement (symbols).Both the measurements were performed at a temperature of T = 5 1C.

Fig. 4 (
Fig.4(left) The temporal behaviour of the measured strain g when a step-stress of s 0 = 500 Pa is applied to a polydimethylsiloxane (PDMS) sample.The inset shows the temporal behaviour of t/g.(Right) Comparison between the PDMS' viscoelastic moduli derived by means of i-Rheo applied to the step-stress data shown on the left (lines) and those obtained via a conventional oscillatory measurement (symbols).Both the measurements were performed at a temperature of T = 5 1C.

Fig. 6
Fig. 6 Left y-axis: storage (G 0 (o), closed symbols) and loss (G 00 (o), open symbols) moduli as a function of strain amplitude g 0 , obtained by dynamic strain sweeps at frequency o = 1 rad s À1 for a sample with f = 0.54.Right y-axis: oscillatory shear s as a function of g 0 .Dashed lines indicate the value of stress applied in creep experiments for which the best agreement between i-Rheo and DFS is obtained.The shaded region indicates the linear response regime of the measurement.

Fig. 7
Fig. 7 Strain g(t) as a function of time t obtained from creep measurements performed with different applied stresses on a sample with f = 0.54.Inset: t/g as a function of t.

Fig. 8
Fig. 8 Linear viscoelastic moduli G 0 (o) (solid lines) and G 00 (o) (dasheddotted lines) obtained from i-Rheo for different applied stresses, as indicated.Results of i-Rheo are compared with results of conventional dynamic frequency sweeps (symbols).Black lines correspond to the selected best result from i-Rheo.The dashed line in the legend separates stress values in the linear and non-linear response regimes.

Fig. 9
Fig. 9 (a and b) Linear viscoelastic moduli G 0 (o) (solid lines) and G 00 (o) (dashed lines) obtained from i-Rheo for samples with different volume fractions f, as a function of frequency o.Moduli obtained from conventional DFS measurements (symbols) are reported for comparison.

Fig. 10
Fig. 10 Strain g(t) as a function of time t obtained in creep experiments for different colloidal volume fractions f, as indicated.The applied stress depends on the sample and corresponds to the value used to obtain the viscoelastic moduli in Fig. 9. Arrows indicate the limit of the response affected by tool inertia effects.