Dissolved organic carbon--contaminant interaction descriptors found by 3D force field calculations

Enthalpies of transfer at 300 K of various partitioning processes were calculated in order to study the suitability of 3D force fields for the calculation of partitioning constants. A 3D fulvic acid (FA) model of dissolved organic carbon (DOC) was built in a MM+ force field using AMI atomic charges and geometrical optimization (GO). 3,5-Dichlorobiphenyl (PCB14), 4,4'-dichlorobiphenyl (PCB15), 1,1,1-trichloro-2,2-bis-(4-chlorophenyl)-ethane (PPDDT) and 2-chloro-4-ethylamino-6-isopropylamino-s-triazine (Atrazine) were inserted into different sites and their interaction energies with FA were calculated. Energies of hydration were calculated and subtracted from FA-contaminant interactions of selected sites. The resulting values for the enthalpies of transfer from water to DOC were 2.8, -1.4, -6.4 and 0.0 kcal/mol for PCB 14, PCB15, PPDDT and Atrazine, respectively. The value of PPDDT compared favorably with the experimental value of -5.0 kcal/mol. Prior to this, the method was studied by the calculation of the enthalpies of vaporization and aqueous solution using various force fields. In the MM + force field GO predicted enthalpies of vaporization deviated by +0.7 (PCB14), +3.6 (PCB15) and -0.7 (PPDDT)kcal/mol from experimental data, whereas enthalpies of aqueous solution deviated by -3.6 (PCB14), +5.8 (PCB15) and +3.7 (PPDDT) kcal/mol. Only for PCB14 the wrong sign of this enthalpy value was predicted. Potential advantages and limitations of the approach were discussed.     

An improved empirical force field for saturated hydrocarbons

A new molecular mechanics force field for alkanes is presented. The force field aims to eliminate some identified failures of the well-known MM2 force field. The new energy function gives an improved prediction of the rotational barriers of highly congested molecules, a better calculation of short nonbonded contacts, and the correct reproduction of bond elongation in small torsion angles. The calculation of sublimation enthalpies is also improved. The standard deviation of the formation enthalpies for a set of 54 compounds is 0.63 kcal/mol; this compares with the reported value of 0.42 calculated with MM2 and MM3 for different sets. The force field parameters were obtained using a least squares method.     

Heats of formation of [2.2]paracyclophane-1-ene and [2.2]paracyclophane-1,9-diene - an experimental study

The enthalpies of formation [Delta(g)] of tricyclo[8.2.2.2(4,7)]hexadeca-1(13),2,4(16),5,7(15),10(14),11-heptaene (2, 1,2-dehydro[2.2]paracyclophane or [2.2]paracyclophane-1-ene) and tricyclo[8.2.2.2(4,7)]hexadeca-1(13),2,4(16),5,7(15),8,10(14),11-octaene (3, 1,2,9,10-dehydro[2.2]paracyclophane or [2.2]paracyclophane-1,9-diene) have been determined by measuring their heats of combustion in a microcalorimeter and their heats of sublimation by the transpiration method. Values of the strain energies (SE) [SE(2) = 34.7 kcal mol(-)(1), SE(3) = 42.0 kcal mol(-)(1)] have been derived from the gas-phase heats of formation and are compared with those from MM3 and PM3 calculations and with the corresponding value SE(1) = 30.1 kcal mol(-)(1) for the parent tricyclo[8.2.2.2(4,7)]hexadeca-1(13),4(16),5,7(15),10(14),11-hexaene (1, [2.2]paracyclophane). The higher strain energies of 2 and 3 (by 4.6 and 11.9 kcal mol(-)(1)) are in accord with the well-known increased reactivities of their aromatic rings as a consequence of their increased bending. As revealed by an X-ray crystal structure analysis, the bending in the monoene 2 corresponds to that of 1 and 3 at one of two bridging corners.     

mezo-和d, l-2, 3-二氯丁烷气态生成焓的测定

本文用直接量热法测量了298.15K时meso-和d, l-2, 3-二氯丁烷的液态蒸发焓值以及其混合物的液态燃烧焓值, 得到了它们的标准气态生成焓值, 考察了二氯烷烃分子内氯原子间的相互影响和相互作用对其气态生成焓值的影响。     

Polar group enhanced gas-phase acidities of carboxylic acids: an investigation of intramolecular electrostatic interaction

We studied the effects of polar groups on the gas-phase acidities of carboxylic acids experimentally and computationally. In this connection, the gas-phase acidities (DeltaH(acid), the enthalpy of deprotonation, and DeltaG(acid), the deprotonation free energy) of borane-complexed methylaminoacetic acid ((CH(3))2N(BH(3))CH(2)CO(2)H) and methylthioacetic acid (CH(3)S(BH(3))CH(2)CO(2)H) were measured using the kinetic method in a flowing afterglow-triple quadrupole mass spectrometer. The values of DeltaH(acid) and DeltaG(acid) of (CH(3))2N(BH(3))CH(2)CO(2)H were determined to be 328.8 +/- 1.9 and 322.1 +/- 1.9 kcal/mol, and those of CH(3)S(BH(3))CH(2)CO(2)H were determined to be 325.8 +/- 1.9 and 319.2 +/- 1.9 kcal/mol, respectively. The theoretical enthalpies of deprotonation of (CH(3))2N(BH(3))CH(2)CO(2)H (329.2 kcal/mol) and CH(3)S(BH(3))CH(2)CO(2)H (325.5 kcal/mol) were calculated at the B3LYP/6-31+G(d) level of theory. The calculated enthalpies of deprotonation of N-oxide-acetic acid (CH(3)NOCH(2)CO(2)H, 329.4 kcal/mol) and S-oxide-acetic acid (CH(3)SOCH(2)CO(2)H, 328.6 kcal/mol) are comparable to the experimental results for borane-complexed methylamino- and methylthioacetic acids. The enthalpy of deprotonation of sulfone-acetic acid (CH(3)SO2CH(2)CO(2)H, 326.1 kcal/mol) is about 2 kcal/mol lower than the S-oxide-acetic acid. The calculated enthalpy of deprotonation of sulfoniumacetic acid, (CH(3))2S+CH(2)CO(2)H, is 243.0 kcal/mol. Compared to the corresponding reference molecules, CH(3)NHCH(2)CO(2)H and CH(3)SCH(2)CO(2)H, the dipolar group and the monopolar group substituted carboxylic acids are stronger acids by 11-14 and 97 kcal/mol, respectively. We correlated the changes of the acidity upon a polar group substitution to the electrostatic free energy within the carboxylate anion. The acidity enhancements in polar group substituted carboxylic acids are the results of the favorable electrostatic interactions between the polar group and the developing charge at the carboxyl group.     

Three-dimensional quantitative structure-property relationship (3D-QSPR) models for prediction of thermodynamic properties of polychlorinated biphenyls (PCBs): enthalpy of vaporization

Three-dimensional Quantitative Structure-Property Relationship (QSPR) models have been derived using Comparative Molecular Field Analysis (CoMFA) to correlate the vaporization enthalpies of a representative set of polychlorinated biphenyls (PCBs) at 298.15 K with their CoMFA-calculated physicochemical properties. Various alignment schemes, such as inertial, as is, and atom fit, were employed in this study. The CoMFA models were also developed using different partial charge formalisms, namely, electrostatic potential (ESP) charges and Gasteiger-Marsili (GM) charges. The most predictive model for vaporization enthalpy (Delta(vap)H(m)(298.15 K)), with atom fit alignment and Gasteiger-Marsili charges, yielded r2 values 0.852 (cross-validated) and 0.996 (conventional). The vaporization enthalpies of PCBs increased with the number of chlorine atoms and were found to be larger for the meta- and para-substituted isomers. This model was used to predict Delta(vap)H(m)(298.15 K) of the entire set of 209 PCB congeners.     

Composition and thermochemistry of silver bromide vapor

The mole fractions of AgBr and Ag3Br3 in the saturated vapor at 840 K have been evaluated from the vapor mass spectrum, by comparison with the corresponding spectrum of AgCl vapor, where the monomer/trimer ratio is known accurately from vapor molecular weight measurements. Combination of these results with new measurements of the vapor pressure of molten AgBr by the torsion-effusion method in the range 805-936 K yielded the third law enthalpies of vaporization and the standard enthalpies of formation DeltafH degrees 298(AgBr, g) = 27.8 +/- 0.3 kcal mol(-1) and DeltafH degrees 298(Ag3Br3, g) = -19.0 +/- 1 kcal mol(-1). The dissociation energy, D degrees 0(AgBr), is found to be 66.4 +/- 0.3 kcal mol(-1), or 2.88 +/- 0.01 eV, some 3.5-5 kcal mol(-1) lower than previous literature values. Approximate thermochemical stabilities of the dimer species Ag2Cl2 and Ag2Br2 have also been evaluated.     

Average sequential water molecule binding enthalpies of M(H2O)(19-124)2+ (M = Co, Fe, Mn, and Cu) measured with ultraviolet photodissociation at 193 and 248 nm

The average sequential water molecule binding enthalpies to large water clusters (between 19 and 124 water molecules) containing divalent ions were obtained by measuring the average number of water molecules lost upon absorption of an UV photon (193 or 248 nm) and using a statistical model to account for the energy released into translations, rotations, and vibrations of the products. These values agree well with the trend established by more conventional methods for obtaining sequential binding enthalpies to much smaller hydrated divalent ions. The average binding enthalpies decrease to a value of ~10.4 kcal/mol for n > ~40 and are insensitive to the ion identity at large cluster size. This value is close to that of the bulk heat of vaporization of water (10.6 kcal/mol) and indicates that the structure of water in these clusters may more closely resemble that of bulk liquid water than ice, owing either to a freezing point depression or rapid evaporative cooling and kinetic trapping of the initial liquid droplet. A discrete implementation of the Thomson equation using parameters for liquid water at 0 °C generally fits the trend in these data but provides values that are ~0.5 kcal/mol too low.     

Conformational study of cis-1,4-di-tert-butylcyclohexane by dynamic NMR spectroscopy and computational methods. Observation of chair and twist-boat conformations

[reaction: see text] Low-temperature 13C NMR spectra of cis-1,4-di-tert-butylcyclohexane (1) showed signals for the twist-boat (1a) and chair (1b) conformations. 13C NMR signals were assigned to specific carbons based on the different populations, different symmetries (time-averaged C(2v) for 1a and time-averaged C(s) for 1b), and calculated chemical shifts (GIAO, HF/6-311+G*). In addition to slow ring inversion and interconversion of the chair and twist-boat conformations, slow rotation of the tert-butyl groups was found. Most of the expected 13C peaks were observed. Free-energy barriers of 6.83 and 6.35 kcal/mol were found for interconversion of 1a (major) and 1b (minor) at -148.1 degrees C. Conformational space was searched with Allinger's MM3 and MM4 programs, and free energies were obtained for several low-energy conformations 1a-c. Calculations were repeated with ab initio methods up to the HF/6-311+G* level. Molecular symmetries, relative free energies, relative enthalpies and entropies, frequencies, and NMR chemical shifts were obtained. A boat conformation (1d; C(2v) symmetry) was generated and optimized as a transition state by ab initio, MM3, and MM4 calculations.     

Comparison of linear-scaling semiempirical methods and combined quantum mechanical/molecular mechanical methods for enzymic reactions. II. An energy decomposition analysis

QM/MM methods have been developed as a computationally feasible solution to QM simulation of chemical processes, such as enzyme-catalyzed reactions, within a more approximate MM representation of the condensed-phase environment. However, there has been no independent method for checking the quality of this representation, especially for highly nonisotropic protein environments such as those surrounding enzyme active sites. Hence, the validity of QM/MM methods is largely untested. Here we use the possibility of performing all-QM calculations at the semiempirical PM3 level with a linear-scaling method (MOZYME) to assess the performance of a QM/MM method (PM3/AMBER94 force field). Using two model pathways for the hydride-ion transfer reaction of the enzyme dihydrofolate reductase studied previously (Titmuss et al., Chem Phys Lett 2000, 320, 169-176), we have analyzed the reaction energy contributions (QM, QM/MM, and MM) from the QM/MM results and compared them with analogous-region components calculated via an energy partitioning scheme implemented into MOZYME. This analysis further divided the MOZYME components into Coulomb, resonance and exchange energy terms. For the model in which the MM coordinates are kept fixed during the reaction, we find that the MOZYME and QM/MM total energy profiles agree very well, but that there are significant differences in the energy components. Most significantly there is a large change (approximately 16 kcal/mol) in the MOZYME MM component due to polarization of the MM region surrounding the active site, and which arises mostly from MM atoms close to (<10 A) the active-site QM region, which is not modelled explicitly by our QM/MM method. However, for the model where the MM coordinates are allowed to vary during the reaction, we find large differences in the MOZYME and QM/MM total energy profiles, with a discrepancy of 52 kcal/mol between the relative reaction (product-reactant) energies. This is largely due to a difference in the MM energies of 58 kcal/mol, of which we can attribute approximately 40 kcal/mol to geometry effects in the MM region and the remainder, as before, to MM region polarization. Contrary to the fixed-geometry model, there is no correlation of the MM energy changes with distance from the QM region, nor are they contributed by only a few residues. Overall, the results suggest that merely extending the size of the QM region in the QM/MM calculation is not a universal solution to the MOZYME- and QM/MM-method differences. They also suggest that attaching physical significance to MOZYME Coulomb, resonance and exchange components is problematic. Although we conclude that it would be possible to reparameterize the QM/MM force field to reproduce MOZYME energies, a better way to account for both the effects of the protein environment and known deficiencies in semiempirical methods would be to parameterize the force field based on data from DFT or ab initio QM linear-scaling calculations. Such a force field could be used efficiently in MD simulations to calculate free energies.     

The thermochemistry of 2,4-pentanedione revisited: observance of a nonzero enthalpy of mixing between tautomers and its effects on enthalpies of formation

The enthalpies of formation of pure liquid and gas-phase (Z)-4-hydroxy-3-penten-2-one and 2,4-pentanedione are examined in the light of some more recent NMR studies on the enthalpy differences between gas-phase enthalpies of the two tautomers. Correlation gas chromatography experiments are used to evaluate the vaporization enthalpies of the pure tautomers. Values of (51.2 +/- 2.2) and (50.8 +/- 0.6) kJ.mol(-1) are measured for pure 2,4-pentanedione and (Z)-4-hydroxy-3-penten-2-one, respectively. The value of (50.8 +/- 0.6) kJ.mol(-1) can be contrasted to a value of (43.2 +/- 0.2) kJ.mol(-1) calculated for pure (Z)-4-hydroxy-3-penten-2-one when the vaporization enthalpy is measured in a mixture of tautomers. The difference is attributed to an endothermic enthalpy of mixing that destabilizes the mixture relative to the pure components. Calculation of new enthalpies of formation for (Z)-4-hydroxy-3-penten-2-one and 2,4-pentanedione in both the gas, Delta(f)H degrees (m)(g) = (-378.2 +/- 1.2) and (-358.9 +/- 2.5) kJ.mol(-1), respectively, and liquid phases, Delta(f)H degrees (m)(l) = (-429.0 +/- 1.0) and (-410.1 +/- 1.2) kJ.mol(-1), respectively, results in enthalpy differences between the two tautomers both in the liquid and gas phases that are identical within experimental error, and in excellent agreement with recent gas-phase NMR studies.     

Alcohols,ethers, carbohydrates,and related compounds. IV. Carbohydrates

Ab initio calculations [B3LYP/6-311++G(2d,2p)] have been carried out on 84 conformations of 12 different sugars (hexoses), in both pyranose and furanose forms, with the idea of generating a data base for carbohydrate structural energies that may be used for developing the predictive value of molecular mechanics calculations for carbohydrates. The average value for the apparent gas phase anomeric effect for a series of 31 pairs of pyranose conformations was found to be 1.83 kcal/mol (vs. 2.67 kcal/mol with a smaller basis set used in earlier calculations). In developing MM4 to reproduce these data, it was necessary first to have good energies for simple alcohols and ethers, together with an adequate treatment of hydrogen bonding, and then to include the anomeric effect, and the ethylene glycol type system, as was previously recognized. It was also found that the so-called delta-2 effect, long recognized in carbohydrates, must be explicitly included, in order to obtain acceptable results. When a force field that included all of these items as developed from the small molecules based on the MM4 hydrocarbon force field was applied without any parameter adjustment to the set of hexopyranose and furanose conformations mentioned earlier, the E(beta) - E(alpha) was found to have an average value of 1.88 kcal/mol, versus 1.74 for the quantum calculations. The signed average and RMS deviations of the MM4 from the QM results were +0.15 and 0.87 kcal/mol.     

Arene-cation interactions of positive quadrupole moment aromatics and arene-anion interactions of negative quadrupole moment aromatics

Intermolecular interactions involving aromatic pi-electron density are widely believed to be governed by the aromatic molecular quadrupole moment, Theta(zz). Arene-cation binding is believed to occur primarily with negative Theta(zz) aromatics, and arene-anion binding is believed to occur largely with positive Theta(zz) aromatics. We have performed quantum mechanical computations that show the cation binding of positive Theta(zz) aromatics and the anion binding of negative Theta(zz) aromatics is quite common in the gas phase. The pi-electron density of hexafluorobenzene, the prototypical positive Theta(zz ) aromatic (experimental Theta(zz) = 9.5 +/- 0.5 DA), has a Li+ binding enthalpy of -4.37 kcal/mol at the MP2(full)/6-311G**level of theory. The RHF/6-311G** calculated Theta(zz) value of 1,4-dicyanobenzene is +11.81 DA, yet it has an MP2(full)/6-311G** Li+ binding enthalpy of -12.65 kcal/mol and a Na+ binding enthalpy of -3.72 kcal/mol. The pi-electron density of benzene, the prototypical negative Theta(zz) aromatic (experimental Theta(zz) = -8.7 +/- 0.5 DA), has a F- binding enthalpy of -5.51 kcal/mol. The RHF/6-311G** calculated Theta(zz) of C6H2I4 is -10.45 DA, yet it has an MP2(full)/6-311++G** calculated F- binding enthalpy of -20.13 kcal/mol. Our results show that as the aromatic Theta(zz) value increases the cation binding enthalpy decreases; a plot of cation binding enthalpies versus aromatic Theta(zz) gives a line of best of fit with R2 = 0.778. No such correlation exists between the aromatic Theta(zz) value and the anion binding enthalpy; the line of best fit has R2 = 0.297. These results are discussed in terms of electrostatic and polarizability contributions to the overall binding enthalpies.     

A multilayered representation,quantum mechanical and molecular mechanics study of the CH3F + OH− reaction in water

The bimolecular nucleophilic substitution (S_N2) reaction of CH₃F + OH^? in aqueous solution was investigated using a combined quantum mechanical and molecular mechanics approach. Reactant complex, transition state, and product complex along the reaction pathway were analyzed in water. The potentials of mean force were calculated using a multilayered representation with the DFT and CCSD(T) level of theory for the reactive region. The obtained free energy activation barrier for this reaction at the CCSD(T)/MM representation is 18.3 kcal/mol which agrees well with the experimental value at ～21.6 kcal/mol. Both the solvation effect and solute polarization effect play key roles on raising the activation barrier height in aqueous solution, with the former raising the barrier height by 3.1 kcal/mol, the latter 1.5 kcal/mol. © 2013 Wiley Periodicals, Inc.     

Modelling of carbohydrate–aromatic interactions: <Emphasis Type="Italic">ab initio </Emphasis>energetics and force field performance

Summary Aromatic amino acid residues are often present in carbohydrate-binding sites of proteins. These binding sites are characterized by a placement of a carbohydrate moiety in a stacking orientation to an aromatic ring. This arrangement is an example of CH/π interactions. Ab initio interaction energies for 20 carbohydrate–aromatic complexes taken from 6 selected ultra-high resolution X-ray structures of glycosidases and carbohydrate-binding proteins were calculated. All interaction energies of a pyranose moiety with a side chain of an aromatic residue were calculated as attractive with interaction energy ranging from −2.8 to −12.3 kcal/mol as calculated at the MP2/6-311+G(d) level. Strong attractive interactions were observed for a wide range of orientations of carbohydrate and aromatic ring as present in selected X-ray structures. The most attractive interaction was associated with apparent combination of CH/π interactions and classical H-bonds. The failure of Hartree–Fock method (interaction energies from +1.0 to −6.9 kcal/mol) can be explained by a dispersion nature of a majority of the studied complexes. We also present a comparison of interaction energies calculated at the MP2 level with those calculated using molecular mechanics force fields (OPLS, GROMOS, CSFF/CHARMM, CHEAT/CHARMM, Glycam/AMBER, MM2 and MM3). For a majority of force fields there was a strong correlation with MP2 values. RMSD between MP2 and force field values were 1.0 for CSFF/CHARMM, 1.2 for Glycam/AMBER, 1.2 for GROMOS, 1.3 for MM3, 1.4 for MM2, 1.5 for OPLS and to 2.3 for CHEAT/CHARMM (in kcal/mol). These results show that molecular mechanics approximates interaction energies very well and support an application of molecular mechanics methods in the area of glycochemistry and glycobiology.     

Thermochemistry of the gaseous vanadium chlorides VCl, VCl2, VCl3, and VCl4

Gaseous equilibria in the V-Ag-Cl system were studied at elevated temperatures by effusion-beam mass spectrometry, where the pertinent species were generated by reaction of Cl 2(g) with V + Ag granules in the effusion cell source. Reaction enthalpies were derived from the equilibrium data, and the standard enthalpies of formation at 298 K of gaseous VCl, VCl2, and VCl3 were found to be +49.7, -34.8, and -85.6 kcal mol(-1), respectively. The corresponding bond dissociation energies at 298 K are D(V-Cl) = 102.9 kcal, D(ClV-Cl) = 113.5 kcal, D(Cl2V-Cl) = 79.8 kcal, and D(Cl3V-Cl) = 69.5 kcal. From these data, the dissociation energy D degrees 0(VCl) = 101.9 kcal mol(-1) or 4.42 eV is obtained. An alternate value, Delta(f)H(o)298(VCl 3,g) = -87.0 kcal mol (-1) was derived from third-law analysis of literature sublimation data for VCl3(s). In addition, literature thermochemical data on VCl4(g) were re-evaluated, leading to Delta(f)H(o)298 = -126.1 kcal mol (-1). The results are compared with various estimates in the literature.     

Effect of geminal substitution at silicon on 1-sila- and 1,3-disilacyclobutanes' strain energies, their 2+2 cycloreversion enthalpies, and Si=C pi-bond energies in silenes.

To evaluate the effect of geminal substitution at silicon on 1-sila- and 1,3-disilacyclobutanes' strain energies, their 2+2 cycloreversion enthalpies, and Si=C pi-bond energies in silenes, an ab initio MO study of silenes, R2Si=CH2 (1), 1-silacyclobutanes, cyclo-R2Si(CH2)3 (2), and 1,3-disilacyclobutanes, cyclo-(R2SiCH2)2 (3), was performed using the level of theory denoted MP4/TZ(d)//MP2/6-31G(d) (TZ means the 6-311G(d) basis set for elements of the second period and hydrogen, and the McLean-Chandler (12s,9p)/[6s,5p](d) basis set for the third period elements). In the series R = H, CH3, SiH3, CH3O, NH2, Cl, F, the growth of the reaction enthalpies and strain energies is proportional to the substituents' electronegativities. 2+2 cycloreversion of 2 is endothermic by 40.6-63.1 kcal/mol, whereas that of 3 is endothermic by 72.7-114.2 kcal/mol. On going from a silicon to a fluorine substituent at the sp2-hybridized silicon atom, the pi-bond energy in 1 weakens by 11.3 kcal/mol, and the Si=C bond length shortens by 0.053 A. The effect of substituents' electronegativities at the double-bonded silicon atom in silenes is formulated as follows: the higher electronegativity, the shorter and the weaker the Si=C pi-bond. The latter is rationalized in terms of more strained geometry resulting from the energetic cost for planarizing the R2SiC moiety. The enthalpies of the ring-opening reaction are 68.0-80.1 kcal/mol (a cleavage of the Si-C bond in 3), 65.0-76.4 kcal/mol (a cleavage of the Si-C bond in 2), and 58.0-64.9 kcal/mol (a cleavage of the C-C bond in 2). The pronounced difference in the enthalpies of 2+2 cycloreversion of 1-sila- and 1,3-disilacyclobutanes is mainly due to the difference in the enthalpies of diradicals' decomposition. The decomposition of diradicals resulting from a cleavage of C-C and Si-C bonds in 2 is exothermic by 24.3-3.3 kcal/mol (apart from the difluoro derivative which is endothermic by 5.1 kcal/mol) and 27.0-13.3 kcal/mol, respectively. The decomposition of a 1,4-diradical resulting from ring opening of 3, apart from the disilyl derivative, is the endothermic process for which the enthalpy varies from 10.6 to 40.4 kcal/mol.     

A semiempirical quantum mechanical study of cationically catalyzed homopolymerization and copolymerization of vinyl ethers and epoxides

The purpose of this study was to use the semiempirical quantum mechanical computational method, AM1, to investigate vinyl ether cationic homopolymerization, epoxide homopolymerization, and copolymerization of selected vinyl ethers with a model epoxide (cyclohexene oxide). Homopolymerization studies of 19 vinyl ethers showed that activation enthalpies ranged between 0.0 and 15 kcal/mol, and that the enthalpies of reaction for homopolymerization were nearly all exothermic. Homopolymerization of three epoxides predicted low activation enthalpies, some of which were virtually activationless. All ring-opening epoxide polymerizations were exothermic. Copolymerization of three vinyl ethers with cyclohexene oxide gave activation enthalpies that varied between 2.7 and 4.0 kcal/mol, and the enthalpies of reaction for copolymerization were all exothermic.