Validation of forcefields in predicting the physical and thermophysical properties of emeraldine base polyaniline

We report a molecular modelling study to validate the forcefields [condensed-phase optimised molecular potentials for atomistic simulation studies (COMPASS) and polymer-consistent forcefield (PCFF)] in predicting the physical and thermophysical properties of polymers. This work comprises of two key steps: (1) generating and validating the molecular model in predicting the material properties of the bulk amorphous emeraldine base polyaniline and (2) modelling the glass–rubber transition of the polymer. From all the molecular dynamics simulation results, it clearly shows that the more recent COMPASS forcefield provides a higher accuracy in predicting the polymer properties than PCFF, and it enables a more accurate prediction of condensed-phase properties (density, glass transition temperature, solubility parameters, etc.) in a broad range of temperature for various applications.     

Dynamics of Small Molecules in a Dense Polymer Matrix: Molecular Dynamics Studies

Abstract

The anomalous diffusion regime appearing in the self-diffusion of small molecules in bulk amorphous polymers has been extensively studied by molecular dynamics simulations. A rather long simulation of duration ～ 10 ^?8 s is performed on a polyethylene-like simple polymer model containing either oxygen molecules or helium atoms as a diffusant. Dynamic properties evaluated for these diffusants are the mean-square displacement, the van Hove self correlation function, and the self part of the density autocorrelation. It is first confirmed that the anomalous diffusion regime appears in a few hundred picoseconds for oxygen molecule, while the Einstein relation adopted beyond this regime results in the self-diffusion coefficient of the order of ～ 10^?5 cm²/s. This anomaly is still observed for helium that diffuses much faster than oxygen. In the anomalous diffusion regime, it is found that the correlation functions for the two diffusants show characteristic features and become essentially the same as time is scaled appropriately. These features allow the estimation of the two characteristic spatial scales which are probably dominated by the microstructure of the polymer matrix, namely, the cage size and the distance between adjacent cages. The time dependence of the mean-square displacements of the two diffusants can be well interpreted by these characteristic spatial scales as time is scaled with the self-diffusion coefficients. It is shown that the anomalous diffusion regime arises from the inhomogeneous microstructure of the polymer matrix.     

Nanoconfined polymers: modelling and simulation approaches

In this article, application of molecular simulation methods for studying molecular pictures of nanoconfined polymers is reviewed and discussed. The simulation methods, covering a range from atomistic to systematically parameterised coarse-grained models, employed in the literature to study nanoconfined polymers are reviewed and their results are compared together. The effect of polymer–surface interactions, surface curvature and surface area on the alteration of polymer structure and dynamics from the unperturbed (bulk) polymer properties are discussed. The length scales over which the surface influences the polymer structure and dynamics and the magnitude of surface effect on dynamics deceleration in the interphase are addressed in terms of different local and global chain properties.     

Structure-property studies on carbohydrate-derived polymers for use as protein-resistant biomaterials

Here we describe structure-property studies on our carbohydrate-derived side-chain ether polymers as protein-resistant biomaterials. A series of side-chain ether polymers, including two polyesters and two polyamides, were prepared by condensation polymerization of monomers derived from simple carbohydrates. The two side-chain permethoxylated polyesters having different stereochemical repeating units demonstrate excellent resistance toward nonspecific protein adsorption as shown by surface plasmon resonance, indicating that the polymer stereochemistry does not have much effect on its protein-resistant properties. The introduction of amide bonds to polymer backbones leads to more pronounced effects. While the polymer degradation stability is significantly enhanced by replacing ester with amide linkages, the protein resistance for the polymer is greatly reduced by introduction of amide bonds. Finally, our results suggest that free hydroxyl and amide groups, while both are hydrogen-bond donors, seem to have different effects on protein resistant properties for polymers. It appears that free amide groups have more detrimental effect on protein resistance than free hydroxyl groups. These results show that the protein-resistant properties of this family of polymers can be tailored by modifying the backbone and side chain functionalities. In combination with the biodegradability and functionalizability, this family of carbohydrate-derived polymers shows promise as versatile biomaterials for biomedical applications.     

Molecular modeling of temperature dependence of solubility parameters for amorphous polymers

A molecular modeling strategy is proposed to describe the temperature (T) dependence of solubility parameter (δ) for the amorphous polymers which exhibit glass-rubber transition behavior. The commercial forcefield "COMPASS" is used to support the atomistic simulations of the polymer. The temperature dependence behavior of δ for the polymer is modeled by running molecular dynamics (MD) simulation at temperatures ranging from 250 up to 650 K. Comparing the MD predicted δ value at 298 K and the glass transition temperature (T(g)) of the polymer determined from δ-T curve with the experimental value confirm the accuracy of our method. The MD modeled relationship between δ and T agrees well with the previous theoretical works. We also observe the specific volume (v), cohesive energy (U(coh)), cohesive energy density (E(CED)) and δ shows a similar temperature dependence characteristics and a drastic change around the T(g). Meanwhile, the applications of δ and its temperature dependence property are addressed and discussed.     

Investigation of moisture diffusion in cross-linked epoxy moulding compound by molecular dynamics simulation

Molecular dynamics (MD) simulation was carried out to study the moisture diffusion in cross-linked epoxy resin, with the influence of temperature, water concentration and polymer conversion taken into account. The simulation results showed that the moisture diffusion coefficients increased with the increase in temperature. And generally, with the increased moisture concentration or decreased polymer conversion, the moisture diffusion coefficients reduced. However, the moisture diffusion was strongly inhibited when the number of epoxy groups in completely reacted epoxy resins was equal to the number of water molecules.     

Molecular dynamics simulation of diffusion behaviour of gas molecules within oil–paper insulation system

The diffusion behaviour of hydrogen, carbon monoxide, carbon dioxide, methane, acetylene, ethylene and ethane in oil and paper medium was examined using molecular dynamics to reveal the diffusion mechanism of gas molecules in transformer oil–paper insulation system at the microscopic level. These compounds are commonly used in the dissolved gas analysis of power transformers and produced during the ageing process of oil–paper composite insulating material. Two groups of models were constructed using molecular dynamics simulation software to simulate the diffusion behaviour of the aforementioned seven types of small gas molecules in oil and paper. The diffusion coefficients, displacement features, free volume characteristics and interaction energies of the gas molecules were investigated. In particular, the diffusion micro-mechanism of the gas molecules was observed. The differences in diffusion features among the gas molecules were discussed, and the factors influencing the diffusion of the gas molecules were compared. Simulation results indicate that the diffusion coefficients of gas molecules in cellulose is an order of magnitude lower than that in oil, and the diffusion coefficients of these gas molecules in the two types of insulation media have different orders. Free volume of gas molecules is the main factor that influences the diffusion behaviour in oil, whereas intermolecular interaction is the main influencing factor of diffusion behaviour in cellulose.     

Temperature dependence of protein-hydration hydrodynamics by molecular dynamics simulations

The dynamics of water molecules near the protein surface are different from those of bulk water and influence the structure and dynamics of the protein itself. To elucidate the temperature dependence hydration dynamics of water molecules, we present results from the molecular dynamic simulation of the water molecules surrounding two proteins (Carboxypeptidase inhibitor and Ovomucoid) at seven different temperatures (T=273 to 303 K, in increments of 5 K). Translational diffusion coefficients of the surface water and bulk water molecules were estimated from 2 ns molecular dynamics simulation trajectories. Temperature dependence of the estimated bulk water diffusion closely reflects the experimental values, while hydration water diffusion is retarded significantly due to the protein. Protein surface induced scaling of translational dynamics of the hydration waters is uniform over the temperature range studied, suggesting the importance protein-water interactions.     

Self-diffusion of polymers in cartilage as studied by pulsed field gradient NMR

Pulsed field gradient (PFG) nuclear magnetic resonance (NMR) was used to investigate the self-diffusion behaviour of polymers in cartilage. Polyethylene glycol and dextran with different molecular weights and in different concentrations were used as model compounds to mimic the diffusion behaviour of metabolites of cartilage. The polymer self-diffusion depends extremely on the observation time: The short-time self-diffusion coefficients (diffusion time Delta approximately 15 ms) are subjected to a rather non-specific obstruction effect that depends mainly on the molecular weights of the applied polymers as well as on the water content of the cartilage. The observed self-diffusion coefficients decrease with increasing molecular weights of the polymers and with a decreasing water content of the cartilage. In contrast, the long-time self-diffusion coefficients of the polymers in cartilage (diffusion time Delta approximately 600 ms) reflect the structural properties of the tissue. Measurements at different water contents, different molecular weights of the polymers and varying observation times suggest that primarily the collagenous network of cartilage but also the entanglements of the polymer chains themselves are responsible for the observed restricted diffusion. Additionally, anomalous restricted diffusion was shown to occur already in concentrated polymer solutions.     

Carbohydrate polymers in amorphous states: an integrated thermodynamic and nanostructural investigation

The effect of water on the structure and physical properties of amorphous polysaccharide matrices is investigated by combining a thermodynamic approach including pressure- and temperature-dependent dilatometry with a nanoscale analysis of the size of intermolecular voids using positron annihilation lifetime spectroscopy. Amorphous polysaccharides are of interest because of a number of unusual properties which are likely to be related to the extensive hydrogen bonding between the carbohydrate chains. Uptake of water by the carbohydrate matrices leads to a strong increase in the size of the holes between the polymer chains in both the glassy and rubbery states while at the same time leading to an increase in matrix free volume. Thermodynamic clustering theory indicates that, in low-moisture carbohydrate matrices, water molecules are closely associated with the carbohydrate chains. Based on these observations, we propose a novel model of plasticization of carbohydrate polymers by water in which the water dynamically disrupts chains the hydrogen bonding between the carbohydrates, leading to an expansion of the matrix originating at the nanolevel and increasing the number of degrees of freedom of the carbohydrate chains. Consequently, even in the glassy state, the uptake of water leads to increased rates of matrix relaxation and mobility of small permeants. In contrast, low-molecular weight sugars plasticize the carbohydrate matrix without appreciably changing the structure and density of the rubbery state, and their role as plasticizer is most likely related to a reduction of the number of molecular entanglements. The improved molecular packing in glassy matrices containing low molecular weight sugars leads to a higher matrix density, explaining, despite the lower glass transition temperature, the reduced mobility of small permeants in such matrices.     

Theoretical study of the surface properties of poly(dimethylsiloxane) and poly(tetrafluoroethylene)

Molecular dynamics (MD) simulations of poly(dimethylsiloxane) (PDMS) and poly(tetrafluoroethylene) (PTFE) were carried out to determine their surface properties and energies. This study helps to gain better insight into the molecular modeling of PDMS and PTFE, in particular how different approaches affect calculations of surface energy. Current experimental and theoretical data were used to further understand the surface properties of PDMS and PTFE as well as to validate and verify results obtained from the combination of density functional theory (DFT) calculations (including periodic boundary conditions) and MD simulations. Detailed analysis of the structure and electronic properties (by calculation of the projected density of states) of the bulk and surface models of PDMS and PTFE was performed. The sensitivity of the surface energy calculation of these two polymers to the chemistry and model preparation was indicated. The balance between the molecular density, weight (which also reflects bond orientation in the surface region), bond flexibility, and intramolecular interactions including bond stretching was revealed to govern the results obtained. In modeling, the structural organization of polymer near a given surface (types and number of end groups and broken bonds due to application of different cut offs of the periodic structure) also significantly affects the final results. Besides the structural differences, certain simulation parameters, such the DFT functionals and simulation boxes utilized, play an important role in determining surface energy. The models used here were shown to be sufficient due to their good agreement with experimental and other theoretical data related to surface properties and surface energies.     

Quantifying the Effect of Polymer Blending through Molecular Modelling of Cyanurate Polymers

Modification of polymer properties by blending is a common practice in the polymer industry. We report here a study of blends of cyanurate polymers by molecular modelling that shows that the final experimentally determined properties can be predicted from first principles modelling to a good degree of accuracy. There is always a compromise between simulation length, accuracy and speed of prediction. A comparison of simulation times shows that 125ps of molecular dynamics simulation at each temperature provides the optimum compromise for models of this size with current technology. This study opens up the possibility of computer aided design of polymer blends with desired physical and mechanical properties.     

Diffusion and electrophoretic mobility of single-stranded RNA from molecular dynamics simulations

Hydrodynamic properties of small single-stranded RNA homopolymers with three and six nucleotides in free solution are determined from molecular dynamics simulations in explicit solvent. We find that the electrophoretic mobility increases with increasing RNA length, consistent with experiment. Diffusion coefficients of RNA, corrected for finite-size effects and solvent viscosity, agree well with those estimated from experiments and hydrodynamic calculations. The diffusion coefficients and electrophoretic mobilities satisfy a Nernst-Einstein relation in which the effective charge of RNA is reduced by the charge of transiently bound counterions. Fluctuations in the counterion atmosphere are shown to enhance the diffusive spread of RNA molecules drifting along the direction of the external electric field. As a consequence, apparent diffusion coefficients measured by capillary zone electrophoresis can be significantly larger than the actual values at certain experimental conditions.     

Dynamical Properties of Condensed Charged Polymer Melts

Abstract

Dynamical properties of condensed charged polymer melts are studied with a two-dimensional model and molecular dynamics simulation. Screened Coulombic interactions are assumed for the interactions between the monomer charges as well as the counterions which were introduced to neutralize the total monomer charges of polymer chains. Through molecular dynamics calculations, we have obtained the radial distribution function and velocity auto-correlation function, and their density dependences. As a structural characteristics in condensed charged polymer melts, we find that the monomers tend to form triangular structures locally. The radial distribution function for the center monomers implies that the polymer chains on our two-dimensional model are not entangled. The diffusion properties for both the counter-ions and innermost monomers are studied in detail. In this paper, we have also computed neutral polymer melts to study the effects of the presence of the long-range Coulombic interactions on the dynamical properties of polymer melts. We find that the Coulombic interactions significantly reduce the self diffusion. Snapshots analyses obtained from the molecular dynamics simulation suggest that the reptation model is not a proper model for two-dimensional polymer-chain melts, in which there are little entanglement effects in agreement with the result of the analyses for the radial distribution function.     

The effect of alumina nanoparticles on the thermal properties of PMMA: a molecular dynamics simulation

Metal oxides, as one of the most promising flame retardant additives, improve the fire retardant and the thermal stability properties of polymers. In the present study, molecular dynamics (MD) simulations based on the united atom model were applied to study the effect of alumina nanoparticles on the density, thermal conductivity, heat capacity, and thermal diffusivity of isotactic poly(methyl methacrylate) (is-PMMA). Thermal diffusivity of PMMA and PMMA/alumina nanocomposite were investigated through calculating thermal conductivity, density and heat capacity in the range of 300–700?K. Heat capacity can be calculated using fluctuations properties of energy. Thermal conductivity was calculated through the nonequilibrium molecular dynamics (NEMD) simulation by Fourier’s law approach. Our results show that the addition of alumina nanoparticles decreases the heat capacity and increases the glass transition temperature (T_g), thermal conductivity and thermal diffusivity of the PMMA. Therefore, the addition of alumina nanoparticles to PMMA improves the fire retardancy of the polymer. In addition, we illustrate the links between the intermolecular and bulk properties of PMMA in the presence of the alumina nanoparticles.     

Polymers of intrinsic microporosity for gas permeation: a molecular simulation study

We report a molecular simulation study for gas permeation in two membranes constructed from polymers of intrinsic microporosity (PIM-1 and PIM-7). With rigid ladder polymer chains, the membranes posses approximately 47.7 and 46.6% fractional free volumes (FFVs) in PIM-1 and PIM-7, respectively. The voids in the membranes have a diameter up to 9 Å and are largely interconnected. The sorption and diffusion of four gases (H₂, O₂, CH₄ and CO₂) were calculated by Monte Carlo and molecular dynamics simulations. The solubility coefficients increase in the order of H₂ < O₂ < CH₄ < CO₂, while the diffusion coefficients increase in the following order: CH₄ < CO₂ < O₂ < H₂. The simulation results agree well with experimental data, particularly for the solubility coefficients. The solubility and diffusion coefficients correlate well separately with the critical temperatures and effective diameters of gases. These molecular-based correlations can be used in the prediction for other gases. As attributed to the microporous structure, PIM-1 and PIM-7 outperform most glassy polymeric membranes in sorption and diffusion. PIM-1 has larger solubility and diffusion coefficients than PIM-7 because the cyano groups in PIM-1 lead to a stronger affinity and a larger FFV. The simulated solubility, diffusivity and permeation selectivities of CO₂/H₂, CO₂/O₂ and CO₂/CH₄ are consistent with experimental data. The quantitative microscopic understanding of gas permeation in the PIM membranes is useful for the new development of high-performance membranes.     

Osmotic properties of poly(ethylene glycols): quantitative features of brush and bulk scaling laws

From glycosylated cell surfaces to sterically stabilized liposomes, polymers attached to membranes attract biological and therapeutic interest. Can the scaling laws of polymer "brushes" describe the physical properties of these coats? We delineate conditions where the Alexander-de Gennes theory of polymer brushes successfully fits the intermembrane distance versus applied osmotic stress data of Kenworthy et al. for poly(ethylene glycol)-grafted multilamellar liposomes. We establish that the polymer density and size in the brush must be high enough that, in a bulk solution of equivalent monomer density, the polymer osmotic pressure is independent of polymer molecular weight (the des Cloizeaux semidilute regime of bulk polymer solutions). The condition that attached polymers behave as semidilute bulk solutions offers a rigorous criterion for brush scaling-law behavior. There is a deep connection between the behaviors of semidilute polymer solutions in bulk and polymers grafted to a surface at a density such that neighbors pack to form a uniform brush. In this regime, two-parameter unconstrained fits of the Alexander-de Gennes brush scaling laws to the Kenworthy et al. data yield effective monomer lengths of 3.3-3.6 A, which agree with structural predictions. The fitted distances between grafting sites are larger than expected from the nominal mole fraction of poly(ethylene glycol)-lipids; the chains apparently saturate the surface. Osmotic stress measurements can be used to estimate the actual densities of membrane-grafted polymers.     

Conformation of a neurokinin antagonist in solution. 2D NMR and restrained molecular dynamics study.

The hexapeptide [cyclo(Leu1 psi(CH2NH2)Leu2-Gln3-Trp4-Phe5-Gly6)]+1 is a potent antagonist of neurokinin A activity in tissues of hamster urinary bladder. The solution conformation of this cyclic hexapeptide has been characterized by the combined use of two dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics. The proton spectrum of the peptide was fully assigned by the sequential assignment procedure. Interproton distances were derived from crosspeak volumes in two dimensional Nuclear Overhauser Effect spectra, and dihedral angles were calculated from appropriate coupling constants. Temperature coefficients of the amide protons were determined. Restrained molecular dynamics simulations were carried out using the backbone interproton distances as constraints. During 210 ps of restrained molecular dynamics the peptide interconverted among three closely related families of conformations. These interconversions occurred at picosecond timescales under the simulation conditions.     

Molecular dynamics simulation of proton transport near the surface of a phospholipid membrane

The structural and dynamical properties of a hydrated proton near the surface of DMPC membrane were studied using a molecular dynamics simulation. The proton transport between water molecules was modeled using the second generation multistate empirical valence bond model. The proton diffusion was found to be inhibited at the membrane surface. The potential of mean force for the proton adsorption to the membrane surface and its release back into the bulk water was also determined, yielding a small barrier in each direction. An efficient algorithm for Ewald summation calculations for the multistate empirical valence bond model is also introduced.