Proton Abstraction from DME n···X+ by OH-, O2-, and XO2-, for X = Li, Na, and K: Implications for Li-O2 Batteries.

The abstraction of a proton by OH-, O2-, and XO2- from DME n···X+, where X is Li, Na, or K, is studied using density functional theory. Both the gas phase and the solution phase are studied. In general, when explicit solvent molecules are added, the difference between the gas-phase and solution results becomes rather small. While the DME n···X+ binding energies differ significantly for various alkali cations, the reaction energies and transition-state energies are far less sensitive to the choice of an alkali cation. XO2- has a lower barrier height than OH-, which, in turn, has a lower barrier height than O2-. The reaction energies follow the same trends.

Products of Criegee intermediate reactions with NO2: experimental measurements and tropospheric implications.

The reactions of Criegee intermediates with NO2 have been proposed as a potentially significant source of the important nighttime oxidant NO3, particularly in urban environments where concentrations of ozone, alkenes and NOx are high. However, previous efforts to characterize the yield of NO3 from these reactions have been inconclusive, with many studies failing to detect NO3. In the present work, the reactions of formaldehyde oxide (CH2OO) and acetaldehyde oxide (CH3CHOO) with NO2 are revisited to further explore the product formation over a pressure range of 4-40 Torr. NO3 is not observed; however, temporally resolved and [NO2]-dependent signal is observed at the mass of the Criegee-NO2 adduct for both formaldehyde- and acetaldehyde-oxide systems, and the structure of this adduct is explored through ab initio calculations. The atmospheric implications of the title reaction are investigated through global modelling.

Pressure-Dependent Competition among Reaction Pathways from First- and Second-O2 Additions in the Low-Temperature Oxidation of Tetrahydrofuran.

We report a combined experimental and quantum chemistry study of the initial reactions in low-temperature oxidation of tetrahydrofuran (THF). Using synchrotron-based time-resolved VUV photoionization mass spectrometry, we probe numerous transient intermediates and products at P = 10-2000 Torr and T = 400-700 K. A key reaction sequence, revealed by our experiments, is the conversion of THF-yl peroxy to hydroperoxy-THF-yl radicals (QOOH), followed by a second O2 addition and subsequent decomposition to dihydrofuranyl hydroperoxide + HO2 or to γ-butyrolactone hydroperoxide + OH. The competition between these two pathways affects the degree of radical chain-branching and is likely of central importance in modeling the autoignition of THF. We interpret our data with the aid of quantum chemical calculations of the THF-yl + O2 and QOOH + O2 potential energy surfaces. On the basis of our results, we propose a simplified THF oxidation mechanism below 700 K, which involves the competition among unimolecular decomposition and oxidation pathways of QOOH.

Carbon radicals. Direct observation and kinetics of a hydroperoxyalkyl radical (QOOH).

Oxidation of organic compounds in combustion and in Earth's troposphere is mediated by reactive species formed by the addition of molecular oxygen (O2) to organic radicals. Among the most crucial and elusive of these intermediates are hydroperoxyalkyl radicals, often denoted "QOOH." These species and their reactions with O2 are responsible for the radical chain branching that sustains autoignition and are implicated in tropospheric autoxidation that can form low-volatility, highly oxygenated organic aerosol precursors. We report direct observation and kinetics measurements of a QOOH intermediate in the oxidation of 1,3-cycloheptadiene, a molecule that offers insight into both resonance-stabilized and nonstabilized radical intermediates. The results establish that resonance stabilization dramatically changes QOOH reactivity and, hence, that oxidation of unsaturated organics can produce exceptionally long-lived QOOH intermediates.

Biofuel combustion. Energetics and kinetics of hydrogen abstraction from carbon-1 in n-butanol by the hydroperoxyl radical calculated by coupled cluster and density functional theories and multistructural variational transition-state theory with multidimensional tunneling.

Multistructural canonical variational transition-state theory with small-curvature multidimensional tunneling (MS-CVT/SCT) is employed to calculate thermal rate constants for hydrogen-atom abstraction from carbon-1 of n-butanol by the hydroperoxyl radical over the temperature range 250-2000 K. The M08-SO hybrid meta-GGA density functional was validated against CCSD(T)-F12a explicitly correlated wave function calculations with the jul-cc-pVTZ basis set. It was then used to compute the properties of all stationary points and the energies and Hessians of a few nonstationary points along the reaction path, which were then used to generate a potential energy surface by the multiconfiguration Shepard interpolation (MCSI) method. The internal rotations in the transition state for this reaction (like those in the reactant alcohol) are strongly coupled to each other and generate multiple stable conformations, which make important contributions to the partition functions. It is shown that neglecting to account for the multiple-structure effects and torsional potential anharmonicity effects that arise from the torsional modes would lead to order-of-magnitude errors in the calculated rate constants at temperatures of interest in combustion.

Multistructural variational transition state theory: kinetics of the hydrogen abstraction from carbon-2 of 2-methyl-1-propanol by hydroperoxyl radical including all structures and torsional anharmonicity.

We calculated the forward and reverse rate constants of the hydrogen abstraction reaction from carbon-2 of 2-methyl-1-propanol by hydroperoxyl radical over the temperature range 250-2400 K by using multistructural canonical variational transition state theory (MS-CVT) including both multiple-structure and torsional potential anharmonicity effects by the multistructural torsional anharmonicity (MS-T) method. In these calculations, multidimensional tunneling (MT) probabilities used to compute the tunneling transmission coefficients were evaluated by the small-curvature tunneling (SCT) approximation. Comparison with the rate constants obtained by the single-structural harmonic oscillator (SS-HO) approximation shows that multistructural anharmonicity increases the forward rate constants for all temperatures, but the reverse rate constants are reduced for temperatures lower than 430 K and increased for higher temperatures. The neglect of multistructural torsional anharmonicity would lead to errors of factors of 1.5, 8.8, and 13 at 300, 1000, and 2400 K, respectively, for the forward reaction, and would lead to errors of factors of 0.76, 3.0, and 6.0, respectively, at these temperatures for the reverse reaction.

Thermochemistry of radicals formed by hydrogen abstraction from 1-butanol, 2-methyl-1-propanol, and butanal.

We calculate the standard state entropy, heat capacity, enthalpy, and Gibbs free energy for 13 radicals important for the combustion chemistry of biofuels. These thermochemical quantities are calculated from recently proposed methods for calculating partition functions of complex molecules by taking into account their multiple conformational structures and torsional anharmonicity. The radicals considered in this study are those obtained by hydrogen abstraction from 1-butanol, 2-methyl-1-propanol, and butanal. Electronic structure calculations for all conformers of the radicals were carried out using both density functional theory and explicitly correlated coupled cluster theory with quasipertubative inclusion of connected triple excitations. The heat capacity and entropy results are compared with sparsely available group additivity data, and trends in enthalpy and free energy as a function of radical center are discussed for the isomeric radicals.

What are the most efficient basis set strategies for correlated wave function calculations of reaction energies and barrier heights?

As electronic structure methods are being used to obtain quantitatively accurate reaction energies and barrier heights for increasingly larger systems, the choice of an efficient basis set is becoming more critical. The optimum strategy for achieving basis set convergence can depend on the way that electron correlation is treated and can take advantage of flexibility in the order in which basis functions are added. Here we study several approaches for estimating accurate reaction energies and barrier heights from post-Hartree-Fock electronic structure calculations. First and second, we evaluate methods of estimating the basis set limit of second order Møller-Plesset perturbation theory and of coupled cluster theory with single and double excitations and a quasiperturbative treatment of connected triple excitations by using explicitly correlated basis functions (in the F12a implementation) along with valence, polarization, and diffuse one-electron basis functions. Third, we test the scheme of adding a higher-order correction to MP2 results (sometimes called MP2∕CBS + ΔCCSD(T)). Finally, we evaluate the basis set requirements of these methods in light of comparisons to Weizmann-3.2, Weizmann-4, and CCSDT(2)(Q)∕CBS+CV+R results.

Multi-structural variational transition state theory: kinetics of the 1,5-hydrogen shift isomerization of the 1-butoxyl radical including all structures and torsional anharmonicity.

We investigate the statistical thermodynamics and kinetics of the 1,5-hydrogen shift isomerization reaction of the 1-butoxyl radical and its reverse isomerization. The partition functions and thermodynamic functions (entropy, enthalpy, heat capacity, and Gibbs free energy) are calculated using the multi-structural torsional (MS-T) anharmonicity method including all structures for three species (reactant, product, and transition state) involved in the reaction. The calculated thermodynamic quantities have been compared to those estimated by the empirical group additivity (GA) method. The kinetics of the unimolecular isomerization reaction was investigated using multi-structural canonical variational transition state theory (MS-CVT) including both multiple-structure and torsional (MS-T) anharmonicity effects. In these calculations, multidimensional tunneling (MT) probabilities were evaluated by the small-curvature tunneling (SCT) approximation and compared to results obtained with the zero-curvature tunneling (ZCT) approximation. The high-pressure-limit rate constants for both the forward and reverse reactions are reported as calculated by MS-CVT/MT, where MT can be ZCT or SCT. Comparison with the rate constants obtained by the single-structural harmonic oscillator (SS-HO) approximation shows the importance of anharmonicity in the rate constants of these reactions, and the effect of multi-structural anharmonicity is found to be very large. Whereas the tunneling effect increases the rate constants, the MS-T anharmonicity decreases them at all temperatures. The two effects counteract each other at temperatures 385 K and 264 K for forward and reverse reactions, respectively, and tunneling dominates at lower temperatures while MS-T anharmonicity has a larger effect at higher temperatures. The multi-structural torsional anharmonicity effect reduces the final reverse reaction rate constants by a much larger factor than it does to the forward ones as a result of the existence of more low-energy structures of the product 4-hydroxy-1-butyl radical than the reactant 1-butoxyl radical. As a consequence there is also a very large effect on the equilibrium constant. The neglect of multi-structural anharmonicity will lead to large errors in the estimation of reverse reaction rate constants.

Statistical thermodynamics of 1-butanol, 2-methyl-1-propanol, and butanal.

The purpose of the present investigation is to calculate partition functions and thermodynamic quantities, viz., entropy, enthalpy, heat capacity, and Gibbs free energies, for 1-butanol, 2-methyl-1-propanol, and butanal in the vapor phase. We employed the multi-structural (MS) anharmonicity method and electronic structure calculations including both explicitly correlated coupled cluster theory and density functional theory. The calculations are performed using all structures for each molecule and employing both the local harmonic approximation (MS-LH) and the inclusion of torsional anharmonicity (MS-T). The results obtained from the MS-T calculations are in excellent agreement with experimental data taken from the Thermodynamics Research Center data series and the CRC Handbook of Chemistry and Physics, where available. They are also compared with Benson's empirical group additivity values, where available; in most cases, the present results are more accurate than the group additivity values. In other cases, where experimental data (but not group additivity values) are available, we also obtain good agreement with experiment. This validates the accuracy of the electronic structure calculations when combined with the MS-T method for estimating the thermodynamic properties of systems with multiple torsions, and it increases our confidence in the predictions made with this method for molecules and temperatures where experimental or empirical data are not available.

Kinetics of the Hydrogen Abstraction from Carbon-3 of 1-Butanol by Hydroperoxyl Radical: Multi-Structural Variational Transition-State Calculations of a Reaction with 262 Conformations of the Transition State.

We estimated rate constants for the hydrogen abstraction from carbon-3 of 1-butanol by hydroperoxyl radical, a critically important reaction in the combustion of biofuel. We employed the recently developed multi-structural variational transition-state theory (MS-VTST), which utilizes a multifaceted dividing surface that allows us to include the contributions of multiple structures for reacting species and transition states. First, multiconfigurational Shepard interpolation-based on molecular-mechanics-guided interpolation of electronic-structure Hessian data obtained by the M08 HX/jun-cc-pVTZ electronic model chemistry-was used to obtain the portion of the potential energy surface needed for single-structure variational transition-state theory rate constants including multidimensional tunneling; then, the M08-HX/MG3S electronic model chemistry was used to calculate multi-structural torsional anharmonicity factors to complete the MS-VTST rate constant calculations. The lowest-energy structures of the transition state have strongly bent hydrogen bonds. Our results indicate that neglect of multi-structural anharmonicity would lead to errors of factors of 0.3, 46, and 171 at 200, 1000, and 2400 K for this reaction.

Practical methods for including torsional anharmonicity in thermochemical calculations on complex molecules: the internal-coordinate multi-structural approximation.

Many methods for correcting harmonic partition functions for the presence of torsional motions employ some form of one-dimensional torsional treatment to replace the harmonic contribution of a specific normal mode. However, torsions are often strongly coupled to other degrees of freedom, especially other torsions and low-frequency bending motions, and this coupling can make assigning torsions to specific normal modes problematic. Here, we present a new class of methods, called multi-structural (MS) methods, that circumvents the need for such assignments by instead adjusting the harmonic results by torsional correction factors that are determined using internal coordinates. We present three versions of the MS method: (i) MS-AS based on including all structures (AS), i.e., all conformers generated by internal rotations; (ii) MS-ASCB based on all structures augmented with explicit conformational barrier (CB) information, i.e., including explicit calculations of all barrier heights for internal-rotation barriers between the conformers; and (iii) MS-RS based on including all conformers generated from a reference structure (RS) by independent torsions. In the MS-AS scheme, one has two options for obtaining the local periodicity parameters, one based on consideration of the nearly separable limit and one based on strongly coupled torsions. The latter involves assigning the local periodicities on the basis of Voronoi volumes. The methods are illustrated with calculations for ethanol, 1-butanol, and 1-pentyl radical as well as two one-dimensional torsional potentials. The MS-AS method is particularly interesting because it does not require any information about conformational barriers or about the paths that connect the various structures.

Perspectives on Basis Sets Beautiful: Seasonal Plantings of Diffuse Basis Functions.

We present a perspective on the use of diffuse basis functions for electronic structure calculations by density functional theory and wave function theory. We especially emphasize minimally augmented basis sets and calendar basis sets. We base our conclusions on our previous experience with commonly computed quantities, such as bond energies, barrier heights, electron affinities, noncovalent (van der Waals and hydrogen bond) interaction energies, and ionization potentials, on Stephens et al.'s results for optical rotation and on our own new calculations (presented here) of polarizabilities and of potential energy curves of van der Waals complexes. We emphasize the benefits of partial augmentation of the higher-zeta basis sets in preference to full augmentation at a lower ζ level. Benefits and limitations of the use of fully, partially, and minimally augmented basis sets are reviewed for different electronic structure methods and molecular properties. We have found that minimal augmentation is almost always enough for density functional theory (DFT) when applied to ionization potentials, electron affinities, atomization energies, barrier heights, and hydrogen-bond energies. For electric dipole polarizabilities, we find that augmentation beyond minimal has an average effect of 8% at the polarized triple-ζ level and 5% at the polarized quadruple-ζ level. The effects are larger for potential energy curves of van der Waals complexes. The effects are also larger for wave function theory (WFT). Even for WFT though, full augmentation is not needed for most purposes, and a level of augmentation between minimal and full is optimal for most problems. The calendar basis sets named after the months provide a convergent sequence of partially augmented basis sets that can be used for such calculations. The jun-cc-pV(T+d)Z basis set is very useful for MP2-F12 calculations of barrier heights and hydrogen bond strengths.

Convergent Partially Augmented Basis Sets for Post-Hartree-Fock Calculations of Molecular Properties and Reaction Barrier Heights.

We present sets of convergent, partially augmented basis set levels corresponding to subsets of the augmented "aug-cc-pV(n+d)Z" basis sets of Dunning and co-workers. We show that for many molecular properties a basis set fully augmented with diffuse functions is computationally expensive and almost always unnecessary. On the other hand, unaugmented cc-pV(n+d)Z basis sets are insufficient for many properties that require diffuse functions. Therefore, we propose using intermediate basis sets. We developed an efficient strategy for partial augmentation, and in this article, we test it and validate it. Sequentially deleting diffuse basis functions from the "aug" basis sets yields the "jul", "jun", "may", "apr", etc. basis sets. Tests of these basis sets for Møller-Plesset second-order perturbation theory (MP2) show the advantages of using these partially augmented basis sets and allow us to recommend which basis sets offer the best accuracy for a given number of basis functions for calculations on large systems. Similar truncations in the diffuse space can be performed for the aug-cc-pVxZ, aug-cc-pCVxZ, etc. basis sets.

Efficient Diffuse Basis Sets for Density Functional Theory.

Eliminating all but the s and p diffuse functions on the non-hydrogenic atoms and all diffuse functions on the hydrogen atoms from the aug-cc-pV(x+d)Z basis sets of Dunning and co-workers, where x = D, T, Q, ..., yields the previously proposed "minimally augmented" basis sets, called maug-cc-pV(x+d)Z. Here, we present extensive and systematic tests of these basis sets for density functional calculations of chemical reaction barrier heights, hydrogen bond energies, electron affinities, ionization potentials, and atomization energies. The tests show that the maug-cc-pV(x+d)Z basis sets are as accurate as the aug-cc-pV(x+d)Z ones for density functional calculations, but the computational cost savings are a factor of about two to seven.

Phase space prediction of product branching ratios: canonical competitive nonstatistical model.

We present a new model for predicting branching ratios of chemical reactions when a branching of the reaction path occurs after the dynamical bottleneck, including the case where it occurs after an intermediate. The model is based on combining nonstatistical phase space theory for the direct component of a reaction with variational transition-state theory for an indirect component of reaction. The competition between direct and indirect processes is treated by an extension of the unified statistical model. This new method provides a way to understand the factors that control this kind of chemical reaction and to perform calculations using high-level electronic structure methods for complex systems. The model is based on quantized energy levels of transition states and products, and it involves the same information as required for calculating transition-state rate constants and equilibrium constants plus a phenomenological relaxation time, which was taken from previous work. For the textbook reaction of the hydroboration of propene by BH(3) it has recently been inferred that the selectivity can only be understood by consideration of dynamical trajectories. However, the calculated branching fraction of this prototype reaction increases from 2%-3% when calculated under the inappropriate assumption of complete equilibration of the intermediate to from 8%-9% when calculated with the new theory, which requires only limited information about the system and does not involve running trajectories. The calculated result is in reasonable agreement with experiment (approximately 10%).

Efficient Diffuse Basis Sets: cc-pVxZ+ and maug-cc-pVxZ.

We combine the diffuse basis functions from the 6-31+G basis set of Pople and co-workers with the correlation-consistent basis sets of Dunning and co-workers. In both wave function and density functional calculations, the resulting basis sets reduce the basis set superposition error almost as much as the augmented correlation-consistent basis sets, although they are much smaller. In addition, in density functional calculations the new basis sets, called cc-pVxZ+ where x = D, T, Q, ..., or x = D+d, T+d, Q+d, ..., give very similar energetic predictions to the much larger aug-cc-pVxZ basis sets. However, energetics calculated from correlated wave function calculations are more slowly convergent with respect to the addition of diffuse functions. We also examined basis sets with the same number and type of functions as the cc-pVxZ+ sets but using the diffuse exponents of the aug-cc-pVxZ basis sets and found very similar performance to cc-pVxZ+; these basis sets are called minimally augmented cc-pVxZ, which we abbreviate as maug-cc-pVxZ.

Mechanism of the reaction catalyzed by DL-2-haloacid dehalogenase as determined from kinetic isotope effects.

dl-2-Haloacid dehalogenase from Pseudomonas sp. 113 is a unique enzyme because it acts on the chiral carbons of both enantiomers, although its amino acid sequence is similar only to that of d-2-haloacid dehalogenase from Pseudomonas putida AJ1 that specifically acts on (R)-(+)-2-haloalkanoic acids. Furthermore, the catalyzed dehalogenation proceeds without formation of an ester intermediate; instead, a water molecule directly attacks the alpha-carbon of the 2-haloalkanoic acid. We have studied solvent deuterium and chlorine kinetic isotope effects for both stereoisomeric reactants. We have found that chlorine kinetic isotope effects are different: 1.0105 +/- 0.0001 for (S)-(-)-2-chloropropionate and 1.0082 +/- 0.0005 for the (R)-(+)-isomer. Together with solvent deuterium isotope effects on V(max)/K(M), 0.78 +/- 0.09 for (S)-(-)-2-chloropropionate and 0.90 +/- 0.13 for the (R)-(+)-isomer, these values indicate that in the case of the (R)-(+)-reactant another step preceding the dehalogenation is partly rate-limiting. Under the V(max) conditions, the corresponding solvent deuterium isotope effects are 1.48 +/- 0.10 and 0.87 +/- 0.27, respectively. These results indicate that the overall reaction rates are controlled by different steps in the catalysis of (S)-(-)- and (R)-(+)-reactants.