<?xml-model href='http://www.tei-c.org/release/xml/tei/custom/schema/relaxng/tei_all.rng' schematypens='http://relaxng.org/ns/structure/1.0'?><TEI xmlns="http://www.tei-c.org/ns/1.0">
	<teiHeader>
		<fileDesc>
			<titleStmt><title level='a'>Effects of shear-induced crystallization on the complex viscosity of lamellar-structured concentrated surfactant solutions</title></titleStmt>
			<publicationStmt>
				<publisher>Royal Society of Chemistry</publisher>
				<date>03/20/2024</date>
			</publicationStmt>
			<sourceDesc>
				<bibl> 
					<idno type="par_id">10526227</idno>
					<idno type="doi">10.1039/D3SM01198D</idno>
					<title level='j'>Soft Matter</title>
<idno>1744-683X</idno>
<biblScope unit="volume">20</biblScope>
<biblScope unit="issue">15</biblScope>					

					<author>Parth U Kelkar</author><author>Matthew Kaboolian</author><author>Ria D Corder</author><author>Marco Caggioni</author><author>Seth Lindberg</author><author>Kendra A Erk</author>
				</bibl>
			</sourceDesc>
		</fileDesc>
		<profileDesc>
			<abstract><ab><![CDATA[Material relationships at low temperatures were determined for concentrated surfactant solutions using a combination of rheological experiments, cross-polarized microscopy, calorimetry, and small angle X-ray scattering. A lamellar structured 70 wt% solution of sodium laureth sulfate in water was used as a model system. At cold temperatures (5 °C and 10 °C), the formation of surfactant crystals resulted in extremely high viscosity. The bulk flow behavior of multi-lamellar vesicles (20 °C) and focal conic defects (90 °C) in the lamellar phase was similar. Shear-induced crystallization at temperatures higher than the equilibrium crystallization temperature range resulted in an unusual complex viscosity peak. The effects of processing-relevant parameters including temperature, cooling time, and applied shear were investigated. Knowledge of key low-temperature structure–property-processing relationships for concentrated feedstocks is essential for the sustainable design and manufacturing of surfactant-based consumer products for applications such as cold-water laundry.]]></ab></abstract>
		</profileDesc>
	</teiHeader>
	<text><body xmlns="http://www.tei-c.org/ns/1.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:xlink="http://www.w3.org/1999/xlink">
<div xmlns="http://www.tei-c.org/ns/1.0"><head>INTRODUCTION</head><p>Consumer cleaning products including shampoos and detergents often contain a large amount of water, 1 a critical component during formulation and production that is typically required at all stages of the product's lifecycle. <ref type="bibr">2</ref> Concentrating products like detergents would lead to reduced water usage, waste generation, and transportation costs. <ref type="bibr">[3]</ref><ref type="bibr">[4]</ref><ref type="bibr">[5]</ref> The average 5 L concentrate bottle can be the equivalent of 500 separate ready-to-use spray bottles, potentially saving up to 29.3 kg of plastic trash over a product's lifetime. <ref type="bibr">[6]</ref><ref type="bibr">[7]</ref><ref type="bibr">[8]</ref> If industries only produced and sold concentrated products, it is projected that up to 20% of the world's disposable plastic packaging by weight could be replaced with reusable packaging. <ref type="bibr">9</ref> Making concentrated formulations, however, is more complicated than simply removing as much water as desired. Surfactants are a key component in consumer cleaning products, and the microstructure and rheology of surfactant solutions radically changes with concentration and exposure to forces during processing. <ref type="bibr">[10]</ref><ref type="bibr">[11]</ref><ref type="bibr">[12]</ref><ref type="bibr">[13]</ref> As the surfactant concentration rises, surfactant solutions typically transform from optically isotropic micellar solutions to optically birefringent infinitely ordered <ref type="bibr">14</ref> liquid crystalline assemblies like the hexagonal and lamellar phases. The lyothermotropic phase transitions of surfactant solutions have traditionally been modeled using the critical packing parameter (CPP) <ref type="bibr">15</ref> or spontaneous curvature considerations. <ref type="bibr">14</ref> Micellar solutions are most often Newtonian fluids while the liquid crystalline phases display non-Newtonian flow behavior. <ref type="bibr">[16]</ref><ref type="bibr">[17]</ref><ref type="bibr">[18]</ref><ref type="bibr">[19]</ref> The influence of supramolecular aggregates on viscosity has been extensively studied in literature. <ref type="bibr">20,</ref><ref type="bibr">21</ref> Interestingly, the apparent viscosity of the more concentrated lamellar phase is often lower than that of the less concentrated hexagonal phase. <ref type="bibr">22,</ref><ref type="bibr">23</ref> Due to its microstructure of closely packed cylindrical micelles, the high-viscosity hexagonal phase can be very difficult to process and is generally avoided by industrial formulators. <ref type="bibr">24</ref> In contrast, the lamellar phase has a microstructure composed of stacked parallel surfactant bilayers which flows easily when exposed to shear forces and requires less energy inputs during processing. <ref type="bibr">22,</ref><ref type="bibr">23</ref> The flow behavior of the lamellar phase in surfactant systems has previously been investigated for a range of systems, as summarized in a review by Berni et al. <ref type="bibr">12</ref> However, the rheology of highly concentrated surfactant systems is rarely studied. Pasch et al. reported high yield stress values and shear-dependent orientation phenomenon at 20 &#176;C in several non-ionic systems at high surfactant concentrations (50 -75 wt.% surfactant in water). <ref type="bibr">25</ref> Robles-Vasquez et al. analyzed room temperature dynamic and steady-state rheological responses of an aqueous ionic surfactant solution as a function of concentration to detect structural defects and shear thinning flow behavior. <ref type="bibr">26</ref> Mongondry et al. used small-angle X-ray scattering and rheology of a highly concentrated ternary mixed surfactant system (70 -90 wt.% surfactant in water) to study the flow behavior (at 25 &#176;C and 45 &#176;C) and temperature-dependent crystal to liquid crystal transition. <ref type="bibr">27</ref> Rheology and scattering (SAXS, SANS, SALS) were used to report the dynamic phase behavior and flow induced microstructural transitions for deuterated solutions of a homologous series of non-ionic polyoxyethylene alkyl ether surfactants (C12E3, C12E4, C12E5) at varying surfactant concentrations (40 -60 wt.%) and temperatures (20 -72 &#176;C). <ref type="bibr">[28]</ref><ref type="bibr">[29]</ref><ref type="bibr">[30]</ref> Veronico et al. recently reported a phase diagram for the nonionic surfactant Brij L4 (10 -90 wt.%) -water binary system at 25 &#176;C. <ref type="bibr">31</ref> Caicedo-Casso et al. used steady-state rheometry and ultrasonic velocimetry at 22 &#176;C to investigate the rheological behavior of aqueous sodium lauryl ether sulfate (SLES) solutions as a function of surfactant concentration (20 -70 wt.%). <ref type="bibr">23</ref> At the highest concentration, the presence of significant flow instabilities such as wall slip, shear banding, and plug flow was detected. <ref type="bibr">23,</ref><ref type="bibr">32</ref> An industrial workhorse used extensively in cleaning product formulations, <ref type="bibr">33</ref> aqueous solutions of anionic SLES were chosen because they can easily recreate the raw feedstocks and microstructures often observed in consumer products. Even though SLES is a common component of many commercial products, there has been limited published research on its flow behavior, temperature dependence, and phase evolution in pure systems. Hendrikse et al. used dissipative particle dynamics to simulate the first full-phase diagram of SLES in water at 25 &#176;C33 and molecular dynamics to investigate conformation changes in SLES molecules at various surfactant concentrations. <ref type="bibr">34</ref> Recently, Ferraro et al. investigated the linear rheological behavior of aqueous SLE3S solutions (25 -72 wt.%) from 30 -60 &#176;C and used polarized optical microscopy to report the presence of four phases (micellar (L1), hexagonal (H), cubic (V1) and lamellar (L&#945;)) at room temperature. <ref type="bibr">35</ref> Due to a lack of knowledge of predictable structure-property-processing relationships, instability initiation and evolution are not well understood and it is currently unclear if flow instabilities can help or hinder the formulation and production of concentrated products. <ref type="bibr">36</ref> Additionally, the type of equipment used for processing and the selected operating conditions can change the structure of the product, which can subsequently affect its physical properties like viscosity and ultimately the product's performance, including shelf-life stability and end-use properties.</p><p>Consumer cleaning products begin as raw feedstock and progress through a series of temperature zones before being fully utilized by the end user. However, most prior work on the rheological behavior of highly concentrated surfactant systems, including the aforementioned studies, has primarily focused on room temperature behavior. Thus, there remains a need to study low-temperature flow behavior and characterize the rheological signature of specimens across multiple temperature zones.</p><p>The impacts of temperature variation on the material relationships of concentrated surfactant solutions are particularly important to investigate as industries are encouraging consumers to use products like detergents at lower temperatures, <ref type="bibr">37,</ref><ref type="bibr">38</ref> and studies have shown that simply lowering the wash temperature can have significant economic and environmental benefits. <ref type="bibr">39,</ref><ref type="bibr">40</ref> For instance, in a European study, heating water required for washing accounted for 60% of the carbon footprint associated with laundry. <ref type="bibr">41</ref> By reducing wash temperatures from 40 &#176;C to 30 &#176;C, annual CO2 emissions in Europe can be reduced by almost 4.9 million metric tonscomparable to taking over a million cars off the road. <ref type="bibr">40,</ref><ref type="bibr">42</ref> Naturally, decreasing the temperature to 20 &#176;C would contribute to an even greater reduction in emissions. <ref type="bibr">43,</ref><ref type="bibr">44</ref> Hot water also hastens the breakdown of clothes and promotes microplastic shedding. It is estimated that the conventional laundering of synthetic clothes is responsible for 35% of primary microplastics in world oceans <ref type="bibr">45</ref> and that a single 6 kg domestic laundry cycle can generate up to 700,000 microplastic particles. <ref type="bibr">46,</ref><ref type="bibr">47</ref> High concentrations of these microplastics are regularly found downstream of water treatment facilities <ref type="bibr">48,</ref><ref type="bibr">49</ref> and they can eventually end up in our potable water sources. <ref type="bibr">[50]</ref><ref type="bibr">[51]</ref><ref type="bibr">[52]</ref> Research has shown that laundering at colder temperatures can have outsized positive impacts on microplastic generation and garment longevity. Lant et al. estimate a 30% reduction in microplastic generation for a washing cycle at 15&#176;C for 30 mins, as opposed to 40&#176;C for 85 mins. <ref type="bibr">53</ref> Cotton et al. concluded that in addition to energy savings, reducing laundry time and temperature results in significantly lesser color loss, dye transfer, and microfiber release. <ref type="bibr">54</ref> Temperature variation -and in particular, temperature reductions -are known to drastically change the material relationships of surfactant solutions. As surfactant solutions cool, crystallization can be a source of concern. <ref type="bibr">55</ref> It can be a desired result in some applications, and crystallization methods are frequently employed for separation, purification, and product definition. <ref type="bibr">56</ref> For instance, in the production of sunscreens or moisturizers, crystalline surfactant aggregations act as emulsifiers and improve product stability. <ref type="bibr">57</ref> It can however also have negative impacts such as pipe blockage or the formation of unwanted precipitates in other areas of the product development process. Considering differences in geography, climatic conditions, and water sources, the annual average tap water temperature in the United States of America varies from 7 &#176;C to 22 &#176;C. <ref type="bibr">58,</ref><ref type="bibr">59</ref> Laundering at these colder temperatures with concentrated products will have significant economic and environmental benefits. However, before that is feasible, it is critical to better understand the impacts of crystallization in concentrated surfactant solutions as well as the factors that can help control and tune crystallization.</p><p>The goal of the present experimental study was to develop low-temperature structureproperty-processing-performance relationships for lamellar structured concentrated sodium lauryl ether sulfate (SLES) solutions. The impacts of three key processing variables -temperature, time, and applied shear forces -were investigated. For this study, three temperature zones were evaluated: (1) low temperature, &lt; 15 &#176;C; (2) room temperature, 15 &#176;C &lt; 40 &#176;C; and (3) high temperature, 40 &#176;C to 90 &#176;C. To characterize the flow behavior, shear, and oscillatory rheometry tests were carried out. The microstructure evolution was investigated by polarized optical microscopy and small angle X-ray scattering (SAXS) measurements in conjunction with temperature ramps. The structure-property-processing relationships identified by this work for anionic surfactants can better enable the sustainable design and manufacturing of low-water concentrated cleaning products with the desired performance.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>EXPERIMENTAL SECTION</head></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Materials.</head><p>A concentrated surfactant solution with anionic 70 wt.% sodium lauryl ether sulfate (SLES) in water (STEOL &#174; CS-170 UB, Stepan Company) was used as received. It has an average hydrophobic chain length of 12 carbon (C) atoms, a range of C10-C16, and a single ethoxy group.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Rheometry. All rheometry experiments were performed on an Anton Paar Modular</head><p>Compact Rheometer (MCR 702) configured with a CC10 concentric cylinder fixture (bob diameter = 10.0 mm, bob length = 14.9 mm, measurement gap = 0.422 mm) with Peltier temperature control.</p><p>The concentric cylinder geometry was chosen to minimize water evaporation. The samples were examined using shear and small amplitude oscillatory shear (SAOS) protocols and new specimens from the same surfactant solution batch were used for all experiments. To ensure that all samples had as close to an identical shear history as possible, they were pre-sheared at 5 s -1 for 1 min and rested for 2 min. The magnitude of applied oscillations (strain amplitude, &#947;0 = 0.1% and angular frequency, &#969; = 10 rad/s) was kept within the sample's linear viscoelastic range (LVER) at all studied temperatures (Supplementary Figure <ref type="figure">S1</ref>). Unless specified otherwise, the temperature was changed at 1 &#176;C/min, with an uncertainty of &#177; 0.1 &#176;C.</p><p>A cycle of forward (0.001 -100 s -1 ) and backward (100 -0.001 s -1 ) rate-controlled flow sweep experiments were performed (7 s/point, 20 points per decade). Data from the second forward shear rate ramp is presented (Supplementary Figure <ref type="figure">S2</ref>). To investigate thermal history effects, starting temperatures for oscillatory temperature ramps were varied from 5 &#176;C to 25 &#176;C. Following the heating ramp, samples were cooled from 35 &#176;C back to the starting temperature. To investigate cooling rate effects, quiescent cooling rates varied from 5 &#176;C/min to 0.5 &#176;C/min, and samples were cooled from 20 &#176;C to 10 &#176;C. In a separate thermal aging study, samples were held at constant temperatures for 20 min. The hold temperature varied from 20 &#176;C to 5 &#176;C and the effect of small oscillations was investigated. All experiments were carried out in triplicate.</p><p>Small Angle X-ray Scattering (SAXS). Scattering experiments were performed using an Anton Paar SAXSPOINT 2.0 with a Cu-&#945; source. Samples were placed in the PasteCell N sample holder inside of the thermally controlled TCStage 150. Samples were cooled from 20 &#8451; to 5 &#8451; at 0.5 &#8451;/min with a thermal equilibration time of 3 minutes. At each temperature of interest, three, 1-minute, 2-dimensional frames were recorded using an Eiger X-ray detector at a sample-detector distance of 575 mm. The three scattering frames were averaged, the shadow of the beam stop was removed, and then integrated using a pie radial integration with a radius of 4 mm, an angle of 60&#176;, and aligned with the positive y-axis.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Dynamic Scanning Calorimetry (DSC).</head><p>All experiments were performed using a TA Instruments Q2000 DSC and hermetically sealed Tzero aluminum pans and lids. Temperature was changed at a rate of 1 &#176;C/min and the temperature range compliments the corresponding rheometry experiments.</p><p>Cross-Polarized Optical Microscopy. Cross-polarized images were captured using a Keyence VHX-F series microscope equipped with a Dual-Objective VH-ZST Zoom Lens (magnification range 20x to 2000x) with polarizers. For in-situ temperature ramping, a digitally controlled Linkam Peltier stage was mounted on the microscope. The initial and end temperatures, as well as the ramp speeds, were pre-programmed. Small volumes (~0.2 ml) of sample were carefully put on a glass slide, followed by a cover slip. Data analysis. Experimental data was fitted to theoretical models using Origin 2022 (OriginLab)</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>RESULTS AND DISCUSSION</head><p>Figure <ref type="figure">1</ref> shows the temperature-dependent microstructure evolution and Figure <ref type="figure">2</ref> displays the shear-dependent flow behavior of the concentrated SLES solution maintained at various temperatures. At 5 &#176;C (Figure <ref type="figure">1(a)</ref>), the lamellar L&#945; structure (the liquid crystalline bilayer phase with non-tilted sheets and disordered hydrocarbon chains <ref type="bibr">60</ref> ) was accompanied by opaque crystalline domains that interacted differently with light. <ref type="bibr">61,</ref><ref type="bibr">62</ref> At room temperature (20 &#176;C) and 35&#176;C, oily streaks and multilamellar vesicles (MLVs) could be seen (Figure <ref type="figure">1(b) -(c</ref>)) which were indicative of the L&#945; structure. <ref type="bibr">16</ref> While it is beyond the practical scope of this study, at the highest investigated temperature (90 &#176;C), Focal Conic Defects (FCDs) were seen within the lamellar structure. Figure <ref type="figure">1(d)</ref> shows FCDs that formed between the oily streaks and Figures <ref type="figure">1(e)</ref> and Figures 1(f) are zoomed-in micrographs. Color variation in the polarized optical microscopy images was orientation-dependent and could be analyzed using a Michel-Levy interference chart. <ref type="bibr">63</ref> However, the quantification of interference colors in these relatively thick (0.1 mm) samples was complex and beyond the scope of this study.</p><p>As shown in Figure <ref type="figure">2</ref>, these solutions were shear thinning at all tested temperatures. The low temperature (5 &#176;C and 10 &#176;C) viscosities were significantly higher than the corresponding values at room temperature and high temperature. Interestingly, the flow behavior and viscosity at 20 &#176;C, 35 &#176;C, and 90 &#176;C across the range of shear rates examined were very similar with slight deviations at the highest applied shear rates that were most likely due to shear-induced flow instabilities. <ref type="bibr">23</ref> The presence of significant flow instabilities in rotational experiments <ref type="bibr">23</ref> was a key driving factor for utilizing oscillatory tests in this study. The apparent viscosity of a polymer or surfactant solution typically will increase with decreasing temperature. <ref type="bibr">64</ref> Here, as temperature &gt; 0 &#176;C, the aqueous component of the solutions was not expected to freeze. Thus, this high viscosity and are made up of many concentric layers of lamellar sheets rolled up like an "onion". <ref type="bibr">23,</ref><ref type="bibr">67,</ref><ref type="bibr">68</ref> The effects of different flow conditions (Figure <ref type="figure">2</ref>  cycled to the predetermined temperature with applied oscillations, and rested for 2 min before shearing. Corresponding Herschel-Bulkley 76,77 fits for shear stress vs shear rate at low shear rates are shown in Supplementary Figure S3.</p><p>FCDs are a common structural defect in lamellar phases (Figure <ref type="figure">1(d)-(f)</ref>) and the flow behavior of FCD rich lamellar phase is shown in Figure <ref type="figure">2</ref>. They were identified about a century ago by G. Friedel 78 and have since been the topic of various studies. Apart from the fundamental papers by Friedel, <ref type="bibr">78</ref> Bragg, 79 and Bouligand, 80 studies by Honglawan et al. <ref type="bibr">81</ref> and Guo et al. <ref type="bibr">82</ref> developed methods to control the generation and morphology of FCDs. Currently, it is unclear whether FCDs in concentrated surfactant solutions are induced by temperature, shear, or a combination of the two, and additional investigations are necessary. Gharbi et al. <ref type="bibr">83</ref> and Ma et al. <ref type="bibr">84</ref> focused on leveraging the hierarchical architecture of FCDs for the development of advanced functional materials. Further, Fujii et al. explored relationships between defect structures in lamellar phases <ref type="bibr">85</ref> and observed similarities between the rheological behavior of FCDs and MLVs <ref type="bibr">86</ref> . This is consistent with the similar flow behaviors at 20 &#176;C and 90 &#176;C shown in Figure <ref type="figure">2</ref>.</p><p>To gain insight into the temperature range of crystallization and the resulting high viscosity at low temperatures, SAXS experiments were performed, and the results displayed in Figure <ref type="figure">3</ref>. At 20 &#8451;, the L&#945; lamellar phase was evident following the characteristic 1, 2 &#8230; q* peak pattern. <ref type="bibr">87</ref> At 20 &#8451; the primary q* peak occurred at 1.55 nm -1 and a secondary peak at 3.06 nm -1 . The location of the primary peak corresponded to a lamellar periodicity spacing of 4.08 nm. <ref type="bibr">23</ref> Due to the similarities in scattering patterns, it was clear that the L&#945; phase's temperature region extended from greater than 20 &#176;C down to 10 &#8451;. <ref type="bibr">87</ref> Within this region there was an increase in the lamellar periodicity spacing from the 4.08 nm at 20 &#8451; to 4.20 nm at 10 &#8451;. At 7.5 &#8451; and below the original second peak became unidentifiable from the background. Interestingly, the primary peak exhibited severe convolution at 7.5 &#8451; with a large broadening of the scattering ring before complete separation of the ring at 5 &#8451; indicative of a more complicated morphology. <ref type="bibr">88</ref> SAXS studies on comparable length cetostearyl alcohols have exhibited similar peak broadening at low temperature prior to complex crystallization. <ref type="bibr">88</ref> The possibility of the observed behavior having been ice crystallization can be disregarded due to the location of primary scattering intensity differing from that of both Cubic (Ic) and Hexagonal (Ih) ice crystals, both occurring at ~1.07 nm -1 . 89 However, the lowest q, highest intensity peak varied from 1.547 nm -1 at 20 &#8451; to 1.292 nm -1 at 5 &#8451;, significantly different from the classic ice peak. with an equilibration time of 3 minutes.</p><p>Complementary rheometry, calorimetry and polarized optical imaging experiments were performed (Figure <ref type="figure">4</ref>) to obtain further insight into the crystallization temperature range, microstructure evolution with cooling, and the rheological behavior. As the specimens were cooled, the growth of an opaque crystalline front was observed (Figure <ref type="figure">4</ref>(a) -(f)). In Figure <ref type="figure">4</ref>(h), cooling from 20 &#176;C showed crystallization beginning at 7.8 &#176;C and a slow corresponding complex viscosity increase below 15 &#176;C (Figure <ref type="figure">4</ref>(g)) -The complex viscosity at 5 &#176;C was 67,610 Pa.s. The subsequent heating ramp from 5 &#176;C displayed a melting between 5.7 &#176;C and 11 &#176;C (Figure <ref type="figure">4(h)</ref>) and a gradual decrease in complex viscosity (152 Pa.s) up to 25 &#176;C (Figure <ref type="figure">4</ref>(g)). The range of melting and crystallization temperatures (i.e., broad peaks) in the DSC data was most likely due to the different tail lengths (C10 -C16) in the as-received SLES solution. Interestingly, the complex viscosity changes in Figure <ref type="figure">4</ref>  20 &#176;C and cooled at 0.5 &#176;C/min. For (g) oscillatory temperature ramp and (h) complementary DSC ramp, samples were loaded at 20 &#176;C, cooled to 5 &#176;C, heated to 35 &#176;C, and immediately cooled back to 20 &#176;C. The temperature was changed at 1 &#176;C/min.</p><p>Significant hysteresis in complex viscosity between the cooling ramp from 20 &#176;C to 5 &#176;C and the subsequent heating ramp from 5 &#176;C to 35 &#176;C was also observed in Figure <ref type="figure">4</ref>(g). For example, the viscosity difference between the cooling and heating ramps at 10 &#176;C was nearly 5 x 10 4 Pa.s. This hysteresis indicated the influence of thermal history on complex viscosity. To investigate this further, the solution's complex viscosity response during temperature ramps was determined for different starting temperatures (Figure <ref type="figure">5</ref>).</p><p>The complex viscosity was indeed influenced by the sample's thermal history. For example, the complex viscosity at 25 &#176;C for a ramp starting at 5 &#176;C was 170 Pa.s -significantly higher than the corresponding complex viscosity values for ramps starting at higher temperatures: 90 Pa.s and 52 Pa.s for 10 &#176;C and 15 &#176;C starting temperatures, respectively. Figure 5 also displays an unusual complex viscosity peak observed for the ramp started at 10 &#176;C. When the ramp was started at 15 &#176;C, 20 &#176;C, and 25 &#176;C, the complex viscosity of these solutions was relatively low and nearly constant, consistent with flow curve behavior in Figure <ref type="figure">2</ref>. pre-determined temperature, pre-sheared and rested before the temperature sweep. The samples were heated up to 35 &#176;C and immediately cooled back to the temperature of interest. The G&#8242; and G&#8243; profiles for ramp starting at 10 &#176;C are shown in Supplementary Figure S4.</p><p>To better understand how the surfactant solutions were impacted by thermal processing history, complementary DSC experiments were performed at starting temperatures of 5 &#176;C and 10 &#176;C (Figure <ref type="figure">6</ref>(a) and (b)). When the temperature ramp was started at 5 &#176;C (Figure <ref type="figure">6</ref>(a)), crystallization and melting peaks were observed between 5 &#176;C and 9 &#176;C, similar to the range in Figure <ref type="figure">4</ref>(h). Heating the specimen from 5 &#176;C resulted in a melting transition peak at 7.8 &#176;C and caused a slow decrease in complex viscosity until to 20 &#176;C. Cooling resulted in crystallization beginning at 6.7 &#176;C and a sharp increase in complex viscosity which was only partially captured in the experiment due to temperature range limitations. In an interesting contrast, for a temperature ramp beginning at 10 &#176;C (Figure <ref type="figure">6</ref>(b)), there were no significant endothermic or exothermic peaks visible in the DSC data that corresponded to the complex viscosity peak observed from 10 &#176;C to 25 &#176;C. different starting temperatures: (a) 5 &#176;C, and (b) 10 &#176;C. The temperature ramp rate was 1&#176;C/min.</p><p>The complex viscosity peak that results during heating from 10 &#176;C (Figure <ref type="figure">5</ref>; Figure <ref type="figure">6b</ref>) was hypothesized to be the result of applied oscillations promoting crystallization of the surfactant molecules at temperatures higher than the crystallization range determined by static DSC experiments (of 5 to 9 &#176;C, Figure <ref type="figure">6(a)</ref>). During additional oscillation experiments, this peak was confirmed to exist in a narrow temperature range from 10 &#176;C to 12.5 &#176;C (Supplementary Figure <ref type="figure">S5</ref>). This small temperature range was consistent with previous findings for some other dilute surfactant systems. <ref type="bibr">90,</ref><ref type="bibr">91</ref> There is a substantial body of research on the effect of applied shear forces on crystallization below the crystallization point in a wide range of complex fluids <ref type="bibr">92,</ref><ref type="bibr">93</ref> such as polymer melts, <ref type="bibr">94,</ref><ref type="bibr">95</ref> colloidal glasses, <ref type="bibr">96,</ref><ref type="bibr">97</ref> and multi-component surfactant systems. <ref type="bibr">98,</ref><ref type="bibr">99</ref> However, the literature on shear-driven crystallization in surfactant solutions above the crystallization temperature range is limited. <ref type="bibr">100,</ref><ref type="bibr">101</ref> Rathee et al. reported a reversible shear-induced crystallization in a cationic-anionic mixed ternary surfactant system. Rheo-optical tests were used to show that the pre-ordering of the isotropic bilayer mesophase 102 in a shear-induced lamellar phase was a precursor to the nucleation of the crystalline phase. <ref type="bibr">100</ref> To further investigate the hypothesized occurrence of shear-induced crystallization, the data displayed in Figure <ref type="figure">7</ref>(a) depicts the specific influence of small applied oscillations on the solution's complex viscosity. Figure <ref type="figure">7</ref>(a) includes data from two different specimens: both cooled from 20 &#176;C to 10 &#176;C at 1 &#176;C/min, one with applied oscillations during cooling (blue circles) and the other was cooled without any applied oscillations (red triangle). At 10 &#176;C, the complex viscosity for the sample oscillated during cooling was 3300 Pa.s (point A) while the complex viscosity for the sample at rest during cooling was 73 Pa.s (point B). Inspired by the theory of the metastable zone in the classical work of Nyvlt et al. <ref type="bibr">103</ref> , specimens were cooled at from 20 &#176;C to 5 &#176;C at cooling rates varying from 0.1 &#176;C/min to 3 &#176;C/min with small oscillations (Figure <ref type="figure">7(b)</ref>).</p><p>Rates higher than 3 &#176;C/min, were too fast for the rheometer and it struggled to achieve a steady thermal state. At the slowest rates (0.1 &#176;C/min and 0.25 &#176;C/min), the effects of shear on nucleation and growth were pronounced and visible significantly above the equilibrium crystallization temperature. The magnitude of complex viscosity at 10 &#176;C for the higher cooling rates (2.5 &#176;C/min and 3 &#176;C/min) was comparable to point B in Figure <ref type="figure">7</ref>(a) where no oscillations were applied while cooling. Thus, as the solution was cooled, while small oscillations within the linear viscoelastic regime affected the microstructure and acted to promote crystallization, this effect was dependent on the cooling rate. loaded at 20 &#176;C and cooled at 1 &#176;C/min and (b) Effects of varying cooling rates as specimens are cooled from 20 &#176;C to 5 &#176;C with small oscillations Figure 8 displays the effect of applied oscillations at a constant temperature. To decouple the effect of applied oscillations from temperature changes, specimens were loaded at 20 &#176;C, cooled to a specified temperature and held there for 20 mins. At each corresponding temperature, Figures 4 (a) -(f) are indicative of the microstructure before pre-shear, rest, and hold in Figure 8.</p><p>On an observable timescale, the complex viscosity plateaus reflect the maximum possible shearinduced crystallization, and the slopes of the curves before the plateau indicate the amount of time required to accomplish maximum crystallization. The time to reach the plateau is a function of temperature (samples at 12.5 &#176;C and 15 &#176;C have not reached a plateau after 1200 seconds).</p><p>Applied oscillations enhanced crystallization at temperatures higher than the crystallization temperature range by primarily enhancing growth. The complex viscosity at t = 0 can be used as an indicator of the effect of the nucleation process. The small applied oscillations increase the rate of mass transfer of surfactant molecules to the crystal surface and promote growth. They can also cause collision breeding, further enhancing growth. <ref type="bibr">103,</ref><ref type="bibr">104</ref> This was seen clearly in the data sets at 7.5 &#176;C, 10 &#176;C, and 12.5 &#176;C in Figure <ref type="figure">8</ref>. The increased complex viscosity was an indicator of increased shear-induced crystallization with time. At 5 &#176;C, the temperature was at the lower end of the crystallization temperature range and the effect of oscillations was not as pronounced. specimens were cooled from 20 &#176;C to a pre-determined temperature (e.g., 5 &#176;C) at 0.5 &#176;C/min without oscillation, pre-sheared, and rested before the aging experiment. G&#8242; and G&#8243; profiles at selected temperatures are shown in Supplementary Figure S6. the samples LVER at 20 &#176;C, 10 &#176;C and 5 &#176;C. The specimens were cooled from 20 &#176;C to a predetermined temperature (e.g., 5 &#176;C) at 0.5 &#176;C/min without oscillation, pre-sheared, and rested before the aging experiment.</p><p>As shown in Figures <ref type="figure">9(a</ref>) and (c), strains within the LVER (0.01% and 0.1%) had the same influence on complex viscosity at temperatures where the solution was at equilibrium (20 &#176;C -L&#945; lamellar phase and 5 &#176;C -crystalline phase). The response is different at 10 &#176;C because the surfactant solution is in a non-equilibrium transient state and small oscillations induce crystallization. The frequency dependence at all temperatures (Figure <ref type="figure">9(b)</ref>, <ref type="figure">(d)</ref> and <ref type="figure">(f)</ref>) is rather more straightforward. The complex viscosity of the solution decreased as the frequency was increased, with a lower baseline complex viscosity at 20 &#176;C. The shear-thinning nature of concentrated lamellar-structured SLES solution (Figure <ref type="figure">2</ref>) was posited to account for this behavior.</p><p>Having investigated the impacts of crystallization temperature range and applied oscillations, the effects of cooling time on rheological behavior were explored. The specimens were cooled from 20 &#176;C to 10 &#176;C at different rates without oscillations and the quiescent cooling time was found to impact the viscosity and microstructure. 10 &#176;C was selected because a complex viscosity peak was previously observed when oscillations were initiated at this temperature (Figure <ref type="figure">5</ref>) and because it was higher than the crystallization range determined by static DSC experiments (Figure <ref type="figure">6</ref>). As shown in Figure <ref type="figure">10</ref>, the slowest cooling rate (0.5 &#176;C/min) resulted in a significant complex viscosity peak during the subsequent heating ramp with oscillations. The slower quiescent cooling rate essentially promoted the development of more crystal nucleation sites which subsequent oscillations then helped to grow. The faster cooling rates provided less time for ordered arrangements and nucleation sites to form at 10 &#176;C and resulted in low complex viscosity peaks. <ref type="bibr">103</ref> The sigmoidal nature of isothermal shear-induced crystallization (Figure <ref type="figure">8</ref>) was analyzed using a modified five-parameter logistic (5PL) model (Eqn. 1). <ref type="bibr">105</ref> Typically used to model immunoassays such as enzyme-linked immunosorbent assays (ELISA) and radioimmunoassay (RIA), the 5PL non-linear regression model builds on the four-parameter logistic model and incorporates an additional parameter (s) to characterize the curve asymmetry. where A is the pre-exponential or Arrhenius factor, Ea is the activation energy, R is the universal gas constant and T is the absolute temperature. In the presence of small oscillations, the activation energy for crystal growth was estimated to be 107 kJ/mol. This value is much lower than that anticipated for micellar aqueous solutions of sodium dodecyl sulfate crystallized without shear (Ea &#8776; 218 &#177; 46 kJ/mol). <ref type="bibr">109</ref> Additional studies evaluating the effects of additives in combination with shear are underway. oscillations. Fit parameters are tabulated in the Supplementary section. Originally developed to model the extent of phase transitions, the Johnson-Mehl-Avrami-Kolmogorov (JMAK) 110-112 formalization, commonly referred to as the Avrami equation (Eqn. 3), was used to model the kinetics of isothermal crystallization (Figure 11(b)). 113 67 ,-67 &gt; --|# * | |# * | $'+ ?0 "#$%#&amp;'()&amp;# = 67 @ + A 67 % Eqn. 3</p><p>where n and K respectively are the intercept and slope. At higher hold temperatures, the crystallization process entered its growth phase after longer inception periods. In addition, it was observed that as hold temperature increased, the gradients were steeper, suggesting a faster approach to the plateau. At 7.5 &#176;C, 10 &#176;C and 12.5 &#176;C, the curves followed classical linear Avrami behavior over a limited timeframe.</p><p>Beyond laundry and personal care, controlled tailoring of microstructures in lamellar phases can have significant applications. For example, polymerizing lamellar phases, such as lyotropic liquid crystal (LLC) templating, has been shown to improve the mechanical and thermal properties of lamellar phases.</p><p>114 Clapper et al. ordered macromolecules within a lamellar phase and developed nanostructured cross-linked biodegradable hydrogels with enhanced swelling and permeability. 115 Bandegi et al. tuned the degree of LLC crystallinity to create ion gels with enhanced mechanical strength and ionic conductivity. 116 Qavi et al. created antimicrobial membranes and ultrafiltration (UF) membranes with superior contamination resistance and water permeability compared to commercial UF membranes. 117</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>CONCLUSIONS</head><p>In this experimental study, low-temperature structure-property relationships for concentrated SLES solutions were developed. At all temperatures and shear rates tested, these solutions were shear thinning (Figure <ref type="figure">2</ref>). The presence of crystals at low temperatures resulted in extremely high viscosity. Despite the presence of FCDs at high temperatures, the viscosity and rheological behavior at 20 &#176;C and 90 &#176;C were very similar. Room temperature viscosity was influenced by the thermal history (Figure <ref type="figure">5</ref>) and a crystallization temperature range was identified (Figure <ref type="figure">3</ref> and Figure <ref type="figure">6</ref>). The peculiar viscosity peaks with increasing temperature were ascribed to a combination of cooling time, shear-induced crystallization, and temperature range.</p><p>Temperature (Figure <ref type="figure">5</ref>), applied oscillations and cooling rates (Figure <ref type="figure">7</ref> -9), and quiescent cooling time (Figure <ref type="figure">10</ref>) are critical processing variables that directly impact the microstructure and rheological properties of concentrated surfactant solutions. Isothermal shear induced crystallization was modelled using a five-parameter logistic model and a relevant activation energy was estimated.</p><p>The presence of crystals and the resulting high viscosity at low temperatures can make cold water laundry problematic. The concentrated product may take longer to dissolve in the washer, lengthening the load cycle and thereby increasing energy use. Oscillations during processing or transit can affect the appearance and performance of finished products. On one hand, oscillating pipelines can cause issues like uneven mixing, blockage, and drippage. On the other hand, if a precise amount of crystallization in a finished product is desired, it could be induced by vibrating the final package for a specific time.</p><p>The workflows developed herein for neat SLES solutions can be applied to study the phase evolution, microstructure, and flow behavior of systems with industrially relevant additives.</p><p>Research into the impact of temperature, shear, and additives on the formation of MLVs and FCDs </p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>ELECTRONIC SUPPLEMENTARY INFORMATION</head><p>Amplitude sweeps at different temperatures, shear stress vs shear rate data at low shear rates with Herschel-Bulkley fits, temperature range for complex viscosity peak, G&#8242; and G&#8243; for selected temperature ramps, cross-polarized micrographs showing the effect of pre-oscillation quiescent cooling rate, 5 Parameter Logistic and Avrami Fits for isothermal shear-induced crystallization.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>CONFLICTS OF INTEREST</head><p>There are no conflicts of interest to declare. high ramp is presented. (b) Time sweep experiments at constant shear rate after low-to-high and high-to-low shear rate sweep. The average of viscosities at each shear rate is shown in (a). Figures (a) and (b) reveal that the procedure used to create the continuous shear rate ramp flow curve achieves steady state conditions. The data is presented as shear stress vs shear rate. Samples were loaded into the fixture at 20 &#176;C, cycled to the predetermined temperature with applied oscillations, and rested for 2 min before shearing. The Herschel -Bulkley model 1,2 is as follows: B = B H + C * (E&#775;:), where &#963; is the shear stress, &#963;y is the yield stress, K is the consistency coefficient, E&#775; is the shear rate and n is the flow behavior index which varies as follows. If n=1, the Herschel -Bulkley model is equivalent to the Bingham Plastic model. If &#963;y = 0, it is equivalent to the Ostwald-de Waele power law model. The yield stresses, consistency coefficients and flow index are summarized in the table below. Herschel -Bulkley Fits: Temperature (&#176;C) Yield Stress (&#963;y) (Pa) Consistency coefficient (K) Flow behavior index (n) R-square (COD)</p><p>5 415 1709 0.268 0.999 10 9.80 123 0.183 0.999 20 5.86 3.83 0.462 0.996 35 4.15 8.54 0.600 0.993 90 3.65 8.26 0.592 0.999</p></div><note xmlns="http://www.tei-c.org/ns/1.0" place="foot" xml:id="foot_0"><p>Procter &amp; Gamble, Leading detergent brands to European consumers: Wash cold, https://www.politico.eu/sponsored-content/leading-detergent-brands-to-europeanconsumers-wash-cold/, (accessed June 4, 2023).</p></note>
		</body>
		</text>
</TEI>
