Computational Design of an Electro-Organocatalyst for Conversion of CO2 into Long Chain Aldehydes.

Density functional theory calculations employing a hybrid implicit/explicit solvation method were used to demonstrate that an electro-organocatalyst designed in our previous work for reducing CO2 to formaldehyde could also be capable of coupling formaldehyde to form long chain aldehydes. The catalytic activity is enabled by an electron-rich vicinal enediamine (>N-C═C-N<) backbone that activates formaldehyde by reversing the polarity on the carbon atom, enabling it to act as a nucleophile in the subsequent aldol addition step. The catalyst then enables reductive dehydroxylation of the aldol addition product by facilitating outer-sphere electron transfer. The optimal pH as well as the limiting potential and formaldehyde concentration are identified and related to the kinetic balance between several rate limiting steps. Finally, the optimal conditions for coupling with the CO2 reduction cycle are discussed, demonstrating that the electro-organocatalyst is capable of efficiently converting CO2 into aldehyde products with a turnover frequency (per carbon atom) on the order of 0.1-1 s-1.

Computational Design of an Electro-Organocatalyst for Conversion of CO2 into Formaldehyde.

Density functional theory calculations employing a hybrid implicit/explicit solvation method were used to explore a new strategy for electrochemical conversion of CO2 using an electro-organocatalyst. A particular structural motif is identified that consists of an electron-rich vicinal enediamine (>N-C═C-N<) backbone, which is capable of activating CO2 by the formation of a C-C bond while subsequently facilitating the transfer of electrons from a chemically inert cathode to ultimately produce formaldehyde. Unlike transition metal-based electrocatalysts, the electro-organocatalyst is not constrained by scaling relations between the formation energies of activated CO2 and adsorbed CO, nor is it expected to be active for the competing hydrogen evolution reaction. The rate-limiting steps are found to occur during two proton-coupled electron transfer (PCET) sequences and are associated with the transfer of a proton from a proton transfer mediator to a carbon atom on the electro-organocatalyst. The difficulty of this step in the second PCET sequence necessitates an electrode potential of -0.85 V vs RHE to achieve the maximum turnover frequency. In addition, it is postulated that the electro-organocatalyst should also be capable of forming long-chain aldehydes by successively carrying out reductive aldol condensation to grow the alkyl chain one carbon at a time.

An implicit electrolyte model for plane wave density functional theory exhibiting nonlinear response and a nonlocal cavity definition.

We have developed and implemented an implicit electrolyte model in the Vienna Ab initio Simulation Package (VASP) that includes nonlinear dielectric and ionic responses as well as a nonlocal definition of the cavities defining the spatial regions where these responses can occur. The implementation into the existing VASPsol code is numerically efficient and exhibits robust convergence, requiring computational effort only slightly higher than the original linear polarizable continuum model. The nonlinear + nonlocal model is able to reproduce the characteristic "double hump" shape observed experimentally for the differential capacitance of an electrified metal interface while preventing "leakage" of the electrolyte into regions of space too small to contain a single water molecule or solvated ion. The model also gives a reasonable prediction of molecular solvation free energies as well as the self-ionization free energy of water and the absolute electron chemical potential of the standard hydrogen electrode. All of this, combined with the additional ability to run constant potential density functional theory calculations, should enable the routine computation of activation barriers for electrocatalytic processes.

Modifying Metastable Sr1-xBO3-δ (B = Nb, Ta, and Mo) Perovskites for Electrode Materials.

The presence of surface/deep defects in 4d- and 5d-perovskite oxide (ABO3, B = Nb, Ta, Mo, etc.) nanoparticles (NPs), originating from multivalent B-site cations, contributes to suppressing their metallic properties. These defect states can be removed using a H2/Ar thermal treatment, enabling the recovery of their electronic properties (i.e., low electrical resistivity, high carrier concentration, etc.) as expected from their electronic structure. Therefore, to engineer the electronic properties of these metastable perovskites, an oxygen-controlled crystallization approach coupled with a subsequent H2/Ar treatment was utilized. A comprehensive study of the effect of the post-treatment time on the electronic properties of these perovskite NPs was performed using a combination of scattering, spectroscopic, and computational techniques. These measurements revealed that a metallic-like state is stabilized in these oxygen-reduced NPs due to the suppression of deep rather than surface defects. Ultimately, this synthetic approach can be employed to synthesize ABO3 perovskite NPs with tunable electronic properties for application into electrochemical devices.

Adsorption of Polarized Molecules for Interfacial Band Engineering of Doped TiO2 Thin Films.

Owing to their chemical and mechanical stability, metal-oxides have emerged as potential alternatives for conventional pure-metal and organic molecule-based solid-state electronic devices. Traditionally, band engineering of these metal-oxides has been performed to improve the efficiency of solar cells and transistors. However, recent advancements in the field of oxide-based electronic devices demand reversible band structure engineering for applications in next-generation adaptive electronics and memory devices. Therefore, this work aims to reversibly engineer the surface band structure of doped metal-oxides using stable organic ligands with weak dipoles. Para-substituted benzoic acid (BZA) ligands with positive and negative dipole moments were adsorbed in situ on the surface of TiO2:Ni2+ thin film to modify the interfacial dipole moment, and the valence band structure was probed using surface-sensitive ultraviolet photoelectron spectroscopy (UPS). UPS, paired with density functional theory (DFT) simulations, demonstrate the ability to selectively tune interfacial electronic/chemical landscapes with ligand-dependent dipole moment. The unique ability to reversibly tune the band bending at the organic-inorganic interface of doped metal-oxide semiconductors using molecular dipoles is expected to play a key role in the development of metal-oxide-based adaptive electronics that outperform the conventional polymer-based and Si-based devices.

Stabilizing the B-site oxidation state in ABO3 perovskite nanoparticles.

The stabilization of the B-site oxidation state in ABO3 perovskites using wet-chemical methods is a synthetic challenge, which is of fundamental and practical interest for energy storage and conversion devices. In this work, defect-controlled (Sr-deficiency and oxygen vacancies) strontium niobium(iv) oxide (Sr1-xNbO3-δ, SNO) metal oxide nanoparticles (NPs) were synthesized for the first time using a low-pressure wet-chemistry synthesis. The experiments were performed under reduced oxygen partial pressure to prevent by-product formation and with varying Sr/Nb molar ratio to favor the formation of Nb4+ pervoskites. At a critical Sr to Nb ratio (Sr/Nb = 1.3), a phase transition is observed forming an oxygen-deficient SrNbO3 phase. Structural refinement on the resultant diffraction pattern shows that the SNO NPs consists of a near equal mixture of SrNbO3 and Sr0.7NbO3-δ crystal phases. A combination of Rietveld refinement and X-ray photoelectron spectroscopy (XPS) confirmed the stabilization of the +4 oxidation state and the formation of oxygen vacancies. The Nb local site symmetry was extracted through Raman spectroscopy and modeled using DFT. As further confirmation, the particles demonstrate the expected absorption highlighting their restored optoelectronic properties. This low-pressure wet-chemical approach for stabilizing the oxidation state of a transition metal has the potential to be extended to other oxygen sensitive, low dimensional perovskite oxides with unique properties.

Lewis-Brønsted Acid Pairs in Ga/H-ZSM-5 To Catalyze Dehydrogenation of Light Alkanes.

The active sites for propane dehydrogenation in Ga/H-ZSM-5 with moderate concentrations of tetrahedral aluminum in the lattice were identified to be Lewis-Brønsted acid pairs. With increasing availability, Ga+ and Brønsted acid site concentrations changed inversely, as protons of Brønsted acid sites were exchanged with Ga+. At a Ga/Al ratio of 1/2, the rate of propane dehydrogenation was 2 orders of magnitude higher than with the parent H-ZSM-5, highlighting the extraordinary activity of the Lewis-Brønsted acid pairs. Density functional theory calculations relate the high activity to a bifunctional mechanism that proceeds via heterolytic activation of the propane C-H bond followed by a monomolecular elimination of H2 and desorption of propene.

Addressing global uncertainty and sensitivity in first-principles based microkinetic models by an adaptive sparse grid approach.

In the last decade, first-principles-based microkinetic modeling has been developed into an important tool for a mechanistic understanding of heterogeneous catalysis. A commonly known, but hitherto barely analyzed issue in this kind of modeling is the presence of sizable errors from the use of approximate Density Functional Theory (DFT). We here address the propagation of these errors to the catalytic turnover frequency (TOF) by global sensitivity and uncertainty analysis. Both analyses require the numerical quadrature of high-dimensional integrals. To achieve this efficiently, we utilize and extend an adaptive sparse grid approach and exploit the confinement of the strongly non-linear behavior of the TOF to local regions of the parameter space. We demonstrate the methodology on a model of the oxygen evolution reaction at the Co3O4 (110)-A surface, using a maximum entropy error model that imposes nothing but reasonable bounds on the errors. For this setting, the DFT errors lead to an absolute uncertainty of several orders of magnitude in the TOF. We nevertheless find that it is still possible to draw conclusions from such uncertain models about the atomistic aspects controlling the reactivity. A comparison with derivative-based local sensitivity analysis instead reveals that this more established approach provides incomplete information. Since the adaptive sparse grids allow for the evaluation of the integrals with only a modest number of function evaluations, this approach opens the way for a global sensitivity analysis of more complex models, for instance, models based on kinetic Monte Carlo simulations.

Assessment of mean-field microkinetic models for CO methanation on stepped metal surfaces using accelerated kinetic Monte Carlo.

First-principles screening studies aimed at predicting the catalytic activity of transition metal (TM) catalysts have traditionally been based on mean-field (MF) microkinetic models, which neglect the effect of spatial correlations in the adsorbate layer. Here we critically assess the accuracy of such models for the specific case of CO methanation over stepped metals by comparing to spatially resolved kinetic Monte Carlo (kMC) simulations. We find that the typical low diffusion barriers offered by metal surfaces can be significantly increased at step sites, which results in persisting correlations in the adsorbate layer. As a consequence, MF models may overestimate the catalytic activity of TM catalysts by several orders of magnitude. The potential higher accuracy of kMC models comes at a higher computational cost, which can be especially challenging for surface reactions on metals due to a large disparity in the time scales of different processes. In order to overcome this issue, we implement and test a recently developed algorithm for achieving temporal acceleration of kMC simulations. While the algorithm overall performs quite well, we identify some challenging cases which may lead to a breakdown of acceleration algorithms and discuss possible directions for future algorithm development.

Constrained-Orbital Density Functional Theory. Computational Method and Applications to Surface Chemical Processes.

We present a method for performing density-functional theory (DFT) calculations in which one or more Kohn-Sham orbitals are constrained to be localized on individual atoms. This constrained-orbital DFT (CO-DFT) approach can be used to tackle two prevalent shortcomings of DFT: the lack of transparency with regard to the governing electronic structure in large (planewave based) DFT calculations and the limitations of semilocal DFT in describing systems with localized electrons or a large degree of static correlation. CO-DFT helps to address the first of these issues by decomposing complex orbital transformations occurring during elementary chemical processes into simpler and more intuitive transformations. The second issue is addressed by using the CO-DFT method to generate configuration states for multiconfiguration Kohn-Sham calculations. We demonstrate both of these applications for elementary reaction steps involved in the oxygen evolution reaction.

Generalized Temporal Acceleration Scheme for Kinetic Monte Carlo Simulations of Surface Catalytic Processes by Scaling the Rates of Fast Reactions.

A novel algorithm is presented that achieves temporal acceleration during kinetic Monte Carlo (KMC) simulations of surface catalytic processes. This algorithm allows for the direct simulation of reaction networks containing kinetic processes occurring on vastly disparate time scales which computationally overburden standard KMC methods. Previously developed methods for temporal acceleration in KMC were designed for specific systems and often require a priori information from the user such as identifying the fast and slow processes. In the approach presented herein, quasi-equilibrated processes are identified automatically based on previous executions of the forward and reverse reactions. Temporal acceleration is achieved by automatically scaling the intrinsic rate constants of the quasi-equilibrated processes, bringing their rates closer to the time scales of the slow kinetically relevant nonequilibrated processes. All reactions are still simulated directly, although with modified rate constants. Abrupt changes in the underlying dynamics of the reaction network are identified during the simulation, and the reaction rate constants are rescaled accordingly. The algorithm was utilized here to model the Fischer-Tropsch synthesis reaction over ruthenium nanoparticles. This reaction network has multiple time-scale-disparate processes which would be intractable to simulate without the aid of temporal acceleration. The accelerated simulations are found to give reaction rates and selectivities indistinguishable from those calculated by an equivalent mean-field kinetic model. The computational savings of the algorithm can span many orders of magnitude in realistic systems, and the computational cost is not limited by the magnitude of the time scale disparity in the system processes. Furthermore, the algorithm has been designed in a generic fashion and can easily be applied to other surface catalytic processes of interest.

Perspective: On the active site model in computational catalyst screening.

First-principles screening approaches exploiting energy trends in surface adsorption represent an unparalleled success story in recent computational catalysis research. Here we argue that our still limited understanding of the structure of active sites is one of the major bottlenecks towards an ever extended and reliable use of such computational screening for catalyst discovery. For low-index transition metal surfaces, the prevalently chosen high-symmetry (terrace and step) sites offered by the nominal bulk-truncated crystal lattice might be justified. For more complex surfaces and composite catalyst materials, computational screening studies will need to actively embrace a considerable uncertainty with respect to what truly are the active sites. By systematically exploring the space of possible active site motifs, such studies might eventually contribute towards a targeted design of optimized sites in future catalysts.

Quantum chemistry of the oxygen evolution reaction on cobalt(ii,iii) oxide - implications for designing the optimal catalyst.

Density functional theory is used to examine the changes in electronic structure that occur during the oxygen evolution reaction (OER) catalyzed by active sites on three different surface terminations of Co3O4. These three active sites have reactive oxo species with differing degrees of coordination by Co cations - a µ(3)-oxo on the (311) surface, a µ(2)-oxo on the (110)-A surface, and an η-oxo on the (110)-B surface. The kinetically relevant step on all surfaces over a wide range of applied potentials is the nucleophilic addition of water to the oxo, which is responsible for formation of the O-O bond. The intrinsic reactivity of a site for this step is found to increase as the coordination of the oxo decreases with the µ(3)-oxo on the (311) surface being the least reactive and the η-oxo on the (110)-B surface being the most reactive. A detailed analysis of the electronic changes occurring during water addition on the three sites reveals that this trend is due to both a decrease in the attractive local Madelung potential on the oxo and a decrease in electron withdrawal from the oxo by Co neighbors. Applying a similar electronic structure analysis to the oxidation steps preceding water addition in the catalytic cycle shows that analogous electronic changes occur during this process, explaining a correlation observed between the oxidation potential of a site and its intrinsic reactivity for water addition. This concept is then used to specify criteria for the design of an optimal OER catalyst at a given applied potential.

Structure Sensitivity of the Oxygen Evolution Reaction Catalyzed by Cobalt(II,III) Oxide.

Quantum chemical calculations and simulated kinetics were used to examine the structure sensitivity of the oxygen evolution reaction on several surface terminations of Co3O4. Active sites consisting of two adjacent Co(IV) cations connected by bridging oxos were identified on both the (001) and (311) surfaces. Formation of the O-O bond proceeds on these sites by nucleophilic attack of water on a bridging oxo. It was found that the relative turnover frequencies for the different sites are highly dependent on the overpotential, with the dual-Co site on the (311) surface being most active at medium overpotentials (0.46-0.77 V), where O-O bond formation by water addition is rate limiting. A similar dual-Co site on the (001) surface is most active at low overpotentials (<0.46 V), where O2 release is rate limiting, and a single-Co site on the (110) surface is most active at overpotentials that are high enough (>0.77 V) to form a significant concentration of highly reactive terminal Co(V)═O species. Two overpotential-dependent Sabatier relationships were identified based on the Brønsted basicity and redox potential of the active site, explaining the change in the active site with overpotential. The (311) dual-Co site that is most active in the medium overpotential range is consistent with recent experimental observations suggesting that a defect site is responsible for the observed oxygen evolution activity and that a modest concentration of superoxo intermediates is present on the surface. Importantly, we find that it is essential to consider the kinetics of the water addition and O2 release steps rather than only the thermodynamics.