The Lowe-Andersen thermostat as an alternative to the dissipative particle dynamics in the mesoscopic simulation of entangled polymers

Slip and momentum transfer mechanisms mediated by Janus rods at polymer interfaces

As an incipient but preeminent technology for multiphase nanomaterials/fluids, exact compatibilizing mechanisms of Janus particles in polymer blends and the consequent morphology remain unknown. The contributions of Janus nanorods to slip suppression and momentum transfer across the interface have been explored through dissipative particle dynamics simulations under shear flow at unentangled polymer-polymer interfaces. Rods have been then grafted with flexible polymer chains to unveil interfacial structure-property relationships at a molecular level when compared with flexible diblock copolymer surfactants. When Janus rods are sparsely grafted with necessarily longer grafts, they favor a greater degree of graft interpenetration with polymer phases. This yields less effective momentum transfer that impacts droplet coalescence processes; dynamic heterogeneities at complex interfaces; and helps map their efficiency as compatibilizers.

Dynamic Interfacial Trapping of Janus Nanorod Aggregates

Taking advantage of both shape and chemical anisotropy on the same nanoparticle offers rich self-assembly possibilities for nanotechnology. Through dissipative particle dynamics calculations, in the present work, the directed assembly of Janus nanorod aggregates and their capability to assemble into metastable novel structures at an interfacial level have been assessed. Symmetric Janus rods become kinetically trapped and exhibit either parallel or antiparallel alignment with respect to their long axis (different compositions). This depends on several factors that have been mapped herein and that can be precisely tuned: Flory-Huggins interaction parameter χ between polymer phases; concentration; shear rate; and even aggregate shape. Ultimately, two different aggregate structures result from rod tumbling that are not observed under quiescent conditions: monolayer-like aggregates exhibiting trapped rods with antiparallel configuration; and stacked nanorod arrays similar to superlattice sheets. These different structures can be controlled by the likelihood with which tumbling Janus rods encounter other aggregate portions showing parallel alignment. Hence, the present study offers fundamental insight into relevant parameters that govern the directed assembly of Janus nanoparticles at an interfacial level. Novel applications may potentially derive from the resulting aggregate structures, such as peculiar displays and sensors.

Interfacial aggregation of Janus rods in binary polymer blends and their effect on phase separation

Janus particles interfacially self-assemble into different structures when incorporated into multiphase systems. Dissipative particle dynamics simulations are employed herein to investigate the interplay between aggregation mechanisms and phase separation in polymer blends. Shorter rods with a standing configuration become increasingly "caged" or trapped in larger aggregates as weight fraction increases, which is reflected in the way that their diffusion is coupled to their aggregation rates. Janus rods of higher aspect ratios that are tilted at the interface aggregate side-by-side and are able to hinder phase separation kinetics. This is due to a combination of individual Janus rod conformations at the interface, their intrinsic aggregation mechanisms, aggregate fractal dimension, and aggregation rates, and can also be traced back to the scaling of the diffusion coefficient of aggregates with their size. Findings presented provide insight into the mechanisms governing two dimensionally growing colloidal aggregates at fluid interfaces, more specifically, those associated with Janus particles, and shed light on the potential of these systems in paving the way for designing new functional materials.

Modeling of Entangled Polymer Diffusion in Melts and Nanocomposites: A Review

This review concerns modeling studies of the fundamental problem of entangled (reptational) homopolymer diffusion in melts and nanocomposite materials in comparison to experiments. In polymer melts, the developed united atom and multibead spring models predict an exponent of the molecular weight dependence to the polymer diffusion very similar to experiments and the tube reptation model. There are rather unexplored parameters that can influence polymer diffusion such as polymer semiflexibility or polydispersity, leading to a different exponent. Models with soft potentials or slip-springs can estimate accurately the tube model predictions in polymer melts enabling us to reach larger length scales and simulate well entangled polymers. However, in polymer nanocomposites, reptational polymer diffusion is more complicated due to nanoparticle fillers size, loading, geometry and polymer-nanoparticle interactions.

Energy Renormalization Method for the Coarse-Graining of Polymer Viscoelasticity

Developing temperature transferable coarse-grained (CG) models is essential for the computational prediction of polymeric glass-forming (GF) material behavior, but their dynamics are often greatly altered from those of all-atom (AA) models mainly because of the reduced fluid configurational entropy under coarse-graining. To address this issue, we have recently introduced an energy renormalization (ER) strategy that corrects the activation free energy of the CG polymer model by renormalizing the cohesive interaction strength ε as a function of temperature T, i.e., ε(T), thus semiempirically preserving the T-dependent dynamics of the AA model. Here we apply our ER method to consider-in addition to T-dependency-the frequency f-dependent polymer viscoelasticity. Through smallamplitude oscillatory shear molecular dynamics simulations, we show that changing the imposed oscillation f on the CG systems requires changes in ε values (i.e., ε(T, f)) to reproduce the AA viscoelasticity. By accounting for the dynamic fragility of polymers as a material parameter, we are able to predict ε(T, f) under coarse-graining in order to capture the AA viscoelasticity, and consequently the activation energy, across a wide range of T and f in the GF regime. Specifically, we showcase our achievements on two representative polymers of distinct fragilities, polybutadiene (PB) and polystyrene (PS), and show that our CG models are able to sample viscoelasticity up to the megahertz regime, which approaches state-of-the-art experimental resolutions, and capture results sampled via AA simulations and prior experiments.

Polymer-mediated nanorod self-assembly predicted by dissipative particle dynamics simulations

Self-assembly of nanoparticles in polymer matrices is an interesting and growing subject in the field of nanoscience and technology. We report herein on modelling studies of the self-assembly and phase behavior of nanorods in a homopolymer matrix, with the specific goal of evaluating the role of deterministic entropic and enthalpic factors that control the aggregation/dispersion in such systems. Grafting polymer brushes from the nanorods is one approach to control/impact their self-assembly capabilities within a polymer matrix. From an energetic point of view, miscible interactions between the brush and the matrix are required for achieving a better dispersibility; however, grafting density and brush length are the two important parameters in dictating the morphology. Unlike in previous computational studies, the present Dissipative Particle Dynamics (DPD) simulation framework is able to both predict dispersion or aggregation of nanorods and determine the self-assembled structure, allowing for the determination of a phase diagram, which takes all of these factors into account. Three types of morphologies are predicted: dispersion, aggregation and partial aggregation. Moreover, favorable enthalpic interactions between the brush and the matrix are found to be essential for expanding the window for achieving a well-dispersed morphology. A three-dimensional phase diagram is mapped on which all the afore-mentioned parameters are taken into account. Additionally, in the case of immiscibility between brushes and the matrix, simulations predict the formation of some new and tunable structures.