Hydrodynamic correlations and diffusion coefficient of star polymers in solution

Self-Diffusion of Star and Linear Polyelectrolytes in Salt-Free and Salt Solutions

This work explored solution properties of linear and star poly(methacrylic acids) with four, six, and eight arms (LPMAA, 4PMAA, PMAA, and 8PMAA, respectively) of matched molecular weights in a wide range of pH, salt, and polymer concentrations. Experimental measurements of self-diffusion were performed by fluorescence correlation spectroscopy (FCS), and the results were interpreted using the scaling theory of polyelectrolyte solutions. While all PMAAs were pH sensitive and showed an increase in hydrodynamic radius (R h) with pH in the dilute regime, the R h of star polymers (measured at basic pH values) was significantly smaller for the star polyacids due to their more compact structure. Fully ionized star PMAAs were also found to be less sensitive to changes in salt concentration and type of the counterion compared to linear PMAA. While R h of fully ionized linear PMAA decreased in the series Li+ > Na+ > K+ > Cs+ in agreement with the Hofmeister series, R h of star PMAAs was virtually independent of type of the counterion for eight-arm PMAA. However, molecular architecture strongly affected interactions of counterions with PMAAs. In particular, 7Li NMR revealed that the spin-lattice relaxation time T 1 of Li+ ions in low-salt solutions of eight-arm PMAA was ∼2-fold smaller than that in the solution of linear PMAA, suggesting slower Li+-ion dynamics within star polymers. An increase in concentration of monovalent chloride salts, c s, above that of the PMAA monomer unit concentration (c m) resulted in shrinking of both linear and star molecules, with the hydrodynamic size R h scaling as R h ∝ c s -0.11±0.01. Self-diffusion of linear and star polyelectrolytes was then studied in a wide range of polyelectrolyte concentrations (10-3 mol/L < c m < 0.5 mol/L) in low-salt (<10-4 mol/L of added salt) and high-salt (1 mol/L) solutions. In both the low-salt and high-salt regimes, diffusion coefficient D was lower for PMAAs with a larger number of arms at a fixed c m. In addition, in both cases, D plateaued at low polymer concentrations and decreased at higher polymer concentrations. However, while in the high-salt conditions, the concentration dependence of D reflected transitions between the dilute to semidilute solution regimes as expected for neutral chains in good and theta solvents, analysis of the diffusion data in the low-salt conditions using the scaling theory revealed a different origin of the concentration dependence of D. Specifically, in the low-salt solutions, both linear and star PMAAs exhibited unentangled (Rouse-like) dynamics in the entire range of polyelectrolyte concentrations.

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.

State-Of-The-Art Quantification of Polymer Solution Viscosity for Plastic Waste Recycling

Solvent-based recycling is a promising approach for closed-loop recovery of plastic-containing waste. It avoids the energy cost to depolymerize the plastic but still allows to clean the polymer of contaminants and additives. However, viscosity plays an important role in handling the polymer solutions at high concentrations and in the cleaning steps. This Review addresses the viscosity behavior of polymer solutions, available data, and (mostly algebraic) models developed. The non-Newtonian viscosity models, such as the Carreau and Yasuda-Cohen-Armstrong models, pragmatically describe the viscosity of polymer solutions at different concentrations and shear rate ranges. This Review also describes how viscosity influences filtration and centrifugation processes, which are crucial steps in the cleaning of the polymer and includes a polystyrene/styrene case study.

Intradimer forces and their implication for conformations of von Willebrand factor multimers

The largest blood glycoprotein von Willebrand factor (VWF) responds to hydrodynamic stresses in the bloodstream with abrupt conformation changes, thus increasing its adhesivity to platelets and collagen. Arterial and microvascular hemostasis relies on mechanical and physicochemical properties of this macromolecule. Recently, it was discovered that the mechanical properties of VWF are controlled by multiple pH-dependent interactions with opposite trends within dimeric subunits. In this work, computer simulations reveal the effect of these intradimer forces on the conformation of VWF multimers in various hydrodynamic conditions. A coarse-grained computer model of VWF has been proposed and parameterized to give a good agreement with experimental data. The simulations suggest that strong attraction between VWF D4 domains increases the resistance to elongation under shear stress, whereas even intermediate attraction between VWF C domains contributes to VWF compaction in nonsheared fluid. It is hypothesized that the detailed subdimer dynamics of VWF concatamers may be one of the biophysical regulators of initial hemostasis and arterial thrombosis.

Hydrodynamics of polymers in an active bath

The conformational and dynamical properties of active polymers in solution are determined by the nature of the activity. Here, the behavior of polymers with self-propelled, active Brownian particle-type monomers differs qualitatively from that of polymers with monomers driven externally by colored-noise forces. We present simulation and theoretical results for polymers in solution in the presence of external active noise. In simulations, a semiflexible bead-spring chain is considered, in analytical calculations, a continuous linear wormlike chain. Activity is taken into account by independent monomer or site velocities, with orientations changing in a diffusive manner. In simulations, hydrodynamic interactions (HIs) are taken into account by the Rotne-Prager-Yamakawa tensor or by an implementation of the active polymer in the multiparticle-collision-dynamics approach for fluids. To arrive at an analytical solution, the preaveraged Oseen tensor is employed. The active process implies a dependence of the stationary-state properties on HIs via the polymer relaxation times. With increasing activity, HIs lead to an enhanced swelling of flexible polymers, and the conformational properties differ substantially from those of polymers with self-propelled monomers in the presence of HIs, or free-draining polymers. The polymer mean-square displacement is enhanced by HIs. Over a wide range of timescales, hydrodynamics leads to a subdiffusive regime of the site mean-square displacement for flexible active polymers, with an exponent of 5/7, larger than that of the Rouse (1/2) and Zimm (2/3) models of passive polymers.

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.

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.

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.

Mesoscopic modelling and simulation of soft matter

The deformability of soft condensed matter often requires modelling of hydrodynamical aspects to gain quantitative understanding. This, however, requires specialised methods that can resolve the multiscale nature of soft matter systems. We review a number of the most popular simulation methods that have emerged, such as Langevin dynamics, dissipative particle dynamics, multi-particle collision dynamics, sometimes also referred to as stochastic rotation dynamics, and the lattice-Boltzmann method. We conclude this review with a short glance at current compute architectures for high-performance computing and community codes for soft matter simulation.

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.

Transient hydrodynamic finite-size effects in simulations under periodic boundary conditions

We use lattice-Boltzmann and analytical calculations to investigate transient hydrodynamic finite-size effects induced by the use of periodic boundary conditions. These effects are inevitable in simulations at the molecular, mesoscopic, or continuum levels of description. We analyze the transient response to a local perturbation in the fluid and obtain the local velocity correlation function via linear response theory. This approach is validated by comparing the finite-size effects on the steady-state velocity with the known results for the diffusion coefficient. We next investigate the full time dependence of the local velocity autocorrelation function. We find at long times a crossover between the expected t^{-3/2} hydrodynamic tail and an oscillatory exponential decay, and study the scaling with the system size of the crossover time, exponential rate and amplitude, and oscillation frequency. We interpret these results from the analytic solution of the compressible Navier-Stokes equation for the slowest modes, which are set by the system size. The present work not only provides a comprehensive analysis of hydrodynamic finite-size effects in bulk fluids, which arise regardless of the level of description and simulation algorithm, but also establishes the lattice-Boltzmann method as a suitable tool to investigate such effects in general.

Self-Organized Velocity Pulses of Dense Colloidal Suspensions in Microchannel Flow

We present a numerical study of dense colloidal suspensions in a pressure-driven microchannel flow in two dimensions. The colloids are modeled as elastic and frictional spheres suspended in a Newtonian fluid, which we simulate using the method of multiparticle collision dynamics. The model reproduces periodic velocity and density pulse trains, traveling upstream in the microchannel, which are found in experiments conducted by Isa et al. [Phys. Rev. Lett. 102, 058302 (2009)PRLTAO0031-900710.1103/PhysRevLett.102.058302]. We show that colloid-wall friction and the resultant force chains are crucial for the formation of these pulses. With an increasing colloid density, first solitary jams occur, which become periodic pulse trains at intermediate densities and unstable solitary pulses at high densities. We formulate a phenomenological continuum model and show how these spatiotemporal flow and density profiles can be understood as homoclinic and periodic orbits in traveling-wave equations.

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.

Self-consistent molecular dynamics calculation of diffusion in higher n-alkanes

Diffusion is one of the key subjects of molecular modeling and simulation studies. However, there is an unresolved lack of consistency between Einstein-Smoluchowski (E-S) and Green-Kubo (G-K) methods for diffusion coefficient calculations in systems of complex molecules. In this paper, we analyze this problem for the case of liquid n-triacontane. The non-conventional long-time tails of the velocity autocorrelation function (VACF) are found for this system. Temperature dependence of the VACF tail decay exponent is defined. The proper inclusion of the long-time tail contributions to the diffusion coefficient calculation results in the consistency between G-K and E-S methods. Having considered the major factors influencing the precision of the diffusion rate calculations in comparison with experimental data (system size effects and force field parameters), we point to hydrogen nuclear quantum effects as, presumably, the last obstacle to fully consistent n-alkane description.

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.

Open boundary molecular dynamics of sheared star-polymer melts

Open boundary molecular dynamics (OBMD) simulations of a sheared star polymer melt under isothermal conditions are performed to study the rheology and molecular structure of the melt under a fixed normal load. Comparison is made with the standard molecular dynamics (MD) in periodic (closed) boxes at a fixed shear rate (using the SLLOD dynamics). The OBMD system exchanges mass and momentum with adjacent reservoirs (buffers) where the external pressure tensor is imposed. Insertion of molecules in the buffers is made feasible by implementing there a low resolution model (blob-molecules with soft effective interactions) and then using the adaptive resolution scheme (AdResS) to connect with the bulk MD. Straining with increasing shear stress induces melt expansion and a significantly different redistribution of pressure compared with the closed case. In the open sample, the shear viscosity is also a bit lowered but more stable against the viscous heating. At a given Weissenberg number, molecular deformations and material properties (recoverable shear strain and normal stress ratio) are found to be similar in both setups. We also study the modelling effect of normal and tangential friction between monomers implemented in a dissipative particle dynamics (DPD) thermostat. Interestingly, the tangential friction substantially enhances the elastic response of the melt due to a reduction of the kinetic stress viscous contribution.

Equilibrium Self-Assembly, Structure, and Dynamics of Clusters of Star-Like Micelles

Hierarchical structure and dynamics of clusters of self-assembled star-like micelles formed by oligocarbonate-fluorene end-functionalized poly(ethylene glycol) triblock copolymers were characterized by small-angle neutron scattering and static and dynamic light scattering at concentrations below the gel point. These micelles persist in equilibrium with concentration-dependent sized hierarchical clusters. When probed at length scales within the clusters by dynamic light scattering, the clusters exhibit Zimm dynamics, reminiscent of dilute mesoscale chains. The ability to form chain-like clusters is attributed to the π-π stacking of the fluorene groups that drives the formation of micelles. This enables a design variable to control the rheology of injectable gels. Further, predictions of the solvent (D₂O) viscosity show deviations consistent with polymers in organic solvents, stressing a need for refinement of molecular theories of polymer dynamics.

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.

Hydrodynamics of discrete-particle models of spherical colloids: a multiparticle collision dynamics simulation study

We investigate the hydrodynamic properties of a spherical colloid model, which is composed of a shell of point particles by hybrid mesoscale simulations, which combine molecular dynamics simulations for the sphere with the multiparticle collision dynamics approach for the fluid. Results are presented for the center-of-mass and angular velocity correlation functions. The simulation results are compared with theoretical results for a rigid colloid obtained as a solution of the Stokes equation with no-slip boundary conditions. Similarly, analytical results of a point-particle model are presented, which account for the finite size of the simulated system. The simulation results agree well with both approaches on appropriative time scales; specifically, the long-time correlations are quantitatively reproduced. Moreover, a procedure is proposed to obtain the infinite-system-size diffusion coefficient based on a combination of simulation results and analytical predictions. In addition, we present the velocity field in the vicinity of the colloid and demonstrate its close agreement with the theoretical prediction. Our studies show that a point-particle model of a sphere is very well suited to describe the hydrodynamic properties of spherical colloids, with a significantly reduced numerical effort.