Large deviations of Rouse polymer chain: first passage problem

The purpose of this paper is to investigate several analytical methods of solving first passage (FP) problem for the Rouse model, a simplest model of a polymer chain. We show that this problem has to be treated as a multi-dimensional Kramers' problem, which presents rich and unexpected behavior. We first perform direct and forward-flux sampling (FFS) simulations, and measure the mean first-passage time $\tau(z)$ for the free end to reach a certain distance $z$ away from the origin. The results show that the mean FP time is getting faster if the Rouse chain is represented by more beads. Two scaling regimes of $\tau(z)$ are observed, with transition between them varying as a function of chain length. We use these simulations results to test two theoretical approaches. One is a well known asymptotic theory valid in the limit of zero temperature. We show that this limit corresponds to fully extended chain when each chain segment is stretched, which is not particularly realistic. A new theory based on the well known Freidlin-Wentzell theory is proposed, where dynamics is projected onto the minimal action path. The new theory predicts both scaling regimes correctly, but fails to get the correct numerical prefactor in the first regime. Combining our theory with the FFS simulations lead us to a simple analytical expression valid for all extensions and chain lengths. One of the applications of polymer FP problem occurs in the context of branched polymer rheology. In this paper, we consider the arm-retraction mechanism in the tube model, which maps exactly on the model we have solved. The results are compared to the Milner-McLeish theory without constraint release, which is found to overestimate FP time by a factor of 10 or more.

Conformational features of cepacian: the exopolysaccharide produced by clinical strains of Burkholderia cepacia

Conformational energy calculations and molecular dynamics investigations, both in water and in dimethyl sulfoxide, were carried out on the exopolysaccharide cepacian produced by the majority of the clinical strains of Burkholderia cepacia, an opportunistic pathogen causing serious lung infection in patients affected by cystic fibrosis, the investigation was aimed at defining the structural and conformational features, which might be relevant for clarification of the structure-function relationships of the polymer. The molecular dynamics calculations were carried out by Ramachandran-type energy plots of the disaccharides that constitute the polymer repeating unit. The dynamics of an oligomer composed of three repeating units were investigated in water and in Me2SO, a non-aggregating solvent. Analysis of the time persistence of hydrogen bonds showed the presence of a large number of favourable interactions in water, which were less evident in Me2SO. The calculations on the cepacian chain indicated that polymer conformational features in water were affected by the lateral chains, but were also largely dictated by the presence of solvent. Moreover, the large number of intra-chain hydrogen bonds in water disappeared in Me2SO solution, increasing the average dimension of the polymer chains. (c) 2005 Elsevier Ltd. All rights reserved.

The physics of baking: rheological and polymer molecular structure-function relationships in breadmaking

Molecular size and structure of the gluten polymers that make up the major structural components of wheat are related to their rheological properties via modem polymer rheology concepts. Interactions between polymer chain entanglements and branching are seen to be the key mechanisms determining the rheology of HMW polymers. Recent work confirms the observation that dynamic shear plateau modulus is essentially independent of variations in MW amongst wheat varieties of varying baking performance and is not related to variations in baking performance, and that it is not the size of the soluble glutenin polymers, but the structural and rheological properties of the insoluble polymer fraction that are mainly responsible for variations in baking performance. The rheological properties of gas cell walls in bread doughs are considered to be important in relation to their stability and gas retention during proof and baking, in particular their extensional strain hardening properties. Large deformation rheological properties of gas cell walls were measured using biaxial extension for a number of doughs of varying breadmaking quality at constant strain rate and elevated temperatures in the range 25-60 degrees C. Strain hardening and failure strain of cell walls were both seen to decrease with temperature, with cell walls in good breadmaking doughs remaining stable and retaining their strain hardening properties to higher temperatures (60 degrees C), whilst the cell walls of poor breadmaking doughs became unstable at lower temperatures (45-50 degrees C) and had lower strain hardening. Strain hardening measured at 50 degrees C gave good correlations with baking volume, with the best correlations achieved between those rheological measurements and baking tests which used similar mixing conditions. As predicted by the Considere failure criterion, a strain hardening value of I defines a region below which gas cell walls become unstable, and discriminates well between the baking quality of a range of commercial flour blends of varying quality. This indicates that the stability of gas cell walls during baking is strongly related to their strain hardening properties, and that extensional rheological measurements can be used as predictors of baking quality. (C) 2004 Elsevier B.V. All rights reserved.

The physics of baking: rheological and polymer molecular structure-function relationships in breadmaking

Molecular size and structure of the gluten polymers that make up the major structural components of wheat are related to their rheological properties via modern polymer rheology concepts. Interactions between polymer chain entanglements and branching are seen to be the key mechanisms determining the rheology of HMW polymers. Recent work confirms the observation that dynamic shear plateau modulus is essentially independent of variations in MW amongst wheat varieties of varying baking performance and is not related to variations in baking performance, and that it is not the size of the soluble glutenin polymers, but the structural and rheological properties of the insoluble polymer fraction that are mainly responsible for variations in baking performance. The rheological properties of gas cell walls in bread doughs are considered to be important in relation to their stability and gas retention during proof and baking, in particular their extensional strain hardening properties. Large deformation rheological properties of gas cell walls were measured using biaxial extension for a number of doughs of varying breadmaking quality at constant strain rate and elevated temperatures in the range 25oC to 60oC. Strain hardening and failure strain of cell walls were both seen to decrease with temperature, with cell walls in good breadmaking doughs remaining stable and retaining their strain hardening properties to higher temperatures (60oC), whilst the cell walls of poor breadmaking doughs became unstable at lower temperatures (45oC to 50oC) and had lower strain hardening. Strain hardening measured at 50oC gave good correlations with baking volume, with the best correlations achieved between those rheological measurements and baking tests which used similar mixing conditions. As predicted by the Considere failure criterion, a strain hardening value of 1 defines a region below which gas cell walls become unstable, and discriminates well between the baking quality of a range of commercial flour blends of varying quality. This indicates that the stability of gas cell walls during baking is strongly related to their strain hardening properties, and that extensional rheological measurements can be used as predictors of baking quality.

Time-dependent orientation coupling in equilibrium polymer melts

The motion in concentrated polymer systems is described by either the Rouse or the reptation model, which both assume that the relaxation of each polymer chain is independent of the surrounding chains. This, however, is in contradiction with several experiments. In this Letter, we propose a universal description of orientation coupling in polymer melts in terms of the time-dependent coupling parameter κ(t). We use molecular dynamics simulations to show that the coupling parameter increases with time, reaching about 50% at long times, independently of the chain length or blend composition. This leads to predictions of component dynamics in mixtures of different molecular weights from the knowledge of monodisperse dynamics for unentangled melts. Finally, we demonstrate that entanglements do not play a significant role in the observed coupling. © 2010 The American Physical Society

Modeling entangled dynamics: comparison between stochastic single-chain and multichain models

To test the effectiveness of stochastic single-chain models in describing the dynamics of entangled polymers, we systematically compare one such model; the slip-spring model; to a multichain model solved using stochastic molecular dynamics(MD) simulations (the Kremer-Grest model). The comparison involves investigating if the single-chain model can adequately describe both a microscopic dynamical and a macroscopic rheological quantity for a range of chain lengths. Choosing a particular chain length in the slip-spring model, the parameter values that best reproduce the mean-square displacement of a group of monomers is determined by fitting toMDdata. Using the same set of parameters we then test if the predictions of the mean-square displacements for other chain lengths agree with the MD calculations. We followed this by a comparison of the time dependent stress relaxation moduli obtained from the two models for a range of chain lengths. After identifying a limitation of the original slip-spring model in describing the static structure of the polymer chain as seen in MD, we remedy this by introducing a pairwise repulsive potential between the monomers in the chains. Poor agreement of the mean-square monomer displacements at short times can be rectified by the use of generalized Langevin equations for the dynamics and resulted in significantly improved agreement.

Effect of headgroup size, charge, and solvent structure on polymer−micelle interactions, studied by molecular dynamics simulations

We performed atomistic molecular dynamics simulations of anionic and cationic micelles in the presence of poly(ethylene oxide) (PEO) to understand why nonionic water-soluble polymers such as PEO interact strongly with anionic micelles but only weakly with cationic micelles. Our micelles include sodium n-dodecyl sulfate (SDS), n-dodecyl trimethylammonium chloride (DTAC), n-dodecyl ammonium chloride (DAC), and micelles in which we artificially reverse the sign of partial charges in SDS and DTAC. We observe that the polymer interacts hydrophobically with anionic SDS but only weakly with cationic DTAC and DAC, in agreement with experiment. However, the polymer also interacts with the artificial anionic DTAC but fails to interact hydrophobically with the artificial cationic SDS, illustrating that large headgroup size does not explain the weak polymer interaction with cationic micelles. In addition, we observe through simulation that this preference for interaction with anionic micelles still exists in a dipolar "dumbbell" solvent, indicating that water structure and hydrogen bonding alone cannot explain this preferential interaction. Our simulations suggest that direct electrostatic interactions between the micelle and polymer explain the preference for interaction with anionic micelles, even though the polymer overall carries no net charge. This is possible given the asymmetric distribution of negative charges on smaller atoms and positive charges oil larger units in the polymer chain.

Molecular dynamics simulation of interactions between a sodium dodecyl sulfate micelle and a poly(ethylene oxide) polymer

We have performed atomistic molecular dynamics simulations of an anionic sodium dodecyl sulfate (SDS) micelle and a nonionic poly(ethylene oxide) (PEO) polymer in aqueous solution. The micelle consisted of 60 surfactant molecules, and the polymer chain lengths varied from 20 to 40 monomers. The force field parameters for PEO were adjusted by using 1,2-dimethoxymethane (DME) as a model compound and matching its hydration enthalpy and conformational behavior to experiment. Excellent agreement with previous experimental and simulation work was obtained through these modifications. The simulated scaling behavior of the PEO radius of gyration was also in close agreement with experimental results. The SDS-PEO simulations show that the polymer resides on the micelle surface and at the hydrocarbon-water interface, leading to a selective reduction in the hydrophobic contribution to the solvent-accessible surface area of the micelle. The association is mainly driven by hydrophobic interactions between the polymer and surfactant tails, while the interaction between the polymer and sulfate headgroups on the micelle surface is weak. The 40-monomer chain is mostly wrapped around the micelle, and nearly 90% of the monomers are adsorbed at low PEO concentration. Simulations were also performed with multiple 20-monomer chains, and gradual addition of polymer indicates that about 120 monomers are required to saturate the micelle surface. The stoichiometry of the resulting complex is in close agreement with experimental results, and the commonly accepted "beaded necklace" structure of the SDS-PEO complex is recovered by our simulations.

The effect of the spacer length on the nature of coupling in side chain liquid crystals polymers and elastomers

Side chain liquid crystal polymers and elastomers exhibit a rich phase behaviour which arises from the antagonistic influences of the entropically disordered polymer chain configuration and the long range orientational ordering of the mesogenic units. This competition arises since the natural macroscopic phase separation is inhibited by the inherent chemical connectivity of the system. At the heart of this connectivity is the spacer link and we consider here its influence on the phase behaviour. In particular we consider a series of elastomers in which the number of alkyl units in the spacer is systematically varied from 2 to 6. The lengthening of the coupling spacer is accompanied by an alternation of the sign of coupling between the polymer chain and the mesogenic unit. These results demonstrate the dominating influence of the so-called hinge effect in determining the phase behaviour. In addition to the alternation of the sign there is some decrease in the magnitude of the coupling with increasing spacer length.

Modelling the structure of a polymer glass poly(methylmethacrylate) through neutron scattering experiments

Determination of the local structure of a polymer glass by scattering methods is complex due to the number of spatial and orientational correlations, both from within the polymer chain (intrachain) and between neighbouring chains (interchain), from which the scattering arises. Recently considerable advances have been made in the structural analysis of relatively simple polymers such as poly(ethylene) through the use of broad Q neutron scattering data tightly coupled to atomistic modelling procedures. This paper presents the results of an investigation into the use of these procedures for the analysis of the local structure of a-PMMA which is chemically more complex with a much greater number of intrachain structural parameters. We have utilised high quality neutron scattering data obtained using SANDALS at ISIS coupled with computer models representing both the single chain and bulk polymer system. Several different modelling approaches have been explored which encompass such techniques as Reverse Monte Carlo refinement and energy minimisation and their relative merits and successes are discussed. These different approaches highlight structural parameters which any realistic model of glassy atactic PMMA must replicate.

Raman spectra of polymer electrolytes based on poly(ethylene glycol) dimethyl ether, lithium perchlorate, and the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate

Raman spectra of polymer electrolytes based on poly(ethylene glycol) dimethyl ether (PEGdME) with LiClO(4), PEGdME/LiClO(4), and the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate, PEGdME/[bmim]PF(6), are compared. Raman spectroscopy suggests stronger interactions in PEGdME/LiClO(4) than PEGdmE/[bmim]PF(6), thus corroborating previous results obtained by molecular dynamics simulations. Quantum Chemistry methods have been used to calculate vibrational frequencies and the equilibrium structure of segments of the polymer chain around the cation. A consistent picture has been obtained from Raman spectroscopy, density functional theory (DFT) calculations, and molecular dynamics simulations for these polymer electrolytes. (C) 2010 Elsevier B.V. All rights reserved.

Electric and magnetic properties of polymer electrolyte/carbon black composites

Measurements of 1H Nuclear Magnetic Resonance (NMR) relaxation times, Electron Paramagnetic Resonance (EPR) and AC Impedance Spectroscopy (IS) are reported for composites based on PEO8:LiClO4 and carbon black (CB), prepared by two methods: solvent and fusion processing. Three nuclear relaxation processes were identified for 1H nuclei: (i) belonging to the polymer chains in the amorphous phase, loosely bound to the CB particles, whose dynamics is almost the same as for unfilled polymer, (ii) belonging to the polymer chains which are tightly attached to the CB particles, and (iii) belonging to the crystalline phase in the loose polymer chain. The paramagnetic electronic susceptibility of the composite samples, measured by EPR, was interpreted by assuming a contribution of localized spin states that follow a Curie law, and a Pauli-like contribution of delocalized spins. A significant change of the EPR linewidth was observed at 40 K, which is the temperature where the Curie and Pauli susceptibilities equally contribute to the paramagnetic electronic susceptibility. The electrical properties are very sensitive to the preparation methods of the composites, which conditions the interaction between carbon particle-carbon particle and carbon particle-polymer chain. Classical statistic models to describe the conductivity in these media were not satisfactory. © 1998 Published by Elsevier Science B.V. All rights reserved.

PWA-diureasils organic-inorganic hybrids: photochromism and effect of the organic chain length

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Self-Assembly of Poly(4-vynil-N-alkylpyridinium) Bromide onto Silica: Effect of Side-Chain Length on Structural Aspects at a Molecular Level

Self-assembly of poly(4-vynil-N-alkyl)pyridinium bromide with alkyl side chains of 2, 5, 7, 10, or 16 carbons from ethanolic solutions onto flat silica surfaces was studied by means of ellipsometry, atomic force microscopy (AFM), contact angle measurements, and sum-frequency generation (SFG) vibrational spectroscopy in the CH3 and CH2 stretch region. Ab initio quantum-chemical calculations on the N-alkylpyridinium side-group with restricted Hartree-Fock (RHF) method and 6-311G (d,p) basis set were C one to estimate the charge distribution along the pyridinium ring and the alkyl side-chain. SFG results showed that longer side chains promote the disorientation of the alkyl groups at the surface, corroborating with the contact angle values. AFM images revealed film homogeneity, regardless the alkyl side group. However, after 24 h contact with water, ringlike structures appeared on the film surfaces, when the polycation alkyl side chain had 7 or less carbons, and as the alkyl chain increased to 10 or 16 carbons, the films dewetted because the hydrophobic interactions prevailed over the electrostatic interactions between the pyridinium charged groups and the negatively charged SiO2 surface. Under acid conditions (HCl 0.1 mol.L-1), the film mean thickness values decreased up to 50% of original values when the alkyl side chains were ethyl or pentyl groups due to ion-pair disruption, but for longer groups they remained unchanged. Quantum-chemical optimization and Mulliken electron population showed that (i) from C2 to C15 the positive charge at the headgroup (HG) decreased 0.025, while the charge at combined HG + alpha-CH2 increased 0.037; and (ii) for C6 or longer, the alkyl side group presents a tilt in the geometry, moving away from the plane. Such effects summed up over the whole polymer chain give support to suggest that when the side chains are longer than 7 carbons, the hydrophobic interaction decreases film stability and increases acid resistance.

Coarse-graining methods and polymer dynamics

This thesis studies molecular dynamics simulations on two levels of resolution: the detailed level of atomistic simulations, where the motion of explicit atoms in a many-particle system is considered, and the coarse-grained level, where the motion of superatoms composed of up to 10 atoms is modeled. While atomistic models are capable of describing material specific effects on small scales, the time and length scales they can cover are limited due to their computational costs. Polymer systems are typically characterized by effects on a broad range of length and time scales. Therefore it is often impossible to atomistically simulate processes, which determine macroscopic properties in polymer systems. Coarse-grained (CG) simulations extend the range of accessible time and length scales by three to four orders of magnitude. However, no standardized coarse-graining procedure has been established yet. Following the ideas of structure-based coarse-graining, a coarse-grained model for polystyrene is presented. Structure-based methods parameterize CG models to reproduce static properties of atomistic melts such as radial distribution functions between superatoms or other probability distributions for coarse-grained degrees of freedom. Two enhancements of the coarse-graining methodology are suggested. Correlations between local degrees of freedom are implicitly taken into account by additional potentials acting between neighboring superatoms in the polymer chain. This improves the reproduction of local chain conformations and allows the study of different tacticities of polystyrene. It also gives better control of the chain stiffness, which agrees perfectly with the atomistic model, and leads to a reproduction of experimental results for overall chain dimensions, such as the characteristic ratio, for all different tacticities. The second new aspect is the computationally cheap development of nonbonded CG potentials based on the sampling of pairs of oligomers in vacuum. Static properties of polymer melts are obtained as predictions of the CG model in contrast to other structure-based CG models, which are iteratively refined to reproduce reference melt structures. The dynamics of simulations at the two levels of resolution are compared. The time scales of dynamical processes in atomistic and coarse-grained simulations can be connected by a time scaling factor, which depends on several specific system properties as molecular weight, density, temperature, and other components in mixtures. In this thesis the influence of molecular weight in systems of oligomers and the situation in two-component mixtures is studied. For a system of small additives in a melt of long polymer chains the temperature dependence of the additive diffusion is predicted and compared to experiments.