Shear-induced order in suspensions of hard spheres

Kinetics of Shear-Induced Structural Ordering in Dense Colloids

The macroscopic rheological response of a colloidal solution is highly correlated with the local microscopic structure, as revealed by an in situ Rheo-SAXS experiment with a high temporal resolution. Oscillatory shear can induce a strain-controlled ordering-to-disorder transition, resulting in a shear-thickening process that is different from the normal shear-thickening behavior that is driven by hydrodynamics and particle friction. We reveal that there is a complex time-dependent kinetics toward structural ordering under different applied strains. When the strain amplitude reaches a critical value that starts to induce disordering in the system, the pathway toward the dynamic equilibrium can also become highly non-monotonic. Within the same oscillatory cycle, there is a strong correlation of ordering with different phases of the oscillation, with the system oscillating between two dynamic metastable states.

Does supercooled water retain its universal nucleation behavior under shear at high pressure?

Understanding the nucleation of homogeneous flow systems at high pressures is vital in protein crystallization and cryopreservation, where high pressure prevents the freezing of biological samples. This study examines the behavior of ice nucleation under shear at various pressures and explores the universal nucleation behavior of the sheared systems applied to supercooled water at higher pressures. In this study, the nucleation rates for the TIP4P/Ice model via a seeding method based on extended classical nucleation theory (CNT) are computed at pressures of 1, 100, 500, 700, and 1000 bar and a constant temperature of 240 K. Using extended CNT with explicitly embodying the shear rate, we analyzed the dependence of pressure on the transport and thermodynamic properties. In line with previous studies, we observed that Δμliq-ice and viscosity decrease while diffusivity increases with an increase in pressure. Furthermore, we showed that the dependence of the nucleation rate on shear at higher pressure is non-monotonic, with the maximum at optimal shear rates between 107 and 108 s-1. Our results demonstrate a non-monotonic pressure dependence of the optimal shear rates, which could originate from a violation of the Stokes-Einstein relation.

Motility-induced shear thickening in dense colloidal suspensions

Phase transitions and collective dynamics of active colloidal suspensions are fascinating topics in soft matter physics, particularly for out-of-equilibrium systems, which can lead to rich rheological behaviours in the presence of steady shear flow. Here the role of self-propulsion in the rheological response of a dense colloidal suspension is investigated by using particle-resolved Brownian dynamics simulations. First, the combined effect of activity and shear in the solid on the disordering transition of the suspension is analyzed. While both self-propulsion and shear destroy order and melt the system if critical values are exceeded, self-propulsion largely lowers the stress barrier needed to be overcome during the transition. We further explore the rheological response of the active sheared system once a steady state is reached. While passive suspensions show a solid-like behaviour, turning on particle motility fluidises the system. At low self-propulsion, the active suspension behaves in the steady state as a shear-thinning fluid. Increasing the self-propulsion changes the behaviour of the liquid from shear-thinning to shear-thickening. We attribute this to clustering in the sheared suspensions induced by motility. This new phenomenon of motility-induced shear thickening (MIST) can be used to tailor the rheological response of colloidal suspensions.

Treating random sequential addition via the replica method

While many physical processes are non-equilibrium in nature, the theory and modeling of such phenomena lag behind theoretical treatments of equilibrium systems. The diversity of powerful theoretical tools available to describe equilibrium systems has inspired strategies that map non-equilibrium systems onto equivalent equilibrium analogs so that interrogation with standard statistical mechanical approaches is possible. In this work, we revisit the mapping from the non-equilibrium random sequential addition process onto an equilibrium multi-component mixture via the replica method, allowing for theoretical predictions of non-equilibrium structural quantities. We validate the above approach by comparing the theoretical predictions to numerical simulations of random sequential addition.

Shear-induced orientational ordering in an active glass former

Dense assemblies of self-propelled particles that can form solid-like states also known as active or living glasses are abundant around us, covering a broad range of length scales and timescales: from the cytoplasm to tissues, from bacterial biofilms to vehicular traffic jams, and from Janus colloids to animal herds. Being structurally disordered as well as strongly out of equilibrium, these systems show fascinating dynamical and mechanical properties. Using extensive molecular dynamics simulation and a number of distinct dynamical and mechanical order parameters, we differentiate three dynamical steady states in a sheared model active glassy system: 1) a disordered state, 2) a propulsion-induced ordered state, and 3) a shear-induced ordered state. We supplement these observations with an analytical theory based on an effective single-particle Fokker-Planck description to rationalize the existence of the shear-induced orientational ordering behavior in an active glassy system without explicit aligning interactions of, for example, Vicsek type. This ordering phenomenon occurs in the large persistence time limit and is made possible only by the applied steady shear. Using a Fokker-Planck description with parameters that can be measured independently, we make testable predictions for the joint distribution of single-particle position and orientation. These predictions match well with the joint distribution measured from direct numerical simulation. Our results are of relevance for experiments exploring the rheological response of dense active colloids and jammed active granular matter systems.

An Experimental and Theoretical Determination of Oscillatory Shear-Induced Crystallization Processes in Viscoelastic Photonic Crystal Media

A study is presented of the oscillatory shear-ordering dynamics of viscoelastic photonic crystal media, using an optical shear cell. The hard-sphere/“sticky”-shell design of these polymeric composite particles produces athermal, quasi-solid rubbery media, with a characteristic viscoelastic ensemble response to applied shear. Monotonic crystallization processes, as directly measured by the photonic stopband transmission, are tracked as a function of strain amplitude, oscillation frequency, and temperature. A complementary generic spatio-temporal model is developed of crystallization due to shear-dependent interlayer viscosity, giving propagating crystalline fronts with increasing applied strain, and a gradual transition from interparticle disorder to order. The introduction of a competing shear-induced flow degradation process, dependent on the global shear rate, gives solutions with both amplitude and frequency dependence. The extracted crystallization timescales show parametric trends which are in good qualitative agreement with experimental observations.

Oscillatory shear-induced bcc-fcc martensitic transformation in a colloidal suspension with long-range repulsive interactions

We perform non-equilibrium Brownian dynamics simulations to investigate the out-of-equilibrium phase behavior of a suspension of charged colloids under external oscillatory shear. We independently vary the frequency f and the maximum strain amplitude γ_max of the oscillations and map out an out-of-equilibrium phase diagram in the f-γ_max plane. Similar to what has been observed in earlier studies on colloidal hard spheres, we find the formation of a twinned face-centered-cubic phase in a specific range of γ_max, which displays a martensitic transition to a body-centered-cubic crystal within half of the oscillation cycle. We provide a comprehensive analysis of these structures and show how the system transforms from one to the other. We also report evidence of a sliding layer phase and a string phase.

Shear-induced glass-to-crystal transition in anisotropic clay-like suspensions

A new numerical framework based on Stokesian dynamics is used to study a shear-induced glass-to-crystal transition in suspensions of clay-like anisotropically charged platelets. The structures obtained in quiescent conditions are in agreement with previous Monte Carlo results: a liquid phase at very short interaction range (high salt concentration), phase separation and a gel without large scale density fluctuations at intermediate interaction ranges, and glassy states at very large interaction ranges. When initially glassy suspensions are sheared, hydrodynamic torques first rotate platelets so they can reach a transient quasi-nematic disordered state. These orientational correlations permit to unlock translational degrees of freedom and the platelets then form strings aligned with the velocity direction and hexagonally packed in the gradient-vorticity plane. Under steady shear, platelet orientations are correlated but the system is not nematic. After flow cessation and relaxation in quiescent conditions, positional and orientational order are further improved as the platelet suspension experiences a transition to a nematic hexagonal crystal. Energy calculations and the existence of residual stress anisotropy after relaxation show that this final structure is not an equilibrium state but rather a new ordered, arrested state. The transient, nematic, disordered state induced by shear immediately after startup and unlocking translational degrees of freedom is thought to be an initial step that may be generic for other suspensions of strongly anisotropic colloids with important translation-orientation coupling induced by long-range interactions.

Wide-angle static and dynamic light scattering under shear

We develop and characterize a wide angle static and dynamic light scattering under shear setup. The apparatus is suitable for the study of the structure and the dynamics of soft materials systems with a sub-micron characteristic length scale. The shear device consists in two parallel plates, and the optical setup allows us to perform light scattering measurements in any plane that contains the gradient of the velocity field direction. We demonstrate several capabilities of our apparatus: a measurement of the evolution with shear of the first peak of the structure factor of a concentrated suspension of spherical particles, both in the compression and extension quadrants of the shear flow, and the measurement of the velocity profile in dynamic light scattering. We present a theoretical treatment of light scattering under flow that takes into account the Gaussian character of the illumination and detection optical paths, in the case where the scattering volume extension is smaller than the gap of the flow cell, and compare with experimental measurements.

Substrate wettability guided oriented self assembly of Janus particles

Self-assembly of Janus particles with spatial inhomogeneous properties is of fundamental importance in diverse areas of sciences and has been extensively observed as a favorably functionalized fluidic interface or in a dilute solution. Interestingly, the unique and non-trivial role of surface wettability on oriented self-assembly of Janus particles has remained largely unexplored. Here, the exclusive role of substrate wettability in directing the orientation of amphiphilic metal-polymer Bifacial spherical Janus particles, obtained by topo-selective metal deposition on colloidal Polymestyere (PS) particles, is explored by drop casting a dilute dispersion of the Janus colloids. While all particles orient with their polymeric (hydrophobic) and metallic (hydrophilic) sides facing upwards on hydrophilic and hydrophobic substrates respectively, they exhibit random orientation on a neutral substrate. The substrate wettability guided orientation of the Janus particles is captured using molecular dynamic simulation, which highlights that the arrangement of water molecules and their local densities near the substrate guide the specific orientation. Finally, it is shown that by spin coating it becomes possible to create a hexagonal close-packed array of the Janus colloids with specific orientation on differential wettability substrates. The results reported here open up new possibilities of substrate-wettability driven functional coatings of Janus particles, which has hitherto remained unexplored.

The first normal stress difference of non-Brownian hard-sphere suspensions in the oscillatory shear flow near the liquid and crystal coexistence region

We carry out a numerical study to investigate the dynamics of non-Brownian hard-sphere suspensions near the liquid and crystal coexistence region in small to large amplitude oscillatory shear flow. The first normal stress difference (N1) and related rheological functions are carefully analyzed, focusing on the strain stiffening phenomenon, which occurs in the large strain amplitude region. Under oscillatory shear, we observe several unique behaviors of N1. A negative nonzero mean value of N1 (N1,0) is observed for the applied strain amplitudes. The change of the sign, from negative to positive, at the maximum value of N1 (N1,max) is observed at a specific point, which is not consistent with the critical strain amplitude (γ0,c) at which the modulus begins to deviate from linear viscoelasticity. The behavior of N1 in the oscillatory shear flow is different from that of N1 in steady shear flow, that is, the characteristics of N1 in strain stiffening and shear thickening are quite distinguished from each other. We also perform structural analysis to confirm the relationship between the rheological properties and microstructure of the suspension. A strong correlation is observed between the global bond order parameter (Ψ6) and the distortions in both nonlinear shear and normal stresses. The most noticeable characteristic is captured through the maximum of the global bond order parameter (Ψ6,max). The strain amplitude at the slope change of Ψ6,max corresponds to the point where a unique behavior of N1 is observed, i.e. the change of the sign in N1,max, but a strong correlation is not captured at γ0,c. This demonstrates that the normal stress responds to particle ordering more sensitively than other rheological functions based on shear stress like dynamic moduli. As far as we are concerned, the behavior of N1 has rarely been fully explored and related with the strain stiffening of non-Brownian suspensions so far. Therefore, this study has significance as the first report to strictly analyze strain stiffening along with the first normal stress difference N1.

Shear thinning and thickening in dispersions of spherical nanoparticles

We present a molecular dynamics study of the flow of rigid spherical nanoparticles in a simple fluid. We evaluate the viscosity of the dispersion as a function of shear rate and nanoparticle volume fraction. We observe shear-thinning behavior at low volume fractions; as the shear rate increases, the shear forces overcome the Brownian forces, resulting in more frequent and more violent collisions between the nanoparticles. This in turn results in more dissipation. We show that in order to be in the shear-thinning regime the nanoparticles have to order themselves into layers longitudinal to the flow to minimize the collisions. As the nanoparticle volume fraction increases there is less room to form the ordered layers; consequently as the shear rate increases the nanoparticles collide more, which results in turn in shear thickening. Most interestingly, we show that at intermediate volume fractions the system exhibits metastability, with successions of ordered and disordered states along the same trajectory. Our results suggest that for nanoparticles in a simple fluid the hydroclustering phenomenon is not present; instead the order-disorder transition is the leading mechanism for the transition from shear thinning to shear thickening.

Memory-induced motion reversal in Brownian liquids

We study the Brownian dynamics of hard spheres under spatially inhomogeneous shear, using event-driven Brownian dynamics simulations and power functional theory. We examine density and current profiles both for steady states and for the transient dynamics after switching on and switching off an external square wave shear force field. We find that a dense hard sphere fluid (volume fraction ≈0.35) undergoes global motion reversal after switching off the shear force field. We use power functional theory with a spatially nonlocal memory kernel to describe the superadiabatic force contributions and obtain good quantitative agreement of the theoretical results with simulation data. The theory provides an explanation for the motion reversal: internal superadiabatic nonequilibrium forces that oppose the externally driven current arise due to memory after switching off. The effect is genuinely viscoelastic: in steady state, viscous forces oppose the current, but they elastically generate an opposing current after switch-off.

Compression-Responsive Photonic Crystals Based on Fluorine-Containing Polymers

Fluoropolymers represent a unique class of functional polymers due to their various interesting and important properties such as thermal stability, resistance toward chemicals, repellent behaviors, and their low refractive indices in comparison to other polymeric materials. Based on the latter optical property, fluoropolymers are particularly of interest for the preparation of photonic crystals for optical sensing application. Within the present study, photonic crystals were prepared based on core-interlayer-shell particles focusing on fluoropolymers. For particle assembly, the melt-shear organization technique was applied. The high order and refractive index contrast of the individual components of the colloidal crystal structure lead to remarkable reflection colors according to Bragg’s law of diffraction. Due to the special architecture of the particles, consisting of a soft core, a comparably hard interlayer, and again a soft shell, the resulting opal films were capable of changing their shape and domain sizes upon applied pressure, which was accompanied with a (reversible) change of the observed reflection colors as well. By the incorporation of adjustable amounts of UV cross-linking agents into the opal film and subsequent treatment with different UV irradiation times, stable and pressure-sensitive opal films were obtained. It is shown that the present strategy led to (i) pressure-sensitive opal films featuring reversibly switchable reflection colors and (ii) that opal films can be prepared, for which the written pattern—resulting from the compressed particles—could be fixed upon subsequent irradiation with UV light. The herein described novel fluoropolymer-containing photonic crystals, with their pressure-tunable reflection color, are promising candidates in the field of sensing devices and as potential candidates for anti-counterfeiting materials.

Classical nucleation theory for the crystallization kinetics in sheared liquids

While statistical mechanics provides a comprehensive framework for the understanding of equilibrium phase behavior, predicting the kinetics of phase transformations remains a challenge. Classical nucleation theory (CNT) provides a thermodynamic framework to relate the nucleation rate to thermodynamic quantities such as pressure difference and interfacial tension through the nucleation work necessary to spawn critical nuclei. However, it remains unclear whether such an approach can be extended to the crystallization of driven melts that are subjected to mechanical stresses and flows. Here, we demonstrate numerically for hard spheres that the impact of simple shear on the crystallization rate can be rationalized within the CNT framework by an additional elastic work proportional to the droplet volume. We extract the local stress and strain inside solid droplets, which yield size-dependent values for the shear modulus that are about half of the bulk value. Finally, we show that for a complete description one also has to take into account the change of interfacial work between the strained droplet and the sheared liquid. From scaling reasons, we expect this extra contribution to dominate the work formation of small nuclei but become negligible compared to the elastic work for droplets composed of a few hundreds of particles.

Boundary-induced inhomogeneity of particle layers in the solidification of suspensions

When a suspension freezes, a compacted particle layer builds up at the solidification front with noticeable implications on the freezing process. In a directional solidification experiment of monodisperse suspensions in thin samples, we evidence a link between the thickness of this layer and the sample depth. We attribute it to an inhomogeneity of particle density that is attested by the evidence of crystallization at the plates and of random close packing far from them. A mechanical model based on the resulting modifications of permeability enables us to relate the layer thickness to this inhomogeneity and to select the distribution of particle density that yields the best fit to our data. This distribution involves an influence length of sample plates of about 11 particle diameters. Altogether, these results clarify the implications of boundaries on suspension freezing. They may be useful to model polydisperse suspensions with large particles playing the role of smooth boundaries with respect to small ones.

Effect of anisotropic diffusion on spinodal decomposition

We study the phase transition dynamics of a fluid system in which the particles diffuse anisotropically in space. The motivation to study such a situation is provided by systems of interacting magnetic colloidal particles subject to the Lorentz force. The Smoluchowski equation for the many-particle probability distribution then acquires an anisotropic diffusion tensor. We show that in comparison to the isotropic case, anisotropic diffusion results in qualitatively different dynamics of spinodal decomposition. Using the method of dynamical density functional theory, we predict that the intermediate-stage decomposition dynamics are slowed down significantly by anisotropy; the coupling between different Fourier modes is strongly reduced. Numerical calculations are performed for a model (Yukawa) fluid that exhibits gas-liquid phase separation.

Confined Polymers as Self-Avoiding Random Walks on Restricted Lattices

Polymers in highly confined geometries can display complex morphologies including ordered phases. A basic component of a theoretical analysis of their phase behavior in confined geometries is the knowledge of the number of possible single-chain conformations compatible with the geometrical restrictions and the established crystalline morphology. While the statistical properties of unrestricted self-avoiding random walks (SAWs) both on and off-lattice are very well known, the same is not true for SAWs in confined geometries. The purpose of this contribution is (a) to enumerate the number of SAWs on the simple cubic (SC) and face-centered cubic (FCC) lattices under confinement for moderate SAW lengths, and (b) to obtain an approximate expression for their behavior as a function of chain length, type of lattice, and degree of confinement. This information is an essential requirement for the understanding and prediction of entropy-driven phase transitions of model polymer chains under confinement. In addition, a simple geometric argument is presented that explains, to first order, the dependence of the number of restricted SAWs on the type of SAW origin.

Assembly of colloidal particles in solution

Advances in both top-down and bottom-up syntheses of a wide variety of complex colloidal building blocks and also in methods of controlling their assembly in solution have led to new and interesting forms of highly controlled soft matter. In particular, top-down lithographic methods of producing monodisperse colloids now provide precise human-designed control over their sub-particle features, opening up a wide range of new possibilities for assembly structures that had been previously limited by the range of shapes available through bottom-up methods. Moreover, an increasing level of control over anisotropic interactions between these colloidal building blocks, which can be tailored through local geometries of sub-particle features as well as site-specific surface modifications, is giving rise to new demonstrations of massively parallel off-chip self-assembly of specific target structures with low defect rates. In particular, new experimental realizations of hierarchical self-assembly and control over the chiral purity of resulting assembly structures have been achieved. Increasingly, shape-dependent, shape-complementary, and roughness-controlled depletion attractions between non-spherical colloids are being used in novel ways to create assemblies that go far beyond early examples, such as fractal clusters formed by diffusion-limited and reaction-limited aggregation of spheres. As self-assembly methods have progressed, a wide variety of advanced directed assembly methods have also been developed; approaches based on microfluidic control and applying structured electromagnetic fields are particularly promising.

Nonequilibrium phase transitions of sheared colloidal microphases: Results from dynamical density functional theory

By means of classical density functional theory and its dynamical extension, we consider a colloidal fluid with spherically symmetric competing interactions, which are well known to exhibit a rich bulk phase behavior. This includes complex three-dimensional periodically ordered cluster phases such as lamellae, two-dimensional hexagonally packed cylinders, gyroid structures, or spherical micelles. While the bulk phase behavior has been studied extensively in earlier work, in this paper we focus on such structures confined between planar repulsive walls under shear flow. For sufficiently high shear rates, we observe that microphase separation can become fully suppressed. For lower shear rates, however, we find that, e.g., the gyroid structure undergoes a kinetic phase transition to a hexagonally packed cylindrical phase, which is found experimentally and theoretically in amphiphilic block copolymer systems. As such, besides the known similarities between the latter and colloidal systems regarding the equilibrium phase behavior, our work reveals further intriguing nonequilibrium relations between copolymer melts and colloidal fluids with competing interactions.

Flow induced crystallisation of penetrable particles

For a system of Brownian particles interacting via a soft exponential potential we investigate the interaction between equilibrium crystallisation and spatially varying shear flow. For thermodynamic state points within the liquid part of the phase diagram, but close to the crystallisation phase boundary, we observe that imposing a Poiseuille flow can induce nonequilibrium crystalline ordering in regions of low shear gradient. The physical mechanism responsible for this phenomenon is shear-induced particle migration, which causes particles to drift preferentially towards the center of the flow channel, thus increasing the local density in the channel center. The method employed is classical dynamical density functional theory.

Effect of particle-size dynamics on properties of dense spongy-particle systems: Approach towards equilibrium

Open-porous deformable particles, often envisaged as sponges, are ubiquitous in biological and industrial systems (e.g., casein micelles in dairy products and microgels in cosmetics). The rich behavior of these suspensions is owing to the elasticity of the supporting network of the particle, and the viscosity of permeating solvent. Therefore, the rate-dependent size change of these particles depends on their structure, i.e., the permeability. This work aims at investigating the effect of the particle-size dynamics and the underlying particle structure, i.e., the particle permeability, on the transient and long-time behavior of suspensions of spongy particles in the absence of applied deformation, using the dynamic two-scale model developed by Hütter et al. [Farad. Discuss. 158, 407 (2012)1359-664010.1039/c2fd20025b]. In the high-density limit, the transient behavior is found to be accelerated by the particle-size dynamics, even at average size changes as small as 1%. The accelerated dynamics is evidenced by (i) the higher short-time diffusion coefficient as compared to elastic-particle systems and (ii) the accelerated formation of the stable fcc crystal structure. Furthermore, after long times, the particle-size dynamics of spongy particles is shown to result in lower stationary values of the energy and normal stresses as compared to elastic-particle systems. This dependence of the long-time behavior of these systems on the permeability, that essentially is a transport coefficient and hence must not affect the equilibrium properties, confirms that full equilibration has not been reached.

Unraveling the Role of Order-to-Disorder Transition in Shear Thickening Suspensions

Using high-resolution in situ small angle x-ray scattering in conjunction with oscillatory shear on highly monodisperse silica suspensions, we demonstrate that an order-to-disorder transition leads to a dynamic shear thickening in a lower stress regime than the standard steady shear thickening. We show that the order-to-disorder transition is controlled by strain, which is distinguishably different from steady shear thickening, which is a stress-related phenomenon. The appearance of this two-step shear thinning and thickening transition is also influenced by the particle size, monodispersity, and measurement conditions (i.e., oscillatory shear versus steady shear). Our results show definitively that the order-to-disorder transition-induced thickening is completely unrelated to the mechanism that drives steady shear thickening.

Chiral symmetry breaking due to impeller size in cooling crystallization of sodium chlorate

We report, for the first time, that the size of impellers is the exception to Denk and Botsaris's work where a stirrer of the smallest surface area led to close to 100% deracemization, while the data from the largest impeller agreed with the results of Denk and Botsaris, displaying a racemic mixture.

Dynamical density functional theory analysis of the laning instability in sheared soft matter

Using dynamical density functional theory (DDFT) methods we investigate the laning instability of a sheared colloidal suspension. The nonequilibrium ordering at the laning transition is driven by nonaffine particle motion arising from interparticle interactions. Starting from a DDFT which incorporates the nonaffine motion, we perform a linear stability analysis that enables identification of the regions of parameter space where lanes form. We illustrate our general approach by applying it to a simple one-component fluid of soft penetrable particles.

Observation of Steady Shear-Induced Nematic Ordering of Selenium Chain Moieties in Arsenic Selenide Liquids

Structural anisotropy induced by steady shear and its mechanistic relation with shear thinning are investigated in As_xSe_100-x glasses (5 ≤ x ≤ 30) quenched from parent liquids subjected to shear rates ranging between 0 and 10⁴ s^-1 using polarized Raman spectroscopy and two-dimensional X-ray diffraction. When taken together, the results demonstrate significant shear-induced partial alignment of -Se-Se-Se- chain moieties in the flow direction of the extruded fibers. This alignment is reminiscent of nematic liquid crystals where orientational order exists without positional order. The degree of this structural alignment in quenched glasses appears to be practically independent of the shear rate, although the parent liquids undergo shear thinning at the highest shear rates. It is conjectured that any causal relationship between structural alignment and shear thinning in the liquid may be masked in the glassy state by the postextrusion structural relaxation of the parent liquid.

Generating Bulk-Scale Ordered Optical Materials Using Shear-Assembly in Viscoelastic Media

We review recent advances in the generation of photonics materials over large areas and volumes, using the paradigm of shear-induced ordering of composite polymer nanoparticles. The hard-core/soft-shell design of these particles produces quasi-solid “gum-like” media, with a viscoelastic ensemble response to applied shear, in marked contrast to the behavior seen in colloidal and granular systems. Applying an oscillatory shearing method to sub-micron spherical nanoparticles gives elastomeric photonic crystals (or “polymer opals”) with intense tunable structural color. The further engineering of this shear-ordering using a controllable “roll-to-roll” process known as Bending Induced Oscillatory Shear (BIOS), together with the interchangeable nature of the base composite particles, opens potentially transformative possibilities for mass manufacture of nano-ordered materials, including advances in optical materials, photonics, and metamaterials/plasmonics.

Molecular Dynamics Simulation of Trimer Self-Assembly Under Shear

The self-assembly of patchy trimer particles consisting of one attractive site and two repulsive sites is investigated with nonequilibrium molecular dynamics simulations in the presence of a velocity gradient, as would be produced by the application of a shear stress on the system. As shear is increased, globular-shaped micellar clusters increase in size and become more elongated. The globular clusters are also more stable at higher temperatures in the presence of shear than at equilibrium. These results help to increase our understanding of the effect of shear on self-assembly for a variety of applications.

Colloidal polycrystalline monolayers under oscillatory shear

In this paper we probe the structural response to oscillatory shear deformations of polycrystalline monolayers of soft repulsive colloids with varying area fraction over a broad range of frequencies and amplitudes. The particles are confined at a fluid interface, sheared using a magnetic microdisk, and imaged through optical microscopy. The structural and mechanical response of soft materials is highly dependent on their microstructure. If crystals are well understood and deform through the creation and mobilization of specific defects, the situation is much more complex for disordered jammed materials, where identifying structural motifs defining plastically rearranging regions remains an elusive task. Our materials fall between these two classes and allow the identification of clear pathways for structural evolution. In particular, we demonstrate that large enough strains are able to fluidize the system, identifying critical strains that fulfill a local Lindemann criterion. Conversely, smaller strains lead to localized and erratic irreversible particle rearrangements due to the motion of structural defects. In this regime, oscillatory shear promotes defect annealing and leads to the growth of large crystalline domains. Numerical simulations help identify the population of rearranging particles with those exhibiting the largest deviatoric stresses and indicate that structural evolution proceeds towards the minimization of the stress stored in the system. The particles showing high deviatoric stresses are localized around grain boundaries and defects, providing a simple criterion to spot regions likely to rearrange plastically under oscillatory shear.

Self-assembly of colloidal micelles in microfluidic channels

The self-assembly of amphiphilic Janus colloids in microfluidic channels is studied using hybrid molecular dynamics simulations with fully resolved hydrodynamic interactions incorporated through the multi-particle collision dynamics algorithm. The simulations are conducted at a density and temperature where the Janus particles spontaneously self-assemble into spherical micelles to minimize the interface between the solvophobic caps and the surrounding solvent. In confined systems, this contact area can also be reduced by aggregation at the channel walls. Indeed, a sizable fraction of free particles and small clusters with three and four members are found at the walls when the microfluidic channel is made up of a comparably solvophobic material as the Janus colloids. When the applied Poiseuille flow is sufficiently strong, the colloidal micelles break up into smaller fragments and isolated particles. However, at intermediate flow rates the shear-induced dissociation and reorganization of aggregates lead to a net growth of the micelles with a sizable amount of particles in icosahedral clusters with 13 particles. Furthermore, the parabolic velocity profile of the flow causes a highly non-uniform cluster size distribution between the channel walls, where the aggregation number decreases close to the walls.

Ultrastretchable Iono-Elastomers with Mechanoelectrical Response

The emerging technologies involving wearable electronics require new materials with high stretchability, resistance to high loads, and high conductivities. We report a facile synthetic strategy based on self-assembly of concentrated solutions of end-functionalized PEO₁₀₆-PPO₇₀-PEO₁₀₆ triblock copolymer in ethylammonium nitrate into face-centered cubic micellar crystals, followed by micelle corona cross-linking to generate elastomeric ion gels (iono-elastomers). These materials exhibit an unprecedented combination of high stretchability, high ionic conductivity, and mechanoelectrical response. The latter consists of a remarkable and counterintuitive increase in ion conductivity with strain during uniaxial extension, which is reversible upon load release. Based on in situ SAXS measurements of reversible crystal structure transformations during deformation, we postulate that the origin of the conductivity increase is a reversible formation of ion nanochannels due to a novel microstructural rearrangement specific to this material.

Large-scale ordering of nanoparticles using viscoelastic shear processing

Despite the availability of elaborate varieties of nanoparticles, their assembly into regular superstructures and photonic materials remains challenging. Here we show how flexible films of stacked polymer nanoparticles can be directly assembled in a roll-to-roll process using a bending-induced oscillatory shear technique. For sub-micron spherical nanoparticles, this gives elastomeric photonic crystals termed polymer opals showing extremely strong tunable structural colour. With oscillatory strain amplitudes of 300%, crystallization initiates at the wall and develops quickly across the bulk within only five oscillations. The resulting structure of random hexagonal close-packed layers is improved by shearing bidirectionally, alternating between two in-plane directions. Our theoretical framework indicates how the reduction in shear viscosity with increasing order of each layer accounts for these results, even when diffusion is totally absent. This general principle of shear ordering in viscoelastic media opens the way to manufacturable photonic materials, and forms a generic tool for ordering nanoparticles.

Dynamic, Directed Self-Assembly of Nanoparticles via Toggled Interactions

Crystals self-assembled from nanoparticles have useful properties such as optical activity and sensing capability. During fabrication, however, gelation and glassification often leave these materials arrested in defective or disordered metastable states. This is a key difficulty preventing adoption of self-assembled nanoparticle materials at scale. Processes which suppress kinetic arrest and defect formation while accelerating growth of ordered materials are essential for bottom-up approaches to creating nanomaterials. Dynamic, directed self-assembly processes in which the interactions between self-assembling components are actuated temporally offer one promising methodology for accelerating and controlling bottom-up growth of nanostructures. In this article, we show through simulation and theory how time-dependent, periodically toggled interparticle attractions can avoid kinetic barriers and yield well-ordered crystalline domains for a dispersion of nanoparticles interacting via a short-ranged, isotropic potential. The growth mechanism and terminal structure of the dispersion are controlled by parameters of the toggling protocol. This control allows for selection of processes that yield rapid self-assembled, low defect crystals. Although self-assembly via periodically toggled attractions is inherently unsteady and out-of-equilibrium, its outcome is predicted by a first-principles theory of nonequilibrium thermodynamics. The theory necessitates equality of the time average of pressure and chemical potential in coexisting phases of the dispersion. These quantities are evaluated using well known equations of state. The phase behavior predicted by this theory agrees well with measurements made in Brownian dynamics simulations of sedimentation equilibrium and homogeneous nucleation. The theory can easily be extended to model dynamic self-assembly directed by other toggled conservative force fields.

Formation of crystal-like structures and branched networks from nonionic spherical micelles

Crystal-like structures at nano and micron scales have promise for purification and confined reactions, and as starting points for fabricating highly ordered crystals for protein engineering and drug discovery applications. However, developing controlled crystallization techniques from batch processes remain challenging. We show that neutrally charged nanoscale spherical micelles from biocompatible nonionic surfactant solutions can evolve into nano- and micro-sized branched networks and crystal-like structures. This occurs under simple combinations of temperature and flow conditions. Our findings not only suggest new opportunities for developing controlled universal crystallization and encapsulation procedures that are sensitive to ionic environments and high temperatures, but also open up new pathways for accelerating drug discovery processes, which are of tremendous interest to pharmaceutical and biotechnological industries.

Nonequilibrium structure of colloidal dumbbells under oscillatory shear

We investigate the nonequilibrium behavior of dense, plastic-crystalline suspensions of mildly anisotropic colloidal hard dumbbells under the action of an oscillatory shear field by employing Brownian dynamics computer simulations. In particular, we extend previous investigations, where we uncovered nonequilibrium phase transitions, to other aspect ratios and to a larger nonequilibrium parameter space, that is, a wider range of strains and shear frequencies. We compare and discuss selected results in the context of scattering and rheological experiments. Both simulations and experiments demonstrate that the previously found transitions from the plastic crystal phase with increasing shear strain also occur at other aspect ratios. We explore the transition behavior in the strain-frequency phase and summarize it in a nonequilibrium phase diagram. Additionally, the experimental rheology results hint at a slowing down of the colloidal dynamics with higher aspect ratio.

Real-time measurements of crystallization processes in viscoelastic polymeric photonic crystals

We present a study of the dynamic shear ordering of viscoelastic photonic crystals, based on core-shell polymeric composite particles. Using an adapted shear-cell arrangement, the crystalline ordering of the material under conditions of oscillatory shear is interrogated in real time, through both video imaging and from the optical transmission spectra of the cell. In order to gain a deeper understanding of the macroscopic influences of shear on the crystallization process in this solvent-free system, the development of bulk ordering is studied as a function of the key parameters including duty cycle and shear-strain magnitude. In particular, optimal ordering is observed from a prerandomized sample at shear strains of around 160%, for 1-Hz oscillations. This ordering reaches completion over time scales of order 10 s. These observations suggest significant local strains are needed to drive nanoparticles through energy barriers, and that local creep is needed to break temporal symmetry in such high-viscosity nanoassemblies. Crystal shear-melting effects are also characterized under conditions of constant shear rate. These quantitative experiments aim to stimulate the development of theoretical models which can deal with the strong local particle interactions in this system.

Tuning colloidal gels by shear

Using a powerful combination of experiments and simulations we demonstrate how the microstructure and its time evolution are linked with mechanical properties in a frustrated, out-of-equilibrium, particle gel under shear. An intermediate volume fraction colloid-polymer gel is used as a model system, allowing quantification of the interplay between interparticle attractions and shear forces. Rheometry, confocal microscopy and Brownian dynamics reveal that high shear rates, fully breaking the structure, lead after shear cessation to more homogeneous and stronger gels, whereas preshear at low rates creates largely heterogeneous weaker gels with reduced elasticity. We find that in comparison, thermal quenching cannot produce structural inhomogeneities under shear. We argue that external shear has strong implications on routes towards metastable equilibrium, and therefore gelation scenarios. Moreover, these results have strong implications for material design and industrial applications, such as mixing, processing and transport protocols coupled to the properties of the final material.

Colloidal Plastic Crystals in a Shear Field

We study the structure and viscoelastic behavior of 3D plastic crystals of colloidal dumbbells in an oscillatory shear field based on a combination of small-angle neutron scattering experiments under shear (rheo-SANS) and Brownian dynamics computer simulations. Sterically stabilized dumbbell-shaped microgels are used as hard dumbbell model systems which consist of dumbbell-shaped polystyrene (PS) cores and thermosensitive poly(N-isopropylacrylamide) (PNIPAM) shells. Under increasing shear strain, a discontinuous transition is found from a twinned-fcc-like crystal to a partially oriented sliding-layer phase with a shear-molten state in between. In the novel partially oriented sliding-layer phase, the hard dumbbells exhibit a small but finite orientational order in the shear direction. We find that this weak correlation is sufficient to perturb the nature of the nonequilibrium phase transition as known for hard sphere systems. The discontinuous transition for hard dumbbells is observed to be accompanied by a novel yielding process with two yielding events in its viscoelastic shear response, while only a single yielding event is observed for sheared hard spheres. Our findings will be useful in interpreting the shear response of anisotropic colloidal systems and in generating novel colloidal crystals from anisotropic systems with applications in colloidal photonics.

Rheology and microstructure of non-Brownian suspensions in the liquid and crystal coexistence region: strain stiffening in large amplitude oscillatory shear

Concentrated hard-sphere suspensions in the liquid and crystal coexistence region show a unique nonlinear behavior under a large amplitude oscillatory shear flow, the so-called strain stiffening, in which the viscosity or modulus suddenly starts to increase near a critical strain amplitude. Even though this phenomenon has been widely reported in experiments, its key mechanism has never been investigated in a systematic way. To have a good understanding of this behavior, a numerical simulation was performed using the lattice Boltzmann method (LBM). Strain stiffening was clearly observed at large strain amplitudes, and the critical strain amplitude showed an angular frequency dependency. The distortion of the shear stress appeared near the critical strain amplitude, and the nonlinear behavior was quantified by the Fourier transformation (FT) and the stress decomposition methods. Above the critical strain amplitude, an increase in the global bond order parameter Ψ(6) was observed at the flow reversal. The maximum of Ψ(6) and the maximum shear stress occurred at the same strain. These results show how strongly the ordered structure of the particles is related to the stress distortion. The ordered particles maintained a bond number of "two" with alignment with the compressive axis, and they were distributed over a narrow range of angular distribution (110°-130°). In addition, the ordered structure was formed near the lowest shear rate region (the flow reversal). The characteristics of the ordered structure were remarkably different from those of the hydroclusters which are regarded as the origin of shear thickening. It is clear that strain stiffening and shear thickening originate from different mechanisms. Our results clearly demonstrate how the ordering of the particles induces strain stiffening in the liquid and crystal coexistence region.

Shear-Triggered Crystallization and Light Emission of a Thermally Stable Organic Supercooled Liquid

Thermodynamics drive crystalline organic molecules to be crystallized at temperatures below their melting point. Even though molecules can form supercooled liquids by rapid cooling, crystalline organic materials readily undergo a phase transformation to an energetically favorable crystalline phase upon subsequent heat treatment. Opposite to this general observation, here, we report molecular design of thermally stable supercooled liquid of diketopyrrolopyrrole (DPP) derivatives and their intriguing shear-triggered crystallization with dramatic optical property changes. Molten DPP8, one of the DPP derivatives, remains as stable supercooled liquid without crystallization through subsequent thermal cycles. More interestingly, under shear conditions, this supercooled liquid DPP8 transforms to its crystal phase accompanied by a 25-fold increase in photoluminescence (PL) quantum efficiency and a color change. By systematic investigation on supercooled liquid formation of crystalline DPP derivatives and their correlation with chemical structures, we reveal that the origin of this thermally stable supercooled liquid is a subtle force balance between aromatic interactions among the core units and van der Waals interactions among the aliphatic side chains acting in opposite directions. Moreover, by applying shear force to a supercooled liquid DPP8 film at different temperatures, we demonstrated direct writing of fluorescent patterns and propagating fluorescence amplification, respectively. Shear-triggered crystallization of DPP8 is further achieved even by living cell attachment and spreading, demonstrating the high sensitivity of the shear-triggered crystallization which is about 6 orders of magnitude more sensitive than typical mechanochromism observed in organic materials.

Kinetics of colloidal deposition, assembly, and crystallization in steady electric fields

We quantify and model the deposition and crystallization kinetics of initially dilute colloidal spheres due to application of a steady, direct current electric field in the thin gap between parallel electrodes. The system studied is poly(12-hydroxystearic acid) (PHSA)-stabilized poly(methyl methacrylate) (PMMA) spheres dispersed in a mixture of cyclohexylbromide (CHB), decalin, and a low concentration of the partially disassociating salt tetrabutylammonium chloride (TBAC). The temporal and spatial evolution of the colloidal volume fraction in the ∼1 mm gap between the electrodes is quantified under conditions of both deposition and relaxation by confocal laser scanning microscopy (CLSM). During deposition assembly, the spatial dependence of the colloid volume fraction approaches steady state at times between hundreds of minutes at the lowest electric field strength (as characterized by a Peclet number, Pe) and at tens of minutes at higher field strengths. During disassembly, the volume fraction relaxes nearly exponentially. The kinetics are modeled by adapting a treatment for sedimentation (Davis and Russel, Phys. Fluids A, 1989, 1, 82) to the case of steady electric fields. The model's predictions show good agreement with the measured kinetics at low Pe; however, agreement progressively deteriorates with increasing Pe. At low Pe the deposits are initially disordered. After an initial delay, 1D crystal growth propagates from the electrode surface at rates of several hundred nm min(-1). The sharp crystal boundary propagates as a characteristic of constant colloidal volume fraction, consistent with an equilibrium crystalline phase transition. The results inform operational ranges for devices that produce active colloidal matter by reversible assembly.

Large-scale high-quality 2D silica crystals: dip-drawing formation and decoration with gold nanorods and nanospheres for SERS analysis

High-quality colloidal crystals (CCs) are important for use in photonic research and as templates for large-scale plasmonic SERS substrates. We investigated how variations in temperature, colloid concentration, and dip-drawing parameters (rate, incubation time, etc) affect the structure of 2D CCs formed by highly monodisperse silica nanoparticles (SiNPs) synthesized in an l-arginine solution and regrown by a modified Stöber method. The best quality 2D CCs were obtained with aqueous 12 wt% colloids at a temperature of 25 °C, an incubation time of 1 min, and a drawing rate of 50 mm min(-1). Assembling of gold nanorods (GNRs) on 2D CCs resulted in the formation of ring-like chains with a preferential tail-to-tail orientation along the hexagonal boundaries. To the best of our knowledge, this is the first time that such nanostructures have been prepared. Owing to the preferential tail-to-tail packing of GNRs, 2D SiNP CC + GNR substrates demonstrated an analytical SERS enhancement of about 8000, which was 10 to 15 times higher than that for self-assembled GNRs on a silicon wafer. In addition, the analytical SERS enhancement was almost 60 times lower after replacing the nanorods in 2D SiNP CC + GNR substrates with 25 nm gold nanospheres.

Rheology of dense suspensions of elastic capsules: normal stresses, yield stress, jamming and confinement effects

We study the shearing rheology of dense suspensions of elastic capsules, taking aggregation-free red blood cells as a physiologically relevant example. Particles are non-Brownian and interact only via hydrodynamics and short-range repulsive forces. An analysis of the different stress mechanisms in the suspension shows that the viscosity is governed by the shear elasticity of the capsules, whereas the repulsive forces are subdominant. Evidence for a dynamic yield stress above a critical volume fraction is provided and related to the elastic properties of the capsules. The shear stress is found to follow a critical jamming scenario and is rather insensitive to the tumbling-to-tank-treading transition. The particle pressure and normal stress differences display some sensitivity to the dynamical state of the cells and exhibit a characteristic scaling, following the behavior of a single particle, in the tank-treading regime. The behavior of the viscosity in the fluid phase is rationalized in terms of effective medium models. Furthermore, the role of confinement effects, which increase the overall magnitude and enhance the shear-thinning of the viscosity, is discussed.