Friction between brush layers of charged and neutral bottle-brush macromolecules. molecular dynamics simulations

Using molecular dynamics simulations, we study the lubricating properties of neutral and charged bottle-brush coatings as a function of the compression and shear stresses and brush grafting density. Our simulations have shown that in charged bottle-brush systems under shear there is a layer with excess counterions located in the middle between brush-bearing surfaces. The main deformation mode of the charged bottle-brush layers is associated with the backbone deformation, resulting in the backbone deformation ratio, α, and shear viscosity, η, being universal functions of the Weissenberg number. In the case of neutral bottle-brush systems, in addition to the backbone deformation there is also side chain deformation. The coupling between backbone and side chain deformation violates universality in the deformation ratio, α, dependence on the Weissenberg number and results in scaling exponents varying with the compression stress and brush grafting density. The existence of different length scales controlling deformation of neutral bottle brushes manifests itself in the shear viscosity, η, dependence on the shear rate, ?γ. Shear viscosity, η, as a function of the shear rate, ?γ, has two plateaus and two shear thinning regimes. The low shear rate plateau and shear thinning regime correspond to the backbone deformation, while the second plateau and shear thinning regime at moderate shear rates are due to side chain deformation. For both systems the value of the friction coefficient increases with increasing shear rate. The values of the friction coefficient for charged bottle-brush systems are about ten times smaller than corresponding values for neutral systems at the same shear rate.     

The Stability of the Micelle Formed by Chain Branch Surfactants and Polymer Under Salt and Shear Force: Insight from Dissipative Particle Dynamics Simulation

The influence of salt and shear force on the stability of the micelle formed by surfactants and polymer are studied using dissipative particle dynamics (DPD) simulation method. The research system mainly includes four types of surfactants with different hydrophilic/hydrophobic chain branches and two kinds of polymers with hydrophilic/hydrophobic properties, respectively. The stability of the micelle is studied based on the analyses of the density peak and root mean square (RMS) of polymer chain under different salt and shear force. The calculated results show that the density peak reduced and RMS increased for all surfactants with the salt concentration and shear force increasing, and then indicate that the micelle has a certain degree of deformation. Whereas, the surfactant chain branch has important influence on the deformation extent of the micelle. For hydrophobic polymer, surfactants containing hydrophobic chain branch (T2H2T2) are beneficial to the stability of the micelle. On the contrary, for hydrophilic polymer, the micelle formed by surfactants with stronger hydrophilic nature such as the hydrophilic groups located in the both ends of the molecule (H1T4H1) have the best salt and shear resistance. The results have certain theoretical significance and can provide theoretical support for the selection of surfactants and polymers in practical application.     

Structure and transport properties of polymer grafted nanoparticles

We perform molecular dynamics simulations on a bead-spring model of pure polymer grafted nanoparticles (PGNs) and of a blend of PGNs with a polymer melt to investigate the correlation between PGN design parameters (such as particle core concentration, polymer grafting density, and polymer length) and properties, such as microstructure, particle mobility, and viscous response. Constant strain-rate simulations were carried out to calculate viscosities and a constant-stress ensemble was used to calculate yield stresses. The PGN systems are found to have less structural order, lower viscosity, and faster diffusivity with increasing length of the grafted chains for a given core concentration or grafting density. Decreasing grafting density causes depletion effects associated with the chains leading to close contacts between some particle cores. All systems were found to shear thin, with the pure PGN systems shear thinning more than the blend; also, the pure systems exhibited a clear yielding behavior that was absent in the blend. Regarding the mechanism of shear thinning at the high shear rates examined, it was found that the shear-induced decrease of Brownian stresses and increase in chain alignment, both correlate with the reduction of viscosity in the system with the latter being more dominant. A coupling between Brownian stresses and chain alignment was also observed wherein the non-equilibrium particle distribution itself promotes chain alignment in the direction of shear.     

Probing the molecular-scale lipid bilayer response to shear flow using nonequilibrium molecular dynamics

A nonequilibrium molecular dynamics simulation of the response of dimyristoylphosphatidylcholine (DMPC) bilayers to a solvent shear flow is presented. Application of shear flow to planar, stationary DMPC bilayers results in a redistribution of the membrane density profile along the bilayer normal due to the alignment of the lipids in the direction of flow and an increase in average lipid chain length. An increase in the intermolecular and intramolecular order of the lipids in response to the shear flow is also observed. This study provides groundwork for understanding the mechanism of the full response of lipid bilayers to externally imposed solvent shear flows, beginning with the response in the absence of collective lipid motions such as undulations and bilayer flow.     

A Langevin dynamics simulation study of the tribology of polymer loop brushes

The tribology of surfaces modified with doubly bound polymer chains (loops) has been investigated in good solvent conditions using Langevin dynamics simulations. The density profiles, brush interpenetration, chain inclination, normal forces, and shear forces for two flat substrates modified by doubly bound bead-necklace polymers and equivalent singly bound polymers (twice as many polymer chains of 12 the molecular weight of the loop chains) were determined and compared as a function of surface separation, grafting density, and shear velocity. The doubly bound polymer layers showed less interpenetration with decreasing separation than the equivalent singly bound layers. Surprisingly, this difference in interpenetration between doubly bound polymer and singly bound polymer did not result in decreased friction at high shear velocity possibly due to the decreased ability of the doubly bound chains to deform in response to the applied shear. However, at lower shear velocity, where deformation of the chains in the flow direction is less pronounced and the difference in interpenetration is greater between the doubly bound and singly bound chains, some reduction in friction was observed.     

Some aspects of theory and simulations of polymers in steady flows

We present a study of various properties of bead spring chains in steady flows. The Langevin equation of the normal modes of the chain is solved by Fourier transformation. From the resulting power spectrum, the autocorrelation functions of all configuration-dependent quantities can be calculated. In equilibrium, the influence of the bead masses on the short-time dynamics is discussed. The influence of different flow fields (shear, elongational and Kramers potential flow) on the mean-square chain dimension is calculated. A comparison with results obtained from non-equilibrium molecular dynamics and Monte Carlo calculations is made. Finally, the influence of shear flow on the configurational and rheological properties of cyclic polymers and on the excluded volume behavior of linear chains is examined.     

Effect of Wax Composition and Shear Force on Wax Aggregation Behavior in Crude Oil: A Molecular Dynamics Simulation Study

To explore the influence of different wax components and the shear effect exerted by the pump and pipe wall in the process of crude oil pipeline transportation on the microbehavior of wax aggregation in crude oil at low temperatures, molecular dynamics models of binary and multivariate systems of crude oil with different wax components are established in this paper. The simulation results are compared with the existing experimental results and the NIST database to verify the rationality and accuracy of the models. By using the established binary model to simulate four crude oil systems containing different wax components, it can be found that the longer the wax molecular chain, the more easily the wax molecules aggregate. The influence of temperature on the aggregation process of wax molecules with different chain lengths is also studied. The lower the temperature, the greater the difference in wax molecular aggregation degree caused by the difference in molecular chain length. Nonequilibrium molecular dynamics is used to simulate the shear process of a multivariate system of crude oil, and the micromechanisms of the shear effect on the aggregation process of wax molecules are studied. Shearing can destroy the stable structure of crude oil, resulting in the orientation and conformational transformation of wax molecules, and obtaining the region of wax molecules sensitive to temperature and shear effects, the temperatures of which are below the wax precipitation point and the shear rate of which is lower than the maximum shear rate to prevent the molecular structure from being destroyed. At the same time, the sensitivity of wax components with different chain lengths to the shear effect is studied. The research results provide theoretical guidance for ensuring the safe and economic operation of waxy crude oil production.     

The Chain Stretching Effect on the Morphology and Rheological Properties of Phase‐Separating Polymer Blends Subjected to Simple Shear Flow

Summary: The spinodal decomposition of a binary mixture subjected to simple shear flow is investigated in the framework of the modified time‐dependent Ginzburg‐Landau (TDGL) equation with an external velocity term. The domain growth and related rheological properties of a binary mixture under shear flow are simulated in three dimensions by means of the cell dynamics scheme (CDS). The simulation results show that the domain growth is anistropic and depends on the terminal relaxation time of the polymer chain. It is found that lamellae‐like domains with the normal parallel to the velocity gradient direction are observed when the terminal relaxation time is long enough. This result has also been confirmed by carefully checking the scattering functions in different incident light directions and the evolution of the domain size in different directions. In addition, when the chain stretching effect is strong, the transients of the excess shear viscosity are much higher than the case without the chain stretching effect. The terminal relaxation time of the chain also has an important effect on the first and second normal stress differences.

The time evolution of the morphology for the case with strong chain stretching effect.     

Conformational and Dynamical Evolution of Block Copolymers in Shear Flow

Conformation and dynamical evolution of block copolymers in shear flow is an important topic in polymer physics that underscores the forming process of various materials. We explored deformation and dynamics of copolymers composed of rigid or flexible blocks in simple shear flow by employing multiparticle collision dynamics integrated with molecular dynamics simulations. We found that compared with the proportion between rigid and flexible blocks, the type of the central blocks plays more important role in the conformational and dynamical evolution of copolymers. That is, if the central block is a coil, the copolymer chain takes end-over-end tumbling motion, while if the central block is a rod, the copolymer chain undergoes U-shape or S-shape deformation at mid shear rate. As the shear strength increases, all copolymers behave similar to flexible polymers at high shear rate. This can be attributed to the fact that shear flow is strong enough to overcome the buckling force of the rigid blocks. These results provide a deeper understanding of the roles played by rod and coil blocks in copolymers for phase interface during forming processing.     

Structural and dynamical properties for confined polymers undergoing planar Poiseuille flow

The authors present the results from nonequilibrium molecular dynamics simulations for the structural and dynamical properties of highly confined linear polymer fluids undergoing planar Poiseuille flow. They study systems confined within pores of several atomic diameters in width and investigate the dependence of the density profiles, the mean squared radius of gyration, the mean squared end-to-end distance, streaming velocity, strain rate, shear stress, and streaming angular velocity as functions of average fluid density and chain length. Their simulation results show that, sufficiently far from the walls, the radius of gyration for molecules under shear in the middle of the pore follows the power law Rg=ANbnu, where Nb is the number of bonds and the exponent has a value of 0.5 which resembles the value for a homogeneous equilibrium fluid. Under the conditions simulated, the authors find the onset of flat velocity profiles but with very little wall slippage. These flat profiles are most likely due to the restricted layering of the fluid into just one or two molecular layers for narrow pore widths compared to chain length, rather than typical plug-flow conditions. The angular velocity is shown to be proportional to half the strain rate in the pore interior when the chain length is sufficiently small compared to the pore width, consistent with the behavior for homogeneous fluids in the linear regime.     

Stochastic dynamics simulation of grafted polymer brushes under shear deformation

We present results of computer simulations of polymer brushes (layers of polymer chains attached at one end onto an impermeable planar surface) under shear deformation at constant shear rate. As the first stage of calculations the behavior of a single brush was studied. The monomer density profile, the distribution of the chain ends, the positions and orientations of different monomers along the chain were calculated. Dimensions of the polymer chains as functions of the shear rate were obtained for different grafting densities. An increase in the brush thickness over the grafting plane with an increase in the shear rate as predicted by the theory of Barrat was observed. However, the magnitude of the effect appears to be small. We explain this by finite extensibility of the grafted chains.     

高分子在受限稀溶液中的结构和动力学性质

利用平衡态及非平衡态耗散粒子动力学模拟方法, 分别研究了平衡态和流场作用下受限高分子在稀溶液中的链结构和动力学. 采用没有滑移和密度涨落的边界条件模拟纳米管道环境, 进而研究了高分子回转半径和扩散系数对受限强度及高分子与溶剂间相互作用的依赖关系. 在非平衡态模拟中, 分别考虑了Poiseuille和Couette两种流场. 研究发现, 在这两种流场作用下, 高分子都随着溶剂与高分子排斥作用的降低而向管道中心迁移. 在强流场下, 在Poiseuille流场中高分子密度呈现出双峰分布, 而在Couette流场中则呈现为单峰分布.     

Chitosan and its derivatives in extensional and shear flows

Allylchitosan and propylchitosan with different degrees of substitution were prepared on the basis of chitosan from shrimp chitin. The dynamics of semidilute electrolytic polymer solutions of chitosan derivatives in acetic acid was studied by measuring birefringence in extensional and shear flows and by means of viscometry. The optical shear coefficient and critical velocity gradients corresponding to the loss of stability of the macromolecular coil in extensional flow were found. The chain relaxation times depending on the polyelectrolyte concentration and ionic strength of solution were determined.     

均聚物链吸附在胶束表面上的非平衡态耗散粒子动力学

利用平衡态耗散粒子动力学（Dissipative Particle Dynamics, DPD）方法,研究了嵌段共聚物自组装形成胶束以及均聚物链在其表面的吸附。进一步采用非平衡态耗散粒子动力学（Non-Equilibrium Dissipative Particle Dynamics, NEDPD）方法,通过加入剪切,研究了均聚物的脱落和胶束的分裂过程。模拟结果揭示了脱落时间随剪切速率的变化规律。     

Divergent Shear Thinning and Shear Thickening Behavior of Supramolecular Polymer Networks in Semidilute Entangled Polymer Solutions

The steady shear behavior of metallo-supramolecular polymer networks formed by bis-Pd(II) cross-linkers and semidilute entangled solutions of poly(4-vinylpyridine) (PVP) in dimethyl sulfoxide (DMSO) or N,N-dimethyl formamide (DMF) is reported. The steady shear behavior of the networks depends on the dissociation rate and association rate of the cross-linkers, the concentration of cross-linkers, and the concentration of the polymer solution. The divergent steady shear behavior-shear thinning versus shear thickening-of samples with identical structure but different cross-linker dynamics (J. Phys. Chem. Lett. 2010, 1, 1683-1686) is further explored in this paper. The divergent steady shear behavior for networks with different cross-linkers is connected to a competition between different time scales: the average time that a cross-linker remains open (τ(1)) and the local relaxation time of a segment of polymer chain (τ(segment)). When τ(1) is larger than τ(segment), shear thickening is observed. When τ(1) is smaller than τ(segment), only shear thinning is observed.     

Static and dynamic properties of the interface between a polymer brush and a melt of identical chains

Molecular-dynamics simulations of a short-chain polymer melt between two brush-covered surfaces under shear have been performed. The end-grafted polymers which constitute the brush have the same chemical properties as the free chains in the melt and provide a soft deformable substrate. Polymer chains are described by a coarse-grained bead-spring model, which includes excluded volume and backbone connectivity of the chains. The grafting density of the brush layer offers a way of controlling the behavior of the surface without altering the molecular interactions. We perform equilibrium and nonequilibrium molecular-dynamics simulations at constant temperature and volume using the dissipative particle dynamics thermostat. The equilibrium density profiles and the behavior under shear are studied as well as the interdigitation of the melt into the brush, the orientation on different length scales (bond vectors, radius of gyration, and end-to-end vector) of free and grafted chains, and velocity profiles. The obtained boundary conditions and slip length show a rich behavior as a function of grafting density and shear velocity.     

Transport properties of Mie(14,7) fluids: molecular dynamics simulation and theory

An extensive computer simulation study is presented for the self-diffusion coefficient, the shear viscosity, and the thermal conductivity of Mie(14,7) fluids. The time-correlation function formalism of Green-Kubo is utilized in conjunction with molecular dynamics (MD) simulations. In addition to molecular simulations, the results of a recent study [A. Eskandari Nasrabad, J. Chem. Phys. 128, 154514 (2008)] for the mean free volume are applied to calculate the self-diffusion coefficients within a free volume theory framework. A detailed comparison between the MD simulation and free volume theory results for the diffusion coefficient is given. The density fluctuation theory of shear viscosity is used to compute the shear viscosity and the results are compared to those from MD simulations. The density and temperature dependences of different time-correlation functions and transport coefficients are studied and discussed.     

Interfacial friction between semiflexible polymers and crystalline surfaces

The results obtained from molecular dynamics simulations of the friction at an interface between polymer melts and weakly attractive crystalline surfaces are reported. We consider a coarse-grained bead-spring model of linear chains with adjustable intrinsic stiffness. The structure and relaxation dynamics of polymer chains near interfaces are quantified by the radius of gyration and decay of the time autocorrelation function of the first normal mode. We found that the friction coefficient at small slip velocities exhibits a distinct maximum which appears due to shear-induced alignment of semiflexible chain segments in contact with solid walls. At large slip velocities, the friction coefficient is independent of the chain stiffness. The data for the friction coefficient and shear viscosity are used to elucidate main trends in the nonlinear shear rate dependence of the slip length. The influence of chain stiffness on the relationship between the friction coefficient and the structure factor in the first fluid layer is discussed.     

Rheological investigation of the influence of molecular structure on natural and accelerated UV degradation of linear low density polyethylene

Metallocene and Ziegler-Natta (ZN) linear low density polyethylenes (LLDPEs) of different branch types and contents as well as linear high density polyethylene (HDPE) were exposed to natural and accelerated weather conditions. The degree of UV degradation of exposed samples was measured by rheological techniques and results were compared with unexposed polymers. Dynamic shear measurements were performed in an ARES rheometer in the linear viscoelastic range. The degree of enhancement or reduction in viscosity and elasticity was used as a measure of the degree of cross-linking or chain scission, respectively. The degradation results of LLDPE suggest that both cross-linking and chain scission are taking place. Chain scission dominated the degradation at high levels of short chain branching (SCB) and long exposure times. The degradation mechanism of m-LLDPE and ZN-LLDPE is similar; however, m-LLDPE showed a higher degradation rate than ZN-LLDPE of similar M_w and average SCB. ZN-LLDPE was found to be more stable than a similar m-LLDPE. Comonomer type had little influence on degradation. Dynamic shear rheology was very useful in revealing the influence of different molecular parameters and it exposed the degradation mechanism.     

Relaxation dynamics of stretched polymer chains traced by shear induced Shish‐Kebab structure

We have developed a morphologic method to investigate the relaxation processing of the stretched polymer chains in melts, in which an atomic force microscope probe was used to shear the surface of an isotactic polypropylene melt to obtain the isolated shish‐kebab structure. We present the results of the time dependence of length of the isolated shish‐kebab structure and the stress dependence of the kebab density along the direction of shish in this paper. Our results demonstrate that the shear‐oriented polymer melts show the relaxation dynamics of worm‐like chain where the length deficit of the isolated shish‐kebab structure is scaled with the relaxation time as a power of 1/3. The melting behavior of shish‐kebab structure was also investigated. © 2013 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2013, 51, 907–914