The authors thank Anton Paar for use of the MCR 702 rheometer through their VIP academic research program. We also thank Dr. Abhishek Shetty for the comments in the instrument setup.

1. Rheometer Setup
<ol>
	<li>With the rheometer configured in the SMT mode (see note), attach the upper and lower drive geometries. To maintain as close to a homogeneous shear field as possible, use a 50 mm plate (PP50) as the lower fixture, and a 2-degree cone (CP50-2) for the upper fixture. 
	Note: The rheometer we use (see the Table of Materials) can be configured in either a combined motor-transducer (CMT) or separate motor transducer (SMT) mode. With only a single motor integrated in the rheometer head, it acts as a traditional CMT stress-controlled rheometer and the data obtained require inertia corrections. With two motors incorporated in a SMT mode, the upper motor operates solely as a torque transducer and the bottom motor acts as a drive unit thus converting the rheometer into a typical strain-controlled rheometer.
	<ol>
		<li>Attach the bottom and top geometries.</li>
		<li>Click the zero-gap button in the control panel.</li>
		<li>Navigate to start service function under the Measuring set tab on top. Run the inertia calibrations for the upper and lower measuring systems, found in the dropdown menu.</li>
		<li>Run adjustments for the upper and lower motors.</li>
		<li>Specify the desired temperature in the control panel. 
		Note: The measurements at which experiments on XG and PEO solutions are performed are 25 &#177; 0.1 &#176;C and 35 &#177; 0.1 &#176;C, respectively.</li>
	</ol>
	</li>
	<li>Load the material of interest on top of the bottom geometry with a spatula or pipette, ensuring no air bubbles are entrained in the sample. 
	Note: Approximate volumes of material required to completely fill a geometry are provided in the rheometry software under Setup | Measuring Systems.
	<ol>
		<li>Load 1.14 mL to fill the cone-and-plate geometry. Load higher viscosity samples with a spatula, and less viscous materials with a pipette. 
		Note: A spatula is used to load the polymer solutions.</li>
		<li>Command the measuring system to trim gap and gently trim the excess material at the edge of geometry with a square-ended spatula, ensuring the spatula remains perpendicular to the axis of the rheometer. 
		Note: The quality of material loading will affect the rheological results significantly and any apparent under- or over-filling should be avoided.</li>
		<li>Press the continue button in the rheometry software to move to the measurement gap. 
		Note: A complete loading process is illustrated in Figure 2.</li>
	</ol>
	</li>
</ol>
2. Running Oscillatory Shear Tests
Note: Two ways of running oscillatory shear tests are introduced. The first approach is designed for sinusoidal stresses and strains only and was used to collect the data we report here. The second method allows for arbitrary stress or strain schedules to be set.
<ol>
	<li>Sinusoidal oscillatory shear 
	<ol>
		<li>Navigate to Large amplitude oscillatory shear-LAOS under My apps in the software. Go to the Measurement box and click strain variable.</li>
		<li>Specify the initial (1%) and final values (4,000%) of a strain amplitude sweep. Specify the imposed frequency of 0.316 rad/s. Define the desired total number of strain amplitudes as 16 in the specified amplitude range, which results in the point density of 5 points per decade.</li>
		<li>Check the Get waveform box at the top to collect transient responses.</li>
		<li>Click the start button at the top to start the experiments and the raw data will be displayed in the rheometry software automatically.</li>
	</ol>
	</li>
	<li>Arbitrary Stress or Strain Schedules
	<ol>
		<li>To impose arbitrary-defined deformation, click Waveform sine generator under My apps in the software.</li>
		<li>Define a list of strain values that correspond to the function that is to be applied (not restricted to sinusoidal waveform). Generate the value list in an external program.</li>
		<li>Click edit under the strain value in the measurement box. Copy and paste these numbers into the value list.</li>
		<li>Specify the number of data points, point duration, and interval time to adjust the imposed frequency. For instance, specify the number of data points and the interval time as 512 points and 6.2832 s, respectively, if a cycle of sinusoidal strain is pasted into strain value list with 512 points and the frequency of 1 rad/s is desired. 
		Note: This approach is not recommended for running sinusoidal oscillatory shear due to the limited number of oscillatory cycles, and also due to the fact that automatic corrections which are enabled in an oscillatory test mode on the rheometer are disabled in this mode. Nonetheless, because there are no assumptions of sinusoidal strain built into the SPP framework, one can arbitrarily define imposed strain functions according to the processing conditions or end use the materials may experience, and the SPP framework remains applicable to analyze the rheological response.</li>
		<li>Check the Get waveform box at the top. Then click the start button at the top to start the experiments.</li>
	</ol>
	</li>
</ol>
3. Performing SPP analysis (SPP-LAOS software)
Note: The SPP analysis software is a MATLAB-based freeware package for analyzing rheological data with the SPP framework and is attached as Supplementary Files 1&#8210;621.
<ol>
	<li>Format the data files to be tab-delimited text (.txt) consisting of four columns in the order of {Time (s), Strain (-), Rate (1/s), Stress (Pa)}. 
	Note: Users may need to modify the number of header lines in the function files to be able to process their data. See sample data files (Supplementary Files 7&#8210;9).</li>
	<li>To run SPP-LAOS software, open the m-file named RunSPPplus_v1.m in MATLAB. 
	Note: While RunSPPplus_v1.m is the main script to run the analysis, the package contains other function files that will be called from the main script, including SPPplus_read_v1.m, SPPplus_fourier_v1.m, SPPplus_numerical_v1.m, SPPplus_print_v1.m and SPPplus_figure_v1.m.</li>
	<li>Navigate to the section labeled User-defined variables, and specify the following variables.
	<ol>
		<li>Filename: Specify the name of .txt file that will be used for the SPP analysis. 
		Note: The file must match the above format requirement.</li>
		<li>Run state: Place the vector as [1, 0] to run Fourier analysis mode for regular oscillatory shear response. 
		Note: The software employs two different methods of calculating the instantaneous SPP moduli,&#160;<img alt="Equation 3" src="/files/ftp_upload/58707/58707eq3.jpg" />&#160;and&#160;<img alt="Equation 4" src="/files/ftp_upload/58707/58707eq4.jpg" />, based on Fourier transformation and numerical differentiation. The Fourier transform approach is designed for periodic input, such as oscillatory shear tests. Arbitrary time-dependent tests, which include, but are not limited to sinusoidal protocols, can be analyzed with the numerical differentiation approach.</li>
		<li>Run state: Input the vector as [0, 1] to run numerical-differentiation analysis mode for arbitrary time-dependent tests.</li>
		<li>Omega (Fourier analysis): Specify the angular frequency of oscillation, with units of rad/s.</li>
		<li>M (Fourier analysis): Define number of higher harmonics to be included in the SPP analysis. Adjust this number to include all the higher harmonics above the noise floor. 
		Note: This number must be a positive odd number and varies with amplitude and material. We include up to the 3rd harmonic in the MAOS regime, and up to the 55th harmonic at the largest amplitude investigated.</li>
		<li>p (Fourier analysis): Specify the total number of periods of measuring time in the input data, which has to be a positive integer. 
		Note: The more periods of data that are collected, the higher the time resolution of the SPP parameters.</li>
		<li>k (numerical differentiation): Define the step size for the numerical differentiation, which has to be a positive integer.</li>
		<li>num_mode (numerical differentiation): Specify num_mode to be either &#8220;0&#8221; (standard differentiation) or &#8220;1&#8221; (looped differentiation). 
		Note: There are two procedures implemented in the numerical differentiation scheme. The &#8220;standard differentiation&#8221; makes no assumptions about the form of the data. It utilizes a forward difference to calculate the derivative for the first 2,000 points of the data, a backward difference for the final 2,000 points, and a centered difference elsewhere. The &#8220;looped differentiation&#8221; assumes that the data is taken under steady-state periodic conditions, and includes an integer number of periods. These assumptions allow a centered difference to be calculated everywhere by looping over the ends of the data.</li>
		<li>Select the run button at the top once all the variables are specified. 
		Note: The software will compute all SPP metrics associated with the data, and then display figures associated with the current analysis run and output a text file containing all the calculated SPP metrics for further analysis.</li>
		<li>Iteratively adjust the number of harmonics to be included in the analysis from the output Fourier spectrum. Include all higher odd harmonics above the noise floor.</li>
	</ol>
	</li>
</ol>
4. Interpreting a LAOS Response
<ol>
	<li>Navigate to the Cole-Cole plot of the instantaneous SPP moduli <img alt="Equation 3" src="/files/ftp_upload/58707/58707eq3.jpg" /> and <img alt="Equation 4" src="/files/ftp_upload/58707/58707eq4.jpg" /> that is automatically generated by the SPP software. 
	Note: A curve in the Cole-Cole plot is considered as the trajectory of the viscoelastic material state, and interpretations can be formed within an oscillation, in intra-cycle processes, or between successive periods, in inter-cycle processes.</li>
	<li>Interpret stiffness by the instantaneous elastic modulus,<img alt="Equation 13" src="/files/ftp_upload/58707/58707eq13.jpg" />, and an increase/decrease of <img alt="Equation 3" src="/files/ftp_upload/58707/58707eq3.jpg" /> that indicates stiffening/softening. See Figure 3.</li>
	<li>Interpret a material&#39;s viscosity based on the instantaneous viscous modulus, <img alt="Equation 4" src="/files/ftp_upload/58707/58707eq4.jpg" />. An increase/decrease in this parameter represents thickening/thinning.</li>
	<li>Transfer the focus to another Cole-Cole plot of the time derivatives of transient moduli <img alt="Equation 14" src="/files/ftp_upload/58707/58707eq14.jpg" /> and <img alt="Equation 15" src="/files/ftp_upload/58707/58707eq15.jpg" />, which provide quantitative information about how much a response is stiffening (<img alt="Equation 16" src="/files/ftp_upload/58707/58707eq16.jpg" />), softening (<img alt="Equation 17" src="/files/ftp_upload/58707/58707eq17.jpg" />), thickening (<img alt="Equation 18" src="/files/ftp_upload/58707/58707eq18.jpg" />), thinning ((<img alt="Equation 19" src="/files/ftp_upload/58707/58707eq19.jpg" />)). See Figure 3. 
	Note: With the values of the derivatives, the rate at which materials undergo stiffening/softening or thickening/thinning can be quantitatively determined.</li>
	<li>Read the center of a trajectory (in a time-weighted average sense) in the Cole-Cole plot of <img alt="Equation 20" src="/files/ftp_upload/58707/58707eq20.jpg" /> as the dynamic moduli, [<img alt="Equation 1" src="/files/ftp_upload/58707/58707eq1.jpg" /> <img alt="Equation 21" src="/files/ftp_upload/58707/58707eq21.jpg" />]. 
	Note: The dynamic moduli are averaged parameters over a cycle of deformation, and are insufficient to provide local information under LAOS.</li>
	<li>Track the relative motion of the trajectory across amplitudes to understand the inter-cycle physics. 
	Note: Focusing on the relative motion of time-weighted average center is equivalent to a traditional strain amplitude sweep of the dynamic moduli. Nonetheless, one can easily analyze the across-amplitude motion of other specific points, for instance, the strain extrema.</li>
	<li>Determine the transient differential viscosity <img alt="Equation 22" src="/files/ftp_upload/58707/58707eq22.jpg" /> and overlay it on top of a steady-shear flow curve. Compare the transient LAOS response with steady-shear conditions.</li>
	<li>Determine the points of maximum <img alt="Equation 12" src="/files/ftp_upload/58707/58707eq12.jpg" /> at the large amplitudes in the Cole-Cole plot of <img alt="Equation 20" src="/files/ftp_upload/58707/58707eq20.jpg" />. See the star labeled in Figure 4c.
	<ol>
		<li>Record the values of <img alt="Equation 23" src="/files/ftp_upload/58707/58707eq23.jpg" /> at those instants.</li>
		<li>Plot them on top of the amplitude sweep of the dynamic moduli. See Figure 4d. 
		Note: Pay attention to any correspondence between the maximum transient elastic modulus and the linear viscoelastic <img alt="Equation 1" src="/files/ftp_upload/58707/58707eq1.jpg" />.</li>
	</ol>
	</li>
	<li>Locate the instants of maximum <img alt="Equation 12" src="/files/ftp_upload/58707/58707eq12.jpg" /> in the elastic Lissajous figure and record the corresponding strain values. See the star labeled in Figure 4a.</li>
	<li>If <img alt="Equation 24" src="/files/ftp_upload/58707/58707eq24.jpg" />, then determine the equilibrium strain <img alt="Equation 25" src="/files/ftp_upload/58707/58707eq25.jpg" /> and the elastic strain <img alt="Equation 26" src="/files/ftp_upload/58707/58707eq26.jpg" />. 
	Note. With the displacement stress <img alt="Equation 27" src="/files/ftp_upload/58707/58707eq27.jpg" />, when <img alt="Equation 28" src="/files/ftp_upload/58707/58707eq28.jpg" /> the equilibrium strain can be determined as <img alt="Equation 29" src="/files/ftp_upload/58707/58707eq29.jpg" /> and the elastic strain can therefore be determined as the difference between strain and equilibrium strain5,13. The requirement of <img alt="Equation 24" src="/files/ftp_upload/58707/58707eq24.jpg" /> is derived and discussed elsewhere15.</li>
	<li>Plot the elastic strain as a function of the imposed strain amplitude. See Figure 4e. If the elastic strain is independent of the strain amplitude, then indicate this critical strain on the amplitude sweep as in Figure 4d.</li>
</ol>

<ol>
	<li>Dolz, M., Hern&#225;ndez, M. J., Delegido, J., Alfaro, M. C., Mu&#241;oz, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Influence+of+xanthan+gum+and+locust+bean+gum+upon+flow+and+thixotropic+behaviour+of+food+emulsions+containing+modified+starch.">Influence of xanthan gum and locust bean gum upon flow and thixotropic behaviour of food emulsions containing modified starch.</a> Journal of Food Engineering. 81 (1), 179-186 (2007).</li><li>Gupta, N., Zeltmann, S. E., Shunmugasamy, V. C., Pinisetty, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Applications+of+Polymer+Matrix+Syntactic+Foams.">Applications of Polymer Matrix Syntactic Foams.</a> JOM. 66 (2), 245-254 (2013).</li><li>Garc&#305;a-Ochoa, F., Santos, V. E., Casas, J. A., G&#243;mez, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Xanthan+gum:+production,+recovery,+and+properties.">Xanthan gum: production, recovery, and properties.</a> Biotechnology Advances. 18 (7), 549-579 (2000).</li><li>Chang, I., Im, J., Prasidhi, A. K., Cho, G. -. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Effects+of+Xanthan+gum+biopolymer+on+soil+strengthening.">Effects of Xanthan gum biopolymer on soil strengthening.</a> Construction and Building Materials. 74, 65-72 (2015).</li><li>Rogers, S. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=In+search+of+physical+meaning:+defining+transient+parameters+for+nonlinear+viscoelasticity.">In search of physical meaning: defining transient parameters for nonlinear viscoelasticity.</a> Rheologica Acta. 56 (5), 501-525 (2017).</li><li>Ferry, J. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Viscoelastic properties of polymers. , (1980).</li><li>Bird, R. B., Armstrong, R. C., Hassager, O. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Dynamics of Polymeric Liquids. Volume 1: Fluid Mechanics. , (1987).</li><li>Park, J. D., Rogers, S. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+transient+behavior+of+soft+glassy+materials+far+from+equilibrium.">The transient behavior of soft glassy materials far from equilibrium.</a> Journal of Rheology. 62 (4), 869-888 (2018).</li><li>Rogers, S., Kohlbrecher, J., Lettinga, M. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+molecular+origin+of+stress+generation+in+worm-like+micelles,+using+a+rheo-SANS+LAOS+approach.">The molecular origin of stress generation in worm-like micelles, using a rheo-SANS LAOS approach.</a> Soft Matter. 8 (30), 7831-7839 (2012).</li><li>Lettinga, M. P., Holmqvist, P., Ballesta, P., Rogers, S., Kleshchanok, D., Struth, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Nonlinear+Behavior+of+Nematic+Platelet+Dispersions+in+Shear+Flow.">Nonlinear Behavior of Nematic Platelet Dispersions in Shear Flow.</a> Phys Rev Lett. 109 (24), 246001 (2012).</li><li>Hyun, K., Wilhelm, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+review+of+nonlinear+oscillatory+shear+tests:+Analysis+and+application+of+large+amplitude+oscillatory+shear+(LAOS).">A review of nonlinear oscillatory shear tests: Analysis and application of large amplitude oscillatory shear (LAOS).</a> Progress in Polymer Science. 36 (12), 1697-1753 (2011).</li><li>Park, J. D., Ahn, K. H., Lee, S. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Structural+change+and+dynamics+of+colloidal+gels+under+oscillatory+shear+flow.">Structural change and dynamics of colloidal gels under oscillatory shear flow.</a> Soft Matter. 11 (48), 9262-9272 (2015).</li><li>Lee, C. -. W., Rogers, S. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+sequence+of+physical+processes+quantified+in+LAOS+by+continuous+local+measures.">A sequence of physical processes quantified in LAOS by continuous local measures.</a> Korea-Australia Rheology Journal. 29 (4), 269-279 (2017).</li><li>Rogers, S. A., Erwin, B. M., Vlassopoulos, D., Cloitre, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+sequence+of+physical+processes+determined+and+quantified+in+LAOS:+Application+to+a+yield+stress+fluid.">A sequence of physical processes determined and quantified in LAOS: Application to a yield stress fluid.</a> Journal of Rheology. 55 (2), 435-458 (2011).</li><li>Wagner, M. H., Rolon-Garrido, V. H., Hyun, K., Wilhelm, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Analysis+of+medium+amplitude+oscillatory+shear+data+of+entangled+linear+and+model+comb+polymers.">Analysis of medium amplitude oscillatory shear data of entangled linear and model comb polymers.</a> Journal of Rheology. 55 (3), 495-516 (2011).</li><li>Radhakrishnan, R., Fielding, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Shear+banding+in+large+amplitude+oscillatory+shear+(LAOStrain+and+LAOStress)+of+soft+glassy+materials.">Shear banding in large amplitude oscillatory shear (LAOStrain and LAOStress) of soft glassy materials.</a> Journal of Rheology. 62 (2), 559-576 (2018).</li><li>Bharadwaj, N. A., Ewoldt, R. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Constitutive+model+fingerprints+in+medium-amplitude+oscillatory+shear.">Constitutive model fingerprints in medium-amplitude oscillatory shear.</a> Journal of Rheology. 59 (2), 557-592 (2015).</li><li>Wilhelm, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Fourier&#8208;Transform+Rheology.">Fourier&#8208;Transform Rheology.</a> Macromolecular Materials and Engineering. 287 (2), 83-105 (2002).</li><li>Ewoldt, R. H., Hosoi, A. E., McKinley, G. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=New+measures+for+characterizing+nonlinear+viscoelasticity+in+large+amplitude+oscillatory+shear.">New measures for characterizing nonlinear viscoelasticity in large amplitude oscillatory shear.</a> Journal of Rheology. 52 (6), 1427-1458 (2008).</li><li>Rogers, S. A., Lettinga, M. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+sequence+of+physical+processes+determined+and+quantified+in+large-amplitude+oscillatory+shear+(LAOS):+Application+to+theoretical+nonlinear+models.">A sequence of physical processes determined and quantified in large-amplitude oscillatory shear (LAOS): Application to theoretical nonlinear models.</a> Journal of Rheology. 56 (1), 1-25 (2011).</li><li>Rogers, S. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+sequence+of+physical+processes+determined+and+quantified+in+LAOS:+An+instantaneous+local+2D/3D+approach.">A sequence of physical processes determined and quantified in LAOS: An instantaneous local 2D/3D approach.</a> Journal of Rheology. 56 (5), 1129-1151 (2012).</li><li>Kim, J., Merger, D., Wilhelm, M., Helgeson, M. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Microstructure+and+nonlinear+signatures+of+yielding+in+a+heterogeneous+colloidal+gel+under+large+amplitude+oscillatory+shear.">Microstructure and nonlinear signatures of yielding in a heterogeneous colloidal gel under large amplitude oscillatory shear.</a> Journal of Rheology. 58 (5), 1359-1390 (2014).</li><li>van der Vaart, K., Rahmani, Y., Zargar, R., Hu, Z., Bonn, D., Schall, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Rheology+of+concentrated+soft+and+hard-sphere+suspensions.">Rheology of concentrated soft and hard-sphere suspensions.</a> Journal of Rheology. 57 (4), 1195-1209 (2013).</li><li>Poulos, A. S., Stellbrink, J., Petekidis, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Flow+of+concentrated+solutions+of+starlike+micelles+under+large-amplitude+oscillatory+shear.">Flow of concentrated solutions of starlike micelles under large-amplitude oscillatory shear.</a> Rheologica Acta. 52 (8-9), 785-800 (2013).</li><li>Armstrong, M. J., Beris, A. N., Rogers, S. A., Wagner, N. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Dynamic+shear+rheology+of+a+thixotropic+suspension:+Comparison+of+an+improved+structure-based+model+with+large+amplitude+oscillatory+shear+experiments.">Dynamic shear rheology of a thixotropic suspension: Comparison of an improved structure-based model with large amplitude oscillatory shear experiments.</a> Journal of Rheology. 60 (3), 433-450 (2016).</li><li>Calabrese, M. A., Wagner, N. J., Rogers, S. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=An+optimized+protocol+for+the+analysis+of+time-resolved+elastic+scattering+experiments.">An optimized protocol for the analysis of time-resolved elastic scattering experiments.</a> Soft Matter. 12 (8), 2301-2308 (2016).</li></ol>

The authors have nothing to disclose.

We have demonstrated how to correctly perform large amplitude oscillatory shear rheometry tests using a commercial rheometer, and to run the SPP analysis freeware to interpret and understand the nonlinear stress responses of two distinct polymer solutions. The SPP framework, which has previously been shown to correlate with structural changes and facilitate understandings of numerous colloidal systems, can be equally applied to polymer systems. The responses of two concentrated polymeric solutions to LAOS have been investigated using the SPP scheme, in which the rheological responses are shown to exhibit complex sequences of processes. These transient intra-cycle interpretations provide essential information on the nonlinear out-of-equilibrium behaviors of polymeric solutions, and provide guidelines for engineers to improve consumer products with desired properties or to transport systems more efficiently.
The gel-like XG solution and the concentrated entangled PEO solution exhibit distinct physical processes that provide clear distinctions between their respective nonlinear behaviors. While the maximum transient elastic modulus of XG remains essentially unchanged across the imposed amplitudes, reminiscent of soft glassy materials that exhibit caging dynamics, the PEO solution displays a local stiffening characteristic that is better described by finite-extensibility concepts typically applied to polymer systems. As a consequence, processes involving each material would be best approximated using glassy and finitely extensible nonlinear elastic (FENE)-type models. In addition to how the maximum elasticity changes with applied strain amplitude, the transient differential viscosity from the two systems show similar behaviors, with apparent overshoots at high shear rates being identified prior to shear thinning. However, the PEO solution displays a lower transient differential viscosity than the steady-state conditions, while the XG solution exhibits no marked difference between steady and dynamic shearing. We therefore identify different pre-yielded processes, but similar post-yield characteristics in the two polymer systems. In both cases, we identify post-yielded conditions that are nearly indistinguishable from steady shearing, showing that it is not necessary to go to the limit of zero frequency in LAOS to obtain reliable information about the flow properties of soft materials.
We identify the nonlinear rheological sequence as containing information about the linear viscoelasticity, the transient flow curves, and the critical strain that is responsible for nonlinear behaviors. This congruence of information obtained via the SPP approach is not possible with any of the FT-based approaches, which treat oscillatory shearing as a special rheological case, with interpretations that are not applicable to other experimental protocols. In contrast, the SPP approach views all materials responses equivalently, providing a clear mechanism for direct comparisons across a range of different tests, such as those made here. We show that the elastic recoverable strain is approximately constant at the point of maximum elasticity for a xanthan gum solution, and this constant elastic strain is indicative of the critical strain of nonlinear regime. We also demonstrate that the transient flow curves can be constructed from the results of the SPP analysis. In a single LAOS test on a concentrated polymeric solution using the SPP approach, we can therefore confidently determine the linear viscoelastic response at that frequency, portions of the steady-state flow curve that correspond to the conditions imposed, and the amplitude above which responses become nonlinear. Overall, this work provides a general approach to performing and understanding nonlinear rheological behaviors of soft matter, with a particular emphasis on polymer solutions. The approach outlined in this work provides an easy-to-implement methodology that provides clear correlation between small and large-amplitude deformation bulk rheology, which can be used to assist in the rational design and optimization of materials under flow.

Concentrated polymeric solutions are used in a variety of industrial applications primarily to increase viscosity, including in foods1 and other consumer products2, enhanced oil recovery3, and soil remediation4. During their processing and use, they are necessarily subjected to large deformations over a range of timescales. Under such processes, they demonstrate rich and complex nonlinear rheological behaviors that depend on the flow or deformation conditions1. Understanding these complex nonlinear rheological behaviors is essential for successfully controlling processes, designing superior products, and maximizing energy efficiency. Aside from the industrial importance, there is a great deal of academic interest in understanding the rheological behaviors of polymeric materials far from equilibrium.
Oscillatory shear tests are a staple component of every thorough rheological characterization because of the orthogonal application of strain and strain rate5, and the ability to independently control the length and time scales probed by tuning the amplitude and frequency. The stress response to small amplitude oscillatory shear strains, which are small enough not to disturb a material&#39;s internal structure, can be decomposed into components in phase with the strain and in phase with strain rate. The coefficients of the components in phase with the strain and the strain rate are collectively referred to as the dynamic moduli6,7, and individually as the storage modulus, <img alt="Equation 1" src="/files/ftp_upload/58707/58707eq1.jpg" />, and loss modulus, <img alt="Equation 2" src="/files/ftp_upload/58707/58707eq2.jpg" />. The dynamic moduli lead to clear elastic and viscous interpretations. However, interpretations based on these dynamic moduli are valid only for small strain amplitudes, where the stress responses to sinusoidal excitations are also sinusoidal. This regime is generally referred to as the small amplitude oscillatory shear (SAOS), or the linear viscoelastic regime. As the imposed deformation becomes larger, changes are induced in the material microstructure, which are reflected in the complexity of the non-sinusoidal transient stress responses8. In this rheologically nonlinear regime, which more closely mimics industrial processing and consumer usage conditions, the dynamic moduli act as poor descriptions of the response. Another way to understand how concentrated soft materials behave out of equilibrium is therefore required.
A number of recent studies9,10,11,12,13,14,15,16 have shown that materials pass through diverse intra-cycle structural and dynamical changes elicited by larger deformations in the medium amplitude oscillatory shear (MAOS)15,17 and large amplitude oscillatory shear (LAOS) regimes. The intra-cycle structural and dynamical changes have different manifestations, such as breakage of microstructure, structural anisotropy, local rearrangements, reformation, and changes in diffusivity. These intra-cycle physical changes in the nonlinear regime lead to the complex nonlinear stress responses that cannot be simply interpreted with the dynamic moduli. As an alternative, several approaches have been suggested for the interpretation of the nonlinear stress responses. Common examples of this are Fourier transform rheology (FT rheology)18, power series expansions11, the Chebyshev description19, and the sequence of physical processes (SPP)5,8,13,14,20 analysis. Although all of these techniques have been shown to be mathematically robust, it is still an unanswered question as to whether any of these techniques can provide clear and reasonable physical explanations of nonlinear oscillatory stress responses. It remains an outstanding challenge to provide concise interpretations of rheological data that correlate to structural and dynamical measures.
In a recent study, the nonlinear stress response of the Soft Glassy Rheology (SGR) model8 and a soft glass made of colloidal star polymers7under oscillatory shear was analyzed through the SPP scheme. Temporal changes in the elastic and viscous properties inherent in nonlinear stress responses were separately quantified by the SPP moduli, <img alt="Equation 3" src="/files/ftp_upload/58707/58707eq3.jpg" /> and <img alt="Equation 4" src="/files/ftp_upload/58707/58707eq4.jpg" />. Furthermore, the rheological transition represented by the transient moduli was accurately correlated to microstructural changes represented by the distribution of mesoscopic elements. In the study of the SGR model8, it was clearly shown that rheological interpretation via the SPP scheme accurately reflects the physical changes under all oscillatory shear conditions in the linear and nonlinear regimes for soft glasses. This unique capability to provide accurate physical interpretation of nonlinear responses of soft glasses makes the SPP method an attractive approach for researchers studying out-of-equilibrium dynamics of polymer solutions and other soft materials.
The SPP scheme is built around viewing rheological behaviors as occurring in a three-dimensional space (<img alt="Equation 5" src="/files/ftp_upload/58707/58707eq5.jpg" />) that consists of the strain (<img alt="Equation 6" src="/files/ftp_upload/58707/58707eq6.jpg" />), strain rate (<img alt="Equation 7" src="/files/ftp_upload/58707/58707eq7.jpg" />), and stress (<img alt="Equation 8" src="/files/ftp_upload/58707/58707eq8.jpg" />)5. In a mathematical sense, the stress responses are treated as multivariable functions of the strain and strain rate (<img alt="Equation 9" src="/files/ftp_upload/58707/58707eq9.jpg" />). As the rheological behavior is regarded as a trajectory in <img alt="Equation 5" src="/files/ftp_upload/58707/58707eq5.jpg" /> (or a multivariable function), a tool for discussing the properties of a trajectory is required. In the SPP approach, the transient moduli <img alt="Equation 3" src="/files/ftp_upload/58707/58707eq3.jpg" /> and <img alt="Equation 4" src="/files/ftp_upload/58707/58707eq4.jpg" /> play such a role. The transient elastic modulus <img alt="Equation 3" src="/files/ftp_upload/58707/58707eq3.jpg" /> and viscous modulus <img alt="Equation 4" src="/files/ftp_upload/58707/58707eq4.jpg" /> are defined as partial derivatives of the stress with respect to the strain (<img alt="Equation 10" src="/files/ftp_upload/58707/58707eq10.jpg" />) and the strain rate (<img alt="Equation 11" src="/files/ftp_upload/58707/58707eq11.jpg" />). Following the physical definition of differential elastic and viscous moduli, the transient moduli quantify the instantaneous influence of strain and strain rate on the stress response respectively, whereas other analysis methods cannot provide any information on elastic and viscous properties separately.
The SPP approach enriches the interpretation of the oscillatory shear tests. With the SPP analysis, the complex nonlinear rheological behaviors of concentrated polymeric solutions in LAOS can be directly related to the linear rheological behaviors in SAOS. We show in this work how the maximum transient elastic modulus (<img alt="Equation 12" src="/files/ftp_upload/58707/58707eq12.jpg" />max) near the strain extrema corresponds to the storage modulus in the linear regime (SAOS). Furthermore, we show how the transient viscous modulus (<img alt="Equation 4" src="/files/ftp_upload/58707/58707eq4.jpg" />) during a LAOS cycle traces the steady state flow curve. In addition to providing details of the complex sequence of processes that concentrated polymer solutions go through under LAOS, the SPP scheme also provides information regarding the recoverable strain in the material. This information, which is not obtainable through other approaches, is a useful measure of how much a material will recoil once stress is removed. Such behavior has impact on the printability of concentrated solutions for 3D printing applications, as well as screen printing, fiber formation, and flow cessation. A number of recent studies5,8,13 clearly indicate that the recoverable strain is not necessarily the same as the strain imposed during LAOS experiments. For instance, a study of soft colloidal glasses under LAOS13 found that the recoverable strain is only 5% when significantly larger total strain (420%) is imposed. Other studies16,21,22,23,24 using the cage modulus21 also conclude that linear elasticity can be observed under LAOS at the point close to the strain maxima, implying that the materials experienced relatively small deformation at those instants. The SPP scheme is the only framework for understanding LAOS that accounts for a shift in the strain equilibrium that leads to a difference between the recoverable and the total strains.
This article aims to facilitate understandings and ease of use of the SPP analysis method by providing a detailed protocol for a LAOS analysis freeware, using two concentrated polymer solutions, a 4 wt% xanthan gum (XG) aqueous solution and a 5 wt% PEO in DMSO solution. These systems are chosen because of their broad range of application and rheologically interesting properties. Xanthan gum, a natural high-molecular-weight polysaccharide, is an exceptionally effective stabilizer for aqueous systems and commonly applied as a food additive to provide desired viscosification or in oil drilling to increase viscosity and yield points of drilling muds. PEO has a unique hydrophilic property and is often used in pharmaceutical products and controlled release systems as well as soil remediation activities. These polymeric systems are tested under various oscillatory shear conditions that are intended to approximate processing, transport, and end-use conditions. Although these practical conditions may not necessarily involve flow reversal as in oscillatory shear, the flow field can be easily approximated and tuned with the independent control of applied amplitude and imposed frequency in an oscillatory test. Furthermore, the SPP scheme can be used as described here to understand a broad range of flow types, including those that do not include flow reversals such as the recently-proposed UD-LAOS25, in which large amplitude oscillations are applied in one direction only (leading to the moniker &#34;uni-directional LAOS&#34;). For simplicity, and for illustrative purposes, we restrict the current study to traditional LAOS, which does include periodic flow reversal. The measured rheological responses are analyzed with the SPP approach. We demonstrate how to use the SPP software with simple explanations on salient calculation steps to improve readers&#39; understanding and usage. A legend for interpreting the SPP analysis results is introduced, according to which the type of rheological transition is identified. Representative SPP analysis results of the two polymers under various oscillatory shear conditions are displayed, in which we clearly identify a sequence of physical processes that contains information on the material&#39;s linear viscoelastic response as well as the steady-state flow properties of the material.
This protocol provides salient details of how to accurately perform nonlinear rheological experiments, as well as a step-by-step guide to analyzing and understanding rheological responses with the SPP framework, as shown in Figure 1. We begin by providing an introduction to the instrument setup and calibrations, followed by specific commands for making a commercially-available rheometer collect high-quality transient data. Once the rheological data have been obtained, we introduce the SPP analysis freeware, with a detailed manual. Further, we discuss how to understand the time-dependent response of the two concentrated polymer solutions within the SPP scheme, by comparing the results obtained from LAOS with the linear-regime frequency sweep and the steady-state flow curve. These results clearly identify that the polymer solutions transition between distinct rheological states within an oscillation, allowing for a more detailed picture of their nonlinear transient rheology to emerge. These data can be used to optimize processing conditions for product formation, transport, and use. These time-dependent responses further provide potential pathways to clearly form structure-property-processing relationships by coupling the rheology with microstructural information obtained from small-angle scattering of neutrons, X-rays, or light (SANS, SAXS, and SALS, respectively), microscopy, or detailed simulations.

<table>
	<tbody>
		<tr>
			<td>SPP analysis software</td>
			<td>Simon Rogers Group (UIUC)</td>
			<td>SPPplus_v1p1</td>
			<td>Attached as supplementary files</td>
		</tr>
		<tr>
			<td>MATLAB</td>
			<td>Mathwork</td>
			<td></td>
			<td></td>
		</tr>
		<tr>
			<td>Rheometer</td>
			<td>Anton Paar</td>
			<td>MCR 702 TwinDrive</td>
			<td></td>
		</tr>
		<tr>
			<td>50mm 2-degree cone</td>
			<td>Anton Paar</td>
			<td>CP50-2</td>
			<td>Upper measuring system</td>
		</tr>
		<tr>
			<td>50mm plate</td>
			<td>Anton Paar</td>
			<td>PP50</td>
			<td>Lower measuring system</td>
		</tr>
		<tr>
			<td>Xanthan gum (XG)</td>
			<td>Sigma-Aldrich</td>
			<td>11138-66-2</td>
			<td></td>
		</tr>
		<tr>
			<td>Polyethylene oxide (PEO)</td>
			<td>Sigma-Aldrich</td>
			<td>25322-68-3</td>
			<td>Mv=1,000,000</td>
		</tr>
		<tr>
			<td>Dimethyl sulfoxide (DMSO)</td>
			<td>Sigma-Aldrich</td>
			<td>67-68-5</td>
			<td></td>
		</tr>
	</tbody>
</table>

studying large amplitude oscillatory shear response of soft materials

We investigate the sequence of physical processes exhibited during large amplitude oscillatory shearing (LAOS) of polyethylene oxide (PEO) in dimethyl sulfoxide (DMSO) and xanthan gum in water &#8212; two concentrated polymer solutions used as viscosifiers in foods, enhanced oil recovery, and soil remediation. Understanding the nonlinear rheological behavior of soft materials is important in the design and controlled manufacturing of many consumer products. It is shown how the response to LAOS of these polymer solutions can be interpreted in terms of a clear transition from linear viscoelasticity to viscoplastic deformation and back again during a period. The LAOS results are analyzed via the fully quantitative Sequence of Physical Processes (SPP) technique, using free MATLAB-based software. A detailed protocol of performing a LAOS measurement with a commercial rheometer, analyzing nonlinear stress responses with the freeware, and interpreting physical processes under LAOS is presented. It is further shown that, within the SPP framework, a LAOS response contains information regarding the linear viscoelasticity, the transient flow curves, and the critical strain responsible for the onset of nonlinearity.

Representative results of the SPP analysis from XG and PEO/DMSO solutions under oscillatory shear tests are presented in Figures 4 and 5. We first present the raw data as elastic (<img alt="Equation 30" src="/files/ftp_upload/58707/58707eq30.jpg" />) and viscous (<img alt="Equation 31" src="/files/ftp_upload/58707/58707eq31.jpg" />) Lissajous-Bowditch curves in Figures 4a, 4b, 5a and 5b. To fully understand the intra-cycle physics, the time-dependent Cole-Cole plots obtained from the SPP freeware are presented in Figures 4c and 5c. Interpretations of the plots are discussed in the manner laid out by the legend in Figure 3 and&#160;protocol steps 4.2-4.7, where the relative motion of the trace quantitatively indicates whether the material undergoes stiffening/softening or thickening/thinning in an intra-cycle sense. The time-weighted centers of these trajectories, which represent the average elastic and viscous moduli, correspond to the dynamic moduli, <img alt="Equation 1" src="/files/ftp_upload/58707/58707eq1.jpg" />&#160;and&#160;<img alt="Equation 32" src="/files/ftp_upload/58707/58707eq32.jpg" />, shown in Figure 4d and 5d. In the case of large deformations, average parameters are insufficient to describe the material response at any particular instant. Forming a bridge between rheological data and microstructural evolutions has proved a difficult task. Microstructural information obtained from either scattering9,26 or simulation12 is often time-resolved and requires a rheological study that matches the temporal resolution. A more complete discussion of the linking the macroscopic SPP analysis and microstructural details can be found in a recent study of soft glassy materials8.
Using the SPP scheme, we are also able to determine the elastic recoverable strain at moments when the material response is predominantly elastic. In particular, the gel-like structure of XG responds in ways that are reminiscent of soft glassy materials, where the responses go through instants of linear-regime viscoelasticity across the large amplitudes as shown in Figure 4d. Indeed, we identify the instantaneous SPP elastic modulus at large amplitudes in the XG solution that is more than three orders of magnitude larger than the traditional storage modulus, showing the clear benefit of the local measures. Similar results have been observed in studies of soft colloidal glasses16,21,22,23,24, where the points of linear-like elasticity also take place at positions near the strain extrema. This indicates that the material equilibrium is well separated from the place where the experiment started, at zero strain. With the SPP analysis, it is shown in Figure 4e that the elastic recoverable strain at the point of maximum elasticity remains nearly constant at 16%, even when the applied strain is as large as 4,000%. This constant recoverable strain of approximately 16% corresponds to the critical strain amplitude, <img alt="Equation 33" src="/files/ftp_upload/58707/58707eq33.jpg" />, above which nonlinear behavior is observed in the strain amplitude sweep of Figure 4d.
In the case of the PEO solution, the maximum transient elastic modulus across different amplitudes is shown in Figure 5d. We identify, using the SPP approach, an increasing stiffness as the amplitude increases, while the storage modulus shows only softening. At the largest amplitudes probed, we identify an instantaneous modulus that is more than an order of magnitude larger than the traditionally-defined storage modulus. The magnitudes of the transient elastic and viscous moduli are comparable at the instants of largest elasticity, meaning that the condition for the SPP to correctly identify the elastic strain is not met.
The major advantage of the quantitative SPP scheme is that elastic and viscous properties can be clearly determined at each point in the cycle. In the previous section, it was established that at instants close to the strain extrema, the XG solution responds as if it is in its linear viscoelastic limit while the PEO solution displays a modulus that is marginally larger than that exhibited in the linear regime. We now turn our attention to the next major component in the sequence of physical processes exhibited by both polymer solutions, the flow condition.
The transient differential viscosity, defined as the transient viscous modulus divided by the frequency,&#160;<img alt="Equation 34" src="/files/ftp_upload/58707/58707eq34.jpg" />, is displayed in Figure 6 on top of the steady-shear flow viscosity, determined from independent steady-shear tests. A similar response is observed from both materials, where the transient differential viscosities initially remain constant at low shear rates, followed by an overshoot, before decreasing rapidly. The transient differential viscosities of both solutions change with shear rate approximately the same as the steady-shear flow viscosity, albeit with transient differential viscosities that are slightly below the steady-state conditions. The steady-shear flow response can be viewed as a LAOS experiment in the limit of zero frequency; nonetheless, with the SPP analysis scheme, the transient flow behaviors at any arbitrary imposed frequency can be quantitatively constructed.
The distinct sequence of physical processes exhibited by XG at a strain amplitude of 4000% is displayed in Figure 7, where the symbols split the Lissajous-Bowditch curve into different processes of interest. We begin in the region labelled as region #1, which we identify as being viscoplastic in nature. In this interval of the response, the SPP analysis scheme shows nearly zero elasticity, as determined by <img alt="Equation 3" src="/files/ftp_upload/58707/58707eq3.jpg" />, which indicates no strain-dependence to the stress. As the shear rate begins to decrease close to the strain extremum, the XG solution stiffens, indicating that the structure responsible for the linear viscoelastic response begins to reform. We term this &#39;restructuring&#39;. The elastic recoverable strain at this point, at around 16%, is much smaller than the total deformation, which is consistent with the linear-regime viscoelasticity of these gel-like and other glassy systems. A rapid transition from elastic to viscous behaviors, reminiscent of yielding or destructuring, takes place once sufficient strain is acquired from reversal, and is followed by a stress overshoot, during which there is a sharp change in the transient moduli. During the portion of the overshoot when the stress is decreasing, the instantaneous viscous modulus, <img alt="Equation 35" src="/files/ftp_upload/58707/58707eq35.jpg" /> is momentarily negative, reflecting the decreasing stress with increasing shear rate. Portions of negative <img alt="Equation 35" src="/files/ftp_upload/58707/58707eq35.jpg" /> are therefore not observed in the PEO solutions because of their lack of any overshoot. Lastly, the system goes back to the viscoplastic deformation regime and experiences the distinct intra-cycle sequence twice over a cycle of oscillation.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/58707/58707fig1.jpg" /> 
Figure 1: A schematic to illustrate a complete process of performing, analyzing and understanding rheological experiments. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/58707fig1large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/58707/58707fig2.jpg" /> 
Figure 2: Detailed procedure of loading materials. (a) Attach the lower (PP50) and upper (CP50-2) geometries followed by setting the zero-gap position. (b) Load the material onto the center of the lower plate with a pipette or spatula while avoiding bubbles. (c) Command the upper geometry to trim gap. Slight overfilling is expected in this step unless pipetting with precise volume. Underfilling should be prevented. (d) Gently trim the overfill at the edge of geometries with a square-ended spatula. (e) Continue to the measurement gap only when the loading and trimming are good, such that no underfilling is observed around the perimeter of the geometry, and the edges show no distinct fractures. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/58707fig2large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/58707/58707fig3.jpg" /> 
Figure 3: Trajectories in time-dependent Cole-Cole plots can be interpreted through these legends. (a) Cole-Cole plot in <img alt="Equation 20" src="/files/ftp_upload/58707/58707eq20.jpg" />-space, (b) in <img alt="Equation 20" src="/files/ftp_upload/58707/58707eq20.jpg" style="font-size: 14px;" />-space. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/58707fig3large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/58707/58707fig4.jpg" /> 
Figure 4: SPP-LAOS analysis from the 4 wt% XG solution at the frequency of 0.316 rad/s. The raw data are presented as elastic (a) and viscous (b) Lissajous-Bowditch curves. (c) Cole-Cole plot of transient moduli <img alt="Equation 37" src="/files/ftp_upload/58707/58707eq37.jpg" />, where the dashed lines represent the linear-regime dynamic moduli. (d) The transient moduli determined at the point of maximum elasticity as a function of strain amplitudes. (e) Elastic recoverable strain at the instant of maximum <img alt="Equation 3" src="/files/ftp_upload/58707/58707eq3.jpg" /> as a function of strain amplitude. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/58707fig4large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/58707/58707fig5.jpg" /> 
Figure 5: SPP-LAOS analysis from 5 wt% PEO in DMSO solution at the frequency of 1.26 rad/s. (a) Elastic and (b) viscous Lissajous-Bowditch curves. (c) Cole-Cole plot of transient moduli <img alt="Equation 37" src="/files/ftp_upload/58707/58707eq37.jpg" />, where the dashed lines represent the linear-regime dynamic moduli. (d) The dynamic moduli as a function of strain amplitudes. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/58707fig5large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 6" class="xfigimg" src="/files/ftp_upload/58707/58707fig6.jpg" /> 
Figure 6: The transient differential viscosity plotted on top of the steady-shear flow curve from the XG (a) and PEO/DMSO (b) systems. Lines show transient differential viscosity <img alt="Equation 22" src="/files/ftp_upload/58707/58707eq22.jpg" /> determined from LAOS tests while star symbols represent steady-shear flow viscosity. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/58707fig6large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 7" class="xfigimg" src="/files/ftp_upload/58707/58707fig7.jpg" /> 
Figure 7: The sequence of physical processes under LAOS from the XG solutions. The symbols shown on elastic Lissajous-Bowditch curves (a) correspond to the ones in the time-dependent Cole-Cole plot of transient moduli (b). <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/58707fig7large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
Supplemental File 1:&#160;<a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/1_RunSPPplus_v1.m" target="_blank">Please click here to download this file.</a>
Supplemental File 2:&#160;<a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/2_SPPplus_print_v1.m" target="_blank">Please click here to download this file.</a>
Supplemental File 3:&#160;<a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/3_SPPplus_numerical_v1.m" target="_blank">Please click here to download this file.</a>
Supplemental File 4:&#160;<a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/4_SPPplus_read_v1.m" target="_blank">Please click here to download this file.</a>
Supplemental File 5:&#160;<a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/5_SPPplus_fourier_v1.m" target="_blank">Please click here to download this file.</a>
Supplemental File 6:&#160;<a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/6_SPPplus_figures_v1.m" target="_blank">Please click here to download this file.</a>
Supplemental File 7:&#160;<a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/7_0100.txt" target="_blank">Please click here to download this file.</a>
Supplemental File 8:&#160;<a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/8_1000.txt" target="_blank">Please click here to download this file.</a>
Supplemental File 9:&#160;<a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/58707/9_0010.txt" target="_blank">Please click here to download this file.</a>

Watch this Scientific Journal Video about Studying Large Amplitude Oscillatory Shear Response of Soft Materials at JoVE.com

Studying Large Amplitude Oscillatory Shear Response of Soft Materials

We present a detailed protocol outlining how to perform nonlinear oscillatory shear rheology on soft materials, and how to run the SPP-LAOS analysis to understand the responses as a sequence of physical processes.

We investigate the sequence of physical processes exhibited during large amplitude oscillatory shearing (LAOS) of polyethylene oxide (PEO) in dimethyl ...

studying-large-amplitude-oscillatory-shear-response-of-soft-materials

Nanyang Technological University

Education

Research

JoVE Core

JoVE Journal

Physics

Engineering

Mechanical Engineering

Electrical Engineering

Unit 1

Transient and Steady-state Response Analysis

Principal Stresses under a Given Loading

Equilibrium and Elasticity

Internal Forces

SCIENTISTS IN ACTION

12.4K Views.  University of Illinois at Urbana-Champaign. We present a detailed protocol outlining how to perform nonlinear oscillatory shear rheology on soft materials, and how to run the SPP-LAOS analysis to understand the responses as a sequence of physical processes.

Video: Studying Large Amplitude Oscillatory Shear Response of Soft Materials

Controlling the Size, Shape and Stability of Supramolecular Polymers in Water

For aqueous based supramolecular polymers, the simultaneous control over shape, size and stability is very difficult1. At the same time, the ability to do so is highly important in view of a number of applications in functional soft matter including electronics, biomedical engineering, and sensors. In the past, successful strategies to control the size and shape of supramolecular polymers typically focused on the use of templates2,3, end cappers4 or selective solvent techniques5. 
Here we disclose a strategy based on self-assembling discotic amphiphiles that leads to the control over stack length and shape of ordered, chiral columnar aggregates. By balancing electrostatic repulsive interactions on the hydrophilic rim and attractive non-covalent forces within the hydrophobic core of the polymerizing building block, we manage to create small and discrete spherical objects6,7. Increasing the salt concentration to screen the charges induces a sphere-to-rod transition. Intriguingly, this transition is expressed in an increase of cooperativity in the temperature-dependent self-assembly mechanism, and more stable aggregates are obtained. 
For our study we select a benzene-1,3,5-tricarboxamide (BTA) core connected to a hydrophilic metal chelate via a hydrophobic, fluorinated L-phenylalanine based spacer (Scheme 1). The metal chelate selected is a Gd(III)-DTPA complex that contains two overall remaining charges per complex and necessarily two counter ions. The one-dimensional growth of the aggregate is directed by &pi;-&pi; stacking and intermolecular hydrogen bonding. However, the electrostatic, repulsive forces that arise from the charges on the Gd(III)-DTPA complex start limiting the one-dimensional growth of the BTA-based discotic once a certain size is reached. At millimolar concentrations the formed aggregate has a spherical shape and a diameter of around 5 nm as inferred from 1H-NMR spectroscopy, small angle X-ray scattering, and cryogenic transmission electron microscopy (cryo-TEM). The strength of the electrostatic repulsive interactions between molecules can be reduced by increasing the salt concentration of the buffered solutions. This screening of the charges induces a transition from spherical aggregates into elongated rods with a length &gt; 25 nm. Cryo-TEM allows to visualise the changes in shape and size. In addition, CD spectroscopy permits to derive the mechanistic details of the self-assembly processes before and after the addition of salt. Importantly, the cooperativity -a key feature that dictates the physical properties of the produced supramolecular polymers- increases dramatically upon screening the electrostatic interactions. This increase in cooperativity results in a significant increase in the molecular weight of the formed supramolecular polymers in water.

The goal of this experiment is to determine and control the size, shape and stability of self-assembled discotic amphiphiles in water. For aqueous based supramolecular polymers such level of control is very difficult. We apply a strategy using both repulsive and attractive interactions. The experimental techniques applied to characterize this system are broadly applicable.

For aqueous based supramolecular polymers, the simultaneous control over shape, size and stability is very difficult1. At the same time, the ability to do ...

Micro 3D Printing Using a Digital Projector and its Application in the Study of Soft Materials Mechanics

Buckling is a classical topic in mechanics. While buckling has long been studied as one of the major structural failure modes1, it has recently drawn new attention as a unique mechanism for pattern transformation. Nature is full of such examples where a wealth of exotic patterns are formed through mechanical instability2-5. Inspired by this elegant mechanism, many studies have demonstrated creation and transformation of patterns using soft materials such as elastomers and hydrogels6-11. Swelling gels are of particular interest because they can spontaneously trigger mechanical instability to create various patterns without the need of external force6-10. Recently, we have reported demonstration of full control over buckling pattern of micro-scaled tubular gels using projection micro-stereolithography (P&mu;SL), a three-dimensional (3D) manufacturing technology capable of rapidly converting computer generated 3D models into physical objects at high resolution12,13. Here we present a simple method to build up a simplified P&mu;SL system using a commercially available digital data projector to study swelling-induced buckling instability for controlled pattern transformation.
A simple desktop 3D printer is built using an off-the-shelf digital data projector and simple optical components such as a convex lens and a mirror14. Cross-sectional images extracted from a 3D solid model is projected on the photosensitive resin surface in sequence, polymerizing liquid resin into a desired 3D solid structure in a layer-by-layer fashion. Even with this simple configuration and easy process, arbitrary 3D objects can be readily fabricated with sub-100 &mu;m resolution.
This desktop 3D printer holds potential in the study of soft material mechanics by offering a great opportunity to explore various 3D geometries. We use this system to fabricate tubular shaped hydrogel structure with different dimensions. Fixed on the bottom to the substrate, the tubular gel develops inhomogeneous stress during swelling, which gives rise to buckling instability. Various wavy patterns appear along the circumference of the tube when the gel structures undergo buckling. Experiment shows that circumferential buckling of desired mode can be created in a controlled manner. Pattern transformation of three-dimensionally structured tubular gels has significant implication not only in mechanics and material science, but also in many other emerging fields such as tunable matamaterials.

We demonstrate controlled pattern transformation of swelling gel tubes by elastic instability. A simple projection micro stereo-lithography setup is built using an off-the-shelf digital data projector to fabricate three-dimensional polymeric structures in a layer-by-layer fashion. Swelling hydrogel tubes under mechanical constraint display various circumferential buckling modes depending on dimension.

Buckling is a classical topic in mechanics. While buckling has long been studied as one of the major structural failure modes1, it has recently drawn new ...

Sub-nanometer Resolution Imaging with Amplitude-modulation Atomic Force Microscopy in Liquid

Atomic force microscopy (AFM) has become a well-established technique for nanoscale imaging of samples in air and in liquid. Recent studies have shown that when operated in amplitude-modulation (tapping) mode, atomic or molecular-level resolution images can be achieved over a wide range of soft and hard samples in liquid. In these situations, small oscillation amplitudes (SAM-AFM) enhance the resolution by exploiting the solvated liquid at the surface of the sample. Although the technique has been successfully applied across fields as diverse as materials science, biology and biophysics and surface chemistry, obtaining high-resolution images in liquid can still remain challenging for novice users. This is partly due to the large number of variables to control and optimize such as the choice of cantilever, the sample preparation, and the correct manipulation of the imaging parameters. Here, we present a protocol for achieving high-resolution images of hard and soft samples in fluid using SAM-AFM on a commercial instrument. Our goal is to provide a step-by-step practical guide to achieving high-resolution images, including the cleaning and preparation of the apparatus and the sample, the choice of cantilever and optimization of the imaging parameters. For each step, we explain the scientific rationale behind our choices to facilitate the adaptation of the methodology to every user's specific system.

We present a method for achieving sub-nanometer resolution images with amplitude-modulation (tapping mode) atomic force microscopy in liquid. The method is demonstrated on commercial atomic force microscopes. We explain the rationale behind our choices of parameters and suggest strategies for resolution optimization.

Atomic force microscopy (AFM) has become a well-established technique for nanoscale imaging of samples in air and in liquid. Recent studies have shown ...

Studying Dynamic Processes of Nano-sized Objects in Liquid using Scanning Transmission Electron Microscopy

Samples fully embedded in liquid can be studied at a nanoscale spatial resolution with Scanning Transmission Electron Microscopy (STEM) using a microfluidic chamber assembled in the specimen holder for Transmission Electron Microscopy (TEM) and STEM. The microfluidic system consists of two silicon microchips supporting thin Silicon Nitride (SiN) membrane windows. This article describes the basic steps of sample loading and data acquisition. Most important of all is to ensure that the liquid compartment is correctly assembled, thus providing a thin liquid layer and a vacuum seal. This protocol also includes a number of tests necessary to perform during sample loading in order to ensure correct assembly. Once the sample is loaded in the electron microscope, the liquid thickness needs to be measured. Incorrect assembly may result in a too-thick liquid, while a too-thin liquid may indicate the absence of liquid, such as when a bubble is formed. Finally, the protocol explains how images are taken and how dynamic processes can be studied. A sample containing AuNPs is imaged both in pure water and in saline.

This protocol describes the operation of a liquid flow specimen holder for scanning transmission electron microscopy of AuNPs in water, as used for the observation of nanoscale dynamic processes.

Samples fully embedded in liquid can be studied at a nanoscale spatial resolution with Scanning Transmission Electron Microscopy (STEM) using a ...

Preparation of Liquid Crystal Networks for Macroscopic Oscillatory Motion Induced by Light

A strategy based on doped liquid crystalline networks is described to create mechanical self-sustained oscillations of plastic films under continuous light irradiation. The photo-excitation of dopants that can quickly dissipate light into heat, coupled with anisotropic thermal expansion and self-shadowing of the film, gives rise to the self-sustained deformation. The oscillations observed are influenced by the dimensions and the modulus of the film, and by the directionality and intensity of the light. The system developed offers applications in energy conversion and harvesting for soft-robotics and automated systems.
The general method described here consists of creating free-standing liquid crystalline films and characterizing the mechanical and thermal effects observed. The molecular alignment is achieved using alignment layers (rubbed polyimide), commonly used in the display manufacturing industry. To obtain actuators with large deformation, the mesogens are aligned and polymerized in a splay/bend configuration, i.e., with the director of the liquid crystals (LCs) going gradually from planar to homeotropic through the film thickness. Upon irradiation, the mechanical and thermal oscillations obtained are monitored with a high-speed camera. The results are further quantified by image analysis using an image processing program.

The goal of the protocol is to create liquid crystalline polymer films that can mechanically oscillate under continuous light irradiation. We describe in great detail the conception of free-standing films, from the liquid crystal alignment method to the photo-actuation. The experimental protocol applied to prepare this material is broadly applicable.

A strategy based on doped liquid crystalline networks is described to create mechanical self-sustained oscillations of plastic films under continuous ...

Atmospheric Pressure Fabrication of Large-Sized Single-Layer Rectangular SnSe Flakes

Tin selenide (SnSe) belongs to the family of layered metal chalcogenide materials with a buckled structure like phosphorene, and has shown potential for applications in two-dimensional nanoelectronics devices. Although many methods to synthesize SnSe nanocrystals have been developed, a simple way to fabricate large-sized single-layer SnSe flakes remains a great challenge. Herein, we show the experimental method to directly grow large-sized single-layer rectangular SnSe flakes on commonly used SiO2/Si insulating substrates using a straightforward two-step fabrication method in an atmospheric pressure quartz tube furnace system. The single-layer rectangular SnSe flakes with an average thickness of ~6.8 &#197; and lateral dimensions of about 30 &#181;m &#215; 50 &#181;m were fabricated by a combination of vapor transport deposition technique and nitrogen etching route. We characterized the morphology, microstructure, and electrical properties of the rectangular SnSe flakes and obtained excellent crystallinity and good electronic properties. This article about the two-step fabrication method can help researchers grow other similar two-dimensional, large-sized, single-layer materials using an atmospheric pressure system.

A protocol is presented demonstrating a two-step fabrication technique to grow large-sized single-layer rectangular shaped SnSe flakes on low-cost SiO2/Si dielectrics wafers in an atmospheric pressure quartz tube furnace system.

Tin selenide (SnSe) belongs to the family of layered metal chalcogenide materials with a buckled structure like phosphorene, and has shown potential for ...

The Diffusion of Passive Tracers in Laminar Shear Flow

A simple method to experimentally observe and measure the dispersion of a passive tracer in a laminar fluid flow is described. The method consists of first injecting fluorescent dye directly into a pipe filled with distilled water and allowing it to diffuse across the cross-section of the pipe to obtain a uniformly distributed initial condition. Following this period, the laminar flow is activated with a programmable syringe pump to observe the competition of advection and diffusion of the tracer through the pipe. Asymmetries in the tracer distribution are studied and correlations between the pipe cross-section and the shape of the distribution is shown: thin channels (aspect ratio &#60;&#60; 1) produce tracers arriving with sharp fronts and tapering tails (front-loaded distributions), while thick channels (aspect ratio ~1) present the opposite behavior (back-loaded distributions). The experimental procedure is applied to capillary tubes of various geometries and is particularly relevant to microfluidic applications by dynamical similarity.

A protocol for the study of the diffusion of passive tracers in laminar pressure-driven flow is presented. The procedure is applicable to various capillary pipe geometries.

A simple method to experimentally observe and measure the dispersion of a passive tracer in a laminar fluid flow is described. The method consists of ...

Fabrication of Soft Pneumatic Network Actuators with Oblique Chambers

Soft pneumatic network actuators have become one of the most promising actuation devices in soft robotics which benefits from their large bending deformation and low input. However, their monotonous bending motion form in two-dimensional (2-D) space keeps them away from wide applications. This paper presents a detailed fabrication method of soft pneumatic network actuators with oblique chambers, to explore their motions in three-dimensional (3-D) space. The design of oblique chambers enables actuators with tunable coupled bending and twisting capabilities, which gives them the possibility to move dexterously in flexible manipulators, to become biologically inspired robots and medical devices. The fabrication process is based on the molding method, including the silicone elastomer preparation, chamber and base parts fabrication, actuator assembly, tubing connections, checks for leaks, and actuator repair. The fabrication method guarantees the rapid manufacturing of a series of actuators with only a few modifications in the molds. The test results show the high quality of the actuators and their prominent bending and twisting capabilities. Experiments of the gripper demonstrate the advantages of the development in adapting to objects with different diameters and providing sufficient friction.

Here we present a fabrication method of soft pneumatic network actuators with oblique chambers. The actuators are capable of generating coupled bending and twisting motions, which broadens their application in soft robotics.

Soft pneumatic network actuators have become one of the most promising actuation devices in soft robotics which benefits from their large bending ...

Rapid Manufacturing of Thin Soft Pneumatic Actuators and Robots

This protocol describes a method for rapid manufacturing of soft pneumatic actuators and robots with an ultrathin form factor using a heat press and a laser cutter machine. The method starts with the lamination of thermoplastic polyurethane (TPU) sheets using a heat press for 10 min at the temperature of ~93 &#176;C. Next, the parameters of the laser cutter machine are optimized to produce a rectangular balloon with maximum burst pressure. Using the optimized parameters, the soft actuators are laser cut/welded three times sequentially. Next, a dispensing needle is attached to the actuator, allowing it to be inflated. The effect of geometrical parameters on the deflection of the actuator are studied systematically by varying the channel width and length. Finally, the performance of the actuator is characterized using an optical camera and a fluid dispenser. Conventional fabrication methods of soft pneumatic actuators based on silicone molding are time consuming (several hours). They also result in strong but bulky actuators, which limits the actuator&#8217;s applications. Moreover, microfabrication of thin pneumatic actuators is both time-consuming and expensive. The proposed manufacturing method in the current work resolves these issues by introducing a fast, simple, and cost-effective fabrication method of ultrathin pneumatic actuators.

This protocol describes a method for rapid manufacturing of soft pneumatic actuators and robots with a thin form factor. The fabrication method starts with lamination of thermoplastic polyurethane (TPU) sheets followed by laser cutting/welding of a two-dimensional pattern to form actuators and robots.

This protocol describes a method for rapid manufacturing of soft pneumatic actuators and robots with an ultrathin form factor using a heat press and a ...

Fabrication of Carbon-Based Ionic Electromechanically Active Soft Actuators

Ionic electromechanically active capacitive laminates are a type of smart material that move in response to electrical stimulation. Due to the soft, compliant and biomimetic nature of this deformation, actuators made of the laminate have received increasing interest in soft robotics and (bio)medical applications. However, methods to easily fabricate the active material in large (even industrial) quantities and with a high batch-to-batch and within-batch repeatability are needed to transfer the knowledge from laboratory to industry. This protocol describes a simple, industrially scalable and reproducible method for the fabrication of ionic carbon-based electromechanically active capacitive laminates and the preparation of actuators made thereof. The inclusion of a passive and chemically inert (insoluble) middle layer (e.g., a textile-reinforced polymer network or microporous Teflon) distinguishes the method from others. The protocol is divided into five steps: membrane preparation, electrode preparation, current collector attachment, cutting and shaping, and actuation. Following the protocol results in an active material that can, for example, compliantly grasp and hold a randomly shaped object as demonstrated in the article.

This article describes a fast and simple manufacturing process of ionic electromechanically active composite materials for actuators in biomedical, biomimetic and soft robotics applications. The key fabrication steps, their importance for the actuators&#39; final properties, and some of the main characterization techniques are described in detail.

Ionic electromechanically active capacitive laminates are a type of smart material that move in response to electrical stimulation. Due to the soft, ...