Equilibrium polymerization of cyclic carbonate oligomers. III. Chain branching and the gel transition

Ring-opening polymerization of cyclic polycarbonate oligomers, where monofunctional active sites act on difunctional monomers to produce an equilibrium distribution of rings and chains, leads to a "living polymer." Monte Carlo simulations [two-dimensional (2D) and three-dimensional (3D)] of the effects of single [J. Chem. Phys. 115, 3895 (2001)] and multiple active sites [J. Chem. Phys. 116, 7724 (2002)] are extended here to trifunctional active sites that lead to branching. Low concentrations of trifunctional particles c(3) reduce the degree of polymerization significantly in 2D, and higher concentrations (up to 32%) lead to further large changes in the phase diagram. Gel formation is observed at high total density and sizable c(3) as a continuous transition similar to percolation. Polymer and gel are much more stable in 3D than in 2D, and both the total density and the value of c(3) required to produce high molecular weight aggregates are reduced significantly. The degree of polymerization in high-density 3D systems is increased by the addition of trifunctional monomers and reduced slightly at low densities and low c(3). The presence of branching makes equilibrium states more sensitive (in 2D and 3D) to changes in temperature T. The stabilities of polymer and gel are enhanced by increasing T, and-for sufficiently high values of c(3)-there is a reversible polymer-gel transformation at a density-dependent floor temperature. (C) 2002 American Institute of Physics.

A reactive force field simulation of liquid-liquid phase transitions in phosphorus

A force field model of phosphorus has been developed based on density functional (DF) computations and experimental results, covering low energy forms of local tetrahedral symmetry and more compact (simple cubic) structures that arise with increasing pressure. Rules tailored to DF data for the addition, deletion, and exchange of covalent bonds allow the system to adapt the bonding configuration to the thermodynamic state. Monte Carlo simulations in the N-P-T ensemble show that the molecular (P-4) liquid phase, stable at low pressure P and relatively low temperature T, transforms to a polymeric (gel) state on increasing either P or T. These phase changes are observed in recent experiments at similar thermodynamic conditions, as shown by the close agreement of computed and measured structure factors in the molecular and polymer phases. The polymeric phase obtained by increasing pressure has a dominant simple cubic character, while the polymer obtained by raising T at moderate pressure is tetrahedral. Comparison with DF results suggests that the latter is a semiconductor, while the cubic form is metallic. The simulations show that the T-induced polymerization is due to the entropy of the configuration of covalent bonds, as in the polymerization transition in sulfur. The transition observed with increasing P is the continuation at high T of the black P to arsenic (A17) structure observed in the solid state, and also corresponds to a semiconductor to metal transition. (C) 2004 American Institute of Physics.

Density functional and Monte Carlo studies of sulfur. II. Equilibrium polymerization of the liquid phase

The equilibrium polymerization of sulfur is investigated by Monte Carlo simulations. The potential energy model is based on density functional results for the cohesive energy, structural, and vibrational properties as well as reactivity of sulfur rings and chains [Part I, J. Chem. Phys. 118, 9257 (2003)]. Liquid samples of 2048 atoms are simulated at temperatures 450less than or equal toTless than or equal to850 K and P=0 starting from monodisperse S-8 molecular compositions. Thermally activated bond breaking processes lead to an equilibrium population of unsaturated atoms that can change the local pattern of covalent bonds and allow the system to approach equilibrium. The concentration of unsaturated atoms and the kinetics of bond interchanges is determined by the energy DeltaE(b) required to break a covalent bond. Equilibrium with respect to the bond distribution is achieved for 15less than or equal toDeltaE(b)less than or equal to21 kcal/mol over a wide temperature range (Tgreater than or equal to450 K), within which polymerization occurs readily, with entropy from the bond distribution overcompensating the increase in enthalpy. There is a maximum in the polymerized fraction at temperature T-max that depends on DeltaE(b). This fraction decreases at higher temperature because broken bonds and short chains proliferate and, for Tless than or equal toT(max), because entropy is less important than enthalpy. The molecular size distribution is described well by a Zimm-Schulz function, plus an isolated peak for S-8. Large molecules are almost exclusively open chains. Rings tend to have fewer than 24 atoms, and only S-8 is present in significant concentrations at all T. The T dependence of the density and the dependence of polymerization fraction and degree on DeltaE(b) give estimates of the polymerization temperature T-f=450+/-20 K. (C) 2003 American Institute of Physics.

Density functional/Monte Carlo study of ring-opening polymerization

Density functional calculations of the structure, potential energy surface and reactivity for organic systems closely related to bisphenol-A-polycarbonate (BPA-PC) provide the basis for a model describing the ring-opening polymerization of its cyclic oligomers by nucleophilic molecules. Monte Carlo simulations using this model show a strong tendency to polymerize that is increased by increasing density and temperature, and is greater in 3D than in 2D. Entropy in the distribution of inter-particle bonds is the driving force for chain formation. (C) 2002 Elsevier Science B.V. All rights reserved.

Equilibrium polymerization of cyclic carbonate oligomers. II. Role of multiple active sites

Ring opening polymerization of bisphenol A polycarbonate is studied by Monte Carlo simulations of a model comprising a fixed number of Lennard-Jones particles and harmonic bonds [J. Chem. Phys. 115, 3895 (2001)]. Bond interchanges produced by a low concentration (0.10%less than or equal toc(a)less than or equal to0.36%) of chemically active particles lead to equilibrium polymerization. There is a continuous transition in both 2D and 3D from unpolymerized cyclic oligomers at low density to a system of linear chains at high density, and the polymeric phase is much more stable in three dimensions than in two. The steepness of the polymerization transition increases rapidly as c(a) decreases, suggesting that it is discontinuous in the limit c(a)-->0. The transition is entropy driven, since the average potential energy increases systematically upon polymerization, and there is a steady decline in the degree of polymerization as the temperature is lowered. The mass distribution functions for open chains and for rings are unimodal, with exponentially decaying tails that can be fitted by Zimm-Schulz functions and simpler exponential forms. (C) 2002 American Institute of Physics.

Equilibrium polymerization of cyclic carbonate oligomers

A model of the polymerization of ring oligomers of bisphenol A polycarbonate (BPA-PC) is used to investigate the influence of dimensionality (2D or 3D), density and temperature on the size distribution of the polymer chains. The polymerization step is catalyzed by a single active particle, conserves the number and type of the chemical bonds, and occurs without a significant gain in either potential energy or configurational entropy. Monte Carlo and molecular dynamics simulations show that polymerization of cyclic oligomers occurs readily at high density and is driven by the entropy associated with the distribution of interparticle bonds. Polymerization competes at lower densities with long range diffusion, which favors small molecular species, and is prevented if the system is sufficiently dilute. Polymerization occurs in 2D via a weakly first order transition as a function of density and is characterized by low hysteresis and large fluctuations in the size of polymer chains. Polymerization occurs more readily in 3D than in 2D, and is favored by increasing temperature, as expected for an entropy-driven process. (C) 2001 American Institute of Physics.