Structures, electron affinities, and harmonic vibrational frequencies of the simplest alkyl peroxyl radicals and their anions

The molecular structures and electron affinities of the R-OO/R-OO(-) (R = CH3, C2H5, n-C3H7, n-C4H9, n-C5H11, i-C3H7, t-C4H9) species have been determined using seven different density functional or hybrid Hartree-Fock density functional methods. The basis set used in this work is of double-zeta plus polarization quality with additional diffuse s-type and p-type functions, denoted DZP++. The geometries are fully optimized with each density functional theory method. Harmonic vibrational frequencies were found to be within 3.1% of available experimental values for most functionals. Two different types of the neutral-anion energy separations reported in this work are the adiabatic electron affinity and the vertical detachment energy. The most reliable adiabatic electron affinities obtained at the DZP++ BP86 level of theory are 1.150 (CH3OO), 1.124 (C2H5OO), 1.146 (n-C3H7OO), 1.173 (n-C4H9OO), 1.184 (n-C5H11OO), 1.145 (i-C3H7OO), and 1.114 eV (t-C4H9OO). Compared with the experimental values, the average absolute error of the BPW91 method is 0.05 eV.     

(+)-儿茶素金属配合物的结构、红外光谱和反应活性

用密度泛函理论(DFT)B3LYP方法, 对非金属原子和金属原子分别采用6-311+G(d,p)基组和LANL2DZ基组, 计算并分析了(+)-儿茶素(Cc, C₁₅H₁₄O₆)及其与金属形成的配合物分子(M-Cc, M=Ca, Zn, Cd, Cu, Al, Cr)的几何构型、红外光谱和反应活性的异同. 计算结果表明, M-Cc的分子结构、红外光谱与反应活性均不同于其前体Cc. 形成金属配合物后, 取代基团及结构的改变使得红外光谱有所差异. 前线分子轨道及概念DFT指数计算结果显示, 一些M-Cc体系的反应活性要强于Cc单体. 金属离子的不同使得配合物的各指数有所差异. 这些结果将为进一步认识(+)-儿茶素及其相关化合物的结构、红外光谱和反应活性提供有益启示.     

Activation energies of pericyclic reactions: performance of DFT, MP2, and CBS-QB3 methods for the prediction of activation barriers and reaction energetics of 1,3-dipolar cycloadditions, and revised activation enthalpies for a standard set of hydrocarbon pericyclic reactions

Activation barriers and reaction energetics for the three main classes of 1,3-dipolar cycloadditions, including nine different reactions, were evaluated with the MPW1K and B3LYP density functional methods, MP2, and the multicomponent CBS-QB3 method. The CBS-QB3 values were used as standards for 1,3-dipolar cycloaddition activation barriers and reaction energetics, and the density functional theory (DFT) and MP2 methods were benchmarked against these values. The MPW1K/6-31G* method and basis set performs best for activation barriers, with a mean absolute deviation (MAD) value of 1.1 kcal/mol. The B3LYP/6-31G* method and basis set performs best for reaction enthalpies, with a MAD value of 2.4 kcal/mol, while the MPW1K method shows large errors for reaction energetics. The MP2 method gives the expected systematic underestimation of barriers. Concerted and nearly synchronous transition structures are predicted by all DFT and MP2 methods. Also reported are revised estimated 0 K experimental activation enthalpies for a standard set of hydrocarbon pericyclic reactions and updated comparisons to experiment for DFT, ab initio, and multicomponent methods. B3LYP and MPW1K methods with MAD values of 1.5 and 2.1 kcal/mol, respectively, fortuitously outperform the multicomponent CBS-QB3 method, which has a MAD value of 2.3. The MAD value of the O3LYP functional improves to 2.4 kcal/mol from the previously reported 3.0 kcal/mol.     

Generalized hybrid orbital for the treatment of boundary atoms in combined quantum mechanical and molecular mechanical calculations using the semiempirical parameterized model 3 method

The application of combined quantum mechanical (QM) and molecular mechanical methods to large molecular systems requires an adequate treatment of the boundary between the two approaches. In this article, we extend the generalized hybrid orbital (GHO) method to the semiempirical parameterized model 3 (PM3) Hamiltonian combined with the CHARMM force field. The GHO method makes use of four hybrid orbitals, one of which is included in the QM region in self-consistent field optimization and three are treated as auxiliary orbitals that do not participate in the QM optimization, but they provide an effective electric field for interactions. An important feature of the GHO method is that the semiempirical parameters for the boundary atom are transferable, and these parameters have been developed for a carbon boundary atom consistent with the PM3 model. The combined GHO-PM3/CHARMM model has been tested on molecular geometry and proton affinity for a series of organic compounds.Acknowledgement We thank the National Institutes of Health for support of this research.Contribution to the Jacopo Tomasi Honorary Issue     

Benchmark of density functional theory methods on the prediction of bond energies and bond distances of noble-gas containing molecules

We have tested three pure density functional theory (DFT) functionals, BLYP, MPWPW91, MPWB95, and ten hybrid DFT functionals, B3LYP, B3P86, B98, MPW1B95, MPW1PW91, BMK, M05-2X, M06-2X, B2GP-PLYP, and DSD-BLYP with a series of commonly used basis sets on the performance of predicting the bond energies and bond distances of 31 small neutral noble-gas containing molecules. The reference structures were obtained using the CCSD(T)∕aug-cc-pVTZ theory and the reference energies were based on the calculation at the CCSD(T)∕CBS level. While in general the hybrid functionals performed significantly better than the pure functionals, our tests showed a range of performance by these hybrid functionals. For the bond energies, the MPW1B95∕6-311+G(2df,2pd), BMK∕aug-cc-pVTZ, B2GP-PLYP∕aug-cc-pVTZ, and DSD-BLYP∕aug-cc-pVTZ methods stood out with mean unsigned errors of 2.0-2.3 kcal∕mol per molecule. For the bond distances, the MPW1B95∕6-311+G(2df,2pd), MPW1PW91∕6-311+G(2df,2pd), and B3P86∕6-311+G(2df,2pd), DSD-BLYP∕6-311+G(2df,2pd), and DSD-BLYP∕aug-cc-pVTZ methods stood out with mean unsigned errors of 0.008-0.013 A? per bond. The current study showed that a careful selection of DFT functionals is very important in the study of noble-gas chemistry, and the most recommended methods are MPW1B95∕6-311+G(2df,2pd) and DSD-BLYP∕aug-cc-pVTZ.     

Electron momentum spectroscopy of CF2Cl2: experimental and theoretical momentum profiles for outer valence orbitals

Electron momentum distributions for outer valence orbitals of CF2Cl2 have been obtained by (e,2e) electron momentum spectroscopy at an incident energy of 1200 eV + binding energy. The experimental electron momentum profiles are compared with Hartree-Fock and density functional theory (DFT) calculations using B3LYP hybrid functional with the 6-31G and 6-311+G* basis sets. Generally, the shapes of the experimental momentum profiles are well reproduced by DFT calculations using larger basis sets 6-311 + G*. An attempt has been made to clarify the ordering of the outer valence orbitals, which have been in controversy, by comparing experimental results with B3LYP/6-311 + G* calculations.     

Quantum mechanical calculations of conformationally relevant 1H and 13C NMR chemical shifts of calixarene systems

[graphs: see text] QM GIAO calculations of 13C and 1H chemical shift values of the ArCH2Ar group have been performed, using the hybrid DFT functional MPW1PW91 and the 6-31G(d,p) basis set, on some representative calixarenes and on a series of simplified calixarene models allowing derivation of chemical shift surfaces versus phi and chi dihedral angles. A good reproduction of experimental data was obtained. The applicability of chemical shift surfaces in the study of calixarene conformational features is illustrated.     

丙烷分子的价轨道(2b2)电子的动量分布测量

电子动量谱学（EMS）是在原子、分子和固体物理中研究电子结构的一种强有力的工具，它基于运动学条件完全确定的（e，2e）碰撞电离反应[1-3]．本文报告用高分辨电子动量谱仪首次测量得到丙烷门3H8）分子的价轨道电子（252）的动量分布·丙烷（C3Hs）价轨道电子的动量分布实验是     

Quantum mechanical calculations of conformationally relevant 1H and 13C NMR chemical shifts of N-, O-, and S-substituted calixarene systems

QM GIAO calculations of (13)C and (1)H chemical shift values of the ArCH(2)Ar group in N-, O-, and S-substituted calixarene systems were performed with a hybrid DFT functional MPW1PW91 and 6-31G(d,p) basis set. A good reproduction of experimental data was obtained for some representative calixarenes and for a series of simplified calixarene models. This allowed the derivation of chemical shift surfaces versus phi and chi dihedral angles. The applicability of chemical shift surfaces in the study of calixarene conformational features is illustrated.     

Optimization of auxiliary basis sets for the LEDO expansion and a projection technique for LEDO-DFT

We present a systematic procedure for the optimization of the expansion basis for the limited expansion of diatomic overlap density functional theory (LEDO-DFT) and report on optimized auxiliary orbitals for the Ahlrichs split valence plus polarization basis set (SVP) for the elements H, Li--F, and Na--Cl. A new method to deal with near-linear dependences in the LEDO expansion basis is introduced, which greatly reduces the computational effort of LEDO-DFT calculations. Numerical results for a test set of small molecules demonstrate the accuracy of electronic energies, structural parameters, dipole moments, and harmonic frequencies. For larger molecular systems the numerical errors introduced by the LEDO approximation can lead to an uncontrollable behavior of the self-consistent field (SCF) process. A projection technique suggested by L?wdin is presented in the framework of LEDO-DFT, which guarantees for SCF convergence. Numerical results on some critical test molecules suggest the general applicability of the auxiliary orbitals presented in combination with this projection technique. Timing results indicate that LEDO-DFT is competitive with conventional density fitting methods.     

The outer valance orbital electron densities of cyclopentane by binary (e,2e) spectroscopy

The binding energy spectra and electron distributions in momentum space of the valence orbitals of cyclopentane (C(5)H(10)) are studied by Electron Momentum Spectroscopy (EMS) in a noncoplanar symmetric geometry. The impact energy was 1200 eV plus binding energy and energy resolution of the EMS spectrometer was 1.2 eV. The experimental momentum profiles of the outer valence orbitals are compared with the theoretical momentum distributions calculated using Hartree-Fock and density functional theory (DFT) methods. The shapes of the experimental momentum distributions are generally quite well described by both the Hartree-Fock and DFT calculations when the large and diffuse basis sets are used.     

Theoretical study of the vibrational frequencies of carbon disulfide

A benchmark comparison for different computational methods and basis sets has been presented. In this study, five computational methods (Hartree–Fock (HF), MP2, B3LYP, MPW1MP91, and PBE1PBE) along with 18 basis sets have been applied to optimize the geometry of carbon disulfide (CS₂), and further calculate the vibrational frequencies of the optimized geometries. The differences between the calculated frequencies and corresponding experimental data are used to evaluate the efficiency of each combination of computational method and basis set. The comparison of frequency difference indicates that B3LYP generally gives the best prediction of frequencies for CS₂, whereas the other two density functional theory (DFT) methods, i.e., MPW1PW91 and PBE1PBE, often give parallel results. Although MP2 predicts the frequencies with accuracy almost as good as those from DFT methods, in a particular case, HF calculation outperforms MP2 as well as MPW1PW91 and PBE1PBE for prediction of the frequency of asymmetrical stretching for CS₂. © 2013 Wiley Periodicals, Inc.     

An assessment of theoretical procedures for predicting the thermochemistry and kinetics of hydrogen abstraction by methyl radical from benzene

The reaction enthalpy (298 K), barrier (0 K), and activation energy and preexponential factor (600-800 K) have been examined computationally for the abstraction of hydrogen from benzene by the methyl radical, to assess their sensitivity to the applied level of theory. The computational methods considered include high-level composite procedures, including W1, G3-RAD, G3(MP2)-RAD, and CBS-QB3, as well as conventional ab initio and density functional theory (DFT) methods, with the latter two classes employing the 6-31G(d), 6-31+G(d,p) and/or 6-311+G(3df,2p) basis sets, and including ZPVE/thermal corrections obtained from 6-31G(d) or 6-31+G(d,p) calculations. Virtually all the theoretical procedures except UMP2 are found to give geometries that are suitable for subsequent calculation of the reaction enthalpy and barrier. For the reaction enthalpy, W1, G3-RAD, and URCCSD(T) give best agreement with experiment, while the large-basis-set DFT procedures slightly underestimate the endothermicity. The reaction barrier is slightly more sensitive to the choice of basis set and/or correlation level, with URCCSD(T) and the low-cost BMK method providing values in close agreement with the benchmark G3-RAD value. Inspection of the theoretically calculated rate parameters reveals a minor dependence on the level of theory for the preexponential factor. There is more sensitivity for the activation energy, with a reasonable agreement with experiment being obtained for the G3 methods and the hybrid functionals BMK, BB1K, and MPW1K, especially in combination with the 6-311+G(3df,2p) basis set. Overall, the high-level G3-RAD composite procedure, URCCSD(T), and the cost-effective DFT methods BMK, BB1K, and MPW1K give the best results among the methods assessed for calculating the thermochemistry and kinetics of hydrogen abstraction by the methyl radical from benzene.     

Calculation of bond dissociation energies for oxygen containing molecules by ab initio and density functional theory methods

Two ab initio (ROHF and MP2), one local (SVWN), four hybrid (BHandH, BHandHLYP, Becke3LYP, and Becke3P86), and two nonlocal (BLYP and BP86) density functional theory (DFT) methods are used for calculating the dissociation energies of molecules that contain H(SINGLE BOND)O, O(SINGLE BOND)O and O(SINGLE BOND)C bonds. The sensitivity to the basis set of the prediction of bond dissociation energies with DFT methods was tested with Becke3LYP on the H(SINGLE BOND)O dissociation energy of water. The 6–31 + G(d) methods are chosen as the smallest basis set which produces reasonable results. The calculated values for all other ab initio and DFT methods were performed with these basis sets and then compared with the experimental data. The suitability of DFT methods for computing reliable bond dissociation energies of oxygen containing molecules is discussed. © 1996 John Wiley & Sons, Inc.     

Bonding and excitation in COCu(001) from a cluster model and density functional treatments

The bonding properties and charge distributions of the COCu(001) system have been studied within density functional theory (DFT) with several density functionals. A Cu18(9,4,5)CO three layer cluster was found to give bond distances and energies in agreement with previous experimental and theoretical results for low coverage systems, provided the atomic basis set includes diffuse orbitals and d-orbitals at the Cu atoms. Charge distributions give insight on the nature of the localized adsorbate bonding. Time-dependent DFT results on excitation energies and on transition and average electric dipoles, relevant to photodesorption, are also presented.     

Combined valence bond-molecular mechanics potential-energy surface and direct dynamics study of rate constants and kinetic isotope effects for the H + C2H6 reaction

This article presents a multifaceted study of the reaction H+C(2)H(6)-->H(2)+C(2)H(5) and three of its deuterium-substituted isotopologs. First we present high-level electronic structure calculations by the W1, G3SX, MCG3-MPWB, CBS-APNO, and MC-QCISD/3 methods that lead to a best estimate of the barrier height of 11.8+/-0.5 kcal/mol. Then we obtain a specific reaction parameter for the MPW density functional in order that it reproduces the best estimate of the barrier height; this yields the MPW54 functional. The MPW54 functional, as well as the MPW60 functional that was previously parametrized for the H+CH(4) reaction, is used with canonical variational theory with small-curvature tunneling to calculate the rate constants for all four ethane reactions from 200 to 2000 K. The final MPW54 calculations are based on curvilinear-coordinate generalized-normal-mode analysis along the reaction path, and they include scaled frequencies and an anharmonic C-C bond torsion. They agree with experiment within 31% for 467-826 K except for a 38% deviation at 748 K; the results for the isotopologs are predictions since these rate constants have never been measured. The kinetic isotope effects (KIEs) are analyzed to reveal the contributions from subsets of vibrational partition functions and from tunneling, which conspire to yield a nonmonotonic temperature dependence for one of the KIEs. The stationary points and reaction-path potential of the MPW54 potential-energy surface are then used to parametrize a new kind of analytical potential-energy surface that combines a semiempirical valence bond formalism for the reactive part of the molecule with a standard molecular mechanics force field for the rest; this may be considered to be either an extension of molecular mechanics to treat a reactive potential-energy surface or a new kind of combined quantum-mechanical/molecular mechanical (QM/MM) method in which the QM part is semiempirical valence bond theory; that is, the new potential-energy surface is a combined valence bond molecular mechanics (CVBMM) surface. Rate constants calculated with the CVBMM surface agree with the MPW54 rate constants within 12% for 534-2000 K and within 23% for 200-491 K. The full CVBMM potential-energy surface is now available for use in variety of dynamics calculations, and it provides a prototype for developing CVBMM potential-energy surfaces for other reactions.     

Ab initio studies of crystalline nitromethane under high pressure

We have studied the mechanical compressibility and band structure of solid nitromethane both in equilibrium and compressed states using Hartree-Fock and density functional theory (DFT) with atom-centered all-electron linear combination of atomic orbitals basis sets. Hartree-Fock calculations with a 6-21G basis set, uncorrected for basis set superposition error, gave the best agreement with experimental compression studies. These results may be due to the cancellation of basis set superposition error with dispersion force errors. The equilibrium DFT band gap is comparable to the lowest-energy feature in electron-impact spectroscopy of nitromethane but underpredicts the optical absorption gap; we interpret these features in terms of the presence of tightly bound excitons. Only minor changes in the gap are observed under hydrostatic compression.     

Ab initio molecular modeling of 13C NMR chemical shifts of polymers. 1. Ethylene–norbornene copolymers

Cycloolefin copolymers (COC) have recently raised much interest because of their excellent thermal and optical properties, largely determined by the chain composition and stereochemistry. Previous force‐field computations allowed us to define the main conformational characteristics of ethylene–norbornene (E–N) copolymers and to contribute to the elucidation of their microstructure on the basis of empirical relationships between conformation and ¹³C nuclear magnetic resonance (NMR) chemical shifts. A thorough test of ab initio ¹³C chemical shifts computations [gauge‐invariant atomic orbitals (GIAO)] on known cases shows that the agreement with experimental data is quite good, especially with the MPW1PW91 density functional theory (DFT), using the 6‐311+G(2d,p) basis set on properly energy‐minimized structures. We applied this method on proper model compounds to confirm the signal assignment of the spectra of E–N copolymers in the presence of norbornene microblocks, where strong effects arising from ring distortions are expected to occur. The results nicely confirm the latest assignment of norbornene signals belonging to ENNE sequences. This shows the great potentialities of GIAO/DFT computations with regard to complex spectra interpretation and polymer microstructural investigations. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002