Hydrodynamic Interactions and Entanglements of Polymer Solutions in Many-Body Dissipative Particle Dynamics

Mesoscale Modeling of Agglomeration of Molecular Bottlebrushes: Focus on Conformations and Clustering Criteria

Using dissipative particle dynamics, we characterize dynamics of aggregation of molecular bottlebrushes in solvents of various qualities by tracking the number of clusters, the size of the largest cluster, and an average aggregation number. We focus on a low volume fraction of bottlebrushes in a range of solvents and probe three different cutoff criteria to identify bottlebrushes belonging to the same cluster. We demonstrate that the cutoff criteria which depend on both the coordination number and the length of the side chain allows one to correlate the agglomeration status with the structural characteristics of bottlebrushes in solvents of various qualities. We characterize conformational changes of the bottlebrush within the agglomerates with respect to those of an isolated bottlebrush in the same solvents. The characterization of bottlebrush conformations within the agglomerates is an important step in understanding the relationship between the bottlebrush architecture and material properties. An analysis of three distinct cutoff criteria to identify bottlebrushes belonging to the same cluster introduces a framework to identify both short-lived transient and long-lived agglomerates; the same approach could be further extended to characterize agglomerates of various macromolecules with complex architectures beyond the specific bottlebrush architecture considered herein.

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.

Nanoparticle assembly modulated by polymer chain conformation in composite materials

Mixing nanoparticles into a strategically selected polymer matrix yields nanocomposites with well-controlled microstructures and unique properties and functions. The modulation of nanoparticle assembly by polymer chain conformation can play a dominant role in determining nanocomposite structures, yet such a physical mechanism remains largely unexplored. We hypothesize that highly ordered microdomains of rigid linear polymers provide a template for nanoparticle assembly into open fractal structures. We conducted mesoscopic computer simulations and physical experiments to elucidate how polymer chain conformation regulates the dynamic evolution of nanoparticle structures during the drying processing of polymer nanocomposite films. The evaporation of polymer-nanoparticle mixtures with varying chain stiffnesses was simulated using dissipative particle dynamics. The formation of distinguished nanoparticle assemblies as a result of matrix selection was further corroborated by probing nanoparticle aggregation in different polymer nanocomposite coatings. The results show that polymer conformation not only influences the dispersion states of individual particles (dispersed vs. aggregated), but also modulates the morphologies of large-scale assembly (globular vs. fractal). The emergence of nematically ordered polymer clusters when the chain rigidity is increased creates local solvent-rich "voids" that promote anisotropic particle aggregates, which then percolate into open fractal structures upon solvent evaporation. The nanoparticle dynamics also exhibits an intriguing non-monotonic behavior attributed to the transitions between the coupling and decoupling with polymer dynamics. The nanoparticle assembly morphologies obtained in simulations match well with the electron microscopy images taken in physical experiments.

Invariance of experimental observables with respect to coarse-graining in standard and many-body dissipative particle dynamics

Dissipative particle dynamics (DPD) is a well-established mesoscale simulation method. However, there have been long-standing ambiguities regarding the dependence of its (purely repulsive) force field parameter on temperature as well as the variation of the resulting experimental observables, such as diffusivity or surface tension, with coarse-graining (CG) degree. Here, we rederive the temperature dependence of DPD interaction parameter and revisit the role of the CG degree in standard DPD simulations. Consequently, we derive a scaling of the input variables that renders the system properties invariant with respect to CG degree and illustrate the versatility of the method by computing the surface tensions of binary solvent mixtures. We then extend this procedure to many-body dissipative particle dynamics and, by computing surface tensions of the same mixtures at a range of CG degrees, demonstrate that this newer method, which has not been widely applied so far, is also capable of simulating complex fluids of practical interest.

Morphology evolution of Janus dumbbell nanoparticles in seeded emulsion polymerization

Emulsion polymerization is a versatile approach to produce different polymeric nanoparticle morphologies, which can be useful in a variety of applications. However, the detailed mechanism of the morphology formation is not entirely clear. We study the kinetics of nanoparticle morphology evolution during a seeded emulsion polymerization using both experimental and computational tools. Lightly crosslinked polystyrene seeds were first synthesized using dispersion polymerization. Then the seed particles were swollen in tert-butyl acrylate and styrene monomers, and subsequently polymerized into nanoparticles of dumbbell and multilobe morphologies. It was discovered that both the seed and final particle morphology were affected by the methanol concentration during the seed synthesis. Systematically adjusting the methanol amount will not only yield spherical seed particles of different size, but also dumbbell particles even without the second monomer polymerization. In addition to methanol concentration, morphology can be controlled by crosslinking density. The kinetics studies revealed an interesting transition from multilobe to dumbbell geometries during the secondary polymerization. Based on the results, a nucleation-growth model has been proposed to describe the morphology evolution and verification was offered by computer simulation. The key discovery is that nanoparticle morphology can be kinetically controlled by diffusion of the protrusions on the seed particles. The condition of seed synthesis and crosslinking density will drastically change the seed and final nanoparticle morphology.

Thinning Approximation for Calculating Two-Dimensional Scattering Patterns in Dissipative Particle Dynamics Simulations under Shear Flow

Modifications to improve thinning approximation (TA) were considered in order to calculate two-dimensional scattering patterns (2DSPs) for dissipative particle dynamics (DPD) simulations of polymer melts under a shear flow. We proposed multipoint TA and adaptive TA because the bond lengths in DPD chains vary widely when compared to those in Kremer⁻Grest (KG) chains, and the effectiveness of these two types of TA for the two major DPD parameter sets were investigated. In this paper, we report our findings on the original DPD model with soft bonds and that with rigid bonds. Based on the behavior of the 2DSPs and the distribution of orientations of the bond vectors, two spot patterns originating from the oriented chain correlations were observed when distinct distributions of the highly oriented bond vectors in the shear direction were obtained. For multipoint TA, we concluded that at least two additional midpoints ( n mid ≥ 2 ) are required to clearly observe the two spot patterns. For adaptive TA, a dividing distance of l ATA ≤ 0.4 is sufficient for clear observation, which is consistent with the requirement of n mid ≥ 2 for multipoint TA.

Multipoint segmental repulsive potential for entangled polymer simulations with dissipative particle dynamics

A model is developed for simulating entangled polymers by dissipative particle dynamics (DPD) using the segmental repulsive potential (SRP). In contrast to previous SRP models that define a single-point interaction on each bond, the proposed SRP model applies a dynamically adjustable multipoint on the bond. Previous SRP models could not reproduce the equilibrium properties of Groot and Warren's original DPD model [R. D. Groot and P. B. Warren, J. Chem. Phys. 107, 4423 (1997)] because the introduction of a single SRP induces a large excluded volume, whereas, the proposed multipoint SRP (MP-SRP) introduces a cylindrical effective excluded bond volume. We demonstrate that our MP-SRP model exhibits equilibrium properties similar to those of the original DPD polymers. The MP-SRP model parameters are determined by monitoring the number of topology violations, thermodynamic properties, and the polymer internal structure. We examine two typical DPD polymers with different bond-length distributions; one of them was used in the modified SRP model by Sirk et al. [J. Chem. Phys. 136, 134903 (2012)], whereas the other was used in the original DPD model. We demonstrate that for both polymers, the proposed MP-SRP model captures the entangled behaviors of a polymer melt naturally, by calculating the slowest relaxation time of a chain in the melt and the shear relaxation modulus. The results indicate that the proposed MP-SRP model can be applied to a variety of DPD polymers.

Self-healing and dewetting dynamics of a polymer nanofilm on a smooth substrate: strategies for dewetting suppression

The self-healing and dewetting dynamics of a polymer nanofilm on a smooth, partial wetting surface are explored by many-body dissipative particle dynamics. Three types of dewetting phenomena are identified, (i) spinodal decomposition, (ii) nucleation and growth, and (iii) metastable self-healing. The outcome depends on the surface wettability (θY), the polymer film thickness (h0), and the radius of the dry hole (R0). The phase diagram of the dewetting mechanism as a function of θY and h0 is obtained for a specified R0. As the surface wettability decreases (increasing θY), the critical film thickness associated with the nucleation/self-healing crossover (hc) grows so that the metastability of the film can be retained by the self-healing process. In addition to θY and R0, hc depends on the polymer length (N) as well. It is found that a longer polymer requires a thicker nanofilm to avoid dewetting by nucleation. Two strategies for dewetting suppression are proposed. The metastability of a film of polymers with a large molecular weight can be promoted either by the addition of short polymers or by employing compact polymers such as star polymers. In the latter approach, the increment of the arm number enhances the nanofilm stability.

An early mechanical coupling of planktonic bacteria in dilute suspensions

It is generally accepted that planktonic bacteria in dilute suspensions are not mechanically coupled and do not show correlated motion. The mechanical coupling of cells is a trait that develops upon transition into a biofilm, a microbial community of self-aggregated bacterial cells. Here we employ optical tweezers to show that bacteria in dilute suspensions are mechanically coupled and show long-range correlated motion. The strength of the coupling increases with the growth of liquid bacterial culture. The matrix responsible for the mechanical coupling is composed of cell debris and extracellular polymer material. The fragile network connecting cells behaves as viscoelastic liquid of entangled extracellular polymers. Our findings point to physical connections between bacteria in dilute bacterial suspensions that may provide a mechanistic framework for understanding of biofilm formation, osmotic flow of nutrients, diffusion of signal molecules in quorum sensing, or different efficacy of antibiotic treatments at low and high bacterial densities.Planktonic bacteria are untethered to surfaces or to each other, and thus are expected to move independently when at low cell densities. Here Sretenovic et al. show, using optical tweezers, that bacteria in dilute suspensions are mechanically coupled and show long-range correlated motion.

Interfacial adsorption of pH-responsive polymers and nanoparticles

Using dissipative particle dynamics (DPD), we model the interfacial adsorption of pH-responsive polyelectrolytes and polyelectrolyte-grafted nanoparticles (PNPs) at a planar water-oil interface. The electrostatic interactions in the presence of the dielectric discontinuity across the interface are modeled by exploiting the Groot method, which uses an iterative method to solve the Poisson equation on a uniform grid with distributed charge. We reveal the effects of the pH and salinity of the aqueous solution and the length of the polyelectrolyte on the adsorption behavior of weak polyelectrolytes. The adsorption kinetics is monitored via the trajectory of the center of mass of the polyelectrolyte in the direction normal to the interface. The residence time at the interface and the pair correlation function between the polyelectrolyte and the oil are measured to quantitatively characterize the adsorption. Similar to the weak polyelectrolytes, the influences of pH, salinity and grafted chain length on the adsorption of an individual PNP are explored. Our results show that by grafting polyelectrolytes, the interfacial behavior of the nanoparticles can be tuned by changing the pH and salinity of the solution, which is dictated by the contact angle, the pair correlation function between the particles and the oil, the desorption energy, and the particle morphology at the interface. We also observe that the electrostatic-driven variations in the interfacial activity and morphology of the PNPs are not sensitive to the length of the grafted polyelectrolytes.