Skip to main content

BRIEF RESEARCH REPORT article

Front. Phys., 15 September 2022
Sec. Soft Matter Physics
This article is part of the Research Topic Probing Out-of-Equilibrium Soft Matter View all 10 articles

Shear-induced vorticity aligned flocs in a temperature responsive colloid-polymer mixture

Ryle RelRyle RelDennis TerwilligerDennis TerwilligerRyan McGorty
Ryan McGorty*
  • Department of Physics and Biophysics, University of San Diego, San Diego, CA, United States

Shear driven patterning is seen in many soft matter systems. We use rheology and optical microscopy to probe the structures formed when we shear a colloid-polymer mixture containing temperature-sensitive microgel particles. By increasing the temperature, we can increase the particle attraction and transition from liquid-like to gel-like behavior. And by applying shear flow to the sample as the temperature and, hence, state of the system changes, we can affect the morphology of mesoscopic colloidal clusters. We can produce gels comprised of fibrous, elongated colloid-dense clusters, or we can form more isotropic clusters. The rheology is measured and shear-induced flocculation observed for colloid-polymer systems with different cluster morphologies. At shear rates high enough to produce elongated clusters but low enough to not break clusters apart, we observe log-like flocs that are aligned with the vorticity direction and roll between the parallel plates of our rheometer.

Introduction

Colloidal suspensions find wide use for their rheological properties and their ability to assemble into complex structures [1]. With a multitude of routes for colloidal synthesis, particle sizes large enough for detection with optical methods, and tunable interparticle interactions, the structure and mechanical properties of colloidal assemblies can be readily observed, measured, and modified. Therefore, the assembly of colloidal particles has been probed to, e.g., understand the kinetics of crystal nucleation and growth [2], the mechanisms of gelation [3], and the fabrication of capsules [4]. Likewise, the rheology of colloidal suspensions has been used to, e.g., develop theoretical models for non-linear flow responses [5] and tailor the properties of 3D printer ink [6]. Whether to test fundamental models of matter or to drive applied research towards product development, these two aspects of colloids—elementary building blocks of complex assemblies and rheological modifiers—are often quite entangled [7]. Flow can assist in the assembly of higher order structures, such as colloidal crystals [8, 9], and the details of a sample’s microstructure can assist in explaining rheological properties, like shear thinning or thickening in dense suspensions [10] or how a gel yields and flows [11, 12].

In many colloidal systems, researchers have observed anisotropic shear-induced structures which align along the vorticity direction. These structures are of great interest due to their implications in determining the (potentially directional dependent) rheology of complex fluids and in revealing routes to fabricate anisotropic materials. Such anisotropic shear-induced structures have been studied experimentally through orthogonal superposition rheology [13, 14], rheo-scattering experiments [15, 16], and real space rheo-optical methods [17, 18]. Of the plethora of anisotropic shear-induced structures observed, a subset feature log-rolling flocs, whose formation strongly depends on the degree of confinement. These vorticity-aligned log-rolling flocculated structures have been observed in emulsions [19]; particle-loaded polymer blends [20]; carbon nanotube suspensions [21, 22]; carbon black suspensions [2325]; mud samples [26]; and other suspensions of attractive particles [27, 28]. A complete understanding of such log-rolling floc formation is crucial for applications where complex fluids might flow through confined geometries (e.g., 3D printing) and for scalable methods of producing anisotropic materials (e.g., films with anisotropic conductivity [29, 30]).

A large advance in the understanding of this type of shear-induced patterning came from the experiments and simulations performed by Varga et al. in 2019 [27]. The authors demonstrated that density fluctuations within a sheared and confined suspension of attractive particles can set up eddies which, due to hydrodynamic coupling between the particle aggregates and the boundaries, promote the stable formation of log-rolling flocs. Further, the authors showed that if the shear forces are sufficiently large, flocs will break apart and stable vorticity-aligned log-rolling structures will not form. Whether the shear force is too large or not depends on the interparticle attraction, a fact which motivates the use of the Mason number, Mn, the ratio of the shear force on a particle to the interparticle force at contact. Stable logs were only observed below a critical Mason number, Mnc, with Mnc decreasing with increasing distance between the confining boundaries.

In this work, we explore the log-rolling floc formation in a colloid-polymer system containing thermoresponsive microgel colloidal particles using an optical microscope attachment on a rheometer with a parallel-plate geometry. Using poly(N-isopropylacrylamide) (pNIPAM) microgel particles allows us to vary the interparticle attraction in situ by changing the temperature, a property used in many prior investigations of colloidal phase behavior and rheology [3136].

We find that, for our system, the morphology of colloidal clusters—assemblies of particles on a mesoscopic scale between that of the particle size (∼1 μm) and system size (>150 μm)—is an important factor influencing log formation. With our thermoresponsive colloid-polymer system, we can create colloidal gels comprised of elongated, highly anisotropic clusters or more isotropic clusters. Gels consisting of more elongated clusters tend to more readily form log-rolling flocs. Our work highlights the importance of investigating the impact of this mesoscopic cluster scale on shear-induced pattern formation. In many previous studies of colloidal gels, a cluster-centric view has proven fruitful. For example, past work has shown that colloidal clusters may determine the length scale for mechanical energy storage in colloidal gels [37] with inter-cluster bonds or cluster-cluster connections governing the elasticity [3840]. Here, we show that a perspective which focuses on the mesoscopic cluster scale may also help elucidate macroscopic pattern formation in attractive colloidal suspensions.

Materials and methods

Materials

We synthesized pNIPAM microgel particles in a manner similar to previous syntheses by our group [41, 42] and following the guidance of previously published protocols [43]. We use ammonium persulfate (APS, purity 98%), N-isopropylacrylamide (NIPAM, 97%), and N,N′- methylenebisacrylamide (BIS, 99%), all used as received from Sigma Aldrich without purification. We synthesize particles in a 250 mL three necked round bottom flask heated to 80°C and under nitrogen atmosphere. In 200 mL of deionised water, we add 3.38 g of NIPAM and 0.15 g of BIS. Once dissolved, this solution is transferred to the three necked flask. We then add 0.11 g of APS dissolved in 3 ml of deionised water. We let the reaction run for 4 h with constant stirring with a magnetic stir bar and continuous bubbling of nitrogen through the solution. We concentrate the pNIPAM particle solution through centrifugation. We centrifuge the suspension at 30,000 g for 90 min. The supernatant is removed and fresh deionized water is added.

We estimate the pNIPAM particle volume fraction by measuring the viscosity of different particle dilutions using a Cannon-Ubbelohde dilution viscometer. According to the Einstein-Bachelor relation, the viscosity of a suspension of spherical particles relative to that of the suspending medium is η/η0=1+2.5ϕ where ϕ is the effective particle volume fraction [44]. We found the volume fraction of our stock solution of concentrated pNIPAM particles to be ϕ 0.6.

We determine the size of our pNIPAM particles around the lower critical solution temperature (LCST) using dynamic light scattering (DLS) measurements taken with a Malvern Zetasizer NanoZS using a 632.8 nm laser with a detection angle of 173°. Above the LCST at 34°C, the diameter is 480 ± 30 nm (mean and standard deviation of four runs on the Zetasizer) and we used a particle concentration of ϕ< 0.01 such that we do not see signs of aggregation in our DLS measurements. Below the LCST at 29°C, the diameter is 840 ± 80 nm. At lower temperatures, the size measurements from the Zetasizer are less consistent perhaps due to large polydispersity though follow up measurements are necessary to quantify this. However, from optical microscopy of the particles at room temperature, we believe that their size at room temperature is approximately 800 nm.

We use xanthan polymer, Ticaxan® Xanthan VI from TIC Gums, to act as a depletant. Xanthan has been used in previous colloid-polymer samples to generate fluid-fluid phase separation or gelation [41, 42, 45]. As reported by the manufacturer, the molecular weight is in the range of 4–12 × 106 g mol−1. We prepare a stock solution of xanthan with a concentration of 0.2% by weight with 0.1 M NaCl (Sigma Aldrich, >99%) and 2 mM NaN3 (Sigma Aldrich, >99.5%).

We made samples to investigate shear-induced pattern formation by combining in equal parts our concentrated pNIPAM stock solution and the 0.2% xanthan solution. This results in a sample with a pNIPAM microgel particle volume fraction of 0.3.

Rheology

Rheological measurements are performed on a Discovery Hybrid Rheometer (DHR-3) from TA Instruments. For all data included here, we use a parallel plate geometry with a 40-mm-diameter upper plate and a bottom 1-mm-thick glass plate. The upper plate is made of stainless steel with a mirror finish to allow for better image quality when using the attached microscope. We use TA Instruments’ Upper Heated Plate accessory to vary the sample temperature between 23°C and 34°C.

We load ∼210 μL of our colloid-polymer mixture onto the rheometer by pipetting the solution onto the center of the bottom glass plate. The upper plate is slowly lowered while rotating at 0.2 rad/s. The upper plate is lowered until the sample completely fills the gap which usually results in a gap height of between 150 and 160 μm. To prevent sample evaporation, we seal the edge of the sample with mineral oil.

Before experimental runs, we try to erase the sample’s shear and loading history by performing a conditioning step where the sample is sheared at 1,000 s−1 for 120 s, sheared at −1,000 s−1 for 120 s, and then allowed to equilibrate at rest for 180 s. This preshearing procedure is similar to ones used in other studies of shear-induced structuring [24, 25, 27, 28]. For the majority of experiments we discuss here, we then proceed to heat the sample from room temperature to 34°C at a rate of 0.5°C/min while steadily shearing the sample at a rate, γ˙T, between 0 and 15 s−1. After the sample reaches 34°C, we usually measure the frequency dependent viscoelastic moduli. We logarithmically sweep the frequency from 100 to 0.1 rad/s with five points per decade and with a strain amplitude of 3%. This is followed by a constant shear rate step where the sample is sheared at a rate γ˙hold for 30 min and kept at 34°C. Typically, we use the same shear rate during this step as was used during the heating stage, i.e., γ˙T=γ˙hold. For our samples, if vorticity-aligned log-rolling flocs form, it is at this step in which they appear. This step with constant shear rate is then followed by another oscillatory test using the same parameters as the previous one. The sample is then returned to 23°C at a rate of 2°C/min while being sheared with rate γ˙T. The steps described above are then repeated. This whole procedure where the sample goes from room temperature to 34°C, is held at 34°C for 30 min at a constant shear rate, and is then brought back to 23°C typically takes ∼75 min. This is usually repeated up to ∼10 times with different shear rates over the course of a day or overnight. Between each run we perform the conditioning step described above (returning the sample to room temperature, shearing at 1,000 s−1 for 120 s, and then at −1,000 s−1 for 120 s) which allows us to conduct multiple runs of the same sample without noticeable effects from the previous run. See the Supplementary Material for a table of the parameters used in our rheological procedures.

Our rheometer is equipped to counter-rotate the bottom glass plate which results in a stagnation plane between the parallel plates. For most of our experimental runs, we do not enable this counterrotation. However, we occasionally use this feature to observe the formation of vorticity-aligned flocs without having such flocs leave the imaging field of view. When this counterrotation is used, we indicate so in the caption.

Microscopy

We use the modular microscopy accessory from TA Instruments to image the shear-induced flocculation. For images included here, we use Nikon objectives with either 2×, 10×, or 40× magnification (2× 0.1 NA Plan Apo; 10× 0.3 NA Plan Fluor; 40× 0.45 NA ELWD SPlan Fluor). Epi-illumination is provided by a 470-nm LED (M470L3, Thorlabs). Images are recorded on an Allied Vision Manta G-033 CCD camera with a resolution of 656 × 492 pixels at a maximum frame rate of 88 frames per second. We use the software Micro-Manager to record image sequences [46].

For macroscopic imaging where we can capture one-quarter of the 40-mm-diameter geometry, we use a macro lens attached to a 1,920 × 1,200 pixel CCD camera (Imaging Source, DMK 33UX174). Images are typically recorded at five frames per second. We illuminate the sample with a dual gooseneck LED light (AmScope, LED-14M).

Results and discussion

We prepare aqueous solutions of pNIPAM microgel particles and xanthan which exhibit temperature-dependent structural and mechanical properties. This colloid-polymer mixture shows gel-like behavior at room temperature as observed with optical microscopy (Figure 1A) and rheology (Figure 1B; Supplementary Figure S1). This behavior is due to the depletion attraction between the pNIPAM particles mediated by the xanthan polymers.

FIGURE 1
www.frontiersin.org

FIGURE 1. (A) At room temperature, the colloid-polymer mixture forms a weak gel as observed with a brightfield microscope using a 60× objective lens. (B) As measured with a DHR-3 instrument using a parallel-plate geometry, the storage (G, solid symbols) and loss (G, open symbols) moduli vary with temperature. The moduli were obtained using a frequency of ω=10 rad/s with an amplitude of γ0=3%. (C) The DHR-3 used for collecting rheology and rheo-microscopy data is shown. The instrument is outfitted with a temperature-controlled upper plate and a bottom glass plate under which a microscope is placed. (D) A view from underneath the parallel-plate geometry shows a 2× microscope objective placed below the bottom glass plate. The sample shows shear-induced patterning as thin white structures oriented radially outward (in the vorticity direction) can be observed. (E) An image of the colloid-polymer mixture in the parallel-plate geometry after the sample was heated from 23°C to 34°C at a rate of 0.5°C/min. During this temperature ramp, the sample was sheared at a constant rate of γ˙T=5s1. This image was taken immediately after the sample first reached 34°C. (F) For the same sample shown in (E), we continue applying a constant shear at a rate γ˙hold=5s1 while the temperature was held at 34°C. This image shows vorticity-aligned flocs have formed at t=540s into this step of constant applied shear rate and constant temperature. The vertical axis of the image is along the flow direction and the horizontal axis is along the vorticity direction. (G) During this period of constant shear described in (F), we measure the stress over time. Increases in the stress coincide with the formation of vorticity-aligned log-like flocs as seen by the uptick in stress at around 500 s which coincides with when the logs seen in (F) begin to form.

When sheared at room temperature (23°C) between parallel plates (40 mm diameter, ∼150 μm gap size), we observe colloid-rich domains aligning along the flow direction. At room temperature, we never observe vorticity-aligned flocs for any of the shear rates we used, 0.5 s−1–15 s−1. This may be due to the fact that even at the lowest shear rate in this range, the shear forces acting on a particle are too large compared to the interparticle depletion attraction for stable vorticity-aligned logs to form.

We next explore if shear-induced patterning is seen at temperatures above pNIPAM’s lower critical solution temperature (LCST) of about 32°C–33°C. Above this temperature, pNIPAM particles become more hydrophobic with an increased interparticle attraction. As the temperature of our sample increases from room temperature to 34°C, the interparticle attraction varies in a non-trivial way. As the temperature increases while remaining below the LCST, the pNIPAM particles decrease in size. This reduces the volume fraction of pNIPAM particles and reduces the strength of the depletion attraction. A quantitative picture of how the pNIPAM interparticle attraction varies as a function of temperature in the presence of a depletant is beyond the scope of this work (though the rich phenomena that similar pNIPAM-containing systems exhibit have been investigated previously by our lab [41] and by others [47]). However, we do see the effect of this change in interparticle attraction in Figure 1B. The storage and loss moduli, G and G, are measured as a function of temperature at a frequency of ω=10 rad/s with a strain amplitude of γ= 3%. We observe that G and G decrease with temperature and that while G>G for lower temperatures, at around 31°C we observe more liquid-like behavior with G<G. We conjecture that this decrease in G with increasing temperature and the transition from gel-like to liquid-like behavior is due to the decreasing volume fraction of pNIPAM particles and decreasing depletion attraction. Above 33°C, we observe a stiffer colloidal gel with G returning to being larger than G. The frequency dependent moduli for our colloid-polymer mixture at both 23°C and 34°C is shown in Supplementary Figure S1. While the pNIPAM particle volume fraction is larger at room temperature, the interparticle attraction is greater above the LCST which results in the elastic moduli being ∼2× larger at 34°C than at 23°C.

At this higher temperature of 34°C, we do observe, with our rheo-optical setup (Figures 1C–F) the formation of vorticity-aligned flocs similar to what has been observed previously in other attractive suspensions or soft matter systems (see also Supplementary Videos S1–S3). The formation of such flocs coincides with an increase in the shear stress along with more variability in the stress as seen in Figure 1G. This increase in stress or viscosity associated with log-rolling floc formation is consistent with previous work on systems of colloidal rods [28], carbon black dispersions [24], and nanotube suspensions [21]. If, after the log-rolling flocs form, we leave the sample undisturbed at 34°C, then the logs seem to remain intact for at least several minutes but we have not investigated this thoroughly. However, if we bring the temperature back down to 23 C, then the pNIPAM particles will disperse away from the flocs.

Interestingly, we only observe the formation of the kind of vorticity-aligned flocs seen in Figures 1D,F in certain circumstances. If, after following the sample conditioning steps described in the Methods section (i.e., applying a shear of 1,000 s−1 for 120 s, −1,000 s−1 for 120 s, and allowing the sample to equilibrate at rest for 180 s), we increase the temperature to 34°C while the sample is at rest (i.e., no applied shear) and then, with the sample temperature held at 34°C, apply a steady shear of anywhere between 0.5 and 15 s−1 we do not observe the formation of flocs. That is, if γ˙T = 0, then we observe no vorticity-aligned flocs for any value of γ˙hold used. Instead, we observe isotropic aggregates of particles in the sheared sample which are smaller than the gap size and which do not appear to flocculate together or change in size or shape over time, at least as observed with our imaging setup. Further, we observe that the measured shear stress is steady with time and does not exhibit large fluctuations as we shear the sample at a constant rate. This is another indicator of little or no change in the sample’s structure as it is sheared.

Instead of heating the sample to 34°C while the sample is quiescent, we also try applying a steady shear to the sample as we heat it from 23°C to 34°C at 0.5°C/min. That is, we set γ˙T to a value greater than zero. In Figure 2A, we show the stress as a function of temperature during this temperature ramp for a few different shear rates. For a given shear rate γ˙T, stress initially decreases with temperature. This is likely due to the decreasing volume fraction of pNIPAM particles as the particles deswell and the concomitant reduction of the depletion attraction. The stress is a minimum at ∼33°C, approximately pNIPAM’s LCST. From 33°C to 34°C, the stress increases, likely due to the stronger interparticle attraction above the LCST. The mesoscopic structures of our samples at 34°C are distinctly different depending on the value of this shear rate, γ˙T. After the sample’s temperature has reached 34°C, colloidal clusters are slightly elliptical in shape for lower shear rates (e.g., Figure 2B where γ˙T=1s1). But for higher γ˙T, clusters take on a fibrous quality where the clusters are highly elongated (e.g., Figure 2C where γ˙T=7s1). See Supplementary Figure S2 for more such images.

FIGURE 2
www.frontiersin.org

FIGURE 2. Shearing the sample as the temperature increases produces distinctly different cluster morphologies. (A) Samples are sheared at a constant rate of γ˙T while the temperature increases from 23°C to 34°C at a rate of 0.5°C/min. The stress is measured during this temperature ramp. For all values of the shear rate γ˙T used, we observe that the stress decreases as the temperature goes from room temperature to ∼33°C. This is likely due to the decreasing volume fraction of the pNIPAM microgel particles as they decrease in size with increasing temperature. Above ∼33°C, the stress increases with temperature as the increased hydrophobicity of the pNIPAM particles leads to stronger interparticle attraction. (B) Image of the sample immediately after this temperature ramp with a shear rate of γ˙T=1s1 taken with a 10× objective lens. (C) Image of the sample immediately after the temperature ramp with a shear rate of γ˙T=7s1. In comparison to (B), the pNIPAM particles have aggregated into more fibrous and higher aspect ratio clusters.

We believe that the shear-rate-dependent shape of clusters can be understood considering the large body of work on how fluid droplets deform in shear flow. Though our colloid-polymer system exhibits gel-like behavior for most of the temperature range, our rheological (Figure 1B) and imaging (Supplementary Video S4) data shows liquid-like behavior at temperatures just below the LCST where colloid-dense and colloid-poor fluids likely coexist. Studies of flow-induced deformation of droplets and soft materials date back at least a century with the seminal work of G. I. Taylor [48, 49]. More recent experimental work on systems more resembling our own demonstrate that droplets sheared in confined geometries take on string-like elongated shapes with aspect ratios that increase with the shear rate [50] and that fluid-fluid phase separated colloid-polymer systems exhibit flow-elongated domains with aspect ratios dependent on the shear rate [51]. In our system, the low contrast between the colloid-dense regions and the continuous phase at temperatures below the LCST (given that the swollen pNIPAM particles are mostly water) make precise measurements of the shape of colloid-dense regions as a function of shear rate difficult. Nevertheless, we can use image autocorrelations to measure the approximate size of colloid-dense domains in the shear and vorticity directions as a function of shear rate at a constant temperature of 31°C. As seen in Supplementary Figure S3, the aspect ratio of these colloid-dense domains increases with increasing shear rate.

We believe that it is the morphology of the colloid-dense regions—a property dependent on the shear rate used as the sample is heated—that largely influences whether vorticity-aligned flocs form. Of course, as previous work has found, other factors like gap height, particle volume fraction, and interparticle attraction also matter, but those factors are not systematically varied in our present work. We contend that a consideration of the mesoscale colloidal-cluster morphology explains why we see the formation of vorticity-aligned flocs for an intermediate range of shear rates.

For 2s1γ˙T=γ˙hold7s1, we observe vorticity-aligned flocs forming within 30 min of steady shear after the sample has reached 34 °C. As seen in Figure 3, vorticity-aligned flocs are easily observed with our rheo-optical setup across a range of scales. Note that while in Figure 3A, vorticity-aligned flocs are seen throughout the sample geometry, in other instances, we observe flocs primarily towards the edges or near the center of the geometry, as one might expect for a parallel plate geometry where the shear rate varies with radial distance. The range of approximately 2 s−1–7 s−1 is where we observe flocs throughout a large fraction of the geometry. Besides visual observation as seen in Figures 3A–C, floc formation can also be inferred from the measured shear stress as a function of time. As seen in Figures 4A–C, for some shear rates we observe that minutes or longer into the steady shear phase of our procedure there is a marked increase and fluctuations in the stress. We observe this for shear rates between 2 and 7 s−1, whereas for shear rates less than 2 s−1 or greater than 7 s−1, the stress varies little with time and shows no large fluctuations or increases. The stress is plotted versus time as logs form in Supplementary Figure S4 as well (see Supplementary Videos S2, S3 for movies of the floc formation seen in Supplementary Figure S4A).

FIGURE 3
www.frontiersin.org

FIGURE 3. Shear-induced vorticity-aligned floc formation across scales. (A) Using a macro lens to capture a quarter of the parallel-plate geometry (diameter of 40 mm), we observe samples having undergone a temperature ramp to 34°C with a constant shear rate of γ˙T=5s1. Images are shown after the sample has reached 34°C and is being sheared at γ˙hold=5s1. The time interval between each image in the set of three (i–iii) is 70 s. The first image (left) was taken seconds after the sample reached 34°C. While here we observe flocs throughout the geometry, we note that in other cases we observe flocs mainly near the center or near the edges of the geometry. (B) The initial formation of a vorticity-aligned floc is observed using a 2× objective while the sample is shared at γ˙hold=4s1 after reading 34°C. The time interval between each image (i–vi) is 11.4 s. Counter-rotating of the bottom glass plate of the parallel-plate geometry was used to keep the floc in the field of view. Supplementary Video S1 captures this floc formation. (C) Image of a vorticity-aligned floc formed at 34°C using a shear rate of 4 s−1 taken with a 10× objective. The width of the floc is approximately equal to the parallel-plate gap height (here, h=157μm).

FIGURE 4
www.frontiersin.org

FIGURE 4. Shear-induced floc formation coincides with increases in the measured stress. (A) The colloid-polymer sample in a parallel-plate geometry is first heated to 34°C as described while subjected to shear at a rate γ˙T. Following this temperature ramp, the sample is subjected to the same shear rate γ˙hold=γ˙T while being held at 34°C. We measure the stress as a function of time during this 30-min period of constant shear rate and constant temperature. Depending on the shear rate used, the stress remains relatively flat for the 30 min or shows an uptick with increased fluctuations. This increase in the stress coincides with floc formation. We observe this indication of floc formation for shear rates of 2–7 s−1. The stress vs. time is shown for shear rates from 0.5 to 15 s−1 with the symbols indicating the mean over typically six runs. The shaded areas represent the mean ± standard error across the different runs. (B) Same as (A), but only the data for shear rates of 2, 3, 4, 5, 6, and 7 s−1 are shown. These are the shear rates for which there is a noticeable increase in the stress with time and where vorticity-aligned logs form. (C) As with (B), we show the same data in (A) but for only shear rates of 0.5, 1, 8, 9, 10, 11, 12, 13, 14, and 15 s−1. (D) Oscillatory measurements are performed after the temperature ramp to 34°C. During this temperature ramp, the sample was sheared at γ˙T. Plotted are the storage (G, solid symbols) and loss (G, open symbols) moduli for a frequency of ω=1rad/s and amplitude γ0=3% as a function of the preceding shear rate γ˙T. (E) Oscillatory measurements are performed after the sample is held at 34°C while subjected to a constant shear rate of γ˙hold (set to the same value as the preceding γ˙T used during the temperature ramp) for 30 min. As in (D), the storage and loss moduli are plotted as a function of γ˙hold. Outside of the range 2–7 s−1, we observe more liquid-like responses (GG). (F) The ratio of the storage modulus measured before and after the 30-min hold at 34°C and constant shear rate γ˙hold (equal to the preceding value for γ˙T). That is, the ratio of G in (D) and (E). The dotted line corresponds to a ratio of 1. Only for shear rates between 2 and 7 s−1 do we observe that the storage modulus increases (ratio > 1) after this 30-min period. The error bars in (D–F), as well as the shaded regions in (A–C), represent the standard error across multiple runs of different samples (of identical composition) loaded onto the rheometer on different days.

Previous studies investigating the formation of vorticity-aligned flocs have noted an upper bound on the shear rate (a critical shear rate, γ˙c) above which such flocs do not form [25, 27, 28]. This critical shear rate depends on the degree of confinement (γ˙ch1.4 from [25, 27]). Our work likewise shows that if our colloidal system is sheared at too high of a rate, vorticity-aligned flocs are not observed. We have not yet explored how this depends on the gap height of our parallel plate geometry, but for the value h150μm we find that γ˙c8s1.

However, our work also shows that there is a minimum shear rate below which vorticity-aligned flocs are not seen. We believe that this is due to how shear is deforming the colloid-dense mesoscale domains as the temperature increases and the interparticle interaction goes from a relatively weak depletion-induced attraction to a stronger hydrophobic attraction. At low shear rates, mesoscale colloidal clusters deform little due to the shear forces. At higher shear rates, colloidal clusters stretch out into high-aspect ratio structures. In comparison to more isotropic colloidal clusters, elongated clusters will form a larger number of connections with neighboring clusters. The fact that anisotropic particles can form a percolating network at a lower concentration than spherical particles has been shown in many studies of rod-shaped particles [28, 52]. In the experiments shown here, the concentration and aspect ratio of our pNIPAM colloidal particles are constants. But, the aspect ratios of the shear-deformed clusters do change. And just as rod-like particles will form space-spanning networks at a lower volume fraction than isotropic particles [52], elongated clusters will form more connections with neighboring clusters than isotropic clusters.

Moreover, the cluster morphology also affects the mechanical properties. After increasing the temperature from 23 to 34 °C during which we apply steady shear at a rate of γ˙T, we perform oscillatory rheological measurements with the sample held at 34°C. The viscoelastic moduli from this data are shown in Figure 4D. We find the largest G values are associated with prior shear rates γ˙T of between 2 and 4 s−1. More work needs to be done to determine why the colloidal gel loses elasticity as γ˙T increases beyond 4 s−1. Though very elongated clusters with a high aspect ratio may form a greater number of connections with surrounding clusters, perhaps the thinness of such clusters leads to a reduced elastic modulus in comparison to networks of thicker clusters with reduced aspect ratios. Or, perhaps at the higher shear rates, elongated clusters are quickly broken apart by shear forces. We also show G and G measured immediately after the sample is sheared at γ˙hold for 30 min at 34°C (before which the sample was sheared at γ˙T=γ˙hold during the temperature ramp to 34°C). As seen in Figures 4A,E, a solid-like mechanical response is observed for the range of shear rates where vorticity-aligned flocs form. Further, we show in Figure 4F the ratio of the elastic moduli after the 30 min of constant shearing and prior (that is, the ratio of the G data in Figure 4E to that in Figure 4D). For shear rates between about 2 and 7 s−1, this ratio is greater than one. Therefore, at these shear rates, the 30-min period of constant shearing increases the elastic moduli. For shear rates lower than 2 s−1 and higher than 7 s−1, the period of constant shear reduced the elastic moduli. This again corroborates our optical imaging data (Figure 3) and stress vs. time data (Figures 4A–C) which shows that shear-driven flocculation is occurring for an intermediate range of shear rates.

Conclusion

Our data suggests that an important factor influencing the propensity for shear-induced macroscopic patterning in our samples is the mesoscopic cluster morphology. Prior studies have shown how the degree of confinement (i.e., gap height), interparticle attraction, and volume fraction dictate whether vorticity-aligned log-rolling flocs will form. Our results show that for systems like ours the structure of the sample on the cluster scale also matters. A sample comprised of elongated, fiber-like colloidal clusters seems more likely to form shear-induced vorticity-aligned flocs than a sample containing isotropic colloidal clusters. This finding dovetails with recent work looking at the shear driven patterns of colloidal rod suspensions [28]. They found that vorticity-aligned flocs formed for rods at lower volume fractions than for spheres, due to the lower volume fraction of rods needed to form a percolating structure. From past work on colloidal rods [52], one can appreciate how shaping the colloidal clusters into elongated structures will more readily lead to a well-connected network if one considers the colloidal clusters as renormalized particles. In our work, we show how the colloidal clusters can be shaped using shear flow applied to the sample when the colloidal particle’s attraction is weak enough to allow colloid-rich domains to stretch out in the direction of flow. However, for other systems, control over the cluster morphology could be achieved by using different external fields or the interparticle potential [53].

More generally, this work highlights the importance of the experimental history of a sample. How a colloidal gel’s rheological properties depend on its shear history (i.e., rheological hysteresis) has been the subject of much research [54, 55]. For our thermoresponsive system, both the shear history and the temperature history are factors in the rheology and structure. This thermokinematic memory has been explored in other industrially relevant materials [56] and relates to recent work probing how temperature changes and the rate of temperature changes can (potentially in combination with shear) structure colloidal materials [47, 5759].

Pattern formation due to flow in colloidal and complex fluids can be used to manufacture small-scale structures. It can also be a hinderance when material needs to flow through confined spaces, such as in 3D printing. Therefore, a complete understanding of shear driven structuration will help streamline current processing and potentially lead to new strategies for tuning soft materials.

Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author contributions

RM conceived and designed the study. RR and DT performed experiments. RM, RR, and DT analysed data. RM and RR wrote the manuscript. All authors helped revise and approve the submitted manuscript.

Funding

Support for this research came from a National Science Foundation Major Research Instrumentation award (CBET 1919429). Support also came from a Cottrell Scholars Award from the Research Corporation for Science Advancement to RM.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2022.955006/full#supplementary-material

References

1. Fernandez-Nieves A, Puertas AM. Fluids, colloids and soft materials: An introduction to soft matter physics. John Wiley & Sons (2016).

Google Scholar

2. Hensley A, Jacobs WM, Rogers WB. Self-assembly of photonic crystals by controlling the nucleation and growth of DNA-coated colloids. Proc Natl Acad Sci U S A (2022) 119(1):e2114050118. doi:10.1073/pnas.2114050118

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Royall CP, Faers M, Fussell S, Hallett J. Real space analysis of colloidal gels: Triumphs, challenges and future directions. J Phys : Condens Matter (2021) 33:453002. doi:10.1088/1361-648x/ac04cb

CrossRef Full Text | Google Scholar

4. Rodriguez AMB, Binks BP. Capsules from Pickering emulsion templates. Curr Opin Colloid Interf Sci (2019) 44:107–29. doi:10.1016/j.cocis.2019.09.006

CrossRef Full Text | Google Scholar

5. Brader JM. Nonlinear rheology of colloidal dispersions. J Phys : Condens Matter (2010) 22(36):363101. doi:10.1088/0953-8984/22/36/363101

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Zhu C, Pascall AJ, Dudukovic N, Worsley MA, Kuntz JD, Duoss EB, et al. Colloidal materials for 3D printing. Annu Rev Chem Biomol Eng (2019) 10(1):17–42. doi:10.1146/annurev-chembioeng-060718-030133

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Vermant J, Solomon MJ. Flow-induced structure in colloidal suspensions. J Phys : Condens Matter (2005) 17(4):R187–216. doi:10.1088/0953-8984/17/4/r02

CrossRef Full Text | Google Scholar

8. Joy M, Muangnapoh T, Snyder MA, Gilchrist JF. Flow-induced alignment of (100) fcc thin film colloidal crystals. Soft Matter (2015) 11(36):7092–100. doi:10.1039/c5sm01076d

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Ackerson BJ, Pusey PN. Shear-induced order in suspensions of hard spheres. Phys Rev Lett (1988) 61(8):1033–6. doi:10.1103/physrevlett.61.1033

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Cheng X, McCoy JH, Israelachvili JN, Cohen I. Imaging the microscopic structure of shear thinning and thickening colloidal suspensions. Science (2011) 333(6047):1276–9. doi:10.1126/science.1207032

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Min Kim J, Eberle APR, Kate Gurnon A, Porcar L, Wagner NJ. The microstructure and rheology of a model, thixotropic nanoparticle gel under steady shear and large amplitude oscillatory shear (Laos). J Rheology (2014) 58(5):1301–28. doi:10.1122/1.4878378

CrossRef Full Text | Google Scholar

12. Boromand A, Jamali S, Maia JM. Structural fingerprints of yielding mechanisms in attractive colloidal gels. Soft Matter (2017) 13(2):458–73. doi:10.1039/c6sm00750c

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Colombo G, Kim S, Schweizer T, Schroyen B, Clasen C, Mewis J, et al. Superposition rheology and anisotropy in rheological properties of sheared colloidal gels. J Rheology (2017) 61(5):1035–48. doi:10.1122/1.4998176

CrossRef Full Text | Google Scholar

14. Wang Y, Ewoldt RH. New insights on carbon black suspension rheology -- anisotropic thixotropy and anti-thixotropy (2022). arXiv:220205772 [cond-mat].

Google Scholar

15. DeGroot JV, Macosko CW, Kume T, Hashimoto T. Flow-induced anisotropic SALS in silica-filled PDMS liquids. J Colloid Interf Sci (1994) 166(2):404–13. doi:10.1006/jcis.1994.1311

CrossRef Full Text | Google Scholar

16. Pignon F, Magnin A, Piau J-M. Butterfly light scattering pattern and rheology of a sheared thixotropic clay gel. Phys Rev Lett (1997) 79(23):4689–92. doi:10.1103/physrevlett.79.4689

CrossRef Full Text | Google Scholar

17. Cheng X, Xu X, Rice SA, Dinner AR, Cohen I. Assembly of vorticity-aligned hard-sphere colloidal strings in a simple shear flow. Proc Natl Acad Sci U S A (2012) 109(1):63–7. doi:10.1073/pnas.1118197108

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Lin NYC, Cheng X, Cohen I. Biaxial shear of confined colloidal hard spheres: The structure and rheology of the vorticity-aligned string phase. Soft Matter (2014) 10(12):1969–76. doi:10.1039/c3sm52880d

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Montesi A, Peña AA, Pasquali M. Vorticity alignment and negative normal stresses in sheared attractive emulsions. Phys Rev Lett (2004) 92(5):058303. doi:10.1103/physrevlett.92.058303

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Mao C, Huang Y, Qiao Y, Zhang J, Kong M, Yang Q, et al. Vorticity-aligned droplet bands in sheared immiscible polymer blends induced by solid particles. Langmuir (2020) 36(16):4383–95. doi:10.1021/acs.langmuir.0c00511

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Lin-Gibson S, Pathak JA, Grulke EA, Wang H, Hobbie EK. Elastic flow instability in nanotube suspensions. Phys Rev Lett (2004) 92(4):048302. doi:10.1103/physrevlett.92.048302

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Ma AWK, Mackley MR, Rahatekar SS. Experimental observation on the flow-induced assembly of Carbon nanotube suspensions to form helical bands. Rheol Acta (2007) 46(7):979–87. doi:10.1007/s00397-007-0183-x

CrossRef Full Text | Google Scholar

23. Osuji CO, Weitz DA. Highly anisotropic vorticity aligned structures in a shear thickening attractive colloidal system. Soft Matter (2008) 4(7):1388–92. doi:10.1039/b716324j

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Negi AS, Osuji CO. New insights on fumed colloidal rheology—Shear thickening and vorticity-aligned structures in flocculating dispersions. Rheol Acta (2009) 48(8):871–81. doi:10.1007/s00397-008-0341-9

CrossRef Full Text | Google Scholar

25. Grenard V, Taberlet N, Manneville S. Shear-induced structuration of confined carbon black gels: Steady-state features of vorticity-aligned flocs. Soft Matter (2011) 7(8):3920–8. doi:10.1039/c0sm01515f

CrossRef Full Text | Google Scholar

26. Shakeel A, MacIver MR, van Kan PJM, Kirichek A, Chassagne C. A rheological and microstructural study of two-step yielding in mud samples from a port area. Colloids Surf A: Physicochemical Eng Aspects (2021) 624:126827. doi:10.1016/j.colsurfa.2021.126827

CrossRef Full Text | Google Scholar

27. Varga Z, Grenard V, Pecorario S, Taberlet N, Dolique V, Manneville S, et al. Hydrodynamics control shear-induced pattern formation in attractive suspensions. Proc Natl Acad Sci U S A (2019) 116:12193–8. doi:10.1073/pnas.1901370116

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Das M, Chambon L, Varga Z, Vamvakaki M, Swan JW, Petekidis G. Shear driven vorticity aligned flocs in a suspension of attractive rigid rods. Soft Matter (2021) 17(5):1232–45. doi:10.1039/d0sm01576h

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Huang J, Xu J, Sheng Y, Zhu Y, Jiang W, Xu D, et al. Fabrication of polymer film with extraordinary conductive anisotropy by forming parallel conductive vorticity-aligned stripes and its formation mechanism. Macromol Mater Eng (2016) 301(6):743–9. doi:10.1002/mame.201600026

CrossRef Full Text | Google Scholar

30. Mao C, Huang J, Zhu Y, Jiang W, Tang Q, Ma X. Tailored parallel graphene stripes in plastic film with conductive anisotropy by shear-induced self-assembly. J Phys Chem Lett (2013) 4(1):43–7. doi:10.1021/jz301811b

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Romeo G, Fernandez‐Nieves A, Wyss HM, Acierno D, Weitz DA. Temperature-controlled transitions between glass, liquid, and gel states in dense p-NIPA suspensions. Adv Mater (2010) 22(31):3441–5. doi:10.1002/adma.200904189

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Yunker PJ, Chen K, Gratale MD, Lohr MA, Still T, Yodh AG. Physics in ordered and disordered colloidal matter composed of poly( N -isopropylacrylamide) microgel particles. Rep Prog Phys (2014) 77(5):056601. doi:10.1088/0034-4885/77/5/056601

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Minami S, Watanabe T, Suzuki D, Urayama K. Rheological properties of suspensions of thermo-responsive poly( N -isopropylacrylamide) microgels undergoing volume phase transition. Polym J (2016) 48(11):1079–86. doi:10.1038/pj.2016.79

CrossRef Full Text | Google Scholar

34. Appel J, Fölker B, Sprakel J. Mechanics at the glass-to-gel transition of thermoresponsive microgel suspensions. Soft Matter (2016) 12(9):2515–22. doi:10.1039/c5sm02940f

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Karg M, Pich A, Hellweg T, Hoare T, Lyon LA, Crassous JJ, et al. Nanogels and microgels: From model colloids to applications, recent developments, and future trends. Langmuir (2019) 35(19):6231–55. doi:10.1021/acs.langmuir.8b04304

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Chaudhary G, Ghosh A, Kang JG, Braun PV, Ewoldt RH, Schweizer KS. Linear and nonlinear viscoelasticity of concentrated thermoresponsive microgel suspensions. J Colloid Interf Sci (2021) 601:886–98. doi:10.1016/j.jcis.2021.05.111

CrossRef Full Text | Google Scholar

37. Ramakrishnan S, Chen Y-L, Schweizer KS, Zukoski CF. Elasticity and clustering in concentrated depletion gels. Phys Rev E (2004) 70(4):040401. doi:10.1103/physreve.70.040401

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Laurati M, Petekidis G, Koumakis N, Cardinaux F, Schofield AB, Brader JM, et al. Structure, dynamics, and rheology of colloid-polymer mixtures: From liquids to gels. J Chem Phys (2009) 130(13):134907. doi:10.1063/1.3103889

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Zaccone A, Wu H, Del Gado E. Elasticity of arrested short-ranged attractive colloids: Homogeneous and heterogeneous glasses. Phys Rev Lett (2009) 103(20):208301. doi:10.1103/physrevlett.103.208301

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Whitaker KA, Varga Z, Hsiao LC, Solomon MJ, Swan JW, Furst EM. Colloidal gel elasticity arises from the packing of locally glassy clusters. Nat Commun (2019) 10(1):2237. doi:10.1038/s41467-019-10039-w

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Dang S, Brady J, Rel R, Surineni S, O’Shaughnessy C, McGorty R. Core–shell droplets and microcapsules formed through liquid–liquid phase separation of a colloid–polymer mixture. Soft Matter (2021) 17:8300–7. doi:10.1039/d1sm01091c

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Wang J, McGorty R. Measuring capillary wave dynamics using differential dynamic microscopy. Soft Matter (2019) 15(37):7412–9. doi:10.1039/c9sm01508f

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Still T, Chen K, Alsayed AM, Aptowicz KB, Yodh AG. Synthesis of micrometer-size poly(N-isopropylacrylamide) microgel particles with homogeneous crosslinker density and diameter control. J Colloid Interf Sci (2013) 405:96–102. doi:10.1016/j.jcis.2013.05.042

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Borrega R, Cloitre M, Betremieux I, Ernst B, Leibler L. Concentration dependence of the low-shear viscosity of polyelectrolyte micro-networks: From hard spheres to soft microgels. Europhys Lett (1999) 47(6):729–35. doi:10.1209/epl/i1999-00451-1

CrossRef Full Text | Google Scholar

45. Jamie EAG, Davies GJ, Howe MD, Dullens RPA, Aarts DGAL. Thermal capillary waves in colloid–polymer mixtures in water. J Phys : Condens Matter (2008) 20(49):494231. doi:10.1088/0953-8984/20/49/494231

CrossRef Full Text | Google Scholar

46. Edelstein AD, Tsuchida MA, Amodaj N, Pinkard H, Vale RD, Stuurman N. Advanced methods of microscope control using μManager software. J Biol Methods (2014) 1(2):10. doi:10.14440/jbm.2014.36

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Fussell SL, Royall CP, van Duijneveldt JS. Controlling phase separation in microgel-polymeric micelle mixtures using variable quench rates (2021). arXiv:210404022 [cond-mat].

Google Scholar

48. Taylor GI. The viscosity of a fluid containing small drops of another fluid. Proc R Soc Lond Ser A, Containing Pap a Math Phys Character (1932) 138(834):41–8.

Google Scholar

49. Taylor GI. the formation of emulsions in definable fields of flow. Proc R Soc Lond A: Math Phys Eng Sci (1934) 146(858):501–23.

Google Scholar

50. Pathak JA, Migler KB. Droplet−String deformation and stability during microconfined shear flow. Langmuir (2003) 19(21):8667–74. doi:10.1021/la0346907

CrossRef Full Text | Google Scholar

51. Derks D, Aarts DGAL, Bonn D, Imhof A. Phase separating colloid polymer mixtures in shear flow. J Phys : Condens Matter (2008) 20(40):404208. doi:10.1088/0953-8984/20/40/404208

CrossRef Full Text | Google Scholar

52. Solomon MJ, Spicer PT. Microstructural regimes of colloidal rod suspensions, gels, and glasses. Soft Matter (2010) 6(7):1391–400. doi:10.1039/b918281k

CrossRef Full Text | Google Scholar

53. Soto-Bustamante F, Valadez-Pérez NE, Liu Y, Castañeda-Priego R, Laurati M. Clusters in colloidal dispersions with a short-range depletion attraction: Thermodynamic identification and morphology. J Colloid Interf Sci (2022) 618:442–50. doi:10.1016/j.jcis.2022.03.061

CrossRef Full Text | Google Scholar

54. Divoux T, Grenard V, Manneville S. Rheological hysteresis in soft glassy materials. Phys Rev Lett (2013) 110(1):018304. doi:10.1103/physrevlett.110.018304

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Jamali S, Armstrong RC, McKinley GH. Multiscale nature of thixotropy and rheological hysteresis in attractive colloidal suspensions under shear. Phys Rev Lett (2019) 123(24):248003. doi:10.1103/physrevlett.123.248003

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Geri M, Venkatesan R, Sambath K, McKinley GH. Thermokinematic memory and the thixotropic elasto-viscoplasticity of waxy crude oils. J Rheology (2017) 61(3):427–54. doi:10.1122/1.4978259

CrossRef Full Text | Google Scholar

57. Smith KM, Hsiao LC. Migration and morphology of colloidal gel clusters in cylindrical channel flow. Langmuir (2021) 37(34):10308–18. doi:10.1021/acs.langmuir.1c01287

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Cardenas-Vasquez ED, Smith KM, Doolan TJ, Hsiao LC. Shear-induced microstructural variations in nanoemulsion-laden organohydrogel fibers. ACS Appl Polym Mater (2020) 2(2):594–603. doi:10.1021/acsapm.9b00979

CrossRef Full Text | Google Scholar

59. Hsiao LC, Doyle PS. Celebrating Soft Matter's 10th Anniversary: Sequential phase transitions in thermoresponsive nanoemulsions. Soft Matter (2015) 11(43):8426–31. doi:10.1039/c5sm01581b

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: colloids, rheology, microgel particles, rheo-optics, shear-induced aggregation

Citation: Rel R, Terwilliger D and McGorty R (2022) Shear-induced vorticity aligned flocs in a temperature responsive colloid-polymer mixture. Front. Phys. 10:955006. doi: 10.3389/fphy.2022.955006

Received: 27 May 2022; Accepted: 22 August 2022;
Published: 15 September 2022.

Edited by:

Emanuela Del Gado, Georgetown University, United States

Reviewed by:

Orhan Özdemir, Istanbul University-Cerrahpasa, Turkey
Ryan Poling-Skutvik, University of Rhode Island, United States
Ali Mohraz, University of California, Irvine, United States

Copyright © 2022 Rel, Terwilliger and McGorty. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ryan McGorty, cm1jZ29ydHlAc2FuZGllZ28uZWR1

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.