Simulation and model development for the equation of state of heteronuclear non‐additive copolymer chains

Molecular dynamics simulations of hard alternating copolymer chains composed of size asymmetric nonadditive segments were performed. Different degrees of polymerization, densities, size ratios and nonadditivities were used. The equation of state for these copolymers was investigated and models based on the first order thermodynamic perturbation theory (TPT1) and the polymeric analog of the Percus‐Yevick approximation (PPY) were developed to predict the compressibility factor of the copolymers. The models predicted the compressibilities of the mixtures accurately at small size ratios, low degrees of polymerization and higher densities. The TPT1 model was generally more accurate in predicting the compressibility factor than the PPY model.     

Molecular-dynamic simulation of nanosized droplet spreading over a continual solid surface

Molecular dynamics methods are employed for a comparative study of the spread of nanosized simple fluid droplets and nanodroplets composed of model chains that reflect molecules of normal alkanes over a solid substrate. The effects of the droplet size, the energy of the molecule-substrate interactions, temperature, and degree of chain rigidity on the spreading kinetics are investigated. The substrate is described in terms of the continual Lennard-Jones model.     

Conformation and diffusion behavior of ring polymers in solution: a comparison between molecular dynamics, multiparticle collision dynamics, and lattice Boltzmann simulations

We have studied the effect of chain topology on the structural properties and diffusion of polymers in a dilute solution in a good solvent. Specifically, we have used three different simulation techniques to compare the chain size and diffusion coefficient of linear and ring polymers in solution. The polymer chain is modeled using a bead-spring representation. The solvent is modeled using three different techniques: molecular dynamics (MD) simulations with a particulate solvent in which hydrodynamic interactions are accounted through the intermolecular interactions, multiparticle collision dynamics (MPCD) with a point particle solvent which has stochastic interactions with the polymer, and the lattice Boltzmann method in which the polymer chains are coupled to the lattice fluid through friction. Our results show that the three methods give quantitatively similar results for the effect of chain topology on the conformation and diffusion behavior of the polymer chain in a good solvent. The ratio of diffusivities of ring and linear polymers is observed to be close to that predicted by perturbation calculations based on the Kirkwood hydrodynamic theory.     

Transitions of tethered polymer chains: a simulation study with the bond fluctuation lattice model

A polymer chain tethered to a surface may be compact or extended, adsorbed or desorbed, depending on interactions with the surface and the surrounding solvent. This leads to a rich phase diagram with a variety of transitions. To investigate these transitions we have performed Monte Carlo simulations of a bond fluctuation model with Wang-Landau and umbrella sampling algorithms in a two-dimensional state space. The simulations' density-of-states results have been evaluated for interaction parameters spanning the range from good- to poor-solvent conditions and from repulsive to strongly attractive surfaces. In this work, we describe the simulation method and present results for the overall phase behavior and for some of the transitions. For adsorption in good solvent, we compare with Metropolis Monte Carlo data for the same model and find good agreement between the results. For the collapse transition, which occurs when the solvent quality changes from good to poor, we consider two situations corresponding to three-dimensional (hard surface) and two-dimensional (very attractive surface) chain conformations, respectively. For the hard surface, we compare tethered chains with free chains and find very similar behavior for both types of chains. For the very attractive surface, we find the two-dimensional chain collapse to be a two-step transition with the same sequence of transitions that is observed for three-dimensional chains: a coil-globule transition that changes the overall chain size is followed by a local rearrangement of chain segments.     

Monte Carlo simulation and self-consistent integral equation theory for polymers in quenched random media

The conformational properties and static structure of freely jointed hard-sphere chains in matrices composed of stationary hard spheres are studied using Monte Carlo simulations and integral equation theory. The simulations show that the chain size is a nonmonotonic function of the matrix density when the matrix spheres are the same size as the monomers. When the matrix spheres are of the order of the chain size the chain size decreases monotonically with increasing matrix volume fraction. The simulations are used to test the replica-symmetric polymer reference interaction site model (RSP) integral equation theory. When the simulation results for the intramolecular correlation functions are input into the theory, the agreement between theoretical predictions and simulation results for the pair-correlation functions is quantitative only at the highest fluid volume fractions and for small matrix sphere sizes. The RSP theory is also implemented in a self-consistent fashion, i.e., the intramolecular and intermolecular correlation functions are calculated self-consistently by combining a field theory with the integral equations. The theory captures qualitative trends observed in the simulations, such as the nonmonotonic dependence of the chain size on media fraction.     

Three-body contribution to the helium interaction potential

Two nonadditive three-body analytic potentials for helium were obtained: one based on three-body symmetry-adapted perturbation theory (SAPT) and the other one on supermolecular coupled-cluster theory with single, double, and noniterative triple excitations [CCSD(T)]. Large basis sets were used, up to the quintuple-zeta doubly augmented size. The fitting functions contain an exponentially decaying component describing the short-range interactions and damped inverse powers expansions for the third- and fourth-order dispersion contributions. The SAPT and CCSD(T) potentials are very close to each other. The largest uncertainty of the potentials comes from the truncation of the level of theory and can be estimated to be about 10 mK or 10% at trimer's minimum configuration. The relative uncertainties for other configurations are also expected to be about 10% except for regions where the nonadditive contribution crosses zero. Such uncertainties are of the same order of magnitude as the current uncertainties of the two-body part of the potential.     

Interchain coupling effects on dynamics of photoexcitations in conjugated polymers

Within an extended Su-Schrieffer-Heeger model including interchain interactions and the extended Hubbard model, the dynamical relaxation of photoexcitations in two coupled conjugated polymer chains is investigated by using a nonadiabatic evolution method. Initially, one of the two chains is photoexcited and the other chain is in the dimerized ground state. Due to the interchain interactions, the electron and/or the hole can be transferred from one chain to the other chain. For weak interchain coupling, the dynamical evolution of the lattice on the photoexcited chain is similar to that found in an isolate single chain case. With interchain interactions increasing, the amplitude of the distortions on the photoexcited chain decreases, and simultaneously, that on the other chain gradually increases. Until stronger interchain coupling, the deformations of the two chains have almost the same amplitude. Besides intrachain polaron-excitons and intrachain oppositely charged polaron pairs as found in single chain case, interchain polaron-excitons and interchain separated charged polaron pairs are obtained. The results show that the yield of interchain products increases and that of intrachain products decreases with interchain coupling increasing. Totally, the yield of charged polarons (including intrachain oppositely charged polaron pairs and interchain oppositely charged polaron pairs) is about 25%, in good agreement with results from experiments.     

Molecular dynamics simulations of capillary rise experiments in nanotubes coated with polymer brushes

The capillary filling of a nanotube coated with a polymer brush is studied by molecular dynamics simulations of a coarse-grained model, assuming various conditions for the fluid-wall and fluid-brush interactions. Whereas the fluid is modeled by simple point particles interacting with Lennard-Jones forces, the (end-grafted, fully flexible) polymers that form the brush coating are described by a standard bead-spring model. Our experiments reveal that capillary filling is observed even for walls that would not be wetted by the fluid, provided the polymer brush coating itself wets. Generally, it is found that the capillary rise always proceeds through a t1/2 law with time t while the underlying molecular mechanism differs for wettable and nonwettable walls. For wettable walls, fluid imbibition is compatible with the Lucas-Washburn mechanism whereby the total influx of matter drops steadily with growing chain length N and the meniscus speed goes through a minimum at intermediate chain lengths. Moreover, because of flow, the polymer brush reorganizes its structure by forming a dense plug of chain segments under the meniscus that follows the meniscus in its motion. When the tube wall does not wet, one observes no meniscus formation for short chains although the fluid seeps through the wet brush. For a brush coating with longer chains, axial segregation between the brush segments and the fluid occurs by a kind of diffusive spreading, reminiscent of invasion percolation transport in a random medium, leading to the formation of a moving meniscus. For even longer chains that reach the tube axis, the rise of a meniscus with vanishing curvature-like imbibition in a porous medium is observed to take place.     

Motion of a branched polymer chain in confinement: a Monte Carlo study

The aim of the study was a theoretical investigation of the polymer molecules located between two parallel and impenetrable surfaces which were also attractive for polymer segments. The chains were constructed of identical segments and were restricted to knots of a simple cubic lattice. Since the chains were at good solvent conditions the only interactions between the segments of the chain were the excluded volume. The properties of the model chains were determined by means of Monte Carlo simulations with a sampling algorithm based on the chain's local changes of conformation. The differences and similarities in the structure for different adsorption regimes and the size of the slit were shown and discussed. It was observed that at certain conditions the polymer chain was adsorbed at one of the confining surfaces, and then after a certain period of time it detached from this surface and approached the opposite wall; this switch was repeated many times. The influence of the strength of the adsorption, the size of the slit, and the chain's length on the frequency of these jumps were determined. The mechanism of the chain's motion during the switch was also shown.     

Dynamics of self-assembled chaining in magnetorheological fluids

The aggregation dynamics of paramagnetic spherical particles embedded in a viscous fluid is investigated via numerical simulations using a fully coupled three-dimensional model. Particles experience simultaneously Brownian motion, dipolar magnetic attraction, and multibody hydrodynamic interactions. When the dipole strength characterizing the ratio of magnetic attraction to random diffusion exceeds a critical value, particles join together forming supraparticle structures. As time evolves, particle/chain and chain/ chain interactions lead to a continuous increase of the cluster size. The mean length of particle chains has a power-law dependence with respect to time, as predicted by the theory of diffusion-limited aggregation. Both the exponent and the characteristic time scale agree very well with the experimental results of Promislow et al.     

THE EFFECT OF POLYMER ON THE CLOUD POINT OF A SERIES OF TRIBLOCK NONIONIC SURFACTANTS IN AQUEOUS SOLUTION

A series of triblock nonionic surfactants with different Propylene oxide and ethylene oxide chain lengths were synthesized. The triblock nonionic surfactants and poly(ethylene glycols) with different molecular weight were used, to find the effects of polymer chain length and size of the micelles on the cloud point of the surfactants. Two possible models are considered on the basis of cloud point changes of the solutions, to describe the polymer- surfactant interactions. One model suggests that flocculation depletion for the polymer chains exist between two regular micelles. This provides the driving force for the neighboring micelles to approach each other and destabilize the colloidal system. The flocculation effect is more important for polymers with a shorter chain block the approach of the micelles, since there is no typical polymer-surfactant association formed but just simple small molecule associations in which the steric and solvation effects of the polymer chains make the inter-micelles interactions repulsive. The other model considers that intra-chain micelles of polysoap are formed among the surfactant monomers and long polymer chains. The bridging attraction between two intra-chain micelles in such structures can enhance the collisions among the micelles, due to the exchange of amphiphilic monomers among the neighboring micelles.     

The Structure of Star-Branched Chains in a Confined Space

Summary. We studied the properties of a simplified model of star-branched polymers confined in a slit formed by two parallel and impenetrable surfaces. The chains were built of identical united atoms (segments) whose positions were restricted to vertices of a simple cubic lattice. The polymer excluded volume and polymer segment-surface contact interactions were also introduced into the model. The properties of the model chains were determined by means of Monte Carlo simulations with a Metropolis-type sampling algorithm based on local changes of chain’s conformation. The structure of star-branched chains was investigated and the influence of the confinement and the temperature on the chain dimensions and structure was studied. It was shown that for chains in the adsorbing slits their sizes do not exhibit a universal behavior contrary to confined athermal polymers. The polymers in narrow slits at higher temperatures still exhibited features of a three-dimensional chain. It was also shown that chains in small slits and at low temperatures were fully adsorbed at one of the surfaces but could also switch the surface rapidly.     

Helical Folding of Conjugated Oligo(phenyleneethynylene): Chain‐Length Dependence,Solvent Effects,and Intermolecular Assembly

As a representative folding system that features a conjugated backbone, a series of monodispersed (o‐phenyleneethynylene)‐alt‐(p‐phenyleneethynylene) (PE) oligomers of varied chain length and different side chains were studied. Molecules with the same backbone but different side‐chain structures were shown to exhibit similar helical conformations in respectively suitable solvents. Specifically, oligomers with dodecyloxy side chains folded into the helical structure in apolar aliphatic solvents, whereas an analogous oligomer with tri(ethylene glycol) (Tg) side chains adopted the same conformation in polar solvents. The fact that the oligomers with the same backbone manifested a similar folded conformation independent of side chains and the nature of the solvent confirmed the concept that the driving force for folding was the intramolecular aromatic stacking and solvophobic interactions. Although all were capable of inducing folding, different solvents were shown to bestow slightly varied folding stability. The chain‐length dependence study revealed a nonlinear correlation between the folding stability with backbone chain length. A critical size of approximately 10 PE units was identified for the system, beyond which folding occurred. This observation corroborated the helical nature of the folded structure. Remarkably, based on the absorption and emission spectra, the effective conjugation length of the system extended more effectively under the folded state than under random conformations. Moreover, as evidenced by the optical spectra and dynamic light‐scattering studies, intermolecular association took place among the helical oligomers with Tg side chains in aqueous solution. The demonstrated ability of such a conjugated foldamer in self‐assembling into hierarchical supramolecular structures promises application potential for the system.     

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.     

Dynamics of polypropylene chains in their binary blends of different stereochemical sequences studied by Monte Carlo simulations

The conformational and dynamic properties of polypropylene （PP） for both pure melts and blends with different chain tacticity were investigated by Monte Carlo simulation of isotactic （iPP）, atactic （aPP） and syndiotactic （sPP） polypropylenes. The simulation of coarse-grained PP models was performed on a high coordination lattice incorporating short- and long-range intramolecular interactions from the rotational isomeric state （RIS） model and Lennard-Jones （LJ） potential function of propane pairs, respectively. The dynamics of chains in binary PP/PP mixture were investigated with the composition of C150H302 with different chain taciticity. The diffusion rates of PP with different stereochemistry are generally in the order as： iPP 〉 aPP 〉〉 sPP. For PP/PP blends with 50：50 wt% binary mixtures, immiscibility was observed when sPP was introduced into the mixtures. The diffusion rate of iPP and aPP became slower after mixing, while sPP diffuses significantly faster in the binary mixtures. The mobility of PP chains depends on both intramolecular （molecular size and chain stiffness） and intermolecular （chain packing） interactions. The effect of intramolecular contribution is greater than that of intermolecular contribution for iPP and aPP chains in binary mixtures. For sPP chain, intermolecular interaction has greater influence on the dynamics than intramolecular contribution.     

类蛋白质分子扩散系数的模拟计算

蛋白质分子是一类比较特殊的高分子，一方面它的性质与氨基酸序列密切相关，另一方面其二级结构（如α螺旋以及β折叠等）极大地影响着它的性质．另外它的内部结构又非常紧密，与一般的高分子链完全不同．同时随着计算机技术的发展，计算机模拟日益作为一种有效的手段广泛应用于高分子科学的研究．但是由于蛋白质分子比一般的高分子结构更加复杂，如氨基酸之间有复杂的互相作用等，因此对蛋白质性质的研究往往建立在简化模型的基础上，如基于格点的HP紧密高分子模型等．虽然建立在格点模型上的蛋白质分子和真实的蛋白质分子存在着一定的差异，但基本上能体现蛋白质分子的主要特征，因此把这一类分子称为类蛋白质分子．     

The interaction between triton X series surfactants and poly (ethylene glycol) in aqueous solutions

 A series of Triton X surfactants with different ethylene oxide chain length and poly(ethylene glycols) with different molecular weight were used, to find the effects of polymer chain length and size of the micelles on the cloud point of the surfactants. Two possible models are considered on the basis of cloud point changes of the solutions, to describe the polymer–surfactant interactions. One model considers that intra-chain micelles of polysoap are formed among the surfactant monomers and long polymer chains. The bridging attraction between two intra-chain micelles in such structures can enhance the collisions among the micelles, due to the exchange of amphiphilic monomers among the neighboring micelles. The other model suggests that flocculation depletion for the polymer chains exists between two regular micelles. This provides the driving force for the neighboring micelles to approach each other and destabilize the colloidal system. The flocculation effect is more significant for polymer with a long chain. Polymers with a shorter chain block the approach of the micelles, since there is no typical polymer–surfactant association formed but just simple small molecule associations in which the steric and solvation effects of the polymer chains make the inter-micelle interactions repulsive. Received: 19 August 1997 Accepted: 11 December 1997     

Structure,conformation and dynamics of polymer chains at solid melt interfaces

We have performed molecular dynamics, and lattice Monte Carlo simulations of polymeric melts in the vicinity of solid surfaces. The structural features of the solid-melt interface were very simple. The interfacial width was comparable to the segment size. Inside this narrow interface the segment density profile was oscillatory. The density oscillations were much less pronounced than those present at solid-atomic liquid interfaces. On a scale much larger than the segment size, chain conformations were found to be identical with those of ideal chains next to a reflective barrier. In particular, the number of surface-segment contacts scaled like the square root of the molecular weight. Extensive molecular dynamics simulations showed that chain desorption times increase with molecular weight but at a rate much slower than the longest relaxation time of Rouse chains. Therefore, sufficiently long chains desorbed almost freely from the surface despite the presence of attractive surface-segment interactions. A study of chain relaxation dynamics confirmed that the Rouse modes constitute an appropriate set of normal coordinates for chains in the melt interacting with a solid surface. The effect of the surface on mode relaxation was significant. All relaxation processes of chains located within a couple of radii of gyration from the surface were slowed down considerably. This effect, however was approximately the same for fast and slow modes and independent of molecular weight for sufficiently long chains.     

Phase behavior of a fluid confined in slitlike pores with walls modified by preadsorbed chain molecules

We have studied the microscopic structure, thermodynamics of adsorption, and phase behavior of Lennard-Jones fluid in slitlike pores with walls modified due to preadsorption of chain molecules. The chain species are grafted at the walls by terminating segments. Our theoretical considerations are based on a density functional approach in the semigrand canonical ensemble. The applied constraint refers to the constant number of grafted chain molecules in the pore without restriction of the number of chains at each of the walls. We have observed capillary condensation of Lennard-Jones fluid combined with the change of the distribution of chains from nonsymmetric to symmetric with respect to the pore walls. The phase diagrams of the model are analyzed in detail, dependent on the pore width, length of chains, and grafted density.