Effects of polar group saturation on physical gelation of amphiphilic polymer solutions

Monte Carlo simulation on the basis of the comblike coarse grained nonpolar/polar (NP) model has been carried out to study the polar group saturation effect on physical gelation of amphiphilic polymer solutions. The effects of polar group saturation due to hydrogen bonding or ion bridging on the sol-gel phase diagram, microstructure of aggregates, and chain conformation of amphiphilic polymer solutions under four different solvent conditions to either the nonpolar backbone or the polar side chain in amphiphilic polymer chains have been investigated. It is found that an increase of polar group saturation results in a monotonically decreased critical concentration of gelation point, which can be qualitatively supported by the dynamic rheological measurements on pectin aqueous solutions. Furthermore, various solvent conditions to either the backbone or the side chain have significant impact on both chain conformation and microstructure of aggregates. When the solvent is repulsive to the nonpolar backbone but attractive to the polar side chain, the polymer chains are collapsed, and the gelation follows the mechanism of colloidal packing; at the other solvent conditions, the gelation follows the mechanism of random aggregation.     

Thermodynamic modelling of polymer solutions with a modified Staverman equation

A model describing the thermodynamic behaviour of polymer solutions is derived which explicitly accounts for the flexibility of the polymer chains. Based on computer simulations on various lattices it is shown that the flexibility of a polymer chain can be modelled by distinguishing different polymer conformations. Here each conformation is characterized by its corresponding number of external contact sites. The equilibrium between the different conformations is then solved for any polymer concentration and any combination of interaction energies utilizing a modified Staverman equation. The model predictions are in good agreement with the results of the computer simulations which were performed using the simple-sampling and the slithering-snake algorithm. Since the knowledge of the distribution of the conformations of a single polymer chain on an empty lattice is a prerequisite to perform the model calculations, Poisson distribution functions are fitted to the results of the corresponding computer simulations. The generalization of these distribution functions not only facilitates the use of the new model but also allows to model polymers of varying chain stiffness.     

Rhodopsin exhibits a preference for solvation by polyunsaturated docosohexaenoic acid

An all-atom molecular dynamics simulation of rhodopsin in a membrane environment has been carried out with lipid composition similar to that of the retinal membrane. The initial conformation of the protein was taken from the X-ray crystallographic structure (1F88), while those of the lipids came from a previous molecular dynamics simulation. During the course of the 12.5 ns simulation, the initially randomly placed lipids adopt an anisotropic solvation structure around the protein. The lipids, having one saturated stearic acid chain and one polyunsaturated docosohexaenoic acid chain with a zwitterionic phosphatidylcholine headgroup, arrange themselves to maximize contact between the polyunsaturated chain and the protein surface. This organization is driven by energetically favorable interactions between the transmembrance helices and the docosohexaenoyl chains that are largely of the van der Waals type. These observations are consistent with various experimental studies on rhodopsin and other G-protein coupled receptors and with the picture of extreme flexibility in polyunsaturated fatty acid chains that has arisen from recent NMR and computational work.     

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.     

How does solvent molecular size affect the microscopic structure in polymer solutions?

Monte Carlo simulation has been used to investigate the effects of linear solvent molecular size on polymer chain conformation in solutions. Increasing the solvent molecular size leads to shrinkage of the polymer chains and increase of the critical overlap concentrations. The root-mean-square radius of gyration of polymer chains (R(g)) is less sensitive to the variation of polymer concentration in solutions of larger solvent molecules. In addition, the dependency of R(g) on polymer concentration under normal solvent conditions and solvent molecular size is in good agreement with scaling laws. When the solvent molecular size approaches the ideal end-to-end distance of the polymer chain, an extra aggregation of polymer chains occurs, and the solvent becomes the so-called medium-sized solvent. When the size of solvent molecules is smaller than the medium size, the polymer chains are swollen or partially swollen. However, when the size of solvent molecules is larger than the medium size, the polymer coils shrink and segregate, enwrapped by the large solvent molecules.     

STUDY ON POLYPHENYLSILSESQUIOXANE''''S SHAPE IN DILUTE SOLUTION

Based on relations between the limiting viscosity number, the translational diffusion coefficient and the molecular weight, the shape of polyphenylsilsesquioxane (PPSQ) in dilute solution was simulated with different models. It is found that the oblate ellipsoid of revolution and free-draining random coil are two suitable models, showing that PPSQ is a polymer with rigid chains, and its shape is closer to sphere due to chain branching rather than flexibility.     

CONFORMATION AND CHAIN PACKING STRUCTURES OF POLY (ρ-PHENYLEN BENZOBISTHIAZOLE) BY MOLECULAR SIMULATION*

The conformation, the chain packing stabilization and the unit cell modeling of poly(p-phenylene benzobisthiazole) have been investigated by using molecular simulation techniquein the present work. A coupling behaviour of σ-bond rotations at either side of the pheny-lene ring or the heterocyclic ring was found surely to exceed a length of the repeat unit ofthe polymer chain. For a single chain model, the stable torsion angle of the repeat resultedat 14°. In the crystalline cell minimization, the dihedral angle along the polymer chaincould even be stabilized in various values. It therefore indicates that the intermolecularinteraction does play a considerable role for this polymer forming the conformation. Ac-cording to cohesive energies calculated for these unit cell models, the torsion angle in themost stable crystalline cell is 0°. When the chains packing together, there exist so manyenergy stable wells along the chain axis 0.35--0.36nm apart from neighbouring chains.Most of the unit cells have quite closed cohesive energies. These factors thus cause thedifficulty of forming an unique perfect crystalline structure of the polymer. The presentstudy suggested a number of reasonable unit cells, and the most stable crystalline structurefor this polymer that is monoclinic, non-primitive unit cell.     

Computer simulation of linkage of two ring chains

We performed off-lattice Monte Carlo simulations of links of two model ring chains with chain length N up to 32,768 in the theta solution or amorphous bulk state by using a random walk model (Model I), and molecular dynamics simulations of two model ring chains in solution with excluded volume interaction (Model II) to investigate topological effects on the geometry of link and ring conformation. In the case of Model I, the mean squared linking number, its distribution, and the size of two chains with fixed linking number are investigated. Our simulation results confirm the previous theoretical prediction that the mean squared linking number decays as pe(-qs(2)) with the distance of centers of chain mass s, where p and q are found to be chain length dependent and q asymptotically approaches to 0.75 as chain length increases. The linking number distribution of two chains has a universal form for long chains, but our simulation results clearly show that the distribution function deviates from the Gaussian distribution, a fact not predicted by any previous theoretical work. A scaling prediction is proposed to predict the link size, and is checked for our simulations for the Model II. The simulation results confirmed the scaling prediction of the blob picture that the link with linking number m occupies a compact volume of m blobs, and the size of the link is asymptotic to R(L) ≈ bN(ν)m(1/3-ν), where N is the chain length, and v is the Flory exponent of polymer in solutions.     

Investigating the Dynamic Viscoelasticity of Single Polymer Chains using Atomic Force Microscopy

In single‐molecule force spectroscopy (SMFS), many studies have focused on the elasticity and conformation of polymer chains, but little attention has been devoted to the dynamic properties of single polymer chains. In this study, we measured the energy dissipation and elastic properties of single polystyrene (PS) chains in toluene, methanol, and N,N‐dimethylformamide using a homemade piezo‐control and data acquisition system externally coupled to a commercial atomic force microscope (AFM), which provided more accurate information regarding the dynamic properties of the PS chains. We quantitatively measured the chain length‐dependent changes in the stiffness and viscosity of a single chain using a phenomenological model consistent with the theory of viscoelasticity for polymer chains in dilute solution. The effective viscosity of a polymer chain can be determined using the Kirkwood model, which is independent of the intrinsic viscosity of the solvent and dependent on the interaction between the polymer and solvent. The results indicated that the viscosity of a single PS chain is dominated by the interaction between the polymer and solvent. © 2019 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2019 , 57, 1736–1743     

Conformation studies on sol-gel transition in triblock copolymer solutions

The gelation of physically associating triblock copolymers in a good solvent was investigated by means of the Monte Carlo simulation and a gelation process based on the conformation transition of the copolymer that was described in detail. In our simulative system, it has been found that the gelation is closely related with chain conformations, and there exist four types of chains defined as free, dangling, loop, and bridge conformations. The copolymer chains with different conformations contribute to the formation of gel in different ways. We proposed a conformational transition model, by which we evaluated the role of these four types of chains in sol-gel transition. It was concluded that the free chains keeping the conformation transition equilibrium and the dangling conformation being the hinge of conformation transition, while the chain with loop conformation enlarges the size of the congeries and the chain with bridge conformations binds the congeries consisted of the copolymer chains. In addition, the effects of temperature and concentration on the physical gelation, the association of the copolymer congeries, and the copolymer chain conformations' distribution were discussed. Furthermore, we employed the structure factor analysis to study the association of copolymer conformations and long-range order of the simulation system and found our results are in agreement with the previous experimental conclusions.     

Conformational behavior of polymers adsorbed on nanotubes

The importance of hydrophobic interactions in determining polymer adsorption and wrapping of carbon nanotubes is still under debate. In this work, we concentrate on the effect of short-ranged weakly attractive hydrophobic interactions between polymers and nanotubes (modeled as an infinitely long and smooth cylindrical surface), neglecting all other interactions apart for chain flexibility. Using coarse-grained Monte Carlo simulation of such simplified systems, we find that uniform adsorption and wrapping of the nanotube occur for all degrees of chain flexibility for tubes with sufficiently large outer radii. However, the adsorbed conformations depend on chain stiffness, ranging from randomly adsorbed conformations of the flexible chain to perfect helical or multihelical conformations (in the case of more concentrated solutions) of the rigid chains. Adsorption appears to occur in a sequential manner, wrapping the nanotube nearly one monomer at a time from the point of contact. Once adsorbed, the chains travel on the surface of the cylinder, retaining their helical conformations for the semiflexible and rigid chains. Our findings may provide additional insight to experimentally observed ordered polymer wrapping of carbon nanotubes.     

Fracture of flexible polymer chains in dilute solution under transient extensional flow

 Applying the technique of Brownian dynamics simulation, we have studied the fracture process of flexible polymer chains when they encounter an extensional flow field of transient character. For this purpose, a mathematical model was made of an experimental device used earlier by other authors to study fracture of polystyrene in dilute solution. The polymer/solvent system studied was a very dilute solution in theta conditions. The polymer molecule was modeled as a FENE bead-spring chain, including a modification to allow for chain fracture. The simulation results showed that the fracture yield depended strongly on flow rate and on molecular weight. We have characterized the molecular-weight dependence of the critical flow rate which is necessary for fracture to occur. The distribution of the result-ing fragments is interpreted in terms of the conformation of the chains prior to fracture. Received: 17 September 1996 Accepted: 29 May 1997     

Probing the pH-dependent chain dynamics of poly(acrylate acid) in concentrated solution by using a cationic AIE fluorophore

We report a novel strategy to study the chain dynamics of poly(acrylic acid) (PAA) in a relative concentrated solution (1.0 g/L). The strategy is based on the fluorescent probe (DCTPE) with unique aggregation-induced emission (AIE) characteristics. Free DCTPE molecules are non-emissive in aqueous solution, but they become highly emissive when trapped in polymer coils. The fluorescence intensity is proportional to the efficiency of trapping DCTPE molecules in polymer coils. By correlation the change of fluorescence intensity with the variation of pH value (from 1.78 to 12.06), the PAA chain’s dynamics in the relatively concentrated solution have been elucidated into three processes. In the pH range from 12.06 to 6.0, PAA chains take an extended and non-folding conformation. Changing pH from 6.0 to 3.86, PAA chains are partially protonated and loosely packed polymer coils are formed. Further lowering the pH value of the solution (from 3.86 to 1.78), protonated segments dominate the PAA chains, and at the same time, the intermolecular hydrogen bonding takes effect, thus the polymer chains posses in the conformation of more compact coils.     

Complexation induced by weak interaction between DNA and PEO-b-P4VP below the CMC of the polymer

The complexation between circular DNA and individual chains of PEO-b-P4VP with a relatively long PEO block and a short P4VP block is highly controllable when the interaction between DNA and the polymer is weak enough. When one circular DNA chain is taken into consideration, and the polymer concentration is far below its critical micelle concentration(CMC), polymer chains are absorbed by DNA chain due to the interaction between the negatively charged DNA chain and the slightly positively charged P4VP block chains. After the adsorption/complexation, the DNA chain is converted into a nanoring(type 1). In the nanoring, the DNA chain is sufficiently wrapped by the polymer and adopts a fully stretched conformation, so that the DNA compact ratio in the nanorings is close to 1. When the polymer concentration is close to but lower than the CMC, the free polymer chains in the solution are adsorbed not only by the DNA chain but also by the polymer chains that have already been adsorbed on the DNA chain. As a result, the circular DNA chain adsorbs more polymer chains, and thus the resultant nanoring(type 2) has a larger width. In the type 2 nanoring, the DNA chain is slightly compressed; the DNA compact ratio is only about 2-3. Therefore, complexation induced by the weak interaction between DNA and PEO-b-P4VP below the CMC can produce narrow-disperse and large nanorings with a perimeter of micrometers, which are difficult to prepare by existing methods.     

Another way to view the chain conformation broadening of the line-width distribution measured in dynamic light scattering

In dynamic laser light scattering (LLS), for a given polydisperse sample, a line-width distribution G(Γ) or the translational diffusion coefficient distribution G(D) can be obtained from the measured time correlation function. For rigid colloid particles, G(Γ) can be directly related to the hydrodynamic size distribution. However, for flexible polymer chains, G(Γ) depends not only on the chain length distribution, but also on the relaxation of the chain conformation; that is, even for a monodisperse polymer sample there still exists a chain conformation distribution. If the time scale of the chain conformation relaxation is comparable to that of the translational diffusion, such as in the case of a very long polymer chain, the conformation relaxation might lead to an additional broadening in G (Γ). This "conformation broadening" has been directly observed for the first time by comparing two G(Γ) s obtained from a poly(N-isopropyl-acrylamide) solution at～25℃ and～32℃ at which the solution is ther     

Polymer-nanoparticle interfacial behavior revisited: a molecular dynamics study

By tuning the polymer-filler interaction, filler size and filler loading, we use a coarse-grained model-based molecular dynamics simulation to study the polymer-filler interfacial structural (the orientations at the bond, segment and chain length scales, chain size and conformation), dynamic and stress-strain properties. Simulated results indicate that the interfacial region is composed of partial segments of different polymer chains, which is consistent with the experimental results presented by Chen et al. (Macromolecules, 2010, 43, 1076). Moreover, it is found that the interfacial region is within one single chain size (R(g)) range, irrespective of the polymer-filler interaction and the filler size, beyond which the bulk behavior appears. In the interfacial region, the orientation and dynamic behaviors are induced by the interfacial enthalpy, while the size and conformation of polymer chains near the filler are controlled by the configurational entropy. In the case of strong polymer-filler interaction (equivalent to the hydrogen bond), the innerest adsorbed polymer segments still undergo adsorption-desorption process, the transport of chain mass center in the interfacial region exhibits away from the glassy behavior, and no plastic-like yielding point appears in the stress-strain curve, which indicates that although the mobility of interfacial polymer chains is restricted, there exist no "polymer glassy layers" surrounding the filler. In addition, it is evidenced that the filler particle prefers selectively adsorbing the long polymer chains for attractive polymer-filler interaction, validating the experimental explanation of the change of the bound rubber (BR). In short, this work provides important information for further experimental and simulation studies of polymer-nanoparticle interfacial behavior.     

Effect of the Solvent on the Conformation of Isolated MEH‐PPV Chains Intercalated Into SnS2

Photophysical processes in conjugated polymers are influenced by two competing effects: the extent of excited state delocalization along a chain, and the electronic interaction between chains. Experimentally, it is often difficult to separate the two because both are controlled by chain conformation. Here we demonstrate that it is possible to modify intra‐chain delocalization without inducing inter‐chain interactions by intercalating polymer monolayers between the sheets of an inorganic layered matrix. The red‐emitting conjugated polymer, MEH‐PPV, is confined to the interlayer space of layered SnS₂. The formation of isolated polymer monolayers between the SnS₂ sheets is confirmed by X‐ray diffraction measurements. Photoluminescence excitation (PLE) and photoluminescence (PL) spectra of the incorporated MEH‐PPV chains reveal that the morphology of the incorporated chains can be varied through the choice of solvent used for chain intercalation. Incorporation from chloroform results in more extended conformations compared to intercalation from xylene. Even highly twisted conformations can be achieved when the incorporation occurs from a methanol:chloroform mixture. The PL spectra of the MEH‐PPV incorporated SnS₂ nanocomposites using the different solvents are in good agreement with the PL spectra of the same solutions, indicating that the conformation of the polymer chains in the solutions is retained upon intercalation into the inorganic host. Therefore, intercalation of conjugated polymer chains into layered hosts enables the study of intra‐chain photophysical processes as a function of chain conformation.     

Modeling of polyethylene and poly (L-lactide) polymer blends and diblock copolymer: chain length and volume fraction effects on structural arrangement

Dissipative particle dynamics (DPD), a mesoscopic simulation approach, has been used to investigate the chain length effect on the structural property of the immiscible polyethylene (PE)/poly(L-lactide) (PLLA) polymer in a polymer blend and in a system with their diblock copolymer. In this work, the interaction parameter in DPD simulation, related to the Flory-Huggins interaction parameter chi, is estimated by the calculation of mixing energy for each pair of components in molecular dynamics simulation. The immiscibility property of PE and PLLA polymers induces the phase separation and exhibits different architectures at different volume fractions. In order to observe the structural property, the radius of gyration is used to observe the detailed arrangement of the polymer chains. It shows that the structure arrangement of a polymer chain is dependent on the phase structure and has a significantly different structural arrangement character for the very short chains in the homopolymer and copolymers. The chain length effect on the degree of stretching or extension of polymers has also been observed. As the chain length increases, the chain exhibits more stretching behavior at lamellae, perforated lamellae, and cylindrical configurations, whereas the chain exhibits a similar degree of stretching or extension at the cluster configuration.     

Conformations produced by interactions of the side chains in poly(silylenemethylene) with the repeating sequence [Si(CH3)R?CH2]x [R = ?O(CH2)3OC6H4C6H5]

A coarse-grained model has been developed for asymmetrically substituted poly(silylenemethylene)s in which the side chain is a flexible spacer terminated by a biphenyl unit. Each monomer unit is represented by four coarse-grained beads that interact via a Lennard–Jones potential and are subject to the first- and second-order interactions deduced from the atomistically detailed model. Metropolis Monte Carlo simulations were performed for isolated syndiotactic, isotactic, and atactic chains. Snapshots from the equilibrated coarse-grained chain on the discrete space of a high coordination lattice were reverse-mapped to atomistically detailed structures in continuous space. At 373 K, the chains were disordered independent of the stereochemical composition. The occupancy of bond pairs depended on the stereochemical composition, with the trans-gauche (tg) sequence being favored by the isotactic chain. When the simulation was performed with the backbone constrained to specific periodic structures, the g helix was the lowest energy structure for either the atactic or isotactic chains. For the syndiotactic chain, the g and gt helices were favored. The appearance of the g helix as the favored periodic structure of the isolated chain was consistent with the chain conformation reported previously for the smectic phase of this polymer in the bulk state. The g helix was disrupted when the backbone was allowed to access nonhelical conformations, even though these conformations may have been slightly higher in energy. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 43: 886–896, 2005     

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.