Conformational and dynamical properties of ultra-soft colloids in semi-dilute solutions under shear flow

Structural transitions of a semi-flexible polyampholyte

Polyampholytes (PAs) are charged polymers composed of positively and negatively charged monomers along their backbone. The sequence of the charged monomers and the bending of the chain significantly influence the conformation and dynamical behavior of the PA. Using coarse-grained molecular dynamics simulations, we comprehensively study the structural and dynamical properties of flexible and semi-flexible PAs. The simulation results demonstrate a flexible PA chain, displaying a transition from a coil to a globule in the parameter space of the charge sequence. In addition, the behavior of the mean-square displacement (MSD), denoted as ⟨(Δr(t))2⟩, reveals distinct dynamics, specifically for the alternating and charge-segregated sequences. The MSD follows a power-law behavior, where ⟨(Δr(t))2⟩ ∼ tβ, with β ≈ 3/5 and β ≈ 1/2 for the alternating sequence and the charge-segregated sequence in the absence of hydrodynamic interactions, respectively. However, when hydrodynamic interactions are incorporated, the exponent β shifts to ∼3/5 for the charge-segregated sequence and 2/3 for the well-mixed alternating sequence. For a semi-flexible PA chain, varying the bending rigidity and electrostatic interaction strength (Γe) leads to distinct, fascinating conformational states, including globule, bundle, and torus-like conformations. We show that PAs acquire circular and hairpin-like conformations in the intermediate bending regime. The transition between various conformations is identified in terms of the shape factor estimated from the ratios of eigenvalues of the gyration tensor.

Explicit characterization of counterion dynamics around a flexible polyelectrolyte

The article presents a comprehensive study of counterion dynamics around a generic linear polyelectrolyte chain with the help of coarse-grained computer simulations. The ion-chain coupling is discussed in the form of binding time, mean-square displacement (MSD) relative to the chain, local ion transport coefficient, and spatiotemporal correlations in the effective charge. We have shown that a counterion exhibits subdiffusive behavior 〈δR^{2}〉∼t^{δ}, δ≈0.9 w.r.t. chain's center of mass. The MSD of ions perpendicularly outward from the chain segment exhibits a smaller subdiffusive exponent compared to the one along the chain backbone. The effective diffusivity of ion is the lowest in chain's close proximity, extending up to the length-scale of radius of gyration R_{g}. Beyond R_{g} at larger distances, they attain diffusivity of free ion with a smooth cross-over from the adsorbed regime to the free ion regime. We have shown that the effective diffusivity drastically decreases for the multivalent ions, while the crossover length scale remains the same. Conversely, with increasing salt concentration the coupling-length scale reduces, while the diffusivity remains unaltered. The effective diffusivity of adsorbed-ion reveals an exponential reduction with electrostatic interaction strength. We further corroborate this from the binding time of ions on the chain, which also grows exponentially with the coupling strength of the ion-polymer duo. Moreover, the binding time of ions exhibits a weak dependence with salt concentration for the monovalent salt, while for multivalent salts the binding time decreases dramatically with concentration. Our work also elucidates fluctuations in the effective charge per site, where it exhibits strong negative correlations at short length-scales.

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.

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.

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.

Network formation and gelation in telechelic star polymers

We investigate the efficiency of gelation and network formation in telechelic star polymer melt, where the tips of polymer arms are dipoles while the rest of the monomers are uncharged. Our work is motivated by the experimental observations [A. Kulkarni et al., Macromolecules 48, 6580 (2015)] in which rheological studies of telechelic star polymers of poly-(L-lactide), a bio-degradable polymer, showed a drastic increase in elastic properties (up to 2000 times) compared to corresponding star polymers without the telechelic arm ends. In contrast to previous studies, we avoid using effective attractive Lennard-Jones potentials or dipolar potentials to model telechelic interactions. Instead we use explicit Coulomb positive and negative charges at the tip of polymer-arms of our bead-spring model of star polymers. By our simulations we show that the dipoles at the tip of star arms aggregate together to form clusters of dipoles. Each cluster has contributions from several stars, and in turn each star contributes to several clusters. Thus the entire polymer melt forms a connected network. Network forming tendencies decrease with a decrease of the value of the effective charge constituting the dipole: this can be experimentally realized by choosing a different ionomer for the star tip. We systematically varied the value of dipole charges, the fraction of star-arms with dipoles at the tip, and the length of the arms. The choice of explicit charges in our calculations enables us to make better quantitative predictions about the onset of gelation; moreover we get qualitatively distinct results about structural organization of dipoles within a dipole-cluster.

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.

Pulling-force-induced elongation and alignment effects on entanglement and knotting characteristics of linear polymers in a melt

We employ a primitive path (PP) algorithm and the Gauss linking integral to study the degree of entanglement and knotting characteristics of linear polymer model chains in a melt under the action of a constant pulling force applied to selected chain ends. Our results for the amount of entanglement, the linking number, the average crossing number, the writhe of the chains and their PPs and the writhe of the entanglement strands all suggest a different response at the length scale of entanglement strands than that of the chains themselves and of the corresponding PPs. Our findings indicate that the chains first stretch at the level of entanglement strands and next the PP (tube) gets oriented with the "flow." These two phases of the extension and alignment of the chains coincide with two phases related to the disentanglement of the chains. Soon after the onset of external force the PPs attain a more entangled conformation, and the number of nontrivially linked end-to-end closed chains increases. Next, the chains disentangle continuously to attain an almost unentangled conformation. Using the linking matrix of the chains in the melt, we furthermore show that these phases are accompanied by a different scaling of the homogeneity of the global entanglement in the system. The homogeneity of the end-to-end closed chains first increases to a maximum and then decreases slowly to a value characterizing a completely unlinked system.

A computational and experimental study of the linear and nonlinear response of a star polymer melt with a moderate number of unentangled arms

We present from simulations and experiments results on the linear and nonlinear rheology of a moderate functionality, low molecular weight unentangled polystyrene (PS) star melt. The PS samples were anionically synthesized and close to monodisperse while their moderate functionality ensures that they do not display a pronounced core effect. We employ a highly coarse-grained model known as Responsive Particle Dynamics where each star polymer is approximated as a point particle. The eliminated degrees of freedom are used in the definition of an appropriate free energy as well as describing the transient pair-wise potential between particles that accounts for the viscoelastic response. First we reproduce very satisfactorily the experimental moduli using simulation. We then consider the nonlinear response of the same polymer melts by implementing a start-up shear protocol for a wide range of shear rates. As in experiments, we observe the development of a stress overshoot with increasing shear rate followed by a steady-state shear stress. We also recover the shear-thinning nature of the melt, although we slightly overestimate the extent of shear-thinning with simulations. In addition, we study relaxations upon the removal of shear where we find encouraging agreement between experiments and simulations, a finding that corroborates our agreement for the linear rheology.