蠕虫状链模型下梳状共聚物分子的均方回转半径的理论分析及其在聚羧酸系减水剂分子中的应用

利用蠕虫状链模型对梳状共聚物分子的均方回转半径进行了理论分析,建立了回转半径与主链轮廓长度、主链持久长度、侧链轮廓长度、侧链持久长度、侧链数目以及侧链沿主链的分布情况(均匀分布和梯度分布)之间的定量数学关系.在此基础上,以被称为"第一代聚羧酸系高性能减水剂"(以下简称为MPEG-type PCE)的甲基丙烯酸(MAA)/烯酸甲酯(MAA-MPEG)梳状共聚物分子为研究对象,结合实验数据,对其聚电解质主链的持久长度进行了分析,并考察了主侧链长度、刚柔性、侧链分布、接枝密度等分子结构参数对PCE回转半径的影响,最后对模型的局限性作简要说明.梳状共聚物分子的蠕虫状链模型物理图像简洁,参数意义明确,应用于PCE分子体系时较之前所报道的柔性链模型要更为合理,能够为分析PCE的分子结构与溶液构象和吸附构象之间的关系提供更科学的视角.     

聚羧酸系高性能减水剂分子在稀溶液中尺寸及其在体积排除色谱表征中表观分子量的模型研究

针对被称为"第一代聚羧酸高性能减水剂"(以下简称为MPEG-type PCE)的甲基丙烯酸(MAA)/烯酸甲酯(MAA-MPEG)梳状共聚物分子,从高分子物理基础理论出发,构建等效自由连接链模型,结合前人的理论结果和实验数据,得到了MPEG-type PCE分子的回转半径、流体力学半径及其相应的支化参数的数学表达式.在此基础上,报道了以下三方面的工作:首先,将计算结果与文献中的实验结果进行比较,检验模型的合理性;其次,利用所建立的数学模型考察主链分子量、侧链分子量和侧链接枝密度对PCE分子的回转半径和流体力学半径的影响;最后,结合近年来发展的体积排除色谱分离理论,对PCE分子的真实分子量与其常规体积排除色谱"表观分子量"(又被称为GPC分子量)两者之间的差异进行了分析.本文所提出的计算模型和数学表达式没有不确定的指前因子,可用来估算MPEG-type PCE分子在稀水溶液中的尺寸以及根据其GPC分子量估算真实分子量.     

Chemomechanics of surface stresses induced by DNA hybridization

When biomolecular reactions occur on one surface of a microcantilever beam, changes in intermolecular forces create surface stresses that bend the cantilever. While this phenomenon has been exploited to create label-free biosensors and biomolecular actuators, the mechanisms through which chemical free energy is transduced to mechanical work in such hybrid systems are not fully understood. To gain insight into these mechanisms, we use DNA hybridization as a model reaction system. We first show that the surface grafting density of single-stranded probe DNA (sspDNA) on a surface is strongly correlated to its radius of gyration in solution, which in turn depends on its persistence length and the DNA chain length. Since the persistence length depends on ionic strength, the grafting density of sspDNA can be controlled by changing a solution's ionic strength. The surface stresses produced by the reaction of complementary single-stranded target DNA (sstDNA) to sspDNA depend on the length of DNA, the grafting density, and the hybridization efficiency. We, however, observe a remarkable trend: regardless of the length and grafting density of sspDNA, the surface stress follows an exponential scaling relation with the density of hybridized sspDNA.     

The variations of the configurational and solvency properties of low molecular weight sodium polyacrylate with ionic strength

The configurational and solvency properties of low molecular weight sodium polyacrylate have been determined for a wide range of ionic strength solutions, from intrinsic viscosity data in the polymer literature.The variations of the polymer properties with ionic strength (I) are described very well by simple mathematical expressions. Thus, a linear relationship was found between the solvency parameter and 1/I ^(1/2), while the variations of the expansion factor and the radius of gyration with 1/I ^(1/2) were described by second order polynomials.LowI solutions (i.e. < 0.01) have a high solvency for sodium polyacrylate. In such solutions the polymer is in a highly expanded configuration. Thus, the radius of gyration of a typical, low molecular weight (ca. 5000 g mol^–1) sodium polyacrylate approaches the limiting value of ca. 4.5 nm atI<0.01.Conversely, high ionic strength solutions (i.e. >0.10) have a low solvency for sodium polyacrylate. In such solutions the polymer is in a virtually unexpanded configuration. Thus, the radius of gyration of a typical, low molecular weight sodium polyacrylate approaches the limiting value of ca. 2.0 nm atI>0.10.     

Structure of flexible and semiflexible polyelectrolyte chains in confined spaces of slit micro/nanochannels

Understanding the behavior of a polyelectrolyte in confined spaces has direct relevance in design and manipulation of microfluidic devices, as well as transport in living organisms. In this paper, a coarse-grained model of anionic semiflexible polyelectrolyte is applied, and its structure and dynamics are fully examined with Brownian dynamics (BD) simulations both in bulk solution and under confinement between two negatively charged parallel plates. The modeling is based on the nonlinear bead-spring discretization of a continuous chain with additional long-range electrostatic, Lennard-Jones, and hydrodynamic interactions between pairs of beads. The authors also consider the steric and electrostatic interactions between the bead and the confining wall. Relevant model parameters are determined from experimental rheology data on the anionic polysaccharide xanthan reported previously. For comparison, both flexible and semiflexible models are developed accompanying zero and finite intrinsic persistence lengths, respectively. The conformational changes of the polyelectrolyte chain induced by confinements and their dependence on the screening effect of the electrolyte solution are faithfully characterized with BD simulations. Depending on the intrinsic rigidity and the medium ionic strength, the polyelectrolyte can be classified as flexible, semiflexible, or rigid. Confined flexible and semiflexible chains exhibit a nonmonotonic variation in size, as measured by the radius of gyration and end-to-end distance, with changing slit width. For the semiflexible chain, this is coupled to the variations in long-range bond vector correlation. The rigid chain, realized at low ionic strength, does not have minima in size but exhibits a sigmoidal transition. The size of confined semiflexible and rigid polyelectrolytes can be well described by the wormlike chain model once the electrostatic effects are taken into account by the persistence length measured at long length scale.     

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.     

The effect of inhomogeneous stickiness on polymer aggregation

The aggregation of polymers is important in the formation of marine aggregates and the vertical transport of material in the ocean. A polymer may be inhomogeneous along its length, with associating groups at some points along its length where bonds are more likely to form. In this paper we investigate the effects of inhomogeneous 'stickiness' along the polymer length. We describe the results of three-dimensional off-lattice simulations of polymer-polymer aggregation for four different types of polymer: polymers which are sticky along their entire length, polymers which are sticky at the ends only and two types of polymer which are slightly sticky along their entire length. We examine the mean radius of gyration and the fractal dimension of the resulting aggregates and the dynamics of aggregation. The slightly sticky polymers and the polymers which are sticky only at the ends form aggregates with a higher fractal dimension than the polymers which are sticky along their entire length. However, the mean radius of gyration of the aggregates formed by polymers which are sticky only at the ends is significantly larger than that of the aggregates formed from slightly sticky polymers. The aggregation dynamics are also different for the polymers which are sticky only at the ends compared to the slightly sticky polymers. A single 'stickiness value' is therefore likely to be inadequate to describe a polymer. We also examine the effect of polymer rigidity; it seems that the effect of inhomogeneous stickiness is greater for almost-straight polymers than for coiled chains.     

Hydrodynamic Radius Characterization of a Vinyl-Type Polynorbornene by Size-Exclusion Chromatography with a Static and Dynamic Laser Light Scattering Detector

Both absolute molecular weight and molecular sizes (radius of gyration and hydrodynamic radius) of a vinyl-type polynorbornene eluting from size-exclusion chromatography columns were determined by combined with a static and dynamic laser light scattering detector. The hydrodynamic radius of polymer fraction eluting from size-exclusion chromatography columns was obtained from dynamic laser light scattering measurements at only a single angle of 90° by introducing a correction factor. According to the scaling relationship between molecular sizes and molecular weight and the ratio between radius of gyration and hydrodynamic radius, the vinyl-type polynorbornene took a random coil conformation in 1,2,4-trichlorobenzene at 150 °C.     

Anomalous nanoparticle diffusion in polymer solutions and melts: a mode-coupling theory study

Mode-coupling theory is employed to study diffusion of nanoparticles in polymer melts and solutions. Theoretical results are directly compared with molecular dynamics simulation data for a similar model. The theory correctly reproduces the effects of the nanoparticle size, mass, particle-polymer interaction strength, and polymer chain length on the nanoparticle diffusion coefficient. In accord with earlier experimental, simulation, and theoretical work, it is found that when the polymer radius of gyration exceeds the nanoparticle radius, the Stokes-Einstein relation underestimates the particle diffusion coefficient by as much as an order of magnitude. Within the mode-coupling theory framework, a microscopic interpretation of this phenomenon is given, whereby the total diffusion coefficient is decomposed into microscopic and hydrodynamic contributions, with the former dominant in the small particle limit, and the latter dominant in the large particle limit. This interpretation is in agreement with previous mode-coupling theory studies of anomalous diffusion of solutes in simple dense fluids.     

Elastic behavior of comb-like polymer chains

<正>Three-dimensional Monte Carlo simulations of comb-like polymer chains with various backbone lengths N_b,arm lengths N_a and arm densities m are carried out to study the elastic behavior of comb-like polymer chains.The radius of gyration,the shape factors and bond length in different cases during elastic process are calculated,and it is found that the comb-like polymer molecules with longer backbone or shorter arm are more close to linear chains.But the arm density m affects the chain conformation non-monotonously.Some thermodynamic properties are also studied.Average Helmholtz free energy and elastic force f all increase with elongation ratioλfor all chains.     

Development of knotting during the collapse transition of polymers

A dynamic Monte Carlo simulation of the collapse transition of polymer chains is presented. The chains are represented as self-avoiding walks on the simple cubic lattice with a nearest-neighbor contact potential to model the effect of solvent quality. The knot state of the chains is determined using the knot group procedure presented in the accompanying paper. The equilibrium knot spectrum and the equilibrium rms radius of gyration as functions of the chain length and the contact potential are reported. The collapse transition was studied following quenches from good-to poor-solvent conditions. Our results confirm the prediction that the newly formed globule is not yet at equilibrium, since it has not yet achieved its equilibrium knot spectrum. For our model system, the relaxation of the knot spectrum is about an order of magnitude slower than that of the radius of gyration. The collapse transition is also studied for a model in which both ends of the chain remain in good-solvent conditions. Over the time scale of these simulations, knot formation is frustrated in this inhomogeneous model, verifying that the mechanism of knotting is the tunneling of chain ends in and out of the globule.     

On the structure and the chain conformation of side-chain liquid crystal polymers

Abstract

Backbone anisotropy and the structure of the mesophases of a series of side-chain liquid crystal polymers have been studied in the bulk by neutron scattering. The backbone conformation is obtained by small-angle neutron scattering on mixtures of hydrogenous polymers with deuteriated backbones. The components of the radius of gyration parallel, R _⊥ and perpendicular, R ∥ to the magnetic field are determined as a function of temperature for both the nematic phase and the smectic phase. It is shown that the polymer backbone is deformed in both phases. When the polymer exhibits only a nematic phase, it adopts a prolate conformation, where the average backbone direction is more or less parallel to the aligned mesogenic groups. Upon transition from the smectic phase to a nematic phase, the backbone in the nematic phase assumes a slightly oblate shape. This tendency towards oblate shape is due to the smectic fluctuations which are always present in such nematic phases. The exentricity of the oblate backbone conformation in the smectic phase is always larger than in the nematic phase. This is attributed to a periodic distribution of the backbone between the mesophase layers. Then, the backbone anisotropy depends not only on the smectic structure (S_A1, S_Ad), but also on the temperature dependence of the density of aligned mesogenic groups in the layers. On the other hand, it is shown that the isotopic mixtures are no longer ideal when polymers deuteriated in the mesogenic moieties are mixed with the corresponding hydrogenous polymers.     

DIMENSIONS OF ATACTIC POLY（α-METHYLSTYRENE） CHAINS WITH SIDE GROUPS

An improved configurational-confomational statistical method is developed and the mean-square radius of gyration for atactic poly(α-methylstyrene)(PαMS)chains is studied,in which the effect of large side groups is considered. The deduced formulas,based on the rotational isomer state theory,are used to investigate the configuration-dependent properties of the atactic polymer chain,and the statistical correlation of the unperturbed polymer chain dimension and structure parameters are calculated.For the fraction of meso dyads w_m=0.4,the dependence of the radius of gyration R_g and the intrinsic viscosity[η]on the molecule mass M are R_g=2.63×10~(-2) M~(0.50) nm and[η]=7.36×10~(-2) M~(0.497),respectively, which are in agreement with the previous experimental data for the PαMS samples.A small hump is detected in the curve of the characteristic ratio of the unperturbed mean-square radius of gyration versus the chain length for short PαMS chains.The R_g increases linearly with the temperature T,and the effects of the chain length and the tacticity on the temperature coefficient are remarkable.These are quite different from the results for PαMS chains not considering side groups or for the monosubstituted polystyrene chain.     

Modelling of θ-conditions for bisphenol—A polycarbonate single chains

Assessing conformational dimensions of macromolecules is a topic of long-standing interest. Because of its simplicity, it is attractive to investigate the chain properties in θ-conditions. Under these special conditions, the effects of excluded volume of the segments of the polymer chain vanish. The molecular chain is only subject to local constraints resulting from the bond structure and the hindrance to rotations about bonds. To model θ-conditions a contour length dependent cutoff is introduced ensuring that only nonbonded interactions of atoms of neighbouring monomeric units are taken into account for energy calculations. Using this energy model we will show that it is possible to model θ-conditions of a single bisphenol-A polycarbonate (BPA-PC) chain in vacuum by two different methods: (i) (Pseudo-) Langevin dynamics simulations and (ii) regular reassignment of randomly generated atom velocities during a molecular dynamics simulation. Both methods can be used to avoid oscillative dynamic behaviour of the chain. Calculations of the end-to-end vector and the radius of gyration of the equilibrium ensembles derived from simulations at different temperatures show good agreement with experimental data. Thus our model techniques are well suited to simulate θ-conditions with small computational expense.     

Noncontinuum effects on the mobility of nanoparticles in unentangled polymer solutions

The results for the diffusivity of nanoparticles in unentangled semidilute polymer solutions obtained using coarse‐grained simulations are presented. The results indicate that for particle sizes smaller than the polymer radius of gyration, the nanoparticle diffusivities deviate from Stokes–Einstein predictions and depend explicitly on the polymer radius of gyration and the polymer solution correlation lengths. Scaling ideas proposed are invoked for rationalizing such noncontinuum effects and demonstrate that the simulation results could be collapsed onto a single universal function of the depletion thickness, the polymer radius of gyration, and the particle radius. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2016 , 54, 2145–2150.     

Influence of Rigid Side Chain Attraction on Stiffness and Conformations of Comb Copolymer Brushes Strongly Adsorbed on a Flat Surface

The influence of side‐chain attraction on the conformational properties of two‐dimensional polymer brushes with rigid side chains is investigated using Monte Carlo simulations. Using a rigid backbone, a characteristic interaction strength is determined by investigating the critical interaction energy for the collapse of the side chains onto the backbone. For a flexible backbone, the persistence length of the backbone is found to decrease with increasing attraction, irrespective of whether side‐chain flipping is allowed or not. This result is in good agreement with the theoretical modeling presented before. If side‐chain flipping is allowed, the attraction between the side chains leads to aggregation of successive side chains at one side of the backbone resulting in a characteristic local spiraling of the backbone.     

Electrostatic interactions in dissipative particle dynamics using the Ewald sums

The electrostatic interactions in dissipative particle dynamics (DPD) simulations are calculated using the standard Ewald [Ann. Phys. 64, 253 (1921)] sum method. Charge distributions on DPD particles are included to prevent artificial ionic pair formation. This proposal is an alternative method to that introduced recently by Groot [J. Chem. Phys. 118, 11265 (2003)] where the electrostatic field was solved locally on a lattice. The Ewald method is applied to study a bulk electrolyte and polyelectrolyte-surfactant solutions. The structure of the fluid is analyzed through the radial distribution function between charged particles. The results are in good agreement with those reported by Groot for the same systems. We also calculated the radius of gyration of a polyelectrolyte in salt solution as a function of solution pH and degree of ionization of the chain. The radius of gyration increases with the net charge of the polymer in agreement with the trend found in static light scattering experiments of polystyrene sulfonate solutions.     

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.     

Influence of Solvent Size on Structural Properties of Bottle-brush Polymers

We proposed a modified dissipative particle dynamics simulation model by which one can effectively avoid the bond-crossing problem, and investigated the effect of solvent size on the structural properties of bottle-brush polymers in dilute solution on the basis of this model. It was found that with the increase of solvent size, the radius of the gyration of the bottle-brush polymer decreases considerably in the athermal solvent but increases in the selective solvent favoring the backbone, respectively.