W4 theory for computational thermochemistry: in pursuit of confident sub-kJ/mol predictions

In an attempt to improve on our earlier W3 theory [J. Chem. Phys. {\bf 120}, 4129 (2004)] we consider such refinements as more accurate estimates for the contribution of connected quadruple excitations ($\hat{T}_4$), inclusion of connected quintuple excitations ($\hat{T}_5$), diagonal Born-Oppenheimer corrections (DBOC), and improved basis set extrapolation procedures. Revised experimental data for validation purposes were obtained from the latest version of the ATcT (Active Thermochemical Tables) Thermochemical Network. We found that the CCSDTQ$-$CCSDT(Q) difference converges quite rapidly with the basis set, and that the formula 1.10[CCSDT(Q)/cc-pVTZ+CCSDTQ/cc-pVDZ$-$CCSDT(Q)/cc-pVDZ] offers a very reliable as well as fairly cost-effective estimate of the basis set limit $\hat{T}_4$ contribution. The largest $\hat{T}_5$ contribution found in the present work is on the order of 0.5 kcal/mol (for ozone). DBOC corrections are significant at the 0.1 kcal/mol level in hydride systems. . Based on the accumulated experience, a new computational thermochemistry protocol for first-and second-row main-group systems, to be known as W4 theory, is proposed. Our W4 atomization energies for a number of key species are in excellent agreement (better than 0.1 kcal/mol on average, 95% confidence intervals narrower than 1 kJ/mol) with the latest experimental data obtained from Active Thermochemical Tables. A simple {\em a priori} estimate for the importance of post-CCSD(T) correlation contributions (and hence a pessimistic estimate for the error in a W2-type calculation) is proposed.Comment: J. Chem. Phys., in press; electronic supporting information available at http://theochem.weizmann.ac.il/web/papers/w4.htm

O(^3P) +CO_2 Collisions at Hyperthermal Energies: Dynamics of Nonreactive Scattering, Oxygen Isotope Exchange, and Oxygen-Atom Abstraction

The dynamics of O(^3P) + CO_2 collisions at hyperthermal energies were investigated experimentally and theoretically. Crossed-molecular-beams experiments at Ecoll = 98.8 kcal mol^(–1) were performed with isotopically labeled ^(12)C^(18)O_2 to distinguish products of nonreactive scattering from those of reactive scattering. The following product channels were observed: elastic and inelastic scattering (^(16)O(^3P) + ^(12)C^(18)O^2), isotope exchange (^(18)O + ^(16)O^(12)C^(18)O), and oxygen-atom abstraction (^(18)O^(16)O + ^(12)C^(18)O). Stationary points on the two lowest triplet potential energy surfaces of the O(^3P) + CO_2 system were characterized at the CCSD(T)/aug-cc-pVTZ level of theory and by means of W4 theory, which represents an approximation to the relativistic basis set limit, full-configuration-interaction (FCI) energy. The calculations predict a planar CO_3(C_(2v),^3A″) intermediate that lies 16.3 kcal mol^(–1) (W4 FCI excluding zero point energy) above reactants and is approached by a C_(2v) transition state with energy 24.08 kcal mol^(–1). Quasi-classical trajectory (QCT) calculations with collision energies in the range 23–150 kcal mol^(–1) were performed at the B3LYP/6-311G(d) and BMK/6-311G(d) levels. Both reactive channels observed in the experiment were predicted by these calculations. In the isotope exchange reaction, the experimental center-of-mass (c.m.) angular distribution, T(θ_(c.m.)), of the ^(16)O^(12)C^(18)O products peaked along the initial CO_2 direction (backward relative to the direction of the reagent O atoms), with a smaller isotropic component. The product translational energy distribution, P(E_T), had a relatively low average of E_T = 35 kcal mol^(–1), indicating that the ^(16)O^(12)C^(18)O products were formed with substantial internal energy. The QCT calculations give c.m. P(E_T) and T(θ_(c.m.)) distributions and a relative product yield that agree qualitatively with the experimental results, and the trajectories indicate that exchange occurs through a short-lived CO_3^* intermediate. A low yield for the abstraction reaction was seen in both the experiment and the theory. Experimentally, a fast and weak ^(16)O^(18)O product signal from an abstraction reaction was observed, which could only be detected in the forward direction. A small number of QCT trajectories leading to abstraction were observed to occur primarily via a transient CO_3 intermediate, albeit only at high collision energies (149 kcal mol^(–1)). The oxygen isotope exchange mechanism for CO_2 in collisions with ground state O atoms is a newly discovered pathway through which oxygen isotopes may be cycled in the upper atmosphere, where O(^3P) atoms with hyperthermal translational energies can be generated by photodissociation of O_3 and O_2

Performance of ab initio and density functional methods for conformational equilibria of CnH2n+2 alkane isomers (n=2-8)

Conformational energies of n-butane, n-pentane, and n-hexane have been calculated at the CCSD(T) level and at or near the basis set limit. Post-CCSD(T) contribution were considered and found to be unimportant. The data thus obtained were used to assess the performance of a variety of density functional methods. Double-hybrid functionals like B2GP-PLYP and B2K-PLYP, especially with a small Grimme-type empirical dispersion correction, are capable of rendering conformational energies of CCSD(T) quality. These were then used as a `secondary standard' for a larger sample of alkanes, including isopentane and the branched hexanes as well as key isomers of heptane and octane. Popular DFT functionals like B3LYP, B3PW91, BLYP, PBE, and PBE0 tend to overestimate conformer energies without dispersion correction, while the M06 family severely underestimates GG interaction energies. Grimme-type dispersion corrections for these overcorrect and lead to qualitatively wrong conformer orderings. All of these functionals also exhibit deficiencies in the conformer geometries, particularly the backbone torsion angles. The PW6B95 and, to a lesser extent, BMK functionals are relatively free of these deficiencies. Performance of these methods is further investigated to derive conformer ensemble corrections to the enthalpy function, $H_{298}-H_0$, and the Gibbs energy function, ${\rm gef}(T)\equiv - [G(T)-H_0]/T$, of these alkanes. While $H_{298}-H_0$ is only moderately sensitive to the level of theory, ${\rm gef}(T)$ exhibits more pronounced sensitivity. Once again, double hybrids acquit themselves very well.Comment: J. Phys. Chem. A, revised [Walter Thiel festschrift

A Halomethane thermochemical network from iPEPICO experiments and quantum chemical calculations

Internal energy selected halomethane cations CH3Cl+, CH2Cl2+, CHCl3+, CH3F+, CH2F2+, CHClF2+ and CBrClF2+ were prepared by vacuum ultraviolet photoionization, and their lowest energy dissociation channel studied using imaging photoelectron photoion coincidence spectroscopy (iPEPICO). This channel involves hydrogen atom loss for CH3F+, CH2F2+ and CH3Cl+, chlorine atom loss for CH2Cl2+, CHCl3+ and CHClF2+, and bromine atom loss for CBrClF2+. Accurate 0 K appearance energies, in conjunction with ab initio isodesmic and halogen exchange reaction energies, establish a thermochemical network, which is optimized to update and confirm the enthalpies of formation of the sample molecules and their dissociative photoionization products. The ground electronic states of CHCl3+, CHClF2+ and CBrClF2+ do not confirm to the deep well assumption, and the experimental breakdown curve deviates from the deep well model at low energies. Breakdown curve analysis of such shallow well systems supplies a satisfactorily succinct route to the adiabatic ionization energy of the parent molecule, particularly if the threshold photoelectron spectrum is not resolved and a purely computational route is unfeasible. The ionization energies have been found to be 11.47 ± 0.01 eV, 12.30 ± 0.02 eV and 11.23 ± 0.03 eV for CHCl3, CHClF2 and CBrClF2, respectively. The updated 0 K enthalpies of formation, ∆fHo0K(g) for the ions CH2F+, CHF2+, CHCl2+, CCl3+, CCl2F+ and CClF2+ have been derived to be 844.4 ± 2.1, 601.6 ± 2.7, 890.3 ± 2.2, 849.8 ± 3.2, 701.2 ± 3.3 and 552.2 ± 3.4 kJ mol–1, respectively. The ∆fHo0K(g) values for the neutrals CCl4, CBrClF2, CClF3, CCl2F2 and CCl3F and have been determined to be –94.0 ± 3.2, –446.6 ± 2.7, –702.1 ± 3.5, –487.8 ± 3.4 and –285.2 ± 3.2 kJ mol–1, respectively

Correct quantum chemistry in a minimal basis from effective Hamiltonians

We describe how to create ab-initio effective Hamiltonians that qualitatively describe correct chemistry even when used with a minimal basis. The Hamiltonians are obtained by folding correlation down from a large parent basis into a small, or minimal, target basis, using the machinery of canonical transformations. We demonstrate the quality of these effective Hamiltonians to correctly capture a wide range of excited states in water, nitrogen, and ethylene, and to describe ground and excited state bond-breaking in nitrogen and the chromium dimer, all in small or minimal basis sets

Benchmark thermochemistry of the C_nH_{2n+2} alkane isomers (n=2--8) and performance of DFT and composite ab initio methods for dispersion-driven isomeric equilibria

The thermochemistry of linear and branched alkanes with up to eight carbons has been reexamined by means of W4, W3.2lite and W1h theories. `Quasi-W4' atomization energies have been obtained via isodesmic and hypohomodesmotic reactions. Our best atomization energies at 0 K (in kcal/mol) are: 1220.04 n-butane, 1497.01 n-pentane, 1774.15 n-hexane, 2051.17 n-heptane, 2328.30 n-octane, 1221.73 isobutane, 1498.27 isopentane, 1501.01 neopentane, 1775.22 isohexane, 1774.61 3-methylpentane, 1775.67 diisopropyl, 1777.27 neohexane, 2052.43 isoheptane, 2054.41 neoheptane, 2330.67 isooctane, and 2330.81 hexamethylethane. Our best estimates for $\Delta H^\circ_{f,298K}$ are: -30.00 n-butane, -34.84 n-pentane, -39.84 n-hexane, -44.74 n-heptane, -49.71 n-octane, -32.01 isobutane, -36.49 isopentane, -39.69 neopentane, -41.42 isohexane, -40.72 3-methylpentane, -42.08 diisopropyl, -43.77 neohexane, -46.43 isoheptane, -48.84 neoheptane, -53.29 isooctane, and -53.68 hexamethylethane. These are in excellent agreement (typically better than 1 kJ/mol) with the experimental heats of formation at 298 K obtained from the CCCBDB and/or NIST Chemistry WebBook databases. However, at 0 K a large discrepancy between theory and experiment (1.1 kcal/mol) is observed for only neopentane. This deviation is mainly due to the erroneous heat content function for neopentane used in calculating the 0 K CCCBDB value. The thermochemistry of these systems, especially of the larger alkanes, is an extremely difficult test for density functional methods. A posteriori corrections for dispersion are essential. Particularly for the atomization energies, the B2GP-PLYP and B2K-PLYP double-hybrids, and the PW6B95 hybrid-meta GGA clearly outperform other DFT functionals.Comment: (J. Phys. Chem. A, in press

Global hybrids from the semiclassical atom theory satisfying the local density linear response

We propose global hybrid approximations of the exchange-correlation (XC) energy functional which reproduce well the modified fourth-order gradient expansion of the exchange energy in the semiclassical limit of many-electron neutral atoms and recover the full local density approximation (LDA) linear response. These XC functionals represent the hybrid versions of the APBE functional [Phys. Rev. Lett. 106, 186406, (2011)] yet employing an additional correlation functional which uses the localization concept of the correlation energy density to improve the compatibility with the Hartree-Fock exchange as well as the coupling-constant-resolved XC potential energy. Broad energetical and structural testings, including thermochemistry and geometry, transition metal complexes, non-covalent interactions, gold clusters and small gold-molecule interfaces, as well as an analysis of the hybrid parameters, show that our construction is quite robust. In particular, our testing shows that the resulting hybrid, including 20\% of Hartree-Fock exchange and named hAPBE, performs remarkably well for a broad palette of systems and properties, being generally better than popular hybrids (PBE0 and B3LYP). Semi-empirical dispersion corrections are also provided.Comment: 12 pages, 4 figure

Dendrimer-Based Fluorescent Indicators: In Vitro and In Vivo Applications

BACKGROUND: The development of fluorescent proteins and synthetic molecules whose fluorescence properties are controlled by the environment makes it possible to monitor physiological and pathological events in living systems with minimal perturbation. A large number of small organic dyes are available and routinely used to measure biologically relevant parameters. Unfortunately their application is hindered by a number of limitations stemming from the use of these small molecules in the biological environment. PRINCIPAL FINDINGS: We present a novel dendrimer-based architecture leading to multifunctional sensing elements that can overcome many of these problems. Applications in vitro, in living cells and in vivo are reported. In particular, we image for the first time extracellular pH in the brain in a mouse epilepsy model. CONCLUSION: We believe that the proposed architecture can represent a useful and novel tool in fluorescence imaging that can be widely applied in conjunction with a broad range of sensing dyes and experimental setups