Solvation effects on equilibria: Triazoles and N-methyl piperidinol

We have performed calculations of the solvation effects on a number of equilibrium constants in water using a recently proposed hybrid quantum classical scheme in which the liquid environment is modelled using classical solvent molecules and the solute electronic structure is computed using modern quantum chemical methods. The liquid phase space is sampled from a fully classical simulation. We find that solvation effects on both triazole tautomeric equilibrium constants and piperidinol conformational equilibrium constants can be interpreted in terms of subtle differences in the local environment which can be seen in probability densities and radial distribution functions. Lower level calculations were performed for comparison and we conclude that the solvation thermodynamics can be predicted from a good classical model of solvent and solute molecules, but the implicit models that we tried are less successful.

Thermochemistry of Ionic Liquid-Catalyzed Reactions. Experimental and Theoretical Study of Chemical Equilibria of Isomerization and Transalkylation of tert-Butylbenzenes

The chemical equilibrium of mutual interconversions of tert-butylbenzenes was studied in the temperature range 286 to 423 K using chloroaluminate ionic liquids as a catalyst. Enthalpies of five reactions of isomerization and transalkylation of tert-butylbenzenes were obtained from temperature dependences of the corresponding equilibrium constants in the liquid phase. Molar enthalpies of vaporization of methyl-tert-butylbenzenes and 1,4-ditert-butylbenzene were obtained by the transpiration method and were used for a recalculation of enthalpies of reactions and equilibrium constants into the gaseous phase. Using these experimental results, ab initio methods (B3LYP and G3MP2) have been tested for prediction thermodynamic functions of the five reactions under study successfully. Thermochemical investigations of tert-butyl benzenes available in the literature combined with experimental results have helped to resolve contradictions in the available thermochemical data for tert-butylbenzene and to recommend consistent and reliable enthalpies of formation for this compound in the liquid and the gaseous state.

Thermochemistry of Ionic Liquid-Catalyzed Reactions: Theoretical and Experimental Study of the Beckmann Rearrangement-Kinetic or Thermodynamic Control?

A thermodynamic analysis of the experimental conditions of the Beckmann rearrangement reaction of oximes into amides has been undertaken to examine whether the reaction is under thermodynamic or kinetic control. To answer this question, the thermodynamic properties of the typical Beckmann rearrangement reactions in the ideal gaseous state-cyclohexanone oxime to caprolactam and 2-butanone oxime to N-methylpropanarnide-were studied by using the quantum mechanical method. Gibbs energy and equilibrium constants of the Beckmann rearrangement have been assessed in the gaseous and the liquid phases. Results of the thermodynamic analysis have shown that Beckmann rearrange ments are kinetically controlled. Thus, a search for possible active ionic liquid based catalysts for the mild reaction conditions has been performed.

Compositional data-driven alternatives to single geochemical component maps

Conventional practice in Regional Geochemistry includes as a final step of any geochemical campaign the generation of a series of maps, to show the spatial distribution of each of the components considered. Such maps, though necessary, do not comply with the compositional, relative nature of the data, which unfortunately make any conclusion based on them sensitive
to spurious correlation problems. This is one of the reasons why these maps are never interpreted isolated. This contribution aims at gathering a series of statistical methods to produce individual maps of multiplicative combinations of components (logcontrasts), much in the flavor of equilibrium constants, which are designed on purpose to capture certain aspects of the data.
We distinguish between supervised and unsupervised methods, where the first require an external, non-compositional variable (besides the compositional geochemical information) available in an analogous training set. This external variable can be a quantity (soil density, collocated magnetics, collocated ratio of Th/U spectral gamma counts, proportion of clay particle fraction, etc) or a category (rock type, land use type, etc). In the supervised methods, a regression-like model between the external variable and the geochemical composition is derived in the training set, and then this model is mapped on the whole region. This case is illustrated with the Tellus dataset, covering Northern Ireland at a density of 1 soil sample per 2 square km, where we map the presence of blanket peat and the underlying geology. The unsupervised methods considered include principal components and principal balances
(Pawlowsky-Glahn et al., CoDaWork2013), i.e. logcontrasts of the data that are devised to capture very large variability or else be quasi-constant. Using the Tellus dataset again, it is found that geological features are highlighted by the quasi-constant ratios Hf/Nb and their ratio against SiO2; Rb/K2O and Zr/Na2O and the balance between these two groups of two variables; the balance of Al2O3 and TiO2 vs. MgO; or the balance of Cr, Ni and Co vs. V and Fe2O3. The largest variability appears to be related to the presence/absence of peat.

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.

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.

Equilibrium structure of erbium-oxygen complexes in crystalline silicon

The equilibrium structure of ErOn (nless than or equal to6) complexes in crystalline silicon has been investigated by density-functional computations. Two different geometries have been considered, corresponding to the substitutional and tetrahedral interstitial site for erbium. All atomic coordinates have been optimized by Car-Parrinello molecular dynamics. The resulting structures have low symmetry, with E-O distances of similar to2.35 Angstrom. The substitutional site is the most stable one for nless than or equal to2, while the tetrahedral interstitial is favored for n>2.

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.

Structure and kinetics of phosphonopyruvate hydrolase from Variovorax sp. Pal2: New insight into the divergence of catalysis within the PEP mutase/isocitrate lyase superfamily .

Phosphonopyruvate (P-pyr) hydrolase (PPH), a member of the phosphoenolpyruvate (PEP) mutase/isocitrate lyase (PEPM/ICL) superfamily, hydrolyzes P-pyr and shares the highest sequence identity and functional similarity with PEPM. Recombinant PPH from Variovorax sp. Pal2 was expressed in Escherichia coli and purified to homogeneity. Analytical gel filtration indicated that the protein exists in solution predominantly as a tetramer. The PPH pH rate profile indicates maximal activity over a broad pH range.The steady-state kinetic constants determined for a rapid equilibrium ordered kinetic mechanism with Mg+2 binding first (Kd =140 ± 40 M), are kcat = 105 ± 2 s-1 and P-pyr Km = 5 ± 1 M. PEP (slow substrate kcat = 2 × 10-4 s-1), oxalate, and sulfopyruvate are competitive inhibitors with Ki values of 2.0 ± 0.1 mM, 17 ± 1 M, and 210 ± 10 M, respectively. Three PPH crystal structures have been determined, that of a ligand-free enzyme, the enzyme bound to Mg2+ and oxalate (inhibitor), and the enzyme bound to Mg2+ and P-pyr (substrate). The complex with the inhibitor was obtained by cocrystallization, whereas that with the substrate was obtained by briefly soaking crystals of the ligand-free enzyme with P-pyr prior to flash cooling. The PPH structure resembles that of the other members of the PEPM/ICL superfamily and is most similar to the functionally related enzyme, PEPM. Each monomer of the dimer of dimers exhibits an (/)8 barrel fold with the eighth helix swapped between two molecules of the dimer. Both P-pyr and oxalate are anchored to the active site by Mg2+. The loop capping the active site is disordered in all three structures, in contrast to PEPM, where the equivalent loop adopts an open or disordered conformation in the unbound state but sequesters the inhibitor from solvent in the bound state. Crystal packing may have favored the open conformation of PPH even when the enzyme was cocrystallized with the oxalate inhibitor. Structure alignment of PPH with other superfamily members revealed two pairs of invariant or conservatively replaced residues that anchor the flexible gating loop. The proposed PPH catalytic mechanism is analogous to that of PEPM but includes activation of a water nucleophile with the loop Thr118 residue.