Henryk A. Witek

Infrared Absorption Spectrum of the Simplest Criegee Intermediate CH2OO

More Criegee Sightings The reaction of ozone with unsaturated hydrocarbons produces short-lived molecules termed Criegee intermediates. The simplest such molecule, H2CO2, was recently detected and monitored in the laboratory. Su et al. (p. 174; see the Perspective by Vereecken) have obtained its vibrational spectrum, which could ultimately enable direct measurements of its reactivity in the atmosphere. Taatjes et al. (p. 177; see the Perspective by Vereecken) report on the laboratory preparation and reactivity of the next heavier Criegee intermediate, which bears a methyl group in place of one of the hydrogen atoms. The vibrational spectrum of an intermediate involved in ozone chemistry may facilitate its direct study in the atmosphere. [Also see Perspective by Vereecken] The Criegee intermediates are carbonyl oxides postulated to play key roles in the reactions of ozone with unsaturated hydrocarbons; these reactions constitute an important mechanism for the removal of unsaturated hydrocarbons and for the production of OH in the atmosphere. Here, we report the transient infrared (IR) absorption spectrum of the simplest Criegee intermediate CH2OO, produced from CH2I + O2 in a flow reactor, using a step-scan Fourier-transform spectrometer. The five observed bands provide definitive identification of this intermediate. The observed vibrational frequencies are more consistent with a zwitterion rather than a diradical structure of CH2OO. The direct IR detection of CH2OO should prove useful for kinetic and mechanistic investigations of the Criegee mechanism.

Parameter calibration of transition-metal elements for the spin-polarized self-consistent-charge density-functional tight-binding (DFTB) method : Sc, Ti, Fe, Co, and Ni

Recently developed parameters for five first-row transition-metal elements (M = Sc, Ti, Fe, Co, and Ni) in combination with H, C, N, and O as well as the same metal (M-M) for the spin-polarized self-consistent-charge density-functional tight-binding (DFTB) method have been calibrated. To test their performance a couple sets of compounds have been selected to represent a variety of interactions and bonding schemes that occur frequently in transition-metal containing systems. The results show that the DFTB method with the present parameters in most cases reproduces structural properties very well, but the bond energies and the relative energies of different spin states only qualitatively compared to the B3LYP/SDD+6-31G(d) density functional (DFT) results. An application to the ONIOM(DFT:DFTB) indicates that DFTB works well as the low level method for the ONIOM calculation.

Systematic study of vibrational frequencies calculated with the self-consistent charge density functional tight-binding method

We present a detailed study of harmonic vibrational frequencies obtained with the self‐consistent charge density functional tight‐binding (SCC‐DFTB) method. Our testing set comprises 66 molecules and 1304 distinct vibrational modes. Harmonic vibrational frequencies are computed using an efficient analytical algorithm developed and coded by the authors. The obtained results are compared to experiment and to other theoretical findings. Scaling factor for the SCC‐DFTB method, determined by minimization of mean absolute deviation of scaled frequencies, is found to be 0.9933. The accuracy of the scaled SCC‐DFTB frequencies is noticeably better than for other semiempirical methods (including standard DFTB method) and approximately twice worse than for other well established scaled ab initio quantum chemistry methods (e.g., HF, BLYP, B3LYP). Mean absolute deviation for the scaled SCC‐DFTB frequencies is 56 cm−1, while standard deviation is 82 cm−1, and maximal absolute deviation is as large as 529 cm−1. Using SCC‐DFTB allows for substantial time savings; computational time is reduced from hours to seconds when compared to standard ab initio techniques.

DFTB Parameters for the Periodic Table: Part 1, Electronic Structure

A parametrization scheme for the electronic part of the density-functional based tight-binding (DFTB) method that covers the periodic table is presented. A semiautomatic parametrization scheme has been developed that uses Kohn-Sham energies and band structure curvatures of real and fictitious homoatomic crystal structures as reference data. A confinement potential is used to tighten the Kohn-Sham orbitals, which includes two free parameters that are used to optimize the performance of the method. The method is tested on more than 100 systems and shows excellent overall performance.

Analytical second-order geometrical derivatives of energy for the self-consistent-charge density-functional tight-binding method

Analytical formulation of the second-order geometrical derivatives of energy for the self-consistent-charge density-functional tight-binding (SCC-DFTB) method is presented. To test its quality and numerical performance, the derived formalism has been coded and applied for calculation of harmonic vibrational frequencies for a set of 17 small and medium size molecules. For this set, the average absolute deviation from experiment is 99 cm(-1) for SCC-DFTB vs 62 cm(-1) for the Møller-Plesset second-order perturbation theory with the cc-pVDZ basis set (MP2/cc-pVDZ) and 32 cm(-1) for the B3LYP density functional method with the same basis set (B3LYP/cc-pVDZ), while the maximal deviation is 465 cm(-1) vs 1,741 cm(-1) for MP2/cc-pVDZ and 112 cm(-1) for B3LYP/cc-pVDZ. The SCC-DFTB results are in reasonable agreement with experiments as well as with ab initio and density-functional results, and are better than other semiempirical methods. The SCC-DFTB method allows for considerable computational time saving when compared to other methods while retaining similar overall accuracy. Data for a series of conjugated polyenes show that an analytical formulation of SCC-DFTB is noticeably faster than its numerical formulation.

Infrared identification of the Criegee intermediates syn- and anti-CH3CHOO, and their distinct conformation-dependent reactivity

The Criegee intermediates are carbonyl oxides that play critical roles in ozonolysis of alkenes in the atmosphere. So far, the mid-infrared spectrum of only the simplest Criegee intermediate CH2OO has been reported. Methyl substitution of CH2OO produces two conformers of CH3CHOO and consequently complicates the infrared spectrum. Here we report the transient infrared spectrum of syn- and anti-CH3CHOO, produced from CH3CHI + O2 in a flow reactor, using a step-scan Fourier-transform spectrometer. Guided and supported by high-level full-dimensional quantum calculations, rotational contours of the four observed bands are simulated successfully and provide definitive identification of both conformers. Furthermore, anti-CH3CHOO shows a reactivity greater than syn-CH3CHOO towards NO/NO2; at the later period of reaction, the spectrum can be simulated with only syn-CH3CHOO. Without NO/NO2, anti-CH3CHOO also decays much faster than syn-CH3CHOO. The direct infrared detection of syn- and anti-CH3CHOO should prove useful for field measurements and laboratory investigations of the Criegee mechanism.

Infrared absorption of methanol clusters (CH3OH)n with n = 2−6 recorded with a time-of-flight mass spectrometer using infrared depletion and vacuum-ultraviolet ionization

We investigated IR spectra in the CH- and OH-stretching regions of size-selected methanol clusters, (CH(3)OH)(n) with n = 2-6, in a pulsed supersonic jet by using the IR-VUV (vacuum-ultraviolet) ionization technique. VUV emission at 118 nm served as the source of ionization in a time-of-flight mass spectrometer. The tunable IR laser emission served as a source of predissociation or excitation before ionization. The variations of intensity of protonated methanol cluster ions (CH(3)OH)(n)H(+) and CH(3)OH(+) and (CH(3)OH)(2)(+) were monitored as the IR laser light was tuned across the range 2650-3750 cm(-1). Careful processing of these action spectra based on photoionization efficiencies and the production and loss of each cluster due to photodissociation yielded IR spectra of the size-selected clusters. Spectra of methanol clusters in the OH region have been extensively investigated; our results are consistent with previous reports, except that the band near 3675 cm(-1) is identified as being associated with the proton acceptor of (CH(3)OH)(2). Spectra in the CH region are new. In the region 2800-3050 cm(-1), bands near 2845, 2956, and 3007 cm(-1) for CH(3)OH split into 2823, 2849, 2934, 2955, 2984, and 3006 cm(-1) for (CH(3)OH)(2) that correspond to proton donor and proton acceptor, indicating that the methanol dimer has a preferred open-chain structure. In contrast, for (CH(3)OH)(3), the splitting diminishes and the bands near 2837, 2954, and 2987 cm(-1) become narrower, indicating a preferred cyclic structure. Anharmonic vibrational wavenumbers predicted for the methanol open-chain dimer and the cyclic trimer with the B3LYP∕VPT2∕ANO1 level of theory are consistent with experimental results. For the tetramer and pentamer, the spectral pattern similar to that of the trimer but with greater widths was observed, indicating that the most stable structures are also cyclic.

Intruder states in multireference perturbation theory: The ground state of manganese dimer

A detailed analysis of a severe intruder state problem in the multistate multireference perturbation theory (MS‐MRPT) calculations on the ground state of manganese dimer is presented. An enormous number of detected intruder states (> 5000) do not permit finding even an approximate shape of the X1Σ  g+ potential energy curve. The intruder states are explicitly demonstrated to originate from quasidegeneracies in the zeroth‐order Hamiltonian spectrum. The electronic configurations responsible for appearance of the quasidegeneracies are identified as single and double excitations from the active orbitals to the external orbitals. It is shown that the quasidegeneracy problem can be completely eliminated using shift techniques despite of its severity. The resultant curves are smooth and continuous. Unfortunately, strong dependence of the spectroscopic parameters of the X1Σ  g+ state on the shift parameter is observed. This finding rises serious controversies regarding validity of employing shift techniques for solving the intruder state problem in MS‐MRPT. Various alternative approaches of removing intruder states (e.g., modification of the basis set or changing the active space) are tested. None of these conventional techniques is able to fully avoid the quasidegeneracies. We believe that the MS‐MRPT calculations on the three lowest Ag states of manganese dimer constitute a perfect benchmark case for studying the behavior of MRPT in extreme situations.

Automatized Parametrization of SCC-DFTB Repulsive Potentials: Application to Hydrocarbons †

In this work, we derive and test a new automatized strategy to construct repulsive potentials for the self-consistent charge density functional tight-binding (SCC-DFTB) method. This approach allows one to explore the parameter space in a systematic fashion in order to find optimal solutions. We find that due to the limited flexibility of the SCC-DFTB electronic part, not all properties can be optimized simultaneously. For example, the optimization of heats of formation is in conflict with the optimization of vibrational frequencies. Therefore, a special parametrization for vibrational frequencies is derived. It is shown that the performance of SCC-DFTB can be significantly improved using a more elaborate fitting strategy. A new fit for C and H is presented, which results in an average error of 2.6 kcal/mol for heats of formations for a large set of hydrocarbons, indicating that the performance of SCC-DFTB can be systematically improved also for other elements.

Modeling carbon nanostructures with the self-consistent charge density-functional tight-binding method : Vibrational spectra and electronic structure of C28, C60, and C70

The self-consistent charge density-functional tight-binding (SCC-DFTB) method is employed for studying various molecular properties of small fullerenes: C(28), C(60), and C(70). The computed bond distances, vibrational infrared and Raman spectra, vibrational densities of states, and electronic densities of states are compared with experiment (where available) and density-functional theory (DFT) calculations using various basis sets. The presented DFT benchmark calculations using the correlation-consistent polarized valence triple zeta basis set are at present the most extensive calculations on harmonic frequencies of these species. Possible limitations of the SCC-DFTB method for the prediction of molecular vibrational and optical properties are discussed. The presented results suggest that SCC-DFTB is a computationally feasible and reliable method for predicting vibrational and electronic properties of such carbon nanostructures comparable in accuracy with small to medium size basis set DFT calculations at the computational cost of standard semiempirical methods.