Hot carrier multiplication in plasmonic photocatalysis.

Light-induced hot carriers derived from the surface plasmons of metal nanostructures have been shown to be highly promising agents for photocatalysis. While both nonthermal and thermalized hot carriers can potentially contribute to this process, their specific role in any given chemical reaction has generally not been identified. Here, we report the observation that the H2-D2 exchange reaction photocatalyzed by Cu nanoparticles is driven primarily by thermalized hot carriers. The external quantum yield shows an intriguing S-shaped intensity dependence and exceeds 100% for high light intensities, suggesting that hot carrier multiplication plays a role. A simplified model for the quantum yield of thermalized hot carriers reproduces the observed kinetic features of the reaction, validating our hypothesis of a thermalized hot carrier mechanism. A quantum mechanical study reveals that vibrational excitations of the surface Cu-H bond is the likely activation mechanism, further supporting the effectiveness of low-energy thermalized hot carriers in photocatalyzing this reaction.

A spur to molecular geometry optimization: Gradient-enhanced universal kriging with on-the-fly adaptive ab initio prior mean functions in curvilinear coordinates.

We present a molecular geometry optimization algorithm based on the gradient-enhanced universal kriging (GEUK) formalism with ab initio prior mean functions, which incorporates prior physical knowledge to surrogate-based optimization. In this formalism, we have demonstrated the advantage of allowing the prior mean functions to be adaptive during geometry optimization over a pre-fixed choice of prior functions. Our implementation is general and flexible in two senses. First, the optimizations on the surrogate surface can be in both Cartesian coordinates and curvilinear coordinates. We explore four representative curvilinear coordinates in this work, including the redundant Coulombic coordinates, the redundant internal coordinates, the non-redundant delocalized internal coordinates, and the non-redundant hybrid delocalized internal Z-matrix coordinates. We show that our GEUK optimizer accelerates geometry optimization as compared to conventional non-surrogate-based optimizers in internal coordinates. We further showcase the power of the GEUK with on-the-fly adaptive priors for efficient optimizations of challenging molecules (Criegee intermediates) with a high-accuracy electronic structure method (the coupled-cluster method). Second, we present the usage of internal coordinates under the complete curvilinear scheme. A complete curvilinear scheme performs both surrogate potential-energy surface (PES) fitting and structure optimization entirely in the curvilinear coordinates. Our benchmark indicates that the complete curvilinear scheme significantly reduces the cost of structure minimization on the surrogate compared to the incomplete curvilinear scheme, which fits the surrogate PES in curvilinear coordinates partially and optimizes a structure in Cartesian coordinates through curvilinear coordinates via the chain rule.

Exploring torsional conformer space with physical prior mean function-driven meta-Gaussian processes.

We present a novel approach for systematically exploring the conformational space of small molecules with multiple internal torsions. Identifying unique conformers through a systematic conformational search is important for obtaining accurate thermodynamic functions (e.g., free energy), encompassing contributions from the ensemble of all local minima. Traditional geometry optimizers focus on one structure at a time, lacking transferability from the local potential-energy surface (PES) around a specific minimum to optimize other conformers. In this work, we introduce a physics-driven meta-Gaussian processes (meta-GPs) method that not only enables efficient exploration of target PES for locating local minima but, critically, incorporates physical surrogates that can be applied universally across the optimization of all conformers of the same molecule. Meta-GPs construct surrogate PESs based on the optimization history of prior conformers, dynamically selecting the most suitable prior mean function (representing prior knowledge in Bayesian learning) as a function of the optimization progress. We systematically benchmarked the performance of multiple GP variants for brute-force conformational search of amino acids. Our findings highlight the superior performance of meta-GPs in terms of efficiency, comprehensiveness of conformer discovery, and the distribution of conformers compared to conventional non-surrogate optimizers and other non-meta-GPs. Furthermore, we demonstrate that by concurrently optimizing, training GPs on the fly, and learning PESs, meta-GPs exhibit the capacity to generate high-quality PESs in the torsional space without extensive training data. This represents a promising avenue for physics-based transfer learning via meta-GPs with adaptive priors in exploring torsional conformer space.

Large Pressure Effects Caused by Internal Rotation in the s-cis-syn-Acrolein Stabilized Criegee Intermediate at Tropospheric Temperature and Pressure.

Criegee intermediates are important atmospheric oxidants, and quantitative kinetics for stabilized Criegee intermediates are key parameters for atmospheric modeling but are still limited. Here we report barriers and rate constants for unimolecular reactions of s-cis-syn-acrolein oxide (scsAO), in which the vinyl group makes it a prototype for Criegee intermediates produced in the ozonolysis of isoprene. We find that the MN15-L and M06-2X density functionals have CCSD(T)/CBS accuracy for the unimolecular cyclization and stereoisomerization of scsAO. We calculated high-pressure-limit rate constants by the dual-level strategy that combines (a) high-level wave function-based conventional transition-state theory (which includes coupled-cluster calculations with quasiperturbative inclusion of quadruple excitations because of the strongly multiconfigurational character of the electronic wave function) and (b) canonical variational transition-state theory with small-curvature tunneling based on a validated density functional. We calculated pressure-dependent rate constants both by system-specific quantum Rice-Ramsperger-Kassel theory and by solving the master equation. We report rate constants for unimolecular reactions of scsAO over the full range of atmospheric temperature and pressure. We found that the unimolecular reaction rates of this larger-than-previously studied Criegee intermediate depend significantly on pressure. Particularly, we found that falloff effects decrease the effective unimolecular cyclization rate constant of scsAO by about a factor of 3, but the unimolecular reaction is still the dominant atmospheric sink for scsAO at low altitudes. The large falloff caused by the inclusion of the stereoisomerization channel in the master equation calculations has broad implications for mechanistic analysis of reactions with competitive internal rotations that can produce stable rotamers.

Reaction of SO3 with HONO2 and Implications for Sulfur Partitioning in the Atmosphere.

Sulfur trioxide is a critical intermediate for the sulfur cycle and the formation of sulfuric acid in the atmosphere. The traditional view is that sulfur trioxide is removed by water vapor in the troposphere. However, the concentration of water vapor decreases significantly with increasing altitude, leading to longer atmospheric lifetimes of sulfur trioxide. Here, we utilize a dual-level strategy that combines transition state theory calculated at the W2X//DF-CCSD(T)-F12b/jun'-cc-pVDZ level, with variational transition state theory with small-curvature tunneling from direct dynamics calculations at the M08-HX/MG3S level. We also report the pressure-dependent rate constants calculated using the system-specific quantum Rice-Ramsperger-Kassel (SS-QRRK) theory. The present findings show that falloff effects in the SO3 + HONO2 reaction are pronounced below 1 bar. The SO3 + HONO2 reaction can be a potential removal reaction for SO3 in the stratosphere and for HONO2 in the troposphere, because the reaction can potentially compete well with the SO3 + 2H2O reaction between 25 and 35 km, as well as the OH + HONO2 reaction. The present findings also suggest an unexpected new product from the SO3 + HONO2 reaction, which, although very short-lived, would have broad implications for understanding the partitioning of sulfur in the stratosphere and the potential for the SO3 reaction with organic acids to generate organosulfates without the need for heterogeneous chemistry.

First-Principles Insights into Plasmon-Induced Catalysis.

The size- and shape-controlled enhanced optical response of metal nanoparticles (NPs) is referred to as a localized surface plasmon resonance (LSPR). LSPRs result in amplified surface and interparticle electric fields, which then enhance light absorption of the molecules or other materials coupled to the metallic NPs and/or generate hot carriers within the NPs themselves. When mediated by metallic NPs, photocatalysis can take advantage of this unique optical phenomenon. This review highlights the contributions of quantum mechanical modeling in understanding and guiding current attempts to incorporate plasmonic excitations to improve the kinetics of heterogeneously catalyzed reactions. A range of first-principles quantum mechanics techniques has offered insights, from ground-state density functional theory (DFT) to excited-state theories such as multireference correlated wavefunction methods. Here we discuss the advantages and limitations of these methods in the context of accurately capturing plasmonic effects, with accompanying examples.

Geometry meta-optimization.

Recent work has demonstrated the promise of using machine-learned surrogates, in particular, Gaussian process (GP) surrogates, in reducing the number of electronic structure calculations (ESCs) needed to perform surrogate model based (SMB) geometry optimization. In this paper, we study geometry meta-optimization with GP surrogates where a SMB optimizer additionally learns from its past "experience" performing geometry optimization. To validate this idea, we start with the simplest setting where a geometry meta-optimizer learns from previous optimizations of the same molecule with different initial-guess geometries. We give empirical evidence that geometry meta-optimization with GP surrogates is effective and requires less tuning compared to SMB optimization with GP surrogates on the ANI-1 dataset of off-equilibrium initial structures of small organic molecules. Unlike SMB optimization where a surrogate should be immediately useful for optimizing a given geometry, a surrogate in geometry meta-optimization has more flexibility because it can distribute its ESC savings across a set of geometries. Indeed, we find that GP surrogates that preserve rotational invariance provide increased marginal ESC savings across geometries. As a more stringent test, we also apply geometry meta-optimization to conformational search on a hand-constructed dataset of hydrocarbons and alcohols. We observe that while SMB optimization and geometry meta-optimization do save on ESCs, they also tend to miss higher energy conformers compared to standard geometry optimization. We believe that further research into characterizing the divergence between GP surrogates and potential energy surfaces is critical not only for advancing geometry meta-optimization but also for exploring the potential of machine-learned surrogates in geometry optimization in general.

Atmospheric Kinetics: Bimolecular Reactions of Carbonyl Oxide by a Triple-Level Strategy.

Criegee intermediates in the atmosphere serve as oxidizing agents to initiate aerosol formation, which are particularly important for atmospheric modeling, and understanding their kinetics is one of the current outstanding challenges in climate change modeling. Because experimental kinetics are still limited, we must rely on theory for the complete picture, but obtaining absolute rates from theory is a formidable task. Here, we report the bimolecular reaction kinetics of carbonyl oxide with ammonia, hydrogen sulfide, formaldehyde, and water dimer by designing a triple-level strategy that combines (i) benchmark results close to the complete-basis limit of coupled-cluster theory with the single, double, triple, and quadruple excitations (CCSDTQ/CBS), (ii) a new hybrid meta density functional (M06CR) specifically optimized for reactions of Criegee intermediates, and (iii) variational transition-state theory with both variable rection coordinates and optimized reaction paths, with multidimensional tunneling, and with pressure effects. For (i) we have found that quadruple excitations are required to obtain quantitative reaction barriers, and we designed new composite methods and strategies to reach CCSDTQ/CBS accuracy. The present findings show that (i) the CH2OO + HCHO reaction can make an important contribution to the sink of HCHO under wide atmospheric conditions in the gas phase and that (ii) CH2OO + (H2O)2 dominates over the CH2OO + H2O reaction below 10 km.

Unimolecular reaction of acetone oxide and its reaction with water in the atmosphere.

Criegee intermediates (i.e., carbonyl oxides with two radical sites) are known to be important atmospheric reagents; however, our knowledge of their reaction kinetics is still limited. Although experimental methods have been developed to directly measure the reaction rate constants of stabilized Criegee intermediates, the experimental results cover limited temperature ranges and do not completely agree well with one another. Here we investigate the unimolecular reaction of acetone oxide [(CH3)2COO] and its bimolecular reaction with H2O to obtain rate constants with quantitative accuracy comparable to experimental accuracy. We do this by using CCSDT(Q)/CBS//CCSD(T)-F12a/DZ-F12 benchmark results to select and validate exchange-correlation functionals, which are then used for direct dynamics calculations by variational transition state theory with small-curvature tunneling and torsional and high-frequency anharmonicity. We find that tunneling is very significant in the unimolecular reaction of (CH3)2COO and its bimolecular reaction with H2O. We show that the atmospheric lifetimes of (CH3)2COO depend on temperature and that the unimolecular reaction of (CH3)2COO is the dominant decay mode above 240 K, while the (CH3)2COO + SO2 reaction can compete with the corresponding unimolecular reaction below 240 K when the SO2 concentration is 9 × 1010 molecules per cubic centimeter. We also find that experimental results may not be sufficiently accurate for the unimolecular reaction of (CH3)2COO above 310 K. Not only does the present investigation provide insights into the decay of (CH3)2COO in the atmosphere, but it also provides an illustration of how to use theoretical methods to predict quantitative rate constants of medium-sized Criegee intermediates.

Enabling Lithium Metal Anode in Nonflammable Phosphate Electrolyte with Electrochemically Induced Chemical Reactions.

Lithium metal anode holds great promises for next-generation battery technologies but is notoriously difficult to work with. The key to solving this challenge is believed to lie in the ability of forming stable solid-electrolyte interphase (SEI) layers. To further address potential safety issues, it is critical to achieve this goal in nonflammable electrolytes. Building upon previous successes in forming stable SEI in conventional carbonate-based electrolytes, here we report that reversible Li stripping/plating could be realized in triethyl phosphate (TEP), a known flame retardant. The critical enabling factor of our approach was the introduction of oxygen, which upon electrochemical reduction induces the initial decomposition of TEP and produces Li3 PO4 and poly-phosphates. Importantly, the reaction was self-limiting, and the resulting material regulated Li plating by limiting dendrite formation. In effect, we obtained a functional SEI on Li metal in a nonflammable electrolyte. When tested in a symmetric Li∥Li cell, more than 300 cycles of stripping/plating were measured at a current density of 0.5 mA cm-2 . Prototypical Li-O2 and Li-ion batteries were also fabricated and tested to further support the effectiveness of this strategy. The mechanism by which the SEI forms was studied by density functional theory (DFT), and the predictions were corroborated by the successful detection of the intermediates and products.

Predicting Bond Dissociation Energies and Bond Lengths of Coordinatively Unsaturated Vanadium-Ligand Bonds.

Understanding the electronic structure of coordinatively unsaturated transition-metal compounds and predicting their physical properties are of great importance for catalyst design. Bond dissociation energy De and bond length re are two of the fundamental quantities for which good predictions are important for a successful design strategy. In the present work, recent experimentally measured bond energies and bond lengths of VX diatomic molecules (X = C, N, S) are used as a gauge to consider the utility of a number of electronic structure methods. Single-reference methods are one focus because of their efficiency and utility in practical calculations, and multireference configuration interaction (MRCISD) methods and a composite coupled cluster (CCC) method are a second focus because of their potential high accuracy. The comparison is especially challenging because of the large multireference M diagnostics of these molecules, in the range 0.15-0.19. For the single-reference methods, Kohn-Sham density functional theory (KS-DFT) has been tested with a variety of approximate exchange-correlation functionals. Of these, MOHLYP provides the bond dissociation energies in best agreement with experiments, and BLYP provides the bond lengths that are in best agreement with experiments; but by requiring good performance for both the De and re of the vanadium compounds, MOHLYP, MN12-L, MGGA_MS1, MGGA_MS0, O3LYP, and M06-L are the most highly recommended functionals. The CCC calculations include up to connected pentuple excitations for the valence electrons and up to connected quadruple excitations for the core-valence terms; this results in highly accurate dissociation energies and good bond lengths. Averaged over the three molecules, the mean unsigned deviation of CCC bond energies from experimental ones is only 0.4 kcal/mol, demonstrating excellent convergence of theory and experiments.

Rationalizing the Hot-Carrier-Mediated Reaction Mechanisms and Kinetics for Ammonia Decomposition on Ruthenium-Doped Copper Nanoparticles.

Localized surface plasmons generated on metallic nanostructures provide an efficient driving force for catalyzing chemical reactions, the kinetics of which cannot be understood properly by means of density functional theory, despite its wide use in simulating heterogeneous catalytic reaction mechanisms. Herein we report reaction pathways for the ammonia decomposition reaction on ruthenium-doped copper studied by the embedded correlated wavefunction method. Our computations provide a qualitative explanation for the experimentally observed change in the reaction order from thermal catalysis to hot-carrier-mediated photocatalysis, as reported very recently in Zhou, L.; et al. Science 2018 , 362 , 69 .

Kinetics of the Strongly Correlated CH3O + O2 Reaction: The Importance of Quadruple Excitations in Atmospheric and Combustion Chemistry.

Kinetics measurements on radical-radical reactions are often unavailable experimentally, and obtaining quantitative rate constants for such reactions by theoretical methods is challenging because the transition states and the reactants are often strongly correlated. Treating strongly correlated systems by coupled cluster theory limited to single, double, and triple connected excitations is often inadequate. We therefore use a new method, called GMM(P), for extrapolation to the complete configuration interaction limit to go beyond triple excitations and in particular to approximate the CCSDTQ(P)/CBS limit. Here, we present this method and use it to investigate the CH3O + O2 reaction. The contribution of connected quadruple excitations to the barrier height energy is found to be -3.13 kcal/mol, and adding a quasiperturbative calculation of the effect of connected pentuple excitations brings the post-connected-triples contributions to -3.44 kcal/mol, which corresponds to Boltzmann factors that increase calculated rate constants by factors of 1.0 × 103, 3.3 × 102, and 18 at 250, 298, and 600 K, respectively. We present rate constants for temperatures from 250 to 2000 K, and we find that the Arrhenius activation energy increases from 0.58 to 9.68 kcal/mol over this range. We also find reasonably good accuracy for the barrier height with the MN15-L exchange-correlation functional, and we calculate rate constants by a combination of GMM(P) and MN15-L electronic structure calculations and conventional and variational transition state theory, in particular canonical variational theory with small-curvature tunneling. The present findings have broad implications for obtaining quantitative rate constants for complex reaction systems in atmospheric and combustion chemistry.

Barrierless association of CF2 and dissociation of C2F4 by variational transition-state theory and system-specific quantum Rice-Ramsperger-Kassel theory.

Bond dissociation is a fundamental chemical reaction, and the first principles modeling of the kinetics of dissociation reactions with a monotonically increasing potential energy along the dissociation coordinate presents a challenge not only for modern electronic structure methods but also for kinetics theory. In this work, we use multifaceted variable-reaction-coordinate variational transition-state theory (VRC-VTST) to compute the high-pressure limit dissociation rate constant of tetrafluoroethylene (C2F4), in which the potential energies are computed by direct dynamics with the M08-HX exchange correlation functional. To treat the pressure dependence of the unimolecular rate constants, we use the recently developed system-specific quantum Rice-Ramsperger-Kassel theory. The calculations are carried out by direct dynamics using an exchange correlation functional validated against calculations that go beyond coupled-cluster theory with single, double, and triple excitations. Our computed dissociation rate constants agree well with the recent experimental measurements.

Relative Rates of Hydrogen Shift Isomerizations Depend Strongly on Multiple-Structure Anharmonicity.

Hydroperoxyalkylperoxy species (OOQOOH) are important intermediates that are generated during the autoignition of transport fuels. A key reaction of hydroperoxyalkylperoxy radicals is a [1,5] hydrogen shift, for which kinetics data are experimentally unavailable. Here we study two typical OOQOOH reactions and compare their kinetics to one another and to a previous study to learn the effect of structural variations of the alkyl group on the competition between alternative [1,5] hydrogen shifts of hydroperoxyalkylperoxy species. We use electronic structure calculations to determine previously missing thermochemical data, and we use variational transition state theory with multidimensional tunneling, multiple structures, torsional potential anharmonicity, and high-frequency anharmonicity to obtain more accurate rate constants than the ones that can be computed by conventional single-structure harmonic transition state theory and than the empirically estimated rate constants that are currently used in combustion modeling. The calculated temperature range is 298-1500 K. The roles of various factors in determining the rates are elucidated, and we find an especially strong effect of multiple structure anharmonicity due to torsions. Thus, even though there is some cancellation between the anharmonicity of the reactant and the anharmonicity of the transition state, and even though the reactants are very similar in structure, differing only by a methyl group, the effect of multiple structure anharmonicity has a large effect on the relative rates, as large as a factor of 17 at room temperature and as large as a factor of 7 at 1500 K. This has broad implications for the estimation of reaction rates in many subfields of chemistry, including combustion chemistry and atmospheric chemistry, where rates of reaction of complex molecules are usually estimated without explicit consideration of this fundamental entropic effect. In addition, the pressure-dependence of the rate constants is modeled by system-specific quantum Rice-Ramsperger-Kassel theory for a reversible isomerization.

Multiconfiguration Pair-Density Functional Theory: A New Way To Treat Strongly Correlated Systems.

The electronic energy of a system provides the Born-Oppenheimer potential energy for internuclear motion and thus determines molecular structure and spectra, bond energies, conformational energies, reaction barrier heights, and vibrational frequencies. The development of more efficient and more accurate ways to calculate the electronic energy of systems with inherently multiconfigurational electronic structure is essential for many applications, including transition metal and actinide chemistry, systems with partially broken bonds, many transition states, and most electronically excited states. Inherently multiconfigurational systems are called strongly correlated systems or multireference systems, where the latter name refers to the need for using more than one ("multiple") configuration state function to provide a good zero-order reference wave function. This Account describes multiconfiguration pair-density functional theory (MC-PDFT), which was developed as a way to combine the advantages of wave function theory (WFT) and density functional theory (DFT) to provide a better treatment of strongly correlated systems. First we review background material: the widely used Kohn-Sham DFT (which uses only a single Slater determinant as reference wave function), multiconfiguration WFT methods that treat inherently multiconfigurational systems based on an active space, and previous attempts to combine multiconfiguration WFT with DFT. Then we review the formulation of MC-PDFT. It is a generalization of Kohn-Sham DFT in that the electron kinetic energy and classical electrostatic energy are calculated from a reference wave function, while the rest of the energy is obtained from a density functional. However, there are two main differences with respent to Kohn-Sham DFT: (i) The reference wave function is multiconfigurational rather than being a single Slater determinant. (ii) The density functional is a function of the total density and the on-top pair density rather than being a function of the spin-up and spin-down densities. In work carried out so far, the multiconfigurational wave function is a multiconfiguration self-consistent-field wave function. The new formulation has the advantage that the reference wave function has the correct spatial and spin symmetry and can describe bond dissociation (of both single and multiple bonds) and electronic excitations in a formally and physically correct way. We then review the formulation of density functionals in terms of the on-top pair density. Finally we review successful applications of the theory to bond energies and bond dissociation potential energy curves of main-group and transition metal bonds, to barrier heights (including pericyclic reactions), to proton affinities, to the hydrogen bond energy of water dimer, to ground- and excited-state charge transfer, to valence and Rydberg excitations of molecules, and to singlet-triplet splittings of radicals. We find that that MC-PDFT can give accurate results not only with complete-active-space multiconfiguration wave functions but also with generalized-active-space multiconfiguration wave functions, which are practical for larger numbers of active electrons and active orbitals than are complete-active-space wave functions. The separated-pair approximation, which is a special case of generalized active space self-consistent-field theory, is especially promising. MC-PDFT, because it requires much less computer time and storage than pure WFT methods, has the potential to open larger and more complex strongly correlated systems to accurate simulation.

Effect of energy dependence of the density of states on pressure-dependent rate constants.

The FE integral for the normalized Boltzmann-weighted number of unimolecular states above the threshold energy is the key quantity for computing the collision efficiency in the pressure-dependent unimolecular rate theory developed by Troe, who calls this the energy dependence factor of the density of states. By using the Whitten-Rabinovitch approximation and assuming that the Whitten-Rabinovitch a(E) function is independent of energy, FE can be approximated by an analytical formula; this approximate formula is widely used because of its convenience and computational efficiency. Here we test its validity by comparing the rate constants computed by using the approximate FE to the ones determined by using the numerically integrated FE. For small-sized molecules and for reactions with high threshold energies E0, the differences are negligible at all temperatures, but in other cases, the approximate formula tends to underestimate FE and thus overestimates the collision efficiency, and this leads to smaller pressure falloff. When a(E) at high energies differs appreciably from a(E0), we find that the underestimation of pressure-dependent rate constants by using the approximate formula can be greater than a factor of 5 at high temperatures. The physical insight we draw from this study is that, for reactions with threshold energies below about 30 kcal mol-1, the rate of collisional energy transfer can be appreciably slowed down by the increase in the density of states at higher energies, and this increases the falloff effect by which finite-pressure rate constants become lower than the high-pressure limit, especially at higher temperatures.

How well can density functional theory and pair-density functional theory predict the correct atomic charges for dissociation and accurate dissociation energetics of ionic bonds?

The accuracy of density functional theory (DFT) is often judged by predicted dissociation energies, but one should also consider charge densities as illustrated here for dissociation of heteronuclear diatomic molecules, including ionic bonds for which local density functionals yield erroneous results. Some hybrid density functionals with 100% exact exchange in Kohn-Sham DFT and the local functionals in multiconfiguration pair-density functional theory give relatively acurate dissociation energies for NaCl, and they correctly yield uncharged dissociated atoms.

Variational transition state theory: theoretical framework and recent developments.

This article reviews the fundamentals of variational transition state theory (VTST), its recent theoretical development, and some modern applications. The theoretical methods reviewed here include multidimensional quantum mechanical tunneling, multistructural VTST (MS-VTST), multi-path VTST (MP-VTST), both reaction-path VTST (RP-VTST) and variable reaction coordinate VTST (VRC-VTST), system-specific quantum Rice-Ramsperger-Kassel theory (SS-QRRK) for predicting pressure-dependent rate constants, and VTST in the solid phase, liquid phase, and enzymes. We also provide some perspectives regarding the general applicability of VTST.

On the Upper Limits of Oxidation States in Chemistry.

The concept of oxidation state (OS) is based on the concept of Lewis electron pairs, in which the bonding electrons are assigned to the more electronegative element. This approach is useful for keeping track of the electrons, predicting chemical trends, and guiding syntheses. Experimental and quantum-chemical results reveal a limit near +8 for the highest OS in stable neutral chemical substances under ambient conditions. OS=+9 was observed for the isolated [IrO4 ]+ cation in vacuum. The prediction of OS=+10 for isolated [PtO4 ]2+ cations is confirmed computationally for low temperatures only, but hasn't yet been experimentally verified. For high OS species, oxidation of the ligands, for example, of O-2 with formation of . O-1 and O-O bonds, and partial reduction of the metal center may be favorable, possibly leading to non-Lewis type structures.