Polycarbonate Simulations with a Density Functional Based Force Field

Full-Atomistic Optimized Potentials for Liquid Simulations and Polymer Consistent Force Field Models for Biocompatible Shape-Memory Poly(ε-caprolactone)

Thermally induced shape memory poly(ε-caprolactone) (PCL)-based polymers are one of the most extensively researched families of biocompatible materials. They are degradable under physiological conditions and have high applicability in general biomedical engineering, with cross-linked PCL networks being particularly useful for tissue engineering. In this study, we used the optimized potentials for liquid simulations (OPLS) force field, which is well suited for describing intermolecular interactions in biomolecules, and the class II polymer consistent force field (PCFF) to investigate the properties of telechelic PCL with diacrylates as reactive functionalities on its end groups. PCFF has been specifically parameterized for simulating synthetic polymeric materials. We compare the findings of all-atom molecular dynamics simulations with known experimental data and theoretical assumptions to verify the applicability of both these force fields. We estimated the melt density, volume, transition temperatures, and mechanical characteristics of two-branched PCL diacrylates with a molecular weight of 2481 Da. Our findings point to the utility of the aforementioned force fields in predicting the properties of PCL-based polymers. It also opens avenues for developing PCL cross-linked polymer models and employing OPLS to investigate the interactions of synthetic polymers with biomolecules.

Effect of chemical substituents on the structure of glassy diphenyl polycarbonates

Polycarbonates offer a wide variety of physical property behavior that is difficult to predict due to complexities at the molecular scale. Here, the physical structure of amorphous glassy polycarbonates having aliphatic and cycloaliphatic chemical groups is explored through atomistic simulations. The influence of chemical structure on solubility parameter, torsion distributions, radial distribution function, scattering structure factor, orientation distributions of phenylene rings and carbonate groups, and free volume distributions, leading to interchain packing effects, are shown. The effect of the cyclohexyl ring at the isopropylidene carbon as compared to the effect of the methyl groups positioned on the phenylene rings results in a larger reduction in the solubility parameter (δ). The interchain distance estimated for polycarbonates in this work is in the range of 5-5.8 Å. The o-methyl groups on the phenylene rings, as compared to a cyclohexyl ring, lead to higher interchain distances. The highest interchain distance is observed with a trimethylcyclohexylidene group at the isopropylidene carbon. Atomistic simulations reveal two different types of packing arrangement of nearest-neighbor chains in the glassy state, one type of which agrees with the NMR experimental data. The fundamental insights provided here can be utilized for design of chemical structures for tailored macroscopic properties.

Simulation of the plastic behavior of amorphous glassy bis-phenol-A-polycarbonate

A protocol for studying the plastic deformation of amorphous glassy polymers is presented. The protocol is based on a viable computational procedure which combines constant-stress molecular dynamics simulations and fixed-cell energy minimizations, followed by kinetic, configurational, and energy analyses. It is shown that the computational results can be accounted for within a "potential energy landscape" theoretical framework, in which the plastic transition is interpreted as a crossing between and a collapse onto each other of "ideal (thermodynamic) structures." The procedure is applied to bis-phenol-A-polycarbonate (BPA-PC), but is equally valid for a wide variety of polymeric species. Allowing for the limited size of the simulation cell, the high strain rate, and the fact that the simulation are conducted at low temperature, the values of the density, Young's modulus, yield strain, yield stress, activation energy, and activation volume are in fair agreement with the experimental data on BPA-PC. The analysis of the results shows that the plastic relaxation for this polymer has both a collective and cooperative character (as in classical percolation theories), involves a significant fraction of the simulation cell, and can be viewed as a "nanoscopic shear band."

BPA-PC on a Ni111 surface: the interplay between adsorption energy and conformational entropy for different chain-end modifications

We extend a previous dual scale modeling approach for the behavior of polymers near a metal surface to a variety of end groups. Our approach combines a coarse-grained polymer model with ab initio DFT calculations. Such a procedure was applied to a melt of phenolic-like terminated Bisphenol A-polycarbonate (BPA-PC) interacting with a (111) nickel surface (Delle Site, L.; Abrams, C. F.; Alavi, A.; Kremer, K. Phys. Rev. Lett. 2002, 89, 156103. Abrams, C. F.; Delle Site, L.; Kremer, K. Phys. Rev. E 2003, 67, 021807). This work extends this study to different chain-end modifications of BPA-PC, p-tert-butylphenolic, p-tetramethylpropylphenolic, and p-cumylphenolic. We show how the interplay between adsorption energies and conformational entropy selects different morphologies for the various melts at the interface. Implications of these results for realistic technical materials are finally discussed.

Dual-resolution coarse-grained simulation of the bisphenol-A-polycarbonate/nickel interface

We present a dual-resolution coarse-graining scheme for efficient molecular dynamics simulations of bisphenol-A-polycarbonate (BP-A-PC) liquids in contact with a (111) nickel surface. The essential feature of this model is the strong adsorption of phenoxy chain ends, and the absence of adsorption by other parts of the chains. Details of how phenoxy chain ends interact with the nickel surface were extracted from Car-Parrinello molecular dynamics calculations of adsorption of phenol on nickel. These calculations show that phenol adsorption on nickel is short ranged (<3 A) and strongly dependent on the C1-C4 orientation of the ring. The structure of BP-A-PC prevents internal phenylene groups from interacting with the surface, due to steric hindrances from the noninteracting isopropylidenes. These dependencies are incorporated in the coarse-grained model of the BP-A-PC chain by resolving chain-terminating carbonate groups with atomistic detail, while the rest of the chain is represented by coarsened "beads." This allows specification of the C1-C4 orientation of the terminal phenoxy groups, while overall allowing for system equilibration with reasonable computer time. We simulate liquids of up to 240 chains of ten chemical repeat units, confined in a slit pore formed by two frozen (111) planes of atoms with the lattice spacing of nickel. We find that the strong adsorption of chain ends has a large effect on the liquid structure through a distance of more than two bulk radii of gyration from the surface. These effects are explained by a competition among single- and double-end adsorption, and dense packing. The structure of the interface less than 10 A from the wall is greatly sensitive to the orientational dependence of the phenoxy adsorption.