Dynamical and rheological properties of soft colloid suspensions

Nonequilibrium relaxation of soft responsive colloids

Stimuli-responsive macromolecules display large conformational changes during their dynamics, sometimes switching between states. Such a multi-stability is useful for the development of soft functional materials. Here, we introduce a mean-field dynamical density functional theory for a model of responsive colloids to study the nonequilibrium dynamics of a colloidal dispersion in time-dependent external fields, with a focus on the coupling of translational and conformational dynamics during their relaxation. Specifically, we consider soft Gaussian particles with a bimodal size distribution between two confining walls with time-dependent (switching-on and off) external gravitational and osmotic fields. We find a rich relaxation behavior of the systems in excellent agreement with particle-based Brownian dynamics computer simulations. In particular, we find time-asymmetric relaxations of integrated observables (wall pressures, mean size, and liquid center-of-mass) for activation/deactivation of external potentials, respectively, which are tunable by the ratio of translational and conformational diffusion time scales. Our work thus paves the way for studying the nonequilibrium relaxation dynamics of complex soft matter with multiple degrees of freedom and hierarchical relaxations.

Novel animal product substitutes: A new category of plant-based alternatives to meat, seafood, egg, and dairy products

Many consumers are adopting plant-centric diets to address the adverse effects of livestock production on the environment, health, and animal welfare. Processed plant-based foods, including animal product analogs (such as meat, seafood, egg, or dairy analogs) and traditional animal product substitutes (such as tofu, seitan, or tempeh), may not be desirable to a broad spectrum of consumers. This article introduces a new category of plant-based foods specifically designed to overcome the limitations of current animal product analogs and substitutes: novel animal product substitutes (NAPS). NAPS are designed to contain high levels of nutrients to be encouraged (such as proteins, omega-3 fatty acids, dietary fibers, vitamins, and minerals) and low levels of nutrients to be discouraged (such as salt, sugar, and saturated fat). Moreover, they may be designed to have a wide range of appearances, textures, mouthfeels, and flavors. For instance, they could be red, orange, green, yellow, blue, or beige; they could be spheres, ovals, cubes, or pyramids; they could be hard/soft or brittle/pliable; and they could be lemon, thyme, curry, or chili flavored. Consequently, there is great flexibility in creating NAPS that could be eaten in situations where animal products are normally consumed, for example, with pasta, rice, potatoes, bread, soups, or salads. This article reviews the science behind the formulation of NAPS, highlights factors impacting their appearance, texture, flavor, and nutritional profile, and discusses methods that can be used to formulate, produce, and characterize them. Finally, it stresses the need for further studies on this new category of foods, especially on their sensory and consumer aspects.

External field-driven property localization in liquids of responsive macromolecules

We explore theoretically the effects of external potentials on the spatial distribution of particle properties in a liquid of explicitly responsive macromolecules. In particular, we focus on the bistable particle size as a coarse-grained internal degree of freedom (DoF, or "property"), σ, that moves in a bimodal energy landscape, in order to model the response of a state-switching (big-to-small) macromolecular liquid to external stimuli. We employ a mean-field density functional theory (DFT) that provides the full inhomogeneous equilibrium distributions of a one-component model system of responsive colloids (RCs) interacting with a Gaussian pair potential. For systems confined between two parallel hard walls, we observe and rationalize a significant localization of the big particle state close to the walls, with pressures described by an exact RC wall theorem. Application of more complex external potentials, such as linear (gravitational), osmotic, and Hamaker potentials, promotes even stronger particle size segregation, in which macromolecules of different size are localized in different spatial regions. Importantly, we demonstrate how the degree of responsiveness of the particle size and its coupling to the external potential tune the position-dependent size distribution. The DFT predictions are corroborated by Brownian dynamics simulations. Our study highlights the fact that particle responsiveness can be used to localize liquid properties and therefore helps to control the property- and position-dependent function of macromolecules, e.g., in biomedical applications.

Loquat seed starch - Emerging source of non-conventional starch: Structure, properties, and novel applications

Recently, non-conventional sources of starch have attracted attention due to their potential to provide cost-effective alternatives to traditional starch. Among non-conventional starches, loquat (Eriobotrya japonica) seed starch is an emerging source of starch consisting of the amount of starch (nearly 20 %). Due to its unique structure, functional properties, and novel applications, it could be utilized as a potential ingredient. Interestingly, this starch has similar properties as commercial starches including high amylose content, small granule size, and high viscosity and heat stability, making it an attractive option for various food applications. Therefore, this review mainly covers the fundamental understanding of the valorization of loquat seeds by extracting the starch using different isolation methods, with preferable structural, morphological, and functional properties. Different isolation and modification methods (wet milling, acid, neutral and alkaline) are effectively used to obtain higher amounts of starch are revealed. Moreover, insight into various analytical techniques including scanning electron microscopy, differential scanning calorimetry, and X-ray diffraction used to characterize the molecular structure of the starch are discussed. In addition, the effect of shear rate and temperature on rheological attributes with solubility index, swelling power, and color is revealed. Besides, this starch contains bioactive compounds that have shown a positive impact on the enhancement of the shelf-life of the fruits. Overall, loquat seed starches have the potential to provide sustainable and cost-effective alternatives to traditional starch sources and can lead to novel applications in the food industry. Further research is needed to optimize processing techniques and develop value-added products that can be produced at a large scale. However, there is relatively limited published scientific evidence on the structural and morphological characteristics of loquat seed starch. Thus, in this review, we focused on different isolation techniques of loquat seed starch, its structural and functional characteristics, along with potential applications.

Rheology of a Nanopolymer Synthesized through Directional Assembly of DNA Nanochambers, for Magnetic Applications

We present a numerical study of the effects of monomer shape and magnetic nature of colloids on the behavior of a single magnetic filament subjected to the simultaneous action of shear flow and a stationary external magnetic field perpendicular to the flow. We find that based on the magnetic nature of monomers, magnetic filaments exhibit a completely different phenomenology. Applying an external magnetic field strongly inhibits tumbling only for filaments with ferromagnetic monomers. Filament orientation with respect to the flow direction is in this case independent of monomer shape. In contrast, reorientational dynamics in filaments with superparamagnetic monomers are not inhibited by applied magnetic fields, but enhanced. We find that the filaments with spherical, superparamagnetic monomers, depending on the flow and external magnetic field strength, assume semipersistent, collapsed, coiled conformations, and their characteristic time of tumbling is a function of field strength. However, external magnetic fields do not affect the characteristic time of tumbling for filaments with cubic, superparamagnetic monomers, but increase how often tumbling occurs.

Active interaction switching controls the dynamic heterogeneity of soft colloidal dispersions

We employ Reactive Dynamical Density Functional Theory (R-DDFT) and Reactive Brownian Dynamics (R-BD) simulations to investigate the dynamics of a suspension of active soft Gaussian colloids with binary interaction switching, i.e., a one-component colloidal system in which every particle stochastically switches at predefined rates between two interaction states with different mobility. Using R-DDFT we extend a theory previously developed to access the dynamics of inhomogeneous liquids [Archer et al., Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2007, 75, 040501] to study the influence of the switching activity on the self and distinct part of the Van Hove function in bulk solution, and determine the corresponding mean squared displacement of the switching particles. Our results demonstrate that, even though the average diffusion coefficient is not affected by the switching activity, it significantly modifies the non-equilibrium dynamics and diffusion coefficients of the individual particles, leading to a crossover from short to long times, with a regime for intermediate times showing anomalous diffusion. In addition, the self-part of the van Hove function has a Gaussian form at short and long times, but becomes non-Gaussian at intermediates ones, having a crossover between short and large displacements. The corresponding self-intermediate scattering function shows the two-step relaxation patters typically observed in soft materials with heterogeneous dynamics such as glasses and gels. We also introduce a phenomenological Continuous Time Random Walk (CTRW) theory to understand the heterogeneous diffusion of this system. R-DDFT results are in excellent agreement with R-BD simulations and the analytical predictions of CTRW theory, thus confirming that R-DDFT constitutes a powerful method to investigate not only the structure and phase behavior, but also the dynamical properties of non-equilibrium active switching colloidal suspensions.

Active binary switching of soft colloids: stability and structural properties

We employ reactive dynamical density functional theory (R-DDFT) and reactive Brownian dynamics (R-BD) simulations to study the non-equilibrium structure and phase behavior of an active dispersion of soft Gaussian colloids with binary interaction switching, i.e., we consider a one-component colloidal system in which every particle can individually switch stochastically between two interaction states (here, sizes 'big' and 'small') at predefined rates. We consider the influence of switching activity on the inhomogeneous density profiles of the colloids confined by various external potentials, as well as on their pair structure and phase behavior in bulk solutions. For the latter, we extend the R-DDFT method to incorporate the Percus test-particle route. Our results demonstrate that switching activity strongly modifies the steady-state density profiles and structural (pair) correlations. In particular, the switching rate interpolates from a near-equilibrium binary colloidal mixture of two states at very low rates to a non-equilibrium, 'one-state liquid' at very high rates characterized by one, average interaction size. The latter limit can be described by an equivalent effective one-component (EOC) equilibrium system, for which the exact analytical expression for the effective pair potential is a diffusion-weighted superposition of the active systems' pair potentials. This leads to the interesting fact that under certain conditions an interacting switching system can behave like a non-interacting (ideal) gas in the limit of high switching rates. Moreover, for colloids that are unstable (i.e., demix) near equilibrium, we demonstrate that phase separation and micro-clustering in both confinement and bulk can be dynamically controlled by the switching rate, and vanish for high rates. All R-DDFT results are in excellent agreement with our R-BD simulations.

Modeling Red Blood Cell Viscosity Contrast Using Inner Soft Particle Suspension

The inner viscosity of a biological red blood cell is about five times larger than the viscosity of the blood plasma. In this work, we use dissipative particles to enable the proper viscosity contrast in a mesh-based red blood cell model. Each soft particle represents a coarse-grained virtual cluster of hemoglobin proteins contained in the cytosol of the red blood cell. The particle interactions are governed by conservative and dissipative forces. The conservative forces have purely repulsive character, whereas the dissipative forces depend on the relative velocity between the particles. We design two computational experiments that mimic the classical viscometers. With these experiments we study the effects of particle suspension parameters on the inner cell viscosity and provide parameter sets that result in the correct viscosity contrast. The results are validated with both static and dynamic biological experiment, showing an improvement in the accuracy of the original model without major increase in computational complexity.

Structure and dynamics of responsive colloids with dynamical polydispersity

Dynamical polydispersity in single-particle properties, for example a fluctuating particle size, shape, charge density, etc, is intrinsic to responsive colloids (RCs), such as biomacromolecules or microgels, but is typically not resolved in coarse-grained mesoscale simulations. Here, we present Brownian dynamics simulations of suspensions of RCs modeling soft hydrogel colloids, for which the size of the individual particles is an explicitly resolved (Gaussian) degree of freedom and dynamically responds to the local interacting environment. We calculate the liquid structure, emergent size distributions, long-time diffusion, and property (size) relaxation kinetics for a wide range of densities and intrinsic property relaxation times in the canonical ensemble. Comparison to interesting reference cases, such as conventional polydisperse suspensions with a frozen parent distribution, or conventional monodisperse systems interacting with an effective pair potential for one fixed size, shows a significant spread in the structure and dynamics. The differences, most apparent in the high density regimes, are due to many-body correlations and the dynamical coupling between property and translation in RC systems, not explicitly accounted for in the conventional treatments. In particular, the translational diffusion in the RC systems is surprisingly close to the free (single RC) diffusion, mainly due to a cancellation of crowding and size compression effects. We show that an effective monodisperse pair potential can be constructed that describes the many-body correlations reasonably well by convoluting the RC pair potential with the density-dependent emergent size distributions and using a mean effective diffusion constant.

Onsager's variational principle in active soft matter

Onsagers variational principle (OVP) was originally proposed by Lars Onsager in 1931 [L. Onsager, Phys. Rev., 1931, 37, 405]. This fundamental principle provides a very powerful tool for formulating thermodynamically consistent models. It can also be employed to find approximate solutions, especially in the study of soft matter dynamics. In this work, OVP is extended and applied to the dynamic modeling of active soft matter such as suspensions of bacteria and aggregates of animal cells. We first extend the general formulation of OVP to active matter dynamics where active forces are included as external non-conservative forces. We then use OVP to analyze the directional motion of individual active units: a molecular motor walking on a stiff biofilament and a toy two-sphere microswimmer. Next we use OVP to formulate a diffuse-interface model for an active polar droplet on a solid substrate. In addition to the generalized hydrodynamic equations for active polar fluids in the bulk region, we have also derived thermodynamically consistent boundary conditions. Finally, we consider the dynamics of a thin active polar droplet under the lubrication approximation. We use OVP to derive a generalized thin film equation and then employ OVP as an approximation tool to find the spreading laws for the thin active polar droplet. By incorporating the activity of biological systems into OVP, we develop a general approach to construct thermodynamically consistent models for better understanding the emergent behaviors of individual animal cells and cell aggregates or tissues.

Structure and position-dependent properties of inhomogeneous suspensions of responsive colloids

Responsive particles, such as biomacromolecules or hydrogels, display a broad and polymodal distribution of conformations and have thus the ability to change their properties (e.g., size, shape, charge density, etc.) substantially in response to external fields or to their local environment (e.g., mediated by cosolutes or pH). Here we discuss the basic statistical mechanics for a model of responsive colloids (RCs) by introducing an additional "property" degree of freedom as a collective variable in a formal coarse-graining procedure. The latter leads to an additional one-body term in the coarse-grained (CG) free energy, defining a single-particle property distribution for an individual polydisperse RC. We argue that in the equilibrium thermodynamic limit such a CG system of RCs behaves like a conventional polydisperse system of nonresponsive particles. We then illustrate the action of external fields, which impose local (position-dependent) property distributions leading to nontrivial effects on the spatial one-body property and density profiles, even for an ideal (noninteracting) gas of RCs. We finally apply density-functional theory in the local density approximation to discuss the effects of particle interactions for specific examples of (i) a suspension of RCs in an external field linear in both position and property, (ii) a suspension of RCs with highly localized properties (sizes) confined between two walls, and (iii) a two-component suspension where an inhomogeneously distributed (nonresponsive) cosolute component, as found, e.g., in the studies of osmolyte- or salt-induced collapse or swelling transitions of thermosensitive polymers, modifies the local properties and density of the RC liquid.

Controlling the Microstructure and Phase Behavior of Confined Soft Colloids by Active Interaction Switching

We explore the microstructure and phase behavior of confined soft colloids which can actively switch their interactions at a predefined kinetic rate. For this, we employ a reactive dynamical density-functional theory and study the effect of a two-state switching of the size of colloids interacting with a Gaussian pair potential in the nonequilibrium steady state. The switching rate interpolates between a near-equilibrium binary mixture at low rates and a nonequilibrium monodisperse liquid for large rates, strongly affecting the one-body density profiles, adsorption, and pressure at confining walls. Importantly, we show that sufficiently fast switching impedes the phase separation of an (in equilibrium) unstable liquid, allowing the control of the degree of mixing and condensation and local microstructuring in a cellular confinement by tuning the switching rate.

Rheology of protein-stabilised emulsion gels envisioned as composite networks 1- Comparison of pure droplet gels and protein gels

HYPOTHESIS

Protein-stabilised emulsion gels can be studied in the theoretical framework of colloidal gels, because both protein assemblies and droplets may be considered as soft colloids. These particles differ in their nature, size and softness, and these differences may have an influence on the rheological properties of the gels they form.

EXPERIMENTS

Pure gels made of milk proteins (sodium caseinate), or of sub-micron protein-stabilised droplets, were prepared by slow acidification of suspensions at various concentrations. Their microstructure was characterised, their viscoelasticity, both in the linear and non-linear regime, and their frequency dependence were measured, and the behaviour of the two types of gels was compared.

FINDINGS

Protein gels and droplet gels were found to have broadly similar microstructure and rheological properties when compared at fixed volume fraction, a parameter derived from the study of the viscosity of the suspensions formed by proteins and by droplets. The viscoelasticity displayed a power law behaviour in concentration, as did the storage modulus in frequency. Additionally, strain hardening was found to occur at low concentration. These behaviours differed slightly between protein gels and droplet gels, showing that some specific properties of the primary colloidal particles play a role in the development of the rheological properties of the gels.

Dissipative particle dynamics with energy conservation: Isoenergetic integration and transport properties

Simulations of nano- to micro-meter scale fluidic systems under thermal gradients require consistent mesoscopic methods accounting for both hydrodynamic interactions and proper transport of energy. One such method is dissipative particle dynamics with energy conservation (DPDE), which has been used for various fluid systems with non-uniform temperature distributions. We propose an easily parallelizable modification of the velocity-Verlet algorithm based on local energy redistribution for each DPDE particle such that the total energy in a simulated system is conserved up to machine precision. Furthermore, transport properties of a DPDE fluid are analyzed in detail. In particular, an analytical approximation for the thermal conductivity coefficient is derived, which allows its a priori estimation for a given parameter set. Finally, we provide approximate expressions for the dimensionless Prandtl and Schmidt numbers, which characterize fluid transport properties and can be adjusted independently by a proper selection of model parameters. In conclusion, our results strengthen the DPDE method as a very robust approach for the investigation of mesoscopic systems with temperature inhomogeneities.

Reticulate collisional structure in boundary-driven granular gases

We report a peculiar head-on collision network between two vibrating boundaries in experiments performed during a parabolic flight and in a laboratory using horizontal vibration. This structure is a new ordering, which is due to an orientation correlation between the relative position and velocity of any particle pair. It weakens the collision frequency and produces a long-range boundary effect. Moreover, we find the molecular chaos assumption is violated in a larger portion of the phase space. Using an anisotropic distribution model, we modify angular integration results and compare them to the results of the kinetic theory.

Hydrodynamic correlations of viscoelastic fluids by multiparticle collision dynamics simulations

The emergent fluctuating hydrodynamics of a viscoelastic fluid modeled by the multiparticle collision dynamics (MPC) approach is studied. The fluid is composed of flexible, Gaussian phantom polymers that interact by local momentum-conserving stochastic MPCs. For comparison, the analytical solution of the linearized Navier-Stokes equation is calculated, where viscoelasticity is taken into account by a time-dependent shear relaxation modulus. The fluid properties are characterized by the transverse velocity autocorrelation function in Fourier space as well as in real space. Various polymer lengths are considered-from dumbbells to (near-)continuous polymers. Viscoelasticity affects the fluid properties and leads to strong correlations, which overall decay exponentially in Fourier space. In real space, the center-of-mass velocity autocorrelation function of individual polymers exhibits a long-time tail, independent of the polymer length, which decays as t^-3/2, similar to a Newtonian fluid, in the asymptotic limit t → ∞. Moreover, for long polymers, an additional power-law decay appears at time scales shorter than the longest polymer relaxation time with the same time dependence, but negative correlations, and the polymer length dependence L^-1/2. Good agreement is found between the analytical and simulation results.

Dynamic Equivalence between Soft Star Polymers and Hard Spheres

Understanding the dynamics of soft colloids, such as star polymers, dendrimers, and microgels, is of scientific and practical importance. It is known that the excluded volume effect plays a key role in colloidal dynamics. Here, we propose a condition of compressibility equivalence that provides a simple method to experimentally evaluate the excluded volume of soft colloids from a thermodynamic view. We apply this condition to survey the dynamics of a series of star polymer dispersions. It is found that, as the concentration increases, the slowing of the long-time self-diffusivity of the star polymer, normalized by the short-time self-diffusivity, can be mapped onto the hard-sphere behavior. This phenomenon reveals the dynamic equivalence between soft colloids and hard spheres, despite the apparent complexity of the interparticle interaction of the soft colloids. The methods for measuring the osmotic compressibility and the self-diffusivities of soft colloidal dispersions are also presented.

Multi-particle collision dynamics for a coarse-grained model of soft colloids

The growing interest in the dynamical properties of colloidal suspensions, both in equilibrium and under an external drive such as shear or pressure flow, requires the development of accurate methods to correctly include hydrodynamic effects due to the suspension in a solvent. In the present work, we generalize Multiparticle Collision Dynamics (MPCD) to be able to deal with soft, polymeric colloids. Our methods build on the knowledge of the monomer density profile that can be obtained from monomer-resolved simulations without hydrodynamics or from theoretical arguments. We hereby propose two different approaches. The first one simply extends the MPCD method by including in the simulations effective monomers with a given density profile, thus neglecting monomer-monomer interactions. The second one considers the macromolecule as a single penetrable soft colloid (PSC), which is permeated by an inhomogeneous distribution of solvent particles. By defining an appropriate set of rules to control the collision events between the solvent and the soft colloid, both linear and angular momenta are exchanged. We apply these methods to the case of linear chains and star polymers for varying monomer lengths and arm number, respectively, and compare the results for the dynamical properties with those obtained within monomer-resolved simulations. We find that the effective monomer method works well for linear chains, while the PSC method provides very good results for stars. These methods pave the way to extend MPCD treatments to complex macromolecular objects such as microgels or dendrimers and to work with soft colloids at finite concentrations.

Understanding conformational and dynamical evolution of semiflexible polymers in shear flow

A clear description of the conformational and dynamical evolution of polymer chains in shear flow is the fundamental basis of microfluidic separations and macroscopic rheological behaviors. We employ graph theory analysis to analyze the local deformation and dynamics of linear polymer chains with different rigidities in shear flow based on the simulation trajectories that record the instantaneous conformations and dynamics. Our results show that all semiflexible chains experience quasi-periodic tumbling motions when the shear strain overwhelms the U-shape (or S-shape) deformation energy barrier. More interestingly, the contact map provides solid evidence for the asymmetric deformation in the whole tumbling motion. In the stretching process: at small and intermediate shear strength, flexible polymers show a quasi-affine deformation while semiflexible ones are initially unfolded from the center of the chains, then both of them follow the extension with half dumbbell- or dumbbell-like ends; at high shear strength, all polymer chains present only a dumbbell-like extension. In the collapse process, all chains prefer to initiate the folding from chain ends. This finding can facilitate our understanding on how semiflexible polymer chains relax and dissipate the stress in shear flow.

Tethered Semiflexible Polymer under Large Amplitude Oscillatory Shear

The properties of a semiflexible polymer with fixed ends exposed to oscillatory shear flow are investigated by simulations. The two-dimensionally confined polymer is modeled as a linear bead-spring chain, and the interaction with the fluid is described by the Brownian multiparticle collision dynamics approach. For small shear rates, the tethering of the ends leads to a more-or-less linear oscillatory response. However, at high shear rates, we found a strongly nonlinear reaction, with a polymer (partially) wrapped around the fixation points. This leads to an overall shrinkage of the polymer. Dynamically, the location probability of the polymer center-of-mass position is largest on a spatial curve resembling a limaçon, although with an inhomogeneous distribution. We found shear-induced modifications of the normal-mode correlation functions, with a frequency doubling at high shear rates. Interestingly, an even-odd asymmetry for the Cartesian components of the correlation functions appears, with rather similar spectra for odd x- and even y-modes and vice versa. Overall, our simulations yielded an intriguing nonlinear behavior of tethered semiflexible polymers under oscillatory shear flow.

Morphology of depletant-induced erythrocyte aggregates

Red blood cells suspended in quiescent plasma tend to aggregate into multicellular assemblages, including linearly stacked columnar rouleaux, which can reversibly form more complex clusters or branching networks. While these aggregates play an essential role in establishing hemorheological and pathological properties, the biophysics behind their self-assembly into dynamic mesoscopic structures remains under-explored. We employ coarse-grained molecular simulations to model low-hematocrit erythrocytes subject to short-range implicit depletion forces, and demonstrate not only that depletion interactions are sufficient to account for a sudden dispersion-aggregate transition, but also that the volume fraction of depletant macromolecules controls small aggregate morphology. We observe a sudden transition from a dispersion to a linear column rouleau, followed by a slow emergence of disorderly amorphous clusters of many short rouleaux at larger volume fractions. This work demonstrates how discocyte topology and short-range, non-specific, physical interactions are sufficient to self-assemble erythrocytes into various aggregate structures, with markedly different morphologies and biomedical consequences.

Unusual polysaccharide rheology of aqueous dispersions of soft phytoglycogen nanoparticles

Phytoglycogen is a natural polysaccharide produced in the form of dense, 35 nm diameter nanoparticles by some varieties of plants such as sweet corn. The highly-branched, dendrimeric structure of phytoglycogen leads to interesting and useful properties such as softness and deformability of the particles, and a strong interaction with water. These properties make the particles ideal for use as unique additives in personal care, nutrition and biomedical formulations. In the present study, we describe rheology measurements of aqueous dispersions of phytoglycogen nanoparticles. The viscosity of the dispersions remained Newtonian up to large concentrations (∼20% w/w). For higher concentrations, the zero-shear viscosity increased dramatically, reaching values that exceeded that of the water solvent by six orders of magnitude at a concentration of 30% w/w and were well described by the Vogel-Fulcher-Tammann relation of glassy dynamics. The very large values of the zero-shear viscosity are coupled with significant deformation of the soft nanoparticles. We quantified the softness of the particles by performing osmotic pressure measurements on concentrated dispersions, obtaining a value of 15 kPa for the compressional modulus. For the most concentrated samples, we observed flow at stresses less than the apparent yield stress value determined by fitting the high strain rate data to the Herschel-Bulkley model. This behavior, similar to that of star polymer glasses, suggests the possibility of a hairy colloid particle geometry. Remarkably, phytoglycogen nanoparticles dispersed in water provide a very simple experimental realization of glass-forming dispersions of soft colloidal particles that can be used to validate theoretical models in detail.

Rotation Dynamics of Star Block Copolymers under Shear Flow

Star block-copolymers (SBCs) are macromolecules formed by a number of diblock copolymers anchored to a common central core, being the internal monomers solvophilic and the end monomers solvophobic. Recent studies have demonstrated that SBCs constitute self-assembling building blocks with specific softness, functionalization, shape and flexibility. Depending on different physical and chemical parameters, the SBCs can behave as flexible patchy particles. In this paper, we study the rotational dynamics of isolated SBCs using a hybrid mesoscale simulation technique. We compare three different approaches to analyze the dynamics: the laboratory frame, the non-inertial Eckart's frame and a geometrical approximation relating the conformation of the SBC to the velocity profile of the solvent. We find that the geometrical approach is adequate when dealing with very soft systems, while in the opposite extreme, the dynamics is best explained using the laboratory frame. On the other hand, the Eckart frame is found to be very general and to reproduced well both extreme cases. We also compare the rotational frequency and the kinetic energy with the definitions of the angular momentum and inertia tensor from recent publications.

Concentration-dependent swelling and structure of ionic microgels: simulation and theory of a coarse-grained model

We study swelling and structural properties of ionic microgel suspensions within a comprehensive coarse-grained model that combines the polymeric and colloidal natures of microgels as permeable, compressible, charged spheres governed by effective interparticle interactions. The model synthesizes the Flory-Rehner theory of cross-linked polymer gels, the Hertz continuum theory of effective elastic interactions, and a theory of density-dependent effective electrostatic interactions. Implementing the model using Monte Carlo simulation and thermodynamic perturbation theory, we compute equilibrium particle size distributions, swelling ratios, volume fractions, net valences, radial distribution functions, and static structure factors as functions of concentration. Trial Monte Carlo moves comprising particle displacements and size variations are accepted or rejected based on the total change in elastic and electrostatic energies. The theory combines first-order thermodynamic perturbation and variational free energy approximations. For illustrative system parameters, theory and simulation agree closely at concentrations ranging from dilute to beyond particle overlap. With increasing concentration, as microgels deswell, we predict a decrease in the net valence and an unusual saturation of pair correlations. Comparison with experimental data for deionized, aqueous suspensions of PNIPAM particles demonstrates the capacity of the coarse-grained model to predict and interpret measured swelling behavior.

Trefoil Knot Hydrodynamic Delocalization on Sheared Ring Polymers

The behavior of unknotted and trefoil-knotted ring polymers under shear flow is here examined by means of mesoscopic simulations. In contrast to most polymers, ring polymers in a hydrodynamic solvent at high shear rates do not get shortened in the vorticity direction. This is a consequence of the backflow produced by the interaction of the sheared solvent with the end-free polymer topology. The extended structures of the ring in the vorticity-flow plane, when they are aligned in a constant velocity plane, favor ring contour fluctuations. This variety of conformations largely suppresses the tank-treading type of rotation with extended conformations in favor of the tumbling type of rotations, where stretched and collapsed conformations alternate. The extension of trefoil knots is also enhanced, so that the knots become delocalized. We anticipate that these effects, which disappear in the absence of hydrodynamic interactions, will have a crucial impact on the rheological properties of concentrated ring solutions, and will also influence the behavior of more complicated systems such as mixtures of polymers with different topologies.

Star Block-Copolymers in Shear Flow

Star block-copolymers (SBCs) have been demonstrated to constitute self-assembling building blocks with specific softness, functionalization, shape, and flexibility. In this work, we study the behavior of an isolated SBC under a shear flow by means of particle-based multiscale simulations. We systematically analyze the conformational properties of low-functionality stars, as well as the formation of attractive patches on their corona as a function of the shear rate. We cover a wide range of system parameters, including functionality, amphiphilicity, and solvent quality. It is shown that SBCs display a richer structural and dynamical behavior than athermal star polymers in a shear flow [ Ripoll Phys. Rev. Lett. , 2006 , 96 , 188302 ], and, therefore, they are also interesting candidates to tune the viscoelastic properties of complex fluids. We identify three factors of patch reorganization under shear that lead to patch numbers and orientations depending on the shear rate, namely, free arms joining existing patches, fusion of medium-sized patches into bigger ones, and fission of large patches into two smaller ones under high shear rates. Because the conformation of single SBC is expected to be preserved in low-density bulk phases, the presented results are a first step in understanding and predicting the rheological properties of semidilute suspensions of this kind of polymers.

Steady state sedimentation of ultrasoft colloids

The structural and dynamical properties of ultra-soft colloids-star polymers-exposed to a uniform external force field are analyzed by applying the multiparticle collision dynamics technique, a hybrid coarse-grain mesoscale simulation approach, which captures thermal fluctuations and long-range hydrodynamic interactions. In the weak-field limit, the structure of the star polymer is nearly unchanged; however, in an intermediate regime, the radius of gyration decreases, in particular transverse to the sedimentation direction. In the limit of a strong field, the radius of gyration increases with field strength. Correspondingly, the sedimentation coefficient increases with increasing field strength, passes through a maximum, and decreases again at high field strengths. The maximum value depends on the functionality of the star polymer. High field strengths lead to symmetry breaking with trailing, strongly stretched polymer arms and a compact star-polymer body. In the weak-field-linear response regime, the sedimentation coefficient follows the scaling relation of a star polymer in terms of functionality and arm length.

Self-assembly of magnetically functionalized star-polymer nano-colloids

We explore the potential of star-polymers that carry super-paramagnetic nano-particles as end-groups with respect to the single-molecule self-assembly process. With the aid of molecular dynamics simulation, the configurations of these macromolecules are analyzed as a function of functionality, magnetic interaction strength, and the length of the polymeric arms. By means of an external magnetic field the nano-particles can be controlled to form static or dynamic dipolar chains, resulting in conformations of isolated stars that can be characterized by the average number of chains and length. The single-molecule conformation diagram in the plane of magnetic interaction strength vs. the star-functionality is obtained. Further, the molecules are characterized by means of various shape and size order parameters.

Modeling shear-induced particle ordering and deformation in a dense soft particle suspension

We apply the lattice Boltzmann method and the bead-spring network model of deformable particles (DPs) to study shear-induced particle ordering and deformation and the corresponding rheological behavior for dense DP suspensions confined in a narrow gap under steady external shear. The particle configuration is characterized with small-angle scattering intensity, the real-space 2D local order parameter, and the particle shape factors including deformation, stretching and tilt angles. We investigate how particle ordering and deformation vary with the particle volume fraction ϕ (=0.45-0.65) and the external shear rate characterized with the capillary number Ca (=0.003-0.191). The degree of particle deformation increases mildly with ϕ but significantly with Ca. Under moderate shear rate (Ca  =  0.105), the inter-particle structure evolves from string-like ordering to layered hexagonal close packing (HCP) as ϕ increases. A long wavelength particle slithering motion emerges for sufficiently large ϕ. For ϕ  =  0.61, the structure maintains layered HCP for Ca  =  0.031-0.143 but gradually becomes disordered for larger and smaller Ca. The correlation in particle zigzag movements depends sensitively on ϕ and particle ordering. Layer-by-layer analysis reveals how the non-slippery hard walls affect particle ordering and deformation. The shear-induced reconfiguration of DPs observed in the simulation agrees qualitatively with experimental results of sheared uniform emulsions. The apparent suspension viscosity increases with ϕ but exhibits much weaker dependence compared to hard-sphere suspensions, indicating that particle deformation and unjamming under shear can significantly reduce the viscous stress. Furthermore, the suspension shear-thins, corresponding to increased inter-DP ordering and particle deformation with Ca. This work provides useful insights into the microstructure-rheology relationship of concentrated deformable particle suspensions.

Microstructure of Flow-Driven Suspension of Hardspheres in Cylindrical Confinement: A Dynamical Density Functional Theory and Monte Carlo Study

We have studied the microstructure of a flow-driven hardsphere suspension inside a cylinder using dynamical density functional theory and Monte Carlo simulations. In order to be representative of various physical conditions that may prevail in experiments, we investigate the problem using both the grand canonical (μVT) ensemble and the canonical (NVT) ensemble. In both ensembles, the hydrodynamic effect on the suspension mediated by the presence of the confining wall is implemented in a mean-field fashion by incorporating the thermodynamic work done by the inertial lift force on the particle given the average flow field. The predicted particle distribution in the μVT ensemble displays strong structural ordering at increasing flow rates due to the correspondingly higher particle concentrations inside the cylinder. In the NVT ensemble, for dilute suspensions we observe a peak in the distribution of density at a location similar to that of the Segré-Silberberg annulus, while for dense suspensions the competing effects of the inertial lift and the hardsphere interaction lead to the formation of several annuli.

Solvent Induced Inversion of Core-Shell Microgels

The morphology of core-shell microgels under different swelling conditions and as a function of the core-shell thickness ratio is systematically characterized by mesoscale hydrodynamic simulations. With increasing hydrophobic interaction of the shell polymers, we observe drastic morphological changes from a core-shell structure to an inverted microgel, where the core is turned to the outside, or a microgel with a patchy surface of core polymers directly exposed to the environment. We establish a phase diagram of the various morphologies. Moreover, we characterize the polymer and microgel conformations. For sufficiently thick shells, the changes of the shell size upon increasing hydrophobic interactions are well described by the Flory-Rehner theory. Additionally, this theory provides a critical line in the phase diagram separating core-shell structures from the distinct two other phases. The appearing new phases provide a novel route to nano- and microscale functionalized materials.

Tuning Adsorption Duration To Control the Diffusion of a Nanoparticle in Adsorbing Polymers

Controlling the nanoparticle (NP) diffusion in polymers is a prerequisite to obtain polymer nanocomposites (PNCs) with desired dynamical and rheological properties and to achieve targeted delivery of nanomedicine in biological systems. Here we determine the suppression mechanism of direct NP-polymer attraction to hamper the NP mobility in adsorbing polymers and then quantify the dependence of the effective viscosity η_eff felt by the NP on the adsorption duration τ_ads of polymers on the NP using scaling theory analysis and molecular dynamics simulations. We propose and confirm that participation of adsorbed chains in the NP motion break up at time intervals beyond τ_ads due to the rearrangement of polymer segments at the NP surface, which accounts for the onset of Fickian NP diffusion on a time scale of t ≈ τ_ads. We develop a power law, η_eff ∼ (τ_ads)^ν, where ν is the scaling exponent of the dependence of polymer coil size on the chain length, which leads to a theoretical basis for the design of PNCs and nanomedicine with desired applications through tuning the polymer adsorption duration.

Internal dynamics of microgels: A mesoscale hydrodynamic simulation study

We analyze the dynamics of polymers in a microgel system under different swelling conditions. A microgel particle consists of coarse-grained linear polymers which are tetra-functionally crosslinked and undergoes conformational changes in response to the external stimuli. Here, a broad range of microgel sizes, extending from tightly collapsed to strongly swollen particles, is considered. In order to account for hydrodynamic interactions, the microgel is embedded in a multiparticle collision dynamics fluid while hydrophobic attraction is modelled by an attractive Lennard-Jones potential and swelling of ionic microgels is described through the Debye-Hückel potential. The polymer dynamics is analyzed in terms of the monomer mean square displacement and the intermediate scattering function S(q, t). The scattering function decays in a stretched-exponential manner, with a decay rate exhibiting a crossover from a collective diffusive dynamics at low magnitudes of the wavevector q to a hydrodynamic-dominated dynamics at larger q. There is little difference between the intermediate scattering functions of microgels under good solvent conditions and strongly swollen gels, but strongly collapsed gels exhibit a faster decay at short times and hydrodynamic interactions become screened. In addition, we present results for the dynamics of the crosslinks, which exhibit an unexpected, semiflexible polymer-like dynamics.

Swelling, structure, and phase stability of compressible microgels

Microgels are soft colloidal particles that, when dispersed in a solvent, swell and deswell in response to changes in environmental conditions, such as temperature, concentration, and pH. Using Monte Carlo simulation, we model bulk suspensions of microgels that interact via Hertzian elastic interparticle forces and can expand or contract via trial moves that allow particles to change size in accordance with the Flory-Rehner free energy of cross-linked polymer gels. We monitor the influence of particle compressibility, size fluctuations, and concentration on bulk structural and thermal properties by computing particle swelling ratios, radial distribution functions, static structure factors, osmotic pressures, and freezing densities. For microgels in the nanoscale size range, particle compressibility and associated size fluctuations suppress crystallization, shifting the freezing transition to a higher density than for the hard-sphere fluid. As densities increase beyond close packing, microgels progressively deswell, while their intrinsic size distribution grows increasingly polydisperse.

Dynamic and rheological properties of soft biological cell suspensions

Quantifying dynamic and rheological properties of suspensions of soft biological particles such as vesicles, capsules, and red blood cells (RBCs) is fundamentally important in computational biology and biomedical engineering. In this review, recent studies on dynamic and rheological behavior of soft biological cell suspensions by computer simulations are presented, considering both unbounded and confined shear flow. Furthermore, the hemodynamic and hemorheological characteristics of RBCs in diseases such as malaria and sickle cell anemia are highlighted.

Hydrodynamic segregation in a bidisperse colloidal suspension in microchannel flow: A theoretical study

Colloids in suspension exhibit shear-induced migration towards regions of low viscous shear. In dense bidisperse colloidal suspensions under pressure driven flow large particles can segregate in the center of a microchannel and the suspension partially demixes. To develop a theoretical understanding of these effects, we formulate a phenomenological model for the particle currents based on the work of Phillips et al. [Phys. Fluids 4, 30 (1992)]. We also simulate hard spheres under pressure-driven flow in two and three dimensions using the mesoscale simulation technique of multi-particle collision dynamics. Using a single fit parameter for the intrinsic diffusivity, our theory accurately reproduces the simulated density profiles across the channel. We present a detailed parameter study on how a monodisperse suspension enriches the channel center and quantitatively confirm the experimental observation that a binary colloidal mixture partially segregates into its two species. In particular, we always find a strong accumulation of large particles in the center. Qualitative differences between two and three dimensions reveal that collective diffusion is more relevant in two dimensions.

Equilibrium and nonequilibrium dynamics of soft sphere fluids

We use computer simulations to test the freezing-point scaling relationship between equilibrium transport coefficients (self-diffusivity, viscosity) and thermodynamic parameters for soft sphere fluids. The fluid particles interact via the inverse-power potential (IPP), and the particle softness is changed by modifying the exponent of the distance-dependent potential term. In the case of IPP fluids, density and temperature are not independent variables and can be combined to obtain a coupling parameter to define the thermodynamic state of the system. We find that the rescaled coupling parameter, based on its value at the freezing point, can approximately collapse the diffusivity and viscosity data for IPP fluids over a wide range of particle softness. Even though the collapse is far from perfect, the freezing-point scaling relationship provides a convenient and effective way to compare the structure and dynamics of fluid systems with different particle softness. We further show that an alternate scaling relationship based on two-body excess entropy can provide an almost perfect collapse of the diffusivity and viscosity data below the freezing transition. Next, we perform nonequilibrium molecular dynamics simulations to calculate the shear-dependent viscosity and to identify the distinct role of particle softness in underlying structural changes associated with rheological properties. Qualitatively, we find a similar shear-thinning behavior for IPP fluids with different particle softness, though softer particles exhibit stronger shear-thinning tendency. By investigating the distance and angle-dependent pair correlation functions in these systems, we find different structural features in the case of IPP fluids with hard-sphere like and softer particle interactions. Interestingly, shear-thinning in hard-sphere like fluids is accompanied by enhanced translational order, whereas softer fluids exhibit loss of order with shear. Our results provide a systematic evaluation of the role of particle softness in equilibrium and nonequilibrium transport properties and their underlying connection with thermodynamic and structural properties.