Limited ionicity in poor protic ionic liquids: Association Gibbs energies.

Protic ionic liquids (PILs), made from anhydrous mixtures of Bronsted acids HA and bases B (HA + B → BH + A), occasionally suffer from limited ionicity. In cases of "poor" PILs (<10% ionicity, e.g.

Semicontinuum (Cluster-Continuum) Modeling of Acid-Catalyzed Aqueous Reactions: Alkene Hydration.

An effort is made to reduce the errors of continuum solvation models (CSMs) with semicontinuum modeling to achieve 3 kcal mol agreement with experiment for acid-catalysis activation Gibbs energies. First, two underappreciated CSM issues are reviewed: errors in the CSM solvation Gibbs energies grow beyond 5 kcal mol (i) as ions are made smaller and (ii) as water clusters grow larger. Second, the computational reproduction of the known Gibbs energies (Δ and Δ) of the paradigmatic reaction ethene + HO + HO → TS → ethanol + HO is attempted.

Improving the Yield and Rate of Acid-Catalyzed Deconstruction of Lignin by Mechanochemical Activation.

Lignin is a potential biomass feedstock from plant material, but it is particularly difficult to economically process. Inspired by recent ball-milling results, state-of-the-art quantum mechanochemistry calculations have been performed to isolate and probe the purely mechanochemical stretching effect alone upon acid-catalyzed deconstruction of lignin. Effects upon cleavage of several exemplary simple ethers are examined first, and with low stretching force they all are predicted to cleave substantially faster, allowing for use of milder acids and lower temperatures.

The origin of the conductivity maximum in molten salts. III. Zinc halides.

In a continuing effort to master the reasons for conductivity maxima vs temperature in semicovalent molten halides, the structure and some transport properties of molten zinc halide are examined with ab initio molecular dynamics. Molten zinc halides are a special class of molten salts, being extremely viscous near their melting point (with a glassy state below it) and low electrical conductivity, and since they are also known (ZnI) or predicted (ZnBr and ZnCl) to exhibit conductivity maxima, they would be useful additional cases to probe, in case the reasons for their maxima are unique. Strong attractive forces in ZnX result in tight tetrahedral coordination, and the known mixture of edge-sharing vs corner-sharing ZnX tetrahedra is observed.

Mechanical Activation Drastically Accelerates Amide Bond Hydrolysis, Matching Enzyme Activity.

Amide bonds, which include peptide bonds connecting amino acids in proteins and polypeptides, give proteins and synthetic polyamides their enormous strength. Although proteins and polyamides sustain mechanical force in nature and technology, how forces affect amide and peptide bond stability is still unknown. Using single-molecule force spectroscopy, we discover that forces of only a few hundred pN accelerate amide hydrolysis 10 -fold, an acceleration hitherto only known from proteolytic enzymes.

The origin of the conductivity maximum vs. mixing ratio in pyridine/acetic acid and water/acetic acid.

Explanations are provided for the first time for the historically known locations of electrical conductivity maxima versus mixing ratio (mole fraction of acid, x) in mixtures of (i) acetic acid with water and (ii) acetic acid with pyridine. To resolve the question for the second system, density-functional-based molecular dynamic simulations were performed, at 1:1, 1:2, 1:3, 1:5, and 1:15 mixing ratios, to gain vital information about speciation. In a zeroth-order picture, the degree of ionization (and hence conductivity) would be maximal at x = 0.

The origin of the conductivity maximum in molten salts. II. SnCl2 and HgBr2.

The phenomenon of electrical conductivity maxima of molten salts versus temperature during orthobaric (closed-vessel) conditions is further examined via ab initio simulations. Previously, in a study of molten BiCl3, a new theory was offered in which the conductivity falloff at high temperatures is due not to traditional ion association, but to a rise in the activation energy for atomic ions hopping from counterion to counterion. Here this theory is further tested on two more inorganic melts which exhibit conductivity maxima: another high-conducting melt (SnCl2, σmax = 2.

The Carbocation Rearrangement Mechanism, Clarified.

The role of protonated cyclopropane (PCP(+)) structures in carbocation rearrangement is a decades-old topic that continues to confound. Here, quantum-chemical computations (PBE molecular dynamics, PBE and CCSD optimizations, CCSD(T) energies) are used to resolve the issue. PCP(+) intermediates are neither edge-protonated nor corner-protonated (normally) but possess "closed" structures mesomeric between these two.

Semicontinuum Solvation Modeling Improves Predictions of Carbamate Stability in the CO2 + Aqueous Amine Reaction.

Quantum chemistry computations with a semicontinuum (cluster + continuum) solvation model have been used to cure long-standing misprediction of aqueous carbamate anion energies in the industrially important CO2 + aqueous amine reaction. Previous errors of over 10 kcal mol(-1) are revealed. Activation energies were also estimated with semicontinuum modeling, and a refined discussion of the competing hypothetical mechanisms for CO2 + monoethanolamine (MEA) is presented.

Challenges in predicting ΔrxnG in solution: the mechanism of ether-catalyzed hydroboration of alkenes.

Ab initio (coupled-cluster and density-functional) calculations of Gibbs reaction energies in solution, with new entropy-of-solvation damping terms, were performed for the ether-catalyzed hydroboration of alkenes. The goal was to test the accuracy of continuum-solvation models for reactions of neutral species in nonaqueous solvents, and the hope was to achieve an accuracy sufficient to address the mechanism in the "Pasto case": B2H6 + alkene in THF solvent. Brown's SN2/SN1 "dissociative" mechanism, of SN2 formation of borane-ether adducts followed by SN1 alkene attack, was at odds with Pasto's original SN2/SN2 hypothesis, and while Brown could prove his mechanism for a variety of cases, he could not perform the experimental test with THF adducts in THF solvent, where the higher THF concentrations might favor an SN2 second step.

Molecular Dynamics Simulations of Proposed Intermediates in the CO2 + Aqueous Amine Reaction.

Ab initio molecular dynamics simulations of up to 210 ps have been performed on various aqueous intermediates postulated in the CO2 + amine reaction, important for CO2 capture. Observations of spontaneous deprotonation of aqueous carbamate zwitterions R1R2NHCOO(±) by bulk water (instead of additional amine, or via umbrella sampling) are reported apparently for the first time. Carbamic acid structures R1R2NCOOH were observed in some simulations, arising from zwitterions not via classical 1,3-H-shifts but via Grotthuss-style multiple-H(+) transfer pathways that involve bulk H2O and require carbamate anions R1R2NCOO(-) as an intermediate stage along the way.

The origin of the conductivity maximum in molten salts. I. Bismuth chloride.

A new theory is presented to explain the conductivity maxima of molten salts (versus temperature and pressure). In the new theory, conductivity is due to ions hopping from counterion to counterion, and its temperature dependence can be explained with an ordinary Arrhenius equation in which the frequency prefactor A (for hopping opportunities) and activation energy E(a) (for hopping) are density dependent. The conductivity maximum is due to competing effects: as density decreases, the frequency of opportunities for hopping increases, but the probability that an opportunity is successfully hopped decreases due to rising E(a) caused by the increased hopping distance.

The mechanism of permanganate oxidation of sulfides and sulfoxides.

Coupled-cluster (CCSD) and density functional computations are used to investigate historically competing mechanisms for the permanganate oxidation of sulfides and sulfoxides. The calculations all lead to a mechanism of 1,3-dipolar cycloaddition of permanganate, as opposed to historical mechanisms of attack of the sulfur atom by one O or by Mn. Such a mechanism, reminiscent of ozonolysis, may prevail in most permanganate oxidations.

Supra-supra, supra-antara, and stepwise-diradical pathways for an observed 16-electron double-[4 + 4] cycloaddition within metal-templated dialkyne dimers (PtX2)2(μ-R2PCCCCPR2)2.

Quantum chemistry calculations are used to provide insight into the cycloaddition of two dialkyne chains in initially monocyclic organoplatinum dimers of the type (PtX(2))(2)(μ-R(2)PC(4)PR(2))(2), where X = Cl or Me and R = Ph or Me. Previous experimental studies showed that the cycloaddition occurs with {X, R} = {Cl, Ph} but not {Me, Ph}. Two concerted pericyclic paths, a D(2h)-symmetry double-[π4s+π4s] "Hückel path" and a D(2)-symmetry double-[π4s+π4a] "Möbius path", were explored via orbital energy correlation diagrams (OECDs) computed using a singly occupied molecular orbital technique developed earlier.

Computational Study of Tungsten(II)-Catalyzed Rearrangements of Norbornadiene.

Rearrangements of norbornadiene (NBD, C7H8) to various alkylidenes, via a hypothetical 7-coordinate tungsten(II) complex W(CO)3I2(NBD), were studied using density-functional theory computations. An extensive search for intermediates and transition states of rearrangement was made. The theoretical method (basis sets and level of DFT) used was justified by new benchmark studies which compare optimized structural parameters to those from crystal structures of several different tungsten complexes.

Photochemistry studied with ab initio orbital-correlation and state-correlation plots: Classic cyclobutene ring opening, and the reaction of N2 with photoexcited O2.

Pericyclic reaction theory arose from ideas presented in 1965, based on orbital-energy correlation diagrams (Woodward and Hoffmann) and state-energy correlation diagrams (Longuet-Higgins and Abrahamson). Here we have used ab initio complete-active-space self-consistent field (CASSCF) calculations to generate such diagrams. First we present diagrams for the classic case of cyclobutene ring opening, to demonstrate agreement between the CASSCF results and the classic diagrams of both Woodward/Hoffmann and Longuet-Higgins/Abrahamson.

Improved results for the excited states of nitric oxide, including the B/C avoided crossing.

The potential energy surfaces of ten electronic states of nitric oxide (NO) have been reexamined computationally, with state energies calculated using ab initio multireference methods. Our wave function expansions of 10x10(6) configurations improve upon the results of de Vivie and Peyerimhoff [J. Chem.

Photodissociation dynamics of the NO dimer. I. Theoretical overview of the ultraviolet singlet excited states.

Molecular orbital theory and calculations are used to describe the ultraviolet singlet excited states of NO dimer. Qualitatively, we derive and catalog the dimer states by correlating them with monomer states, and provide illustrative complete active space self-consistent field calculations. Quantitatively, we provide computational estimates of vertical transition energies and absorption intensities with multireference configuration interaction and equations-of-motion coupled-cluster methods, and examine an important avoided crossing between a Rydberg and a valence state along the intermonomer and intramonomer stretching coordinates.

A theoretical comparison of Lewis acid vs bronsted acid catalysis for n-hexane --> propane + propene.

Cracking of an all-trans n-alkane, via idealized Lewis acid and Bronsted acid catalysis, was examined using density functional theory. Optimized geometries and transitions states were determined for catalyst-reactant complexes, using AlCl3 and HCl.AlCl3 as the Lewis and Bronsted acids.

Competing isomerizations: a combined experimental/theoretical study of phenylpentenone isomerism.

The possible competition of Z/E versus hydrogen-shift isomerization in (E)-5-phenyl-3-penten-2-one (E-1) and (E)-5-phenyl-4-penten-2-one (E-2) was studied, both experimentally and theoretically. Iodine-catalyzed isomerization experiments and computational modeling studies show that the equilibrated system consists predominantly of E-1 and E-2, with E-2 in moderate excess, and with no detectable amounts of the Z (cis) diastereoisomers. Density functional theory (DFT) calculations corroborated the free energy difference (Delta(r) and Delta(r) were -0.

Length and substituent-scrambling energies of parent and halogen-substituted conjugated polyynes.

Conjugated polyynes are a class of species of diverse and increasing interest. Length-scrambling and substituent scrambling reaction energies were examined using ab initio quantum chemistry calculations to investigate issues concerning the energetic effects of the molecular ends (substituent communication). Computations were performed for the parent, monohalogenated, and dihalogenated (F, Cl, Br, I) polyynes of up to 60 carbon atoms.

Mechanism of cis/trans equilibration of alkenes via iodine catalysis.

Iodine is commonly used to speed the equilibration of Wittig cis/trans alkene products. This study uses computational chemistry to study the catalyzed isomerization mechanism in detail for seven different examples of 1,2-disubstituted alkenes. We find that the iodo intermediates of the conventional three-step reaction path are weakly stable, bound by less than 7 kJ mol(-1) in five cases and nonexistent in the other two.