Mesoscopic Simulation of Entangled Polymer Brushes under Shear: Compression and Rheological Properties

All-atom molecular dynamics simulations of polymer and polyelectrolyte brushes

Densely grafted polymer and polyelectrolyte (PE) brushes, owing to their significant abilities to functionalize surfaces for a plethora of applications in sensing, diagnostics, current rectification, surface wettability modification, drug delivery, and oil recovery, have attracted significant attention over the past several decades. Unfortunately, most of the attention has primarily focused on understanding the properties of the grafted polymer and the PE chains with little attention devoted to studying the behavior of the brush-supported ions (counterions needed to screen the PE chains) and water molecules. Over the past few years, our group has been at the forefront of addressing this gap: we have employed all-atom molecular dynamics (MD) simulations for studying a wide variety of polymer and PE brush systems with specific attention to unraveling the properties and behavior of the brush-supported water molecules and ions. Our findings have revealed some of the most fascinating properties of such brush-supported ions and water molecules, including the most remarkable control of nanofluidic transport afforded by the specific ion and water responses induced by the PE brushes grafted on the inner walls of the nanochannel. This feature article aims to summarize some of our key contributions associated with such atomistic simulations of polymer and PE brushes and brush-supported water molecules and counterions.

Future Trends in Obolodiplosis robiniae Distribution across Eurasian Continent under Global Climate Change

Obolodiplosis robiniae was discovered in Eurasia at the beginning of the 21st century. In this study, we explore the present and future (in the years 2050 and 2070) trends in the potential distribution of O. robiniae in Eurasia under diverse climate change scenarios based on a maximum entropy model. Our findings indicated that the current potential distribution area of O. robiniae is within the range of 21°34' and 65°39' N in the Eurasian continent. The primary factor controlling the distribution of O. robiniae is temperature. The highly and moderately suitable areas are mainly distributed in the semi-humid and semi-arid regions, which also happen to be the locations where the host black locust (Robinia pseudoacacia L.) grows at its fastest rate. The forecast of the potential distribution area of O. robiniae revealed that the species would benefit from global warming. The region suitable for the habitat of O. robiniae is characterized by a large-scale northward expansion trend and an increase in temperature. This information would help the forestry quarantine departments of Asian and European countries provide early warnings on the probable distribution areas of O. robiniae and provide a scientific basis for the prevention and control of O. robiniae spread and outbreaks.

Dissipative particle dynamics simulations in colloid and Interface science: a review

Dissipative particle dynamics (DPD) is one of the most efficient mesoscale coarse-grained methodologies for modeling soft matter systems. Here, we comprehensively review the progress in theoretical formulations, parametrization strategies, and applications of DPD over the last two decades. DPD bridges the gap between the microscopic atomistic and macroscopic continuum length and time scales. Numerous efforts have been performed to improve the computational efficiency and to develop advanced versions and modifications of the original DPD framework. The progress in the parametrization techniques that can reproduce the engineering properties of experimental systems attracted a lot of interest from the industrial community longing to use DPD to characterize, help design and optimize the practical products. While there are still areas for improvements, DPD has been efficiently applied to numerous colloidal and interfacial phenomena involving phase separations, self-assembly, and transport in polymeric, surfactant, nanoparticle, and biomolecules systems.

Backbone oriented anisotropic coarse grains for efficient simulations of polymers

Despite the fact that anisotropic particles have been introduced to describe molecular interactions for decades, they have been poorly used for polymers because of their computing time overhead and the absence of a relevant proof of their impact in this field. We first report a method using anisotropic beads for polymers, which solves the computing time issue by considering that beads keep their principal orientation alongside the mean local backbone vector of the polymer chain, avoiding the computation of torques during the dynamics. Applying this method to a polymer bulk, we study the effect of anisotropic interactions vs isotropic ones for various properties such as density, pressure, topology of the chain network, local structure, and orientational order. We show that for different classes of potentials traditionally used in molecular simulations, those backbone oriented anisotropic beads can solve numerous issues usually encountered with isotropic interactions. We conclude that the use of backbone oriented anisotropic beads is a promising approach for the development of realistic coarse-grained potentials for polymers.

Shear-induced microstructures and dynamics processes of phospholipid cylinders in solutions

Shear-induced microstructures and their corresponding dynamic processes are investigated for phospholipid cylinders in aqueous solution by dissipative particle dynamic simulation. Various phospholipid cylinders with cross-sections, which are formed under shear-free flow, are selected to examine the effects of shear flow on their structures and dynamic processes. Shear flow induces the transition from cylinders into vesicles at weak rate and the transition into vesicle–lamella mixtures with increased shear rate and lamella structures at the strong shear rate. Then, the average radius of gyration and shape factors of the polymer chains in the dynamic processes are discussed in detail. Results show that shear flow causes the structure of the polymer chains to be elongated along the shear direction, and the configuration of the polymer chain can be rapidly transformed into an ellipsoid structure under strong shear.

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.

Tribological Behavior of Grafted Nanoparticle on Polymer-Brushed Walls: A Dissipative Particle Dynamics Study

Two contacting surfaces grafted with polymer brushes have potential applications due to their extraordinary lubricating behavior. However, the polymer brushes may have poor mechanical stability under high normal and shear stresses, which is a challenge for practical usage of polymer brush systems. In this study, we propose the use of grafted nanoparticles as nanobearings on polymer-brush-coated surfaces to alleviate the harsh working conditions of polymer brushes and to improve their mechanical stability. We have performed dissipative particle dynamics (DPD) simulations to investigate the tribological interaction between grafted nanoparticle and parallel walls with noncharged polymer brushes in the presence of explicit solvent. The influences of several parameters (solvent quality, brush miscibility, etc.) on the tribological behavior of the system are investigated. The grafted nanoparticle obviously acts as a nanobearing that partially replaces the sliding contact between two brushed walls with rolling contact between the grafted nanoparticle and two brushed walls and reduces the number of DPD particles withstanding high force. Although the introduction of the grafted nanoparticle into polymer-brushed walls increases the friction coefficient by 20-30%, it does not greatly decrease lubrication of the brushed walls, while still helping in stabilizing the system of polymer brushes to be used with liquids with low viscosity, such as water. The DPD simulation results and analysis performed in this study would be beneficial in designing systems with polymer-brushed surfaces and grafted nanoparticles.

Self-assembly of phospholipid molecules in solutions under shear flows: Microstructures and phase diagrams

Shear-induced microstructures and their phase diagrams were investigated for phospholipid molecules in aqueous solution by dissipative particle dynamic simulation. Self-assembled microstructures, including spherical and cylindrical micelles, spherical vesicles, lamellae, undulated lamellae, perforated lamellae, and continuous networks, were observed under various shear flows and phospholipid concentrations, where the spatial inhomogeneity and symmetry were analysed. A series of phase diagrams were constructed based on the chain lengths under various phospholipid concentrations. The phase distributions showed that the structures with spherical symmetry could be shear-induced to structures with cylindrical symmetry in the dilute solutions. In the semi-concentrated solutions, the lamellae were located in most spaces under zero shear flows, which could be shear-induced into undulated lamellae and then into cylindrical micelles. For the concentrated solutions, the strong shear flows oriented the directions of multilayer lamellae and phase transitions appeared between several cylindrical network structures. These observations on shear-induced microstructures and their distributions revealed a promising approach that could be used to design bio-microstructures based on phospholipid molecules under shear flows.

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.

Poly(ethylene oxide) Mushrooms Adsorbed at Silica-Ionic Liquid Interfaces Reduce Friction

The adsorbed layer conformation and lubricity of 35, 100, and 300 kDa PEO adsorbed to ionic liquid (IL)-silica interfaces from 0.01 wt % solutions have been investigated using colloid probe atomic force microscopy. The ILs used were propylammonium nitrate (PAN) and 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF4]), which are protic and aprotic ILs, respectively. Normal force curves reveal steric interactions consistent with adsorbed polymer layers which are best fit using the mushroom model. Friction measurements show that the adsorbed polymer layer markedly reduces friction compared to surfaces sliding in the pure ILs and that lubricity increases with polymer length. When polymer is adsorbed to the sliding surfaces, friction is controlled by the creation and disruption of intermolecular interactions between entangled chains and the dragging of polymer chains through the interpenetration region. These experiments show that added polymer can reduce friction while maintaining the useful properties of ILs as lubricants.

Structure and rheology of star polymers in confined geometries: a mesoscopic simulation study

Mesoscopic simulations of star polymer melts adsorbed onto solid surfaces are performed using the dissipative particle dynamics (DPD) method. A set of parameters is developed to study the low functionality star polymers under shear. The use of a new bond-angle potential between the arms of the star creates more rigid chains and discriminates between different functionalities at equilibrium, but still allows the polymers to deform appropriately under shear. The rheology of the polymer melts is studied by calculating the kinetic friction and viscosity and there is good agreement with experimental properties of these systems. The study is completed with predictive simulations of star polymer solutions in an athermal solvent.

Prediction of structural and thermomechanical properties of polymers from multiscale simulations

We report mesoscale simulations of polymer melts and crosslinked polymer networks by using realistic coarse-grained (CG) models that are developed from atomistic simulations of polymer melts.

Mesoscopic modeling of structural and thermodynamic properties of fluids confined by rough surfaces

Simulations show that the ordering of particles confined by rough surfaces induces a structural phase transition while the interfacial tension is insensitive to it.

Impact of chain morphology on the lubricity of surface-grafted polysaccharides

The impact of brush-like structure and disorder on the lubricating ability of dextran chains at low and high loads has been investigated by means of AFM.

Load-induced transitions in the lubricity of adsorbed poly(L-lysine)-g-dextran as a function of polysaccharide chain density

Chain-density gradients of poly(l-lysine)-graft-dextran (PLL-g-dex), a synthetic comblike copolymer with a poly(l-lysine) backbone grafted with dextran side chains, were fabricated on an oxidized silicon substrate. The influence of the changing dextran chain density along the gradient on the local coefficient of friction was investigated via colloidal-probe lateral force microscopy. Both in composition and structure, PLL-g-dex shares many similarities with bottlebrush biomolecules present in natural lubricating systems, while having the advantage of being well-characterized in terms of both architecture and adsorption behavior on negatively charged oxide surfaces. The results indicate that the transition of the dextran chain density from the mushroom into the brush regime coincides with a sharp reduction in friction at low loads. Above a critical load, the friction increases by more than an order of magnitude, likely signaling a pressure-induced change in the brush conformation at the contact area and a corresponding change in the mechanism of sliding. The onset of this higher-friction regime is moved to higher loads as the chain density of the film is increased. While in the low-load (and low-friction) regime, increased chain density leads to lower friction, in the high-load (high-friction) regime, increased chain density was found to lead to higher friction.

Solvent-induced self-assembly of polymer-tethered nanorods

Self-assembly behaviors of polymer-tethered nanorods in the selective solvent are systematically investigated via a dissipative particle dynamics (DPD) simulation method. Three types of polymer-tethered nanorods are considered: one end tethered, both ends tethered, and middle tethered. The solvent-induced diverse morphologies and morphological transitions depend on the topology, rod/tether length ratio, solvent selectivity, and mixed solvent content. In the pure rod-selective solvent (solvent I) or the pure tether-selective solvent (solvent II), the ordered micellar structures include: cylinders, hexagonal cylinders, bilayer lamellae, lamellae/cylinder mixed phases, inverted hollow cylinders, and nematic bundles. These micelles are formed by the competition among the stretching of tethers, liquid crystalline of rods, interfacial energy, and solvent selectivity. In the I/II mixed solvent, with varying mixed solvent content in sequence (i.e., changing the solvent quality for the blocks), the reversible morphological transitions and fantastic intermediate phases (e.g., liquid crystalline phase) are observed, which correspond directly to the case of that induced by varying the rod/tether length ratio in the pure solvent. It is concluded that improving the selective solvent content is equivalent to increasing the soluble block ratio. The present study reveals that the morphology and morphological transition of polymer-tethered nanorods could be significantly manipulated through topology, block length, and solvent, especially the selectivity.

Calculation of the surface tension from multibody dissipative particle dynamics and Monte Carlo methods

We report the calculation of the coexisting densities and surface tensions of the liquid-vapor equilibrium using the multibody dissipative particle dynamics and Monte Carlo (MMC) methods. We focus on the calculation of the surface tension by using the thermodynamic and mechanical routes. It is the first time that the test-area method is applied on the many-body conservative potential. We discuss the mechanical equilibrium of these two-phase systems by analyzing the profiles of the normal and tangential components of the pressure tensor using the Irving-Kirkwood and Kirkwood-Buff approaches. The profile of the configurational temperature is shown to establish the thermal equilibrium of these two-phase simulations carried out with large time steps. We complete this study to show the impact of the range of the many-body repulsive term of the conservative force on the surface tension. We conclude that the MMC method is an efficient sampling scheme to compute the interfacial properties of liquid-vapor interfaces using the multibody soft potential.