Showing papers by "Martin Head-Gordon published in 2018"

Mechanism of CO2 Reduction at Copper Surfaces: Pathways to C2 Products

Abstract: On the basis of constraints from reported experimental observations and density functional theory simulations, we propose a mechanism for the reduction of CO2 to C2 products on copper electrodes. To model the effects of an applied potential bias on the reactions, calculations are carried out with a variable, fractional number of electrons on the unit cell, which is optimized so that the Fermi level matches the actual chemical potential of electrons (i.e., the applied bias); an implicit electrolyte model allows for compensation of the surface charge so that neutrality is maintained in the overall simulation cell. Our mechanism explains the presence of the seven C2 species that have been detected in the reaction, as well as other notable experimental observations. Furthermore, our results shed light on the difference in activities toward C2 products between the (100) and (111) facets of copper. We compare our methodologies and findings with those in other recent mechanistic studies of the copper-catalyzed C...

Resonance-stabilized hydrocarbon-radical chain reactions may explain soot inception and growth

Abstract: Mystery surrounds the transition from gas-phase hydrocarbon precursors to terrestrial soot and interstellar dust, which are carbonaceous particles formed under similar conditions. Although polycyclic aromatic hydrocarbons (PAHs) are known precursors to high-temperature carbonaceous-particle formation, the molecular pathways that initiate particle formation are unknown. We present experimental and theoretical evidence for rapid molecular clustering–reaction pathways involving radicals with extended conjugation. These radicals react with other hydrocarbon species to form covalently bound complexes that promote further growth and clustering by regenerating resonance-stabilized radicals through low-barrier hydrogen-abstraction and hydrogen-ejection reactions. Such radical–chain reaction pathways may lead to covalently bound clusters of PAHs and other hydrocarbons that would otherwise be too small to condense at high temperatures, thus providing the key mechanistic steps for rapid particle formation and surface growth by hydrocarbon chemisorption.

How Accurate Is Density Functional Theory at Predicting Dipole Moments? An Assessment Using a New Database of 200 Benchmark Values

Abstract: Dipole moments are a simple, global measure of the accuracy of the electron density of a polar molecule. Dipole moments also affect the interactions of a molecule with other molecules as well as electric fields. To directly assess the accuracy of modern density functionals for calculating dipole moments, we have developed a database of 200 benchmark dipole moments, using coupled cluster theory through triple excitations, extrapolated to the complete basis set limit. This new database is used to assess the performance of 88 popular or recently developed density functionals. The results suggest that double hybrid functionals perform the best, yielding dipole moments within about 3.6–4.5% regularized RMS error versus the reference values—which is not very different from the 4% regularized RMS error produced by coupled cluster singles and doubles. Many hybrid functionals also perform quite well, generating regularized RMS errors in the 5–6% range. Some functionals, however, exhibit large outliers, and local f...

Survival of the most transferable at the top of Jacob’s ladder: Defining and testing the ωB97M(2) double hybrid density functional

Abstract: A meta-generalized gradient approximation, range-separated double hybrid (DH) density functional with VV10 non-local correlation is presented. The final 14-parameter functional form is determined by screening trillions of candidate fits through a combination of best subset selection, forward stepwise selection, and random sample consensus (RANSAC) outlier detection. The MGCDB84 database of 4986 data points is employed in this work, containing a training set of 870 data points, a validation set of 2964 data points, and a test set of 1152 data points. Following an xDH approach, orbitals from the ωB97M-V density functional are used to compute the second-order perturbation theory correction. The resulting functional, ωB97M(2), is benchmarked against a variety of leading double hybrid density functionals, including B2PLYP-D3(BJ), B2GPPLYP-D3(BJ), ωB97X-2(TQZ), XYG3, PTPSS-D3(0), XYGJ-OS, DSD-PBEP86-D3(BJ), and DSD-PBEPBE-D3(BJ). Encouragingly, the overall performance of ωB97M(2) on nearly 5000 data points clearly surpasses that of all of the tested density functionals. As a Rung 5 density functional, ωB97M(2) completes our family of combinatorially optimized functionals, complementing B97M-V on Rung 3, and ωB97X-V and ωB97M-V on Rung 4. The results suggest that ωB97M(2) has the potential to serve as a powerful predictive tool for accurate and efficient electronic structure calculations of main-group chemistry.

An assessment of strategies for the development of solid-state adsorbents for vehicular hydrogen storage

Abstract: Nanoporous adsorbents are a diverse category of solid-state materials that hold considerable promise for vehicular hydrogen storage. Although impressive storage capacities have been demonstrated for several materials, particularly at cryogenic temperatures, materials meeting all of the targets established by the U.S. Department of Energy have yet to be identified. In this Perspective, we provide an overview of the major known and proposed strategies for hydrogen adsorbents, with the aim of guiding ongoing research as well as future new storage concepts. The discussion of each strategy includes current relevant literature, strengths and weaknesses, and outstanding challenges that preclude implementation. We consider in particular metal–organic frameworks (MOFs), including surface area/volume tailoring, open metal sites, and the binding of multiple H2 molecules to a single metal site. Two related classes of porous framework materials, covalent organic frameworks (COFs) and porous aromatic frameworks (PAFs), are also discussed, as are graphene and graphene oxide and doped porous carbons. We additionally introduce criteria for evaluating the merits of a particular materials design strategy. Computation has become an important tool in the discovery of new storage materials, and a brief introduction to the benefits and limitations of computational predictions of H2 physisorption is therefore presented. Finally, considerations for the synthesis and characterization of hydrogen storage adsorbents are discussed.

Characterizing the interplay of Pauli repulsion, electrostatics, dispersion and charge transfer in halogen bonding with energy decomposition analysis.

Abstract: The halogen bond is a class of non-covalent interaction that has attracted considerable attention recently. A widespread theory for describing them is the σ-hole concept, which predicts that the strength of the interaction is proportional to the size of the σ-hole, a region of positive electrostatic potential opposite a σ bond. Previous work shows that in the case of CX3I, with X equal to F, Cl, Br, and I, the σ-hole trend is exactly opposite to the trend in binding energy with common electron pair donors. Using energy decomposition analysis (EDA) applied to a potential energy scan as well as the recent adiabatic EDA technique, we show that the observed trend is a result of charge transfer. Therefore a picture of the halogen bond that excludes charge transfer cannot be complete, and permanent and induced electrostatics do not always provide the dominant stabilizing contributions to halogen bonds. Overall, three universally attractive factors, polarization, dispersion and charge transfer, together with permanent electrostatics, which is usually attractive, drive halogen bonding, against Pauli repulsion.

Is Subsurface Oxygen Necessary for the Electrochemical Reduction of CO2 on Copper

Abstract: It has recently been proposed that subsurface oxygen is crucial for the adsorption and subsequent electroreduction of CO2 on copper. Using density functional theory, we have studied the stability and diffusion of subsurface oxygen in single crystals of copper exposing (111) and (100) facets. Oxygen is at least 1.5 eV more stable on the surface than beneath it for both crystal orientations; interstitial sites are too small to accommodate oxygen. The rate of atomic oxygen diffusion from one layer below a Cu(111) surface to the surface is 5 × 103 s–1. Oxygen can survive longer in deeper layers, but it does not promote CO2 adsorption there. Diffusion of subsurface oxygen is easier to the less-dense Cu(100) surface, even from lower layers (rate ≈ 1 × 107 s–1). Once the applied voltage and dispersion forces are properly modeled, we find that subsurface oxygen is unnecessary for CO2 adsorption on copper.

Regularized Orbital-Optimized Second-Order Møller–Plesset Perturbation Theory: A Reliable Fifth-Order-Scaling Electron Correlation Model with Orbital Energy Dependent Regularizers

Abstract: We derive and assess two new classes of regularizers that cope with offending denominators in the single-reference second-order Moller-Plesset perturbation theory (MP2). In particular, we discuss the use of two types of orbital energy dependent regularizers, κ and σ, in conjunction with orbital-optimized MP2 (OOMP2). The resulting fifth-order-scaling methods, κ-OOMP2 and σ-OOMP2, have been examined for bond-breaking, thermochemistry, nonbonded interactions, and biradical problems. Both methods with strong enough regularization restore restricted to unrestricted instability (i.e., Coulson-Fischer points) that unregularized OOMP2 lacks when breaking bonds in H2, C2H6, C2H4, and C2H2. The training of the κ and σ regularization parameters was performed with the W4-11 set. We further developed scaled correlation energy variants, κ-S-OOMP2 and σ-S-OOMP2, by training on the TAE140 subset of the W4-11 set. Those new OOMP2 methods were tested on the RSE43 set and the TA13 set where unmodified OOMP2 itself performs very well. The modifications we made were found insignificant in these data sets. Furthermore, we tested the new OOMP2 methods on singlet biradicaloids using Yamaguchi's approximate spin-projection. Unlike the unregularized OOMP2, which fails to converge these systems due to the singularity, we show that regularized OOMP2 methods successfully capture strong biradicaloid characters. While further assessment on larger data sets is desirable, κ-OOMP2 with κ = 1.45 E h-1 appears to combine favorable recovery of Coulson-Fischer points with good numerical performance.

How accurate are static polarizability predictions from density functional theory? An assessment over 132 species at equilibrium geometry

Abstract: Static polarizabilities are the first response of the electron density to electric fields, and are therefore important for predicting intermolecular and molecule-field interactions. They also offer a global measure of the accuracy of the treatment of excited states by density functionals in a formally exact manner. We have developed a database of benchmark static polarizabilities for 132 small species at equilibrium geometry, using coupled cluster theory through triple excitations (extrapolated to the complete basis set limit), for the purpose of developing and assessing density functionals. The performance of 60 popular and recent functionals are also assessed, which indicates that double hybrid functionals perform the best, having RMS relative errors in the range of 2.5-3.8%. Many hybrid functionals also give quite reasonable estimates with 4-5% RMS relative error. A few meta-GGAs like mBEEF and MVS yield performance comparable to hybrids, indicating potential for improved excited state predictions relative to typical local functionals. Some recent functionals however are found to be prone to catastrophic failure (possibly as a consequence of overparameterization), indicating a need for caution in applying these.

Computational Modeling of the Nature and Role of Ga Species for Light Alkane Dehydrogenation Catalyzed by Ga/H-MFI

Abstract: Ga-exchanged H-MFI zeolites are highly active for the dehydrogenation of light alkanes; however, both the nature of the active gallium species and the associated dehydrogenation mechanism have been difficult to establish. In this study, we examine the activity of Ga species in Ga/H-MFI by calculating the free energy landscapes on which all reactions occur. To this end, we use a hybrid quantum mechanics/molecular mechanics model for all electronic structure calculations. Quantum chemical calculations were carried out with a range-corrected functional and a good representation of dispersive interactions. The molecular mechanics part of our approach captures the long-range effects of Coulombic and dispersive interactions due to atoms in the extended framework. The rate-determining TS (RDTS) is identified by analysis of the free energy landscape for each mechanism, using the energetic span model. Our analysis reveals that, for reduced Ga/H-MFI, univalent and divalent gallium hydrides, [GaH2]+ and [GaH]2+, res...

Characterization of Isolated Ga3+ Cations in Ga/H-MFI Prepared by Vapor-Phase Exchange of H-MFI Zeolite with GaCl3

Abstract: Ga/H-MFI was prepared by vapor-phase reaction of GaCl3 with Bronsted acid O–H groups in dehydrated H-MFI zeolite. The resulting [GaCl2]+ cations in the as-exchanged zeolite are treated in H2 at 823...

Delocalization Errors in Density Functional Theory Are Essentially Quadratic in Fractional Occupation Number.

Abstract: Approximate functionals used in practical density functional theory (DFT) deviate from the piecewise linear behavior of the exact functional for fractional charges. This deviation causes excess charge delocalization, which leads to incorrect densities, molecular properties, barrier heights, band gaps, and excitation energies. We present a simple delocalization function for characterizing this error and find it to be almost perfectly linear vs the fractional electron number for systems spanning in size from the H atom to the C12H14 polyene. This causes the delocalization energy error to be a quadratic polynomial in the fractional electron number, which permits us to assess the comparative performance of 47 popular and recent functionals through the curvature. The quadratic form further suggests that information about a single fractional charge is sufficient to eliminate the principal source of delocalization error. Generalizing traditional two-point information like ionization potentials or electron affinities to account for a third, fractional charge-based data point could therefore permit fitting/tuning of functionals with lower delocalization error.

Non-orthogonal configuration interaction with single substitutions for the calculation of core-excited states

Abstract: In this paper, we present the non-orthogonal configuration interaction singles (NOCIS) method for calculating core-excited states of closed-shell molecules. NOCIS is a black-box variant of NOCI, which uses A different core-ionized determinants for a molecule with A atoms of a given element to form single substitutions. NOCIS is a variational, spin-pure, size-consistent ab initio method that dramatically improves on standard CIS by capturing essential orbital relaxation effects, in addition to essential configuration interaction. We apply it to the calculation of core-excitations for several smaller molecules and demonstrate that it performs competitively with other Hartree-Fock and DFT-based methods. We also benchmark it in several basis sets.

On the Computational Characterization of Charge-Transfer Effects in Noncovalently Bound Molecular Complexes.

Abstract: Charge-transfer (CT) is an important binding force in the formation of intermolecular complexes, and there have been a variety of theoretical models proposed to quantify this effect. These approaches, which typically rely on a definition of a “CT-free” state based on a partition of the system, sometimes yield significantly different results for a given intermolecular complex. Two widely used definitions of the “CT-free” state, the absolutely localized molecular orbitals (ALMO) method (where only on-fragment orbital mixings are permitted) and the constrained density functional theory (CDFT) approach (where fragment electron populations are fixed), are carefully examined in this work. Natural bond orbital (NBO) and the regularized symmetry-adapted perturbation theory (SAPT) are also briefly considered. Results for the ALMO and CDFT definitions of CT are compared on a broad range of model systems, including hydrogen-bonding systems, borane complexes, metal–carbonyl complexes, and complexes formed by water an...

Unraveling Substituent Effects on Frontier Orbitals of Conjugated Molecules Using an Absolutely Localized Molecular Orbital Based Analysis

Abstract: It is common to introduce electron-donating or electron-withdrawing substituent groups into functional conjugated molecules (such as dyes) to tune their electronic structure properties (such as frontier orbital energy levels) and photophysical properties (such as absorption and emission wavelengths). However, there lacks a generally applicable tool that can unravel the underlying interactions between orbitals from a substrate molecule and those from its substituents in modern electronic structure calculations, despite the long history of qualitative molecular orbital theory. In this work, the absolutely localized molecular orbitals (ALMO) based analysis is extended to analyze the effects of substituent groups on the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of a given system. This provides a bottom-up avenue towards quantification of effects from distinct physical origins (e.g. permanent electrostatics/Pauli repulsion, mutual polarization, inter-fragment orbital mixing). For the example case of prodan (a typical dye molecule), it is found that inter-fragment orbital mixing plays a key role in narrowing the HOMO-LUMO gap of the naphthalene core. Specifically, an out-of-phase mixing of high-lying occupied orbitals on the naphthalene core and the dimethylamino group leads to an elevated HOMO, whereas an in-phase combination of LUMOs on the naphthalene core and the propionyl group lowers the LUMO energy of the entire molecule. We expect this ALMO-based analysis to bridge the gap between concepts from qualitative orbital interaction analysis and quantitative electronic structure calculations.

Energy decomposition analysis for exciplexes using absolutely localized molecular orbitals.

Abstract: An energy decomposition analysis (EDA) scheme is developed for understanding the intermolecular interaction involving molecules in their excited states. The EDA utilizes absolutely localized molecular orbitals to define intermediate states and is compatible with excited state methods based on linear response theory such as configuration interaction singles and time-dependent density functional theory. The shift in excitation energy when an excited molecule interacts with the environment is decomposed into frozen, polarization, and charge transfer contributions, and the frozen term can be further separated into Pauli repulsion and electrostatics. These terms can be added to their counterparts obtained from the ground state EDA to form a decomposition of the total interaction energy. The EDA scheme is applied to study a variety of systems, including some model systems to demonstrate the correct behavior of all the proposed energy components as well as more realistic systems such as hydrogen-bonding complexes (e.g., formamide-water, pyridine/pyrimidine-water) and halide (F-, Cl-)-water clusters that involve charge-transfer-to-solvent excitations.

Communication: xDH double hybrid functionals can be qualitatively incorrect for non-equilibrium geometries: Dipole moment inversion and barriers to radical-radical association using XYG3 and XYGJ-OS.

Abstract: Double hybrid (DH) density functionals are amongst the most accurate density functional approximations developed so far, largely due to the incorporation of correlation effects from unoccupied orbitals via second order perturbation theory (PT2). The xDH family of DH functionals calculate energy directly from orbitals optimized by a lower level approach like B3LYP, without self-consistent optimization. XYG3 and XYGJ-OS are two widely used xDH functionals that are known to be quite accurate at equilibrium geometries. Here, we show that the XYG3 and XYGJ-OS functionals can be ill behaved for stretched bonds well beyond the Coulson-Fischer point, predicting unphysical dipole moments and humps in potential energy curves for some simple systems like the hydrogen fluoride molecule. Numerical experiments and analysis show that these failures are not due to PT2. Instead, a large mismatch at stretched bond-lengths between the reference B3LYP orbitals and the optimized orbitals associated with the non-PT2 part of XYG3 leads to an unphysically large non-Hellman-Feynman contribution to first order properties like forces and electron densities.

Nonempirical Meta-Generalized Gradient Approximations for Modeling Chemisorption at Metal Surfaces.

Abstract: We assess the accuracy of popular nonempirical GGAs (PBE, PBEsol, RPBE) and meta-GGAs (TPSS, revTPSS, and SCAN) for describing chemisorption reactions at metal surfaces. Except for RPBE, all the functionals tend to overbind the adsorbate significantly. We then propose a nonempirical meta-GGA, denoted as RTPSS, that is based on RPBE in the same way that TPSS is based on PBE. The RTPSS functional remedies the overbinding problem and improves the description of chemisorption energies. As an example of an application of RTPSS, we study the adsorption of CO on Cu surfaces (a notably difficult problem for semilocal functionals) and find that RTPSS is the only tested functional that predicts accurate chemisorption energies and the preferred adsorption site of CO. Although RTPSS gives an accurate description of chemisorption, nonlocal correlation may be necessary to describe physisorption if long-range van der Waals interactions are involved (however, this is true for semilocal functionals in general). We suggest...

Understanding Brønsted-Acid Catalyzed Monomolecular Reactions of Alkanes in Zeolite Pores by Combining Insights from Experiment and Theory

Abstract: Acidic zeolites are effective catalysts for the cracking of large hydrocarbon molecules into lower molecular weight products required for transportation fuels. However, the ways in which the zeolite structure affects the catalytic activity at Bronsted protons are not fully understood. One way to characterize the influence of the zeolite structure on the catalysis is to study alkane cracking and dehydrogenation at very low conversion, conditions for which the kinetics are well defined. To understand the effects of zeolite structure on the measured rate coefficient (kapp ), it is necessary to identify the equilibrium constant for adsorption into the reactant state (Kads-H+ ) and the intrinsic rate coefficient of the reaction (kint ) at reaction temperatures, since kapp is proportional to the product of Kads-H+ and kint . We show that Kads-H+ cannot be calculated from experimental adsorption data collected near ambient temperature, but can, however, be estimated accurately from configurational-bias Monte Carlo (CBMC) simulations. Using monomolecular cracking and dehydrogenation of C3 -C6 alkanes as an example, we review recent efforts aimed at elucidating the influence of the acid site location and the zeolite framework structure on the observed values of kapp and its components, Kads-H+ and kint .

Efficient Implementation of NOCI-MP2 Using the Resolution of the Identity Approximation with Application to Charged Dimers and Long C-C Bonds in Ethane Derivatives.

Abstract: An efficient implementation of the perturb-then-diagonalize nonorthogonal configuration interaction method with second-order Moller-Plesset perturbation theory (NOCI-MP2) is presented. Relative to other low scaling multireference perturbation theories, NOCI-MP2 often requires a much smaller active space because of the use of nonorthogonal reference configurations. Reworking the NOCI-MP2 equations with the resolution of the identity (RI) approximation enables the method to have the same memory requirements and computational scaling as single reference RI-MP2. The working equations are extended to include single substitutions as required when the reference determinants do not satisfy the Hartree-Fock equations. A detailed computational algorithm is presented along with timings to establish the performance of the implementation. NOCI-MP2 is applied to the binding energy and charge resonance energy in dication and monocation π dimers, as well as didiamantane ethane, and hexaphenylethane. A well-defined set of nonorthogonal determinants are obtained using absolutely localized molecular orbitals (ALMOs), as solutions to the self-consistent field for molecular interactions (SCF-MI) equations corresponding to covalent and ionic determinants. Agreement with experimental information where available, and other multireference methods, is satisfactory, with the use of an 0.3 au level shift to guard against large MP2 amplitudes. For didiamantane ethane and hexaphenylethane, large dispersion forces help stabilize the molecules despite the steric repulsion. By contrast, in the case of hexaphenylethane, the energy penalty from the geometric distortion of the fragments significantly weakens the bond.

Energy Decomposition Analysis for Excimers Using Absolutely Localized Molecular Orbitals within Time-Dependent Density Functional Theory and Configuration Interaction with Single Excitations.

Abstract: We present an improved energy decomposition analysis (EDA) scheme for understanding intermolecular interactions in delocalized excited states, especially in excimers. In the EDA procedure, excited states are treated with linear response theory such as configuration interaction singles (CIS) or time-dependent density functional theory (TDDFT), and absolutely localized molecular orbitals (ALMOs) are used to define the intermediate (frozen, excitonic coupling, and polarized) states. The intermolecular interaction energy is thereby separated into frozen, excitonic splitting, polarization, and charge transfer contributions. The excitonic splitting term describes the delocalization effect as two or more degenerate local excitations coupled with each other, which is often an important binding force in excimers. A maximum overlap state-tracking procedure is introduced to connect the initial fragment excitations to the constrained intermediate states and finally to the unconstrained delocalized states of the complex. The EDA scheme is applied to several excimer systems, including the He2* and Ne2* noble gas excimers, the doubly hydrogen-bonded 2-pyridone dimer, and the aromatic benzene and perylene excimers. We are able to gain some useful insights into the role each term is playing in the formation of these excimers, and the resulting method may also be useful for understanding a range of other complexes in excited states.

Open-shell coupled-cluster valence-bond theory augmented with an independent amplitude approximation for three-pair correlations: Application to a model oxygen-evolving complex and single molecular magnet

Abstract: We report the failure of coupled-cluster valence-bond (CCVB) theory with two-pair configurations [D. W. Small and M. Head-Gordon, J. Chem. Phys. 130, 084103 (2009)] for open-shell (OS) spin-frustrated systems where including three-pair configurations is necessary to properly describe strong spin-correlations. We extend OS-CCVB by augmenting the model with three-pair configurations within the independent amplitude approximation. The resulting new electronic structure model, OS-CCVB+i3, involves only a quadratic number of independent wavefunction parameters. It includes the recently reported closed-shell CCVB+i3 as a special case. Its cost is dominated by integral transformations, and it is capable of breaking multiple bonds exactly for all systems examined so far. The strength of OS-CCVB+i3 is highlighted in realistic systems including the [CaMn3O4] cubane subunit of the oxygen-evolving complex and a molecular magnet with the [Cr9] core unit as well as model systems such as N3, V3O3, and P5. We show that OS-CCVB+i3 is only slightly dependent on the underlying perfect-pairing reference, while OS-CCVB shows a stronger dependence. We also emphasize the compactness of the OS-CCVB+i3 wavefunction compared to the heat-bath configuration interaction wavefunction, a recently introduced soft exponential-scaling approach.

Orbital optimisation in the perfect pairing hierarchy: applications to full-valence calculations on linear polyacenes

Abstract: We describe the implementation of orbital optimisation for the models in the perfect pairing hierarchy. Orbital optimisation, which is generally necessary to obtain reliable results, is pursued at perfect pairing (PP) and perfect quadruples (PQ) levels of theory for applications on linear polyacenes, which are believed to exhibit strong correlation in the π space. While local minima and σ-π symmetry breaking solutions were found for PP orbitals, no such problems were encountered for PQ orbitals. The PQ orbitals are used for single-point calculations at PP, PQ and perfect hextuples (PH) levels of theory, both only in the π subspace, as well as in the full σπ valence space. It is numerically demonstrated that the inclusion of single excitations is necessary also when optimised orbitals are used. PH is found to yield good agreement with previously published density matrix renormalisation group data in the π space, capturing over 95% of the correlation energy. Full-valence calculations made possible b...

Reaction mechanism of the selective reduction of CO2 to CO by a tetraaza [CoIIN4H]2+ complex in the presence of protons

Abstract: The tetraaza [CoIIN4H]2+ complex (1) is remarkable for its ability to selectively reduce CO2 to CO with 45% Faradaic efficiency and a CO to H2 ratio of 3 : 2. We employ density functional theory (DFT) to determine the reasons behind the unusual catalytic properties of 1 and the most likely mechanism for CO2 reduction. The selectivity for CO2 over proton reduction is explained by analyzing the catalyst's affinity for the possible ligands present under typical reaction conditions: acetonitrile, water, CO2, and bicarbonate. After reduction of the catalyst by two electrons, formation of [CoIN4H]+-CO2- is strongly favored. Based on thermodynamic and kinetic data, we establish that the only likely route for producing CO from here consists of a protonation step to yield [CoIN4H]+-CO2H, followed by reaction with CO2 to form [CoIIN4H]2+-CO and bicarbonate. This conclusion corroborates the idea of a direct role of CO2 as a Lewis acid to assist in C-O bond dissociation, a conjecture put forward by other authors to explain recent experimental observations. The pathway to formic acid is predicted to be forbidden by high activation barriers, in accordance with the products that are known to be generated by 1. Calculated physical observables such as standard reduction potentials and the turnover frequency for our proposed catalytic cycle are in agreement with available experimental data reported in the literature. The mechanism also makes a prediction that may be experimentally verified: that the rate of CO formation should increase linearly with the partial pressure of CO2.

Independent amplitude approximations in coupled cluster valence bond theory: Incorporation of 3-electron-pair correlation and application to spin frustration in the low-lying excited states of a ferredoxin-type tetrametallic iron-sulfur cluster.

Abstract: Coupled cluster valence bond (CCVB) is a simple electronic structure method based on a perfect pairing (PP) reference with 2-pair recouplings for strong electron correlation problems. CCVB is spin-pure, size-consistent, and can exactly (in its active space) separate any molecule into atoms for which unrestricted Hartree-Fock (UHF) at dissociation is the sum of the ground state UHF energies of the atoms. However CCVB is far from a complete description of strong correlations. Its first failure to exactly describe spin-recouplings arises at the level of 3 electron pairs, such as the recoupling of 3 triplet oxygen atoms in the dissociation of singlet ozone. Such situations are often associated with spin frustration. To address this limitation, an extension of CCVB, termed CCVB+i3, is reported here that includes an independent (i) amplitude approximation to the 3-pair recouplings. CCVB+i3 thereby has the same basic computational requirements as those of CCVB, which has previously been shown to be an efficient method. CCVB+i3 correctly separates molecules that CCVB cannot. As a by-product, an independent 2-pair amplitude approximation to CCVB, called PP+i2, is also defined. Remarkably, PP+i2 can also correctly separate systems that CCVB cannot. CCVB+i3 is validated on the symmetric dissociation of D3h ozone. CCVB+i3 is then used to explore the role of 3-pair recouplings in an [Fe4S4(SCH3)4]2- cluster that has been used to model the iron-sulfur core of [Fe4S4] ferredoxins. Using localized PP orbitals, such recouplings are demonstrated to be large in some low-lying singlet excited states of the cluster. Significant 3 pair recoupling amplitudes include the usual triangular motif associated with spin frustration and other geometric arrangements of the 3 entangled pairs across the 4 iron centers.

Open-Shell Coupled-Cluster Valence-Bond Theory Augmented with an Independent Amplitude Approximation for Three-Pair Correlations: Application to a Model Oxygen-Evolving Complex and Single Molecular Magnet

Abstract: We report the failure of coupled-cluster valence-bond (CCVB) theory with two-pair configurations [J. Chem. Phys. 2009, 130, 084103 (2009)] for open-shell (OS) spin-frustrated systems where including three-pair configurations is necessary to properly describe strong spin-correlations. We extend OS-CCVB by augmenting the model with three-pair configurations within the independent amplitude approximation (IAA). The resulting new electronic structure model, OS-CCVB+i3, involves only a quadratic number of independent wavefunction parameters. It includes the recently reported closed-shell CCVB+i3 as a special case. Its cost is dominated by integral transformations and it is capable of breaking multiple bonds exactly for all systems examined so far. The strength of OS-CCVB+i3 is highlighted in realistic systems including the [CaMn$_3$O$_4$] cubane subunit of the oxygen-evolving complex and a molecular magnet with the [Cr$_9$] core unit as well as model systems such as N$_3$, V$_3$O$_3$, and P$_5$. We show that OS-CCVB+i3 is only slightly dependent on the underlying perfect-pairing reference while OS-CCVB shows a stronger dependence. We also emphasize the compactness of the OS-CCVB+i3 wavefunction compared to the heat-bath configuration interaction wavefunction, a recently introduced soft exponential-scaling approach.

Understanding Non-Covalent Interactions: Correlated Energy Decomposition Analysis and Applications to Halogen Bonding.

Abstract: Non-covalent interactions play a primordial role in chemistry. Beyond their quantification, the detailed understanding of their physical processes is necessary to rationalize chemical trends and improve designs of chemical systems. Energy decomposition analyses allow detailed insight into non-covalent interactions by extracting electrostatics, Pauli repulsion, polarization, dispersion and charge transfer components from interaction energies. Recent work has demonstrated that electronic correlation influenced significantly all of these energy components, whereas previous decompositions only partitioned correlation between dispersion and charge transfer. The MP2 energy decomposition analysis with Absolutely Localized Molecular Orbitals (MP2 ALMO-EDA) takes these results fully into account and offers a correlation correction for each extracted component. A recent detailed investigation of the CCSD dispersion energy showed that a small number of virtual orbitals is sufficient to describe dispersion interactions accurately in the long-range, which potentially offers a basis-set independent definition of dispersion. Finally, we present an application of MP2 ALMO-EDA to a series of unusual halogen bonding complexes where charge transfer dominates over the electrostatic σ-hole interaction.

Bimolecular Reaction Dynamics in the Phenyl-Silane System: Exploring the Prototype of a Radical Substitution Mechanism

Abstract: We present a combined experimental and theoretical investigation of the bimolecular gas-phase reaction of the phenyl radical (C6H5) with silane (SiH4) under single collision conditions to investigate the chemical dynamics of forming phenylsilane (C6H5SiH3) via a bimolecular radical substitution mechanism at a tetracoordinated silicon atom. Verified by electronic structure and quasiclassical trajectory calculations, the replacement of a single carbon atom in methane by silicon lowers the barrier to substitution, thus defying conventional wisdom that tetracoordinated hydrides undergo preferentially hydrogen abstraction. This reaction mechanism provides fundamental insights into the hitherto unexplored gas-phase chemical dynamics of radical substitution reactions of mononuclear main group hydrides under single collision conditions and highlights the distinct reactivity of silicon compared to its isovalent carbon. This mechanism might be also involved in the synthesis of cyanosilane (SiH3CN) and methylsilane (CH3SiH3) probed in the circumstellar envelope of the carbon star IRC+10216.

Erratum: "Non-orthogonal configuration interaction with single substitutions for the calculation of core-excited states" [J. Chem. Phys. 149, 044116 (2018)].

Abstract: Author(s): Oosterbaan, Katherine J; White, Alec F; Head-Gordon, Martin | Abstract: In our recent publication,1 we incorrectly stated some of the CIS and TDDFT k-edge excitations. Specifically the errors were in the CIS k-edge for all molecules except C2N2 and C2H6 and the TDDFT k-edge for only C2H2, N2, CO2 O, F2, and C2H2, all of which are shown in Table II of the original paper. A corrected version of that table can be found below. This changes our results slightly in that the RMSEs of the CIS and the TDDFT results are slightly reduced, but it does not change the primary conclusion that NOCIS is a promising method for calculating core excitations though it lacks dynamical correlation. (Table Presented).