Dynamics of polymers in a particle-based mesoscopic solvent

Multi-particle collision dynamics for a coarse-grained model of soft colloids applied to ring polymers

The simulation of polymer solutions often requires the development of methods that accurately include hydrodynamic interactions. Resolution on the atomistic scale is too computationally expensive to cover mesoscopic time and length scales on which the interesting polymer phenomena are observed. Therefore, coarse-graining methods have to be applied. In this work, the solvent is simulated using the well-established multi-particle collision dynamics scheme, and for the polymer, different coarse-graining methods are employed and compared against the monomer resolved Kremer-Grest model by their resulting diffusion coefficients. This research builds on previous work [Ruiz-Franco et al., J. Chem. Phys. 151, 074902 (2019)], in which star polymers and linear chains in a solvent were simulated and two different coarse-graining methods were developed, in order to increase computational efficiency. The present work extends this approach to ring polymers and seeks to refine one of the authors' proposed model: the penetrable soft colloid model. It was found that both proposed models are not well suited to ring polymers; however, the introduction of a factor to the PSC model delivers satisfying results for the diffusion behavior by regulating the interaction intensity with the solvent.

Particle-based mesoscopic model for phase separation in a binary fluid mixture

A mesoscopic simulation model to study the phase separation in a binary fluid mixture in three dimensions (3D) is presented here by augmenting the existing particle-based multiparticle collision dynamics (MPCD) algorithm. The approach describes the nonideal equation of the fluid state by incorporating the excluded-volume interaction between the two components within the framework of stochastic collision, which depends on the local fluid composition and velocity. Calculating the nonideal contribution to the pressure both from simulation and analytics shows the model to be thermodynamically consistent. A phase diagram to explore the range of parameters that give rise to phase separation in the model is investigated. The interfacial width and phase growth obtained from the model agree with the literature for a wide range of temperatures and parameters.

Single-polymer dynamics of starch-like branched ring polymers in steady shear flow

The stretching dynamics and dynamical behaviors of individual branched ring polymer (BRP), a coarse-grained model for some types of the starch, in steady shear flow are studied by using a hybrid mesoscale simulation approach that combines multiparticle collision dynamics with standard molecular dynamics. By analyzing the stretched configuration of BRPs, we find the polymer size increases nonmonotonically with increasing branch length. Meanwhile, the decrease of the alignment angle of the stretched configuration of BRPs follows a universal power law during the first downward phase as the shear rate increases. Constructing the three-dimensional surface of the polymer's ring backbone and tracing the temporal fluctuations of the surface's normal vector along the simulation trajectory, the tumbling and tank-treading motion are clearly reflected by periodic and non-periodic changes of the normal vector. Interestingly, these temporal changes are much more regular than that of the gyration tensor. Thus, a novel cross-correlation function, which is the correlation between fluctuations of the normal vector along the flow direction and the velocity-gradient direction, is proposed to analyze the tumbling motion that usually coexists with the tank-treading motion. This function can naturally address the fails of traditional method that analyzing the tumbling motion by determining the correlation of temporal fluctuations of the gyration tensor G_αα. By analyzing the dynamical behaviors of BRPs, diverse dependences of the tumbling frequency ω_TB and tank-treading frequency ω_TT on the shear rate γ̇ are observed at a wide range of shear rates and polymer sizes. Furthermore, our simulations also reveal that the tank-treading motion is more stable than the tumbling motion for small-branch-size BRPs but the tumbling motion is more stable than the tank-treading motion for large-branch-size BRPs.

Effect of Charge Distribution on the Dynamics of Polyampholytic Disordered Proteins

The stability and physiological function of many biomolecular coacervates depend on the structure and dynamics of intrinsically disordered proteins (IDPs) that typically contain a significant fraction of charged residues. Although the effect of relative arrangement of charged residues on IDP conformation is a well-studied problem, the associated changes in dynamics are far less understood. In this work, we systematically interrogate the effects of charge distribution on the chain-level and segmental dynamics of polyampholytic IDPs in dilute solutions. We study a coarse-grained model polyampholyte consisting of an equal fraction of two oppositely charged residues (glutamic acid and lysine) that undergoes a transition from an ideal chain-like conformation for uniformly charge-patterned sequences to a semi-compact conformation for highly charge-segregated sequences. Changes in the chain-level dynamics with increasing charge segregation correlate with changes in conformation. The chain-level and segmental dynamics conform to simple homopolymer models for uniformly charge-patterned sequences but deviate with increasing charge segregation, both in the presence and absence of hydrodynamic interactions. We discuss the significance of these findings, obtained for a model polyampholyte, in the context of a charge-rich intrinsically disordered region of the naturally occurring protein LAF-1. Our findings have important implications for understanding the effects of charge patterning on the dynamics of polyampholytic IDPs in dilute conditions using polymer scaling theories.

Clustering of self-thermophilic asymmetric dimers: the relevance of hydrodynamics

Self-thermophilic dimers are characterized by a net phoretic attraction which, in combination with hydrodynamic interactions, results in the formation of crystalline-like aggregates. To distinguish the effect of the different contributions is frequently an important challenge. We present a simulation investigation done in parallel in the presence and the absence of hydrodynamic interactions for the case of asymmetric self-thermophoretic dimers. In the absence of hydrodynamics, the clusters have the standard heads-in configurations. In contrast, in the presence of hydrodynamics, clusters with heads-in conformation are being formed, in which dimers with their propulsion velocity pointing out of the cluster are assembled and stabilized by strong hydrodynamic osmotic flows. Significant variation in the material properties is to be expected from such differences in the collective behavior, whose understanding and control is of great relevance for the development of new synthetic active materials.

Dynamical and Structural Properties of Comb Long-Chain Branched Polymer in Shear Flow

Using hybrid multi-particle collision dynamics (MPCD) and a molecular dynamics (MD) method, we investigate the effect of arms and shear flow on dynamical and structural properties of the comb long-chain branched (LCB) polymer with dense arms. Firstly, we analyze dynamical properties of the LCB polymer by tracking the temporal changes on the end-to-end distance of both backbones and arms as well as the orientations of the backbone in the flow-gradient plane. Simultaneously, the rotation and tumbling behaviors with stable frequencies are observed. In other words, the LCB polymer undergoes a process of periodic stretched–folded–stretched state transition and rotation, whose period is obtained by fitting temporal changes on the orientation to a periodic function. In addition, the impact induced by random and fast motions of arms and the backbone will descend as the shear rate increases. By analyzing the period of rotation behavior of LCB polymers, we find that arms have a function in keeping the LCB polymer’s motion stable. Meanwhile, we find that the rotation period of the LCB polymer is mainly determined by the conformational distribution and the non-shrinkable state of the structure along the velocity-gradient direction. Secondly, structural properties are numerically characterized by the average gyration tensor of the LCB polymer. The changes in gyration are in accordance with the LCB polymer rolling when varying the shear rate. By analyzing the alignment of the LCB polymer and comparing with its linear and star counterparts, we find that the LCB polymer with very long arms, like the corresponding linear chain, has a high speed to reach its configuration expansion limit in the flow direction. However, the comb polymer with shorter arms has stronger resistance on configuration expansion against the imposed flow field. Moreover, with increasing arm length, the comb polymer in shear flow follows change from linear-polymer-like to capsule-like behavior.

Electroosmotic Flow Induced Lift Forces on Polymer Chains in Nanochannels

A major objective of research in nanofluidics is to achieve better selectivity in manipulating the fluxes of nano-objects and in particular of biopolymers. Numerical simulations allow one to better understand the physical mechanisms at play in such situations. We performed hybrid mesoscale simulations to investigate the properties of polymers under flows in slit pores at the nanoscale. We use multiparticle collision dynamics, an algorithm that includes hydrodynamics and thermal fluctuations, to investigate the properties of fully flexible and stiff polymers under several types of flow, showing that Poiseuille flows and electroosmotic flows can lead to quantitatively and qualitatively different behaviors of the chain. In particular, a counterintuitive phenomenon occurs in the presence of an electroosmotic flow: When the monomers are attracted by the solid surfaces through van der Waals forces, shear-induced forces lead to a stronger repulsion of the polymers from these surfaces. Such focusing of the chain in the middle of the channel increases its flowing velocity, a phenomenon that may be exploited to separate different types of polymers.

Simulating wet active polymers by multiparticle collision dynamics

The conformational and dynamical properties of active Brownian polymers embedded in a fluid depend on the nature of the driving mechanism, e.g., self-propulsion or external actuation of the monomers. Implementations of self-propelled and actuated active Brownian polymers in a multiparticle collision (MPC) dynamics fluid are presented, which capture the distinct differences between the two driving mechanisms. The active force-free nature of self-propelled monomers requires adaptations of the MPC simulation scheme, with its streaming and collision steps, where the monomer self-propulsion velocity has to be omitted in the collision step. Comparison of MPC simulation results for active polymers in dilute solution with results of Brownian dynamics simulations accounting for hydrodynamics via the Rotne-Prager-Yamakawa tensor confirm the suitability of the implementation. The polymer conformational and dynamical properties are analyzed by the static and dynamic structure factor, and the scaling behavior of the latter with respect to the wave number and time dependence are discussed. The dynamic structure factor displays various activity-induced temporal regimes, depending on the considered wave number, which reflect the persistent diffusive motion of the whole polymer at small wave numbers, and the activity-enhanced internal dynamics at large wave numbers. The obtained simulation results are compared with theoretical predictions.

Hydrodynamics of immiscible binary fluids with viscosity contrast: a multiparticle collision dynamics approach

We present a multiparticle collision dynamics (MPC) implementation of layered immiscible fluids A and B of different shear viscosities separated by planar interfaces. The simulated flow profile for imposed steady shear motion and the time-dependent shear stress functions are in excellent agreement with our continuum hydrodynamics results for the composite fluid. The wave-vector dependent transverse velocity auto-correlation functions (TVAF) in the bulk-fluid regions of the layers decay exponentially, and agree with those of single-phase isotropic MPC fluids. In addition, we determine the hydrodynamic mobilities of an embedded colloidal sphere moving steadily parallel or transverse to a fluid-fluid interface, as functions of the distance from the interface. The obtained mobilities are in good agreement with hydrodynamic force multipoles calculations, for a no-slip sphere moving under creeping flow conditions near a clean, ideally flat interface. The proposed MPC fluid-layer model can be straightforwardly implemented, and it is computationally very efficient. Yet, owing to the spatial discretization inherent to the MPC method, the model can not reproduce all hydrodynamic features of an ideally flat interface between immiscible fluids.

Crowded solutions of single-chain nanoparticles under shear flow

Single-chain nanoparticles (SCNPs) are ultrasoft objects obtained through purely intramolecular cross-linking of single polymer chains. By means of computer simulations with implemented hydrodynamic interactions, we investigate for the first time the effect of the shear flow on the structural and dynamic properties of SCNPs in semidilute and concentrated solutions. We characterize the dependence of several conformational and dynamic observables on the shear rate and the concentration, obtaining a set of power-law scaling laws. The concentration has a very different effect on the shear rate dependence of the former observables in SCNPs than in simple linear chains. Whereas for the latter the scaling behaviour is marginally dependent on the concentration, two clearly different scaling regimes are found for the SCNPs below and above the overlap concentration. At fixed shear rate SCNPs and linear chains also respond very differently to crowding. Whereas, at moderate and high Weissenberg numbers the linear chains swell, the SCNPs exhibit a complex non-monotonic behaviour. We suggest that these findings are inherently related to the topological interactions preventing concatenation of the SCNPs, which lead to less interpenetration than for linear chains, and to the limitation to stretching imposed by the permanent cross-links in the SCNPs, which itself limits the ways to spatially arrange in the shear flow.

Stratification of polymer mixtures in drying droplets: Hydrodynamics and diffusion

We study the evaporation-induced stratification of a mixture of short and long polymer chains in a drying droplet using molecular simulations. We systematically investigate the effects of hydrodynamic interactions (HI) on this process by comparing hybrid simulations accounting for HI between polymers through the multiparticle collision dynamics technique with free-draining Langevin dynamics simulations neglecting the same. We find that the dried supraparticle morphologies are homogeneous when HI are included but are stratified in core-shell structures (with the short polymers forming the shell) when HI are neglected. The simulation methodology unambiguously attributes this difference to the treatment of the solvent in the two models. We rationalize the presence (or absence) of stratification by measuring phenomenological multicomponent diffusion coefficients for the polymer mixtures. The diffusion coefficients show the importance of not only solvent backflow but also HI between polymers in controlling the dried supraparticle morphology.

Wall entrapment of peritrichous bacteria: a mesoscale hydrodynamics simulation study

Microswimmers such as E. coli bacteria accumulate and exhibit an intriguing dynamics near walls, governed by hydrodynamic and steric interactions. Insight into the underlying mechanisms and predominant interactions demand a detailed characterization of the entrapment process. We employ a mesoscale hydrodynamics simulation approach to study entrapment of an E. coli-type cell at a no-slip wall. The cell is modeled by a spherocylindrical body with several explicit helical flagella. Three stages of the entrapment process can be distinguished: the approaching regime, where a cell swims toward a wall on a nearly straight trajectory; a scattering regime, where the cell touches the wall and reorients; and a surface-swimming regime. Our simulations show that steric interactions may dominate the entrapment process, yet, hydrodynamic interactions slow down the adsorption dynamics close to the boundary and imply a circular motion on the wall. The locomotion of the cell is characterized by a strong wobbling dynamics, with cells preferentially pointing toward the wall during surface swimming.

Coarsening Kinetics of Complex Macromolecular Architectures in Bad Solvent

This study reports a general scenario for the out-of-equilibrium features of collapsing polymeric architectures. We use molecular dynamics simulations to characterize the coarsening kinetics, in bad solvent, for several macromolecular systems with an increasing degree of structural complexity. In particular, we focus on: flexible and semiflexible polymer chains, star polymers with 3 and 12 arms, and microgels with both ordered and disordered networks. Starting from a powerful analogy with critical phenomena, we construct a density field representation that removes fast fluctuations and provides a consistent characterization of the domain growth. Our results indicate that the coarsening kinetics presents a scaling behaviour that is independent of the solvent quality parameter, in analogy to the time-temperature superposition principle. Interestingly, the domain growth in time follows a power-law behaviour that is approximately independent of the architecture for all the flexible systems; while it is steeper for the semiflexible chains. Nevertheless, the fractal nature of the dense regions emerging during the collapse exhibits the same scaling behaviour for all the macromolecules. This suggests that the faster growing length scale in the semiflexible chains originates just from a faster mass diffusion along the chain contour, induced by the local stiffness. The decay of the dynamic correlations displays scaling behavior with the growing length scale of the system, which is a characteristic signature in coarsening phenomena.

Enhanced Rotational Motion of Spherical Squirmer in Polymer Solutions

The rotational diffusive motion of a self-propelled, attractive spherical colloid immersed in a solution of self-avoiding polymers is studied by mesoscale hydrodynamic simulations. A drastic enhancement of the rotational diffusion by more than an order of magnitude in the presence of activity is obtained. The amplification is a consequence of two effects, a decrease of the amount of adsorbed polymers by active motion and an asymmetric encounter with polymers on the squirmer surface, which yields an additional torque and random noise for the rotational motion. Our simulations suggest a way to control the rotational dynamics of squirmer-type microswimmers by the degree of polymer adsorption and system heterogeneity.

Dynamics control of an in-plane-switching liquid crystal cell using heterogeneous substrates

The control of surface anchoring strength can be achieved by using heterogeneous substrates. In contrast to conventional substrates that control the anchoring strength by using temperature or chemical processes, heterogeneous substrates provide surface anchoring to liquid crystal molecules by using a mixed composition of (1) a zero anchoring surface and (2) planar-anchoring patches. To study the dynamics of in-plane-switching liquid crystal displays (IPS-LCDs) under external fields, a new particle-based numerical algorithm is developed to simulate both nematic liquid crystals and heterogeneous surfaces. This new method allows us to create different heterogeneous surfaces easily by adopting predefined distributions of numerical particles. The generated effective anchoring strength from the heterogeneous surface is thus calculated, and the dynamical response is found to be similar to that of conventional homogeneously processed substrates. The results suggest that the use of a heterogeneous LCD cell is a suitable alternative for creating desirable LCD substrates, for which chemical/temperature dependence can be transferred to a more controllable configurational dependence. Interestingly, we found master curves in the peak transmittance/recovery time phase space, and they appeared to be dependent solely on the cell thickness. This discovery clarifies the fundamental optical dynamics of IPS-LCD cells.

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.

Nanoparticle transport phenomena in confined flows

Nanoparticles submerged in confined flow fields occur in several technological applications involving heat and mass transfer in nanoscale systems. Describing the transport with nanoparticles in confined flows poses additional challenges due to the coupling between the thermal effects and fluid forces. Here, we focus on the relevant literature related to Brownian motion, hydrodynamic interactions and transport associated with nanoparticles in confined flows. We review the literature on the several techniques that are based on the principles of non-equilibrium statistical mechanics and computational fluid dynamics in order to simultaneously preserve the fluctuation-dissipation relationship and the prevailing hydrodynamic correlations. Through a review of select examples, we discuss the treatments of the temporal dynamics from the colloidal scales to the molecular scales pertaining to nanoscale fluid dynamics and heat transfer. As evident from this review, there, indeed has been little progress made in regard to the accurate modeling of heat transport in nanofluids flowing in confined geometries such as tubes. Therefore the associated mechanisms with such processes remain unexplained. This review has revealed that the information available in open literature on the transport properties of nanofluids is often contradictory and confusing. It has been very difficult to draw definitive conclusions. The quality of work reported on this topic is non-uniform. A significant portion of this review pertains to the treatment of the fluid dynamic aspects of the nanoparticle transport problem. By simultaneously treating the energy transport in ways discussed in this review as related to momentum transport, the ultimate goal of understanding nanoscale heat transport in confined flows may be achieved.

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.

A hybrid method for micro-mesoscopic stochastic simulation of reaction-diffusion systems

The present paper introduces a new micro-meso hybrid algorithm based on the Ghost Cell Method concept in which the microscopic subdomain is governed by the Reactive Multi-Particle Collision (RMPC) dynamics. The mesoscopic subdomain is modeled using the Reaction-Diffusion Master Equation (RDME). The RDME is solved by means of the Inhomogeneous Stochastic Simulation Algorithm. No hybrid algorithm has hitherto used the RMPC dynamics for modeling reactions and the trajectories of each individual particle. The RMPC is faster than other molecular based methods and has the advantage of conserving mass, energy and momentum in the collision and free streaming steps. The new algorithm is tested on three reaction-diffusion systems. In all the systems studied, very good agreement with the deterministic solutions of the corresponding differential equations is obtained. In addition, it has been shown that proper discretization of the computational domain results in significant speed-ups in comparison with the full RMPC algorithm.

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.

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.

Guidance and Self-Sorting of Active Swimmers: 3D Periodic Arrays Increase Persistence Length of Human Sperm Selecting for the Fittest

Male infertility is a reproductive disease, and existing clinical solutions for this condition often involve long and cumbersome sperm sorting methods, including preprocessing and centrifugation-based steps. These methods also fall short when sorting for sperm free of reactive oxygen species, DNA damage, and epigenetic aberrations. Although several microfluidic platforms exist, they suffer from structural complexities, i.e., pumps or chemoattractants, setting insurmountable barriers to clinical adoption. Inspired by the natural filter-like capabilities of the female reproductive tract for sperm selection, a model-driven design, featuring pillar arrays that efficiently and noninvasively isolate highly motile and morphologically normal sperm, with lower epigenetic global methylation, from raw semen, is presented. The Simple Periodic ARray for Trapping And isolatioN (SPARTAN) created here modulates the directional persistence of sperm, increasing the spatial separation between progressive and nonprogressive motile sperm populations within an unprecedentedly short 10 min assay time. With over 99% motility of sorted sperm, a 5-fold improvement in morphology, 3-fold increase in nuclear maturity, and 2-4-fold enhancement in DNA integrity, SPARTAN offers to standardize sperm selection while eliminating operator-to-operator variations, centrifugation, and flow. SPARTAN can also be applied in other areas, including conservation ecology, breeding of farm animals, and design of flagellar microrobots for diagnostics.

Clustering and dynamics of particles in dispersions with competing interactions: theory and simulation

Dispersions of particles with short-range attractive and long-range repulsive interactions exhibit rich equilibrium microstructures and a complex phase behavior. We present theoretical and simulation results for structural and, in particular, short-time diffusion properties of a colloidal model system with such interactions, both in the dispersed-fluid and equilibrium-cluster phase regions. The particle interactions are described by a generalized Lennard-Jones-Yukawa pair potential. For the theoretical-analytical description, we apply the hybrid Beenakker-Mazur pairwise additivity (BM-PA) scheme. The static structure factor input to this scheme is calculated self-consistently using the Zerah-Hansen integral equation theory approach. In the simulations, a hybrid simulation method is adopted, combing molecular dynamics simulations of colloids with the multiparticle collision dynamics approach for the fluid, which fully captures hydrodynamic interactions. The comparison of our theoretical and simulation results confirms the high accuracy of the BM-PA scheme for dispersed-fluid phase systems. For particle attraction strengths exceeding a critical value, our simulations yield an equilibrium cluster phase. Calculations of the mean lifetime of the appearing clusters and the comparison with the analytical prediction of the dissociation time of an isolated particle pair reveal quantitative differences pointing to the importance of many-particle hydrodynamic interactions for the cluster dynamics. The cluster lifetime in the equilibrium-cluster phase increases far stronger with increasing attraction strength than that in the dispersed-fluid phase. Moreover, significant changes in the cluster shapes are observed in the course of time. Hence, an equilibrium-cluster dispersion cannot be treated dynamically as a system of permanent rigid bodies.

The effect of hydrodynamic interactions on nanoparticle diffusion in polymer solutions: a multiparticle collision dynamics study

The diffusion of nanoparticles (NPs) in polymer solutions is studied by a combination of a mesoscale simulation method, multiparticle collision dynamics (MPCD), and molecular dynamics (MD) simulations. We investigate the long-time diffusion coefficient D as well as the subdiffusive behavior in the intermediate time region. The dependencies of both D and subdiffusion factor α on NP size and polymer concentration, respectively, are explicitly calculated. Particular attention is paid to the role of hydrodynamic interaction (HI) in the NP diffusion dynamics. Our simulation results show that the long-time diffusion coefficients satisfy perfectly the scaling relation found by experimental observations. Meanwhile, the subdiffusive factor decreases with the increase in polymer concentration but is of little relevance to the NP size. By parallel simulations with and without HI, we reveal that HI will generally enhance D, while the enhancement effect is non-monotonous with increasing polymer concentration, and it becomes most pronounced at semidilute concentrations. With the aid of a scaling law based on the diffusive activation energy model, we understand that HI affects diffusion through decreasing the diffusive activation energy on the one hand while increasing the effective diffusion size on the other. In addition, HI will certainly influence the subdiffusive behavior of the NP, leading to a larger subdiffusion exponent.

Anisotropic thermophoresis

Colloidal migration in a temperature gradient is referred to as thermophoresis. In contrast to particles with a spherical shape, we show that elongated colloids may have a thermophoretic response that varies with the colloid orientation. Remarkably, this can translate into a non-vanishing thermophoretic force in the direction perpendicular to the temperature gradient. Opposite to the friction force, the thermophoretic force of a rod oriented with the temperature gradient can be larger or smaller than when oriented perpendicular to it. The precise anisotropic thermophoretic behavior clearly depends on the colloidal rod aspect ratio, and also on its surface details, which provides an interesting tunability to the devices constructed based on this principle. By means of mesoscale hydrodynamic simulations, we characterize this effect for different types of rod-like colloids.

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.

Shear viscosity in hard-sphere and adhesive colloidal suspensions with reverse non-equilibrium molecular dynamics

We employ the reverse non-equilibrium molecular dynamics method (RNEMD) of Müller-Plathe [Phys. Rev. E, 1999, 59, 4894] to calculate the shear viscosity of colloidal suspensions within the stochastic rotation dynamics-molecular dynamics (SRD-MD) simulation method. We examine the influence of different coupling schemes in SRD-MD on the colloidal volume fraction ϕ_c dependent viscosity from the dilute limit up to ϕ_c = 0.3. Our results demonstrate that the RNEMD method is a robust and reliable method for calculating rheological properties of colloidal suspensions. To obtain quantitatively accurate results beyond the dilute regime, the hydrodynamic interactions between the effective fluid particles in the SRD and the MD colloidal particles must be carefully considered in the coupling scheme. We benchmark the method by comparing with the hard sphere suspension case, and then calculate relative viscosities for colloids with mutually attractive interactions. We show that the viscosity displays a sharp increase at the onset of aggregation of the colloidal particles with increasing volume fraction and attraction.

Diffusion in systems crowded by active force-dipole molecules

Experimental studies of systems containing active proteins that undergo conformational changes driven by catalytic chemical reactions have shown that the diffusion coefficients of passive tracer particles and active molecules are larger than the corresponding values when chemical activity is absent. Various mechanisms have been proposed for such behavior, including, among others, force dipole interactions of molecular motors moving on filaments and collective hydrodynamic effects arising from active proteins. Simulations of a multi-component system containing active dumbbell molecules that cycle between open and closed states, a passive tracer particle and solvent molecules are carried out. Consistent with experiments, it is shown that the diffusion coefficients of both passive particles and the dumbbells themselves are enhanced when the dumbbells are active. The dependence of the diffusion enhancement on the volume fraction of dumbbells is determined, and the effects of crowding by active dumbbell molecules are shown to differ from those due to inactive molecules.

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.

Polymer Conformations in Ionic Microgels in the Presence of Salt: Theoretical and Mesoscale Simulation Results

We investigate the conformational properties of polymers in ionic microgels in the presence of salt ions by molecular dynamics simulations and analytical theory. A microgel particle consists of coarse-grained linear polymers, which are tetra-functionally crosslinked. Counterions and salt ions are taken into account explicitly, and charge-charge interactions are described by the Coulomb potential. By varying the charge interaction strength and salt concentration, we characterize the swelling of the polyelectrolytes and the charge distribution. In particular, we determine the amount of trapped mobile charges inside the microgel and the Debye screening length. Moreover, we analyze the polymer extension theoretically in terms of the tension blob model taking into account counterions and salt ions implicitly by the Debye–Hückel model. Our studies reveal a strong dependence of the amount of ions absorbed in the interior of the microgel on the electrostatic interaction strength, which is related to the degree of the gel swelling. This implies a dependence of the inverse Debye screening length κ on the ion concentration; we find a power-law increase of κ with the Coulomb interaction strength with the exponent 3/5 for a salt-free microgel and an exponent 1/2 for moderate salt concentrations. Additionally, the radial dependence of polymer conformations and ion distributions is addressed.

Tunable slow dynamics in a new class of soft colloids

By means of extensive simulations, we investigate concentrated solutions of globular single-chain nanoparticles (SCNPs), an emergent class of synthetic soft nano-objects. By increasing the concentration, the SCNPs show flattening, and even reentrant behaviour, in the density dependence of their structural and dynamic correlations, as well as a soft caging regime and weak dynamic heterogeneity. The latter is confirmed by validation of the Stokes-Einstein relation up to concentrations far beyond the overlap density. Therefore SCNPs arise as a new class of soft colloids, exhibiting slow dynamics and actualizing in a real system the structural and dynamic anomalies proposed by models of ultrasoft particles. Quantitative differences in the dynamic behaviour depend on the SCNP deformability, which can be tuned through the degree of internal cross-linking.

Bacterial swarmer cells in confinement: a mesoscale hydrodynamic simulation study

A wide spectrum of Peritrichous bacteria undergo considerable physiological changes when they are inoculated onto nutrition-rich surfaces and exhibit a rapid and collective migration denoted as swarming. Thereby, the length of such swarmer cells and their number of flagella increases substantially. In this article, we investigated the properties of individual E. coli-type swarmer cells confined between two parallel walls via mesoscale hydrodynamic simulations, combining molecular dynamics simulations of the swarmer cell with the multiparticle particle collision dynamics approach for the embedding fluid. E. coli-type swarmer cells are three-times longer than their planktonic counter parts, but their flagella density is comparable. By varying the wall separation, we analyze the confinement effect on the flagella arrangement, on the distribution of cells in the gap between the walls, and on the cell dynamics. We find only a weak dependence of confinement on the bundle structure and dynamics. The distribution of cells in the gap changes from a geometry-dominated behavior for very narrow to fluid-dominated behavior for wider gaps, where cells are preferentially located in the gap center for narrower gaps and stay preferentially next to one of the walls for wider gaps. Dynamically, the cells exhibit a wide spectrum of migration behaviors, depending on their flagella bundle arrangement, and ranges from straight swimming to wall rolling.

Viscoelastic properties of marginal networks in a solvent

Polymer networks at the margins of mechanical stability are known to be highly sensitive to applied forces and fields and to exhibit an anomalously large resistance to deformation. In this paper, we study the effects of hydrodynamic interactions on the behavior of marginal networks using a hybrid molecular dynamics and multiparticle collision dynamics simulation technique. We examine how the filament and solvent properties affect the response of marginal networks to shear. We find that the stiffening of the network shows a stronger dependence on the shear frequency when hydrodynamic interactions are present than when they are not. The network shear modulus scales as G'∼ω(α(c)), with a critical stiffening exponent α(c) that can be controlled by varying the relative concentrations of the network and the solvent. Our results show that this arises due to the solvent aiding the relaxation of the network and suppressing the network nonaffinity, with the system deforming more affinely when hydrodynamic interactions are maximized.

Self-organization in suspensions of end-functionalized semiflexible polymers under shear flow

The nonequilibrium dynamical behavior and structure formation of end-functionalized semiflexible polymer suspensions under flow are investigated by mesoscale hydrodynamic simulations. The hybrid simulation approach combines the multiparticle collision dynamics method for the fluid, which accounts for hydrodynamic interactions, with molecular dynamics simulations for the semiflexible polymers. In equilibrium, various kinds of scaffold-like network structures are observed, depending on polymer flexibility and end-attraction strength. We investigate the flow behavior of the polymer networks under shear and analyze their nonequilibrium structural and rheological properties. The scaffold structure breaks up and densified aggregates are formed at low shear rates, while the structural integrity is completely lost at high shear rates. We provide a detailed analysis of the shear- rate-dependent flow-induced structures. The studies provide a deeper understanding of the formation and deformation of network structures in complex materials.

Hydrodynamic correlations and diffusion coefficient of star polymers in solution

The center-of-mass dynamics of star polymers in dilute solution is analyzed by hybrid mesoscale simulations. The fluid is modeled by the multiparticle collision dynamics approach, a particle-based hydrodynamic simulation technique, which is combined with molecular dynamics simulations for the polymers. Star polymers of various functionalities are considered. We determine the center-of-mass velocity correlation functions, the corresponding mean square displacements, and diffusion coefficients. The velocity correlation functions exhibit a functionality-dependent and structure-specific intermediate time regime, with a slow decay. It is followed by the long-time tail t(-3/2), which is solely determined by the fluid. Infinite-system-size diffusion coefficients are determined from the velocity correlation function by a combination of simulation and analytical results, as well as from the center-of-mass mean square displacement for various systems sizes and extrapolation. In terms of the hydrodynamic radius, the star polymer hydrodynamic diffusion coefficient exhibits the same universal system-size dependence as a spherical colloid. The functionality dependence of the ratio of hydrodynamic radii and the radii of gyration agrees well with experimental predictions.