Advances in quantum chemical methods in combination with exponential growth in the computational speed of computers have enabled researchers in the field of catalysis to apply electronic structure calculations to a wide variety of increasingly complex problems. Such calculations provide insights into why and how changes in the composition and structure of catalytically active sites affect their activity and selectivity for targeted reactions. The aim of this review is to survey the recent advances in the methods used to make quantum chemical calculations and to define transition states as well as to illustrate the application of these methods to a selected series of examples taken from the authors' recent work.
There is increasing interest in the possibility of photoelectrochemical (PEC) reduction of CO<sub>2</sub> to C<sub>2+</sub> products; however, the criteria for maximizing PEC solar-to-C<sub>2+</sub> (STC<sub>2+</sub>) rates are not well understood. We report here a continuum-scale model of PEC CO<sub>2</sub> reduction (CO<sub>2</sub>R) on Cu in 0.1 M CsHCO<sub>3</sub> and use it to optimize the design and operating conditions for generating C<sub>2+</sub>products. Furthermore, we demonstrate that the potential-dependent product distribution of CO<sub>2</sub>R on Cu requires operating near the potential that maximizes C<sub>2+</sub> generation rates ($V$<sub>id</sub>), unlike PEC water splitting, which desires operation at the maximum photocurrent density. Because of this requirement, the criterion for a high STC<sub>2+</sub> rate includes high-photocurrent semiconductors with photovoltages near $V$<sub>id</sub> and low series resistance. The STC<sub>2+</sub> rate in these systems is enhanced by optimal CO<sub>2</sub> transport and exhibits low sensitivity to dirunal solar irradiance variations.
The characterization of native point defects in ZnO is still a question of debate. For example, experimental evidence for ZnO with an excess of Zn is inconclusive as to whether the dominant defects are metal interstitials or oxygen vacancies. This information is essential to understand the behavior of the material and to tailor its numerous technological applications. We use the first-principles pseudopotential method to determine the electronic structure, atomic geometry, and formation energy of native point defects in ZnO. Interstitials, vacancies, and antisites in their relevant charge states are considered and the effects of dopants are also discussed. The results show that both the Zn and O vacancies are the relevant defects in ZnO. We also propose a possible transition mechanism and defect center responsible for the experimentally observed green luminescence.
A major challenge in the modeling of electrochemical phenomena is the accurate description of the interface between an electrolyte and a charged conductor. Polarizable continuum models (PCM) have been gaining popularity because they offer a computationally inexpensive method of modeling the electrolyte. In this Perspective, we discuss challenges from using one such model which treats the ions using a linearized Poisson-Boltzmann (LPB) distribution. From a physical perspective, this model places charge unphysically close to the surface and adsorbates, and it includes excessively steep ramping of the dielectric constant from the surface to the bulk solvent. Both of these issues can be somewhat mitigated by adjusting parameters built into the model, but in doing so, the resultant capacitance deviates from experimental values. Likewise, hybrid explicit-implicit approaches to the solvent may offer a more realistic description of hydrogen bonding and solvation to reaction intermediates, but the corresponding capacitances also deviate from experimental values. These deviations highlight the need for a careful adjustment of parameters in order to reproduce not only solvation energies but also other physical properties of solid-liquid interfaces. Continuum approaches alone also necessarily do not capture local variations in the electric field from cations at the interface, which can affect the energetics of intermediates with substantial dipoles or polarizability. Finally, since the double-layer charge can be varied continuously, LPB/PCM models provide a way to determine electrochemical barriers at constant potential. However, double-layer charging and the atomic motion associated with reaction events occur on significantly different timescales. We suggest that more detailed approaches, such as the modified Poisson-Boltzmann model and/or the addition of a Stern layer, may be able to mitigate some but not all of the challenges discussed.
Establishing how Cu facilitates the electrochemical CO<sub>2</sub> reduction reaction (CO2RR) to C<sub>2+</sub> products remains a critical challenge. Under typical reaction conditions, the pH near the electrode is considerably more alkaline than that in the bulk due to mass transport limitations. Challenges with probing alkaline pathways using computational methods have limited understanding of the CO2RR under experimentally relevant conditions. In this work, using the Volmer reaction on Cu (100), we demonstrate that predicted activation barriers can substantially differ between acidic and alkaline pathways. We compute reaction energetics for alkaline *CO protonation and find that, while the formation of *CHO is preferred thermodynamically, the formation of *COH is favored kinetically at high overpotential. However, we find formation of *CHO via reaction of *H and *CO feasible at room temperature. Further, we report potential-dependent energetics for forming the first C-C bond in CO2RR and find that CO dimerization likely dominates. Finally, we investigate how long-range van der Waals interactions impact our results by comparing to the meta-GGA B97M-rV.
Predicting and characterizing the crystal structure of materials is a key\nproblem in materials research and development. We report the results of ab\ninitio LDA/GGA computations for the following systems: AgAu, AgCd, AgMg, AgMo*,\nAgNa, AgNb*, AgPd, AgRh*, AgRu*, AgTc*, AgTi, AgY, AgZr, AlSc, AuCd, AuMo*,\nAuNb, AuPd, AuPt*, AuRh*, AuRu*, AuSc, AuTc*, AuTi, AuY, AuZr, CdMo*, CdNb*,\nCdPd, CdPt, CdRh, CdRu*, CdTc*, CdTi, CdY, CdZr, CrMg*, MoNb, MoPd, MoPt, MoRh,\nMoRu, MoTc*, MoTi, MoY*, MoZr, NbPd, NbPt, NbRh, NbRu, NbTc, NbY*, NbZr*, PdPt,\nPdRh*, PdRu*, PdTc, PdTi, PdY, PdZr, PtRh, PtRu, PtY, PtTc, PtTi, PtZr, RhRu,\nRhTc, RhTi, RhY, RhZr, RuTi, RuTc, RuY, RuZr, TcTi, TcY, TcZr, TiZr*, YZr* (*=\nsystems in which the ab initio method predicts that no compounds are stable). A\ndetailed comparison to experimental data confirms the high accuracy with which\nab initio methods can predict ground states.\n Keywords: Binary Alloys, Ab initio, Intermetallics, Transition Metals,\nStructureAluminum, Cadmium, Gold, Magnesium, Molybdenum, Niobium, Palladium,\nPlatinum, Rhodium, Ruthenium, Scandium, Silver, Sodium, Titanium, Technetium,\nYttrium, Zirconium.\n
This project aimed to develop fundamental understanding of the chemistry of NO adsorption and reaction in Pd/zeolites so as to facilitate the rational design of passive NOx adsorber catalysts. The approach adopted combined both experimental and computational methods, which together allow a deeper understanding of the governing chemistry than the use of either method alone. The workflow began with Pd/H-CHA and Pd/H-BEA catalyst synthesis and characterization, in which the Si/Al ratio and Al siting were systematically varied. This was followed by catalyst evaluation using temperature-programed adsorption/desorption methods, as well as in situ spectroscopic measurements to probe the chemistry of NO adsorption. In parallel, the adsorption of NO and other relevant species (H<sub>2</sub>O, CO, HCs) was studied by means of quantum chemical calculations in order to rationalize the experimental data and provide additional insights. Catalyst aging studies were also performed with the aim of elucidating the mechanism of catalyst degradation. Finally, the insights gained in this project were applied to the preparation of an optimized HC/NOx adsorber catalyst, the performance of which was studied using exhaust gas from an engine dynamometer.
Read moreMembrane-electrode assemblies utilize ionomer-coated electrocatalysts to achieve facile ion transport. Consequently, isolation of intrinsic catalyst kinetics from measured polarization curves is challenging, as the properties of the catalyst and ionomer both affect the measurements. Here, we employ a Pt microelectrode coated by a thin perfluorosulfonic acid (PFSA) layer to measure polarization curves for the hydrogen oxidation reaction/hydrogen evolution reaction (HER/HOR). Intrinsic electrode kinetics are isolated by theoretical analysis of the local catalyst microenvironment, accounting for mass transport and thermodynamics. The observed enhancements in HER and HOR rates with increasing relative humidity (RH) at the working electrode are attributable to two competing factors: the decrease in activity of H+ in the ionomer and the dominant decrease in the water reorganization energy in the Marcus–Hush–Chidsey representation of HER/HOR kinetics. The increase in intrinsic rate with increasing RH is attributed to increased H+ transfer dynamics resulting from reduced confinement of water within the subnanometer water layer between the catalyst and ionomer as RH increases.
Read moreA decade ago, the U.S. chemical industry was in decline. Of the more than 40 chemical manufacturing plants being built worldwide in the mid-2000s with more than $1 billion in capitalization, none were under construction in the United States. Today, as a result of abundant domestic supplies of affordable natural gas and natural gas liquids resulting from the dramatic rise in shale gas production, the U.S. chemical industry has gone from the world’s highest-cost producer in 2005 to among the lowest-cost producers today. The low cost and increased supply of natural gas and natural gas liquids provides an opportunity to discover and develop new catalysts and processes to enable the direct conversion of natural gas and natural gas liquids into value-added chemicals with a lower carbon footprint. The economic implications of developing advanced technologies to utilize and process natural gas and natural gas liquids for chemical production could be significant, as commodity, intermediate, and fine chemicals represent a higher-economic-value use of shale gas compared with its use as a fuel. To better understand the opportunities for catalysis research in an era of shifting feedstocks for chemical production and to identify the gaps in the current research portfolio, the National Academies of Sciences, Engineering, and Medicine conducted an interactive, multidisciplinary workshop in March 2016. The goal of this workshop was to identify advances in catalysis that can enable the United States to fully realize the potential of the shale gas revolution for the U.S. chemical industry and, as a result, to help target the efforts of U.S. researchers and funding agencies on those areas of science and technology development that are most critical to achieving these advances. This publication summarizes the presentations and discussions from the workshop.
Read moreWe describe catalytic sequences for converting biomass-derived carboxylic acids, to fuels and lubricants that are compatible with the existing energy infrastructure.
Read moreAbstract Synthesis of a pentasil‐type zeolite with ultra‐small few‐unit‐cell crystalline domains, which we call FDP (few‐unit‐cell crystalline domain pentasil), is reported. FDP is made using bis‐1,5(tributyl ammonium) pentamethylene cations as structure directing agent (SDA). This di‐quaternary ammonium SDA combines butyl ammonium, in place of the one commonly used for MFI synthesis, propyl ammonium, and a five‐carbon nitrogen‐connecting chain, in place of the six‐carbon connecting chain SDAs that are known to fit well within the MFI pores. X‐ray diffraction analysis and electron microscopy imaging of FDP indicate ca. 10 nm crystalline domains organized in hierarchical micro‐/meso‐porous aggregates exhibiting mesoscopic order with an aggregate particle size up to ca. 5 μm. Al and Sn can be incorporated into the FDP zeolite framework to produce active and selective methanol‐to‐hydrocarbon and glucose isomerization catalysts, respectively.
Read moresolid solutions of structure type have been synthesized from homogeneous hydroxide precursors. X‐ray powder diffraction and scanning transmission electron microscopy show that single‐phase solid solutions are formed up to y ≈ 0.5 at 800°C in air. Electrochemical tests using nonaqueous liquid electrolyte cells shows that the open‐circuit voltage and working voltage increase with Al content as predicted by ab initio calculations [Ceder et al., Nature, 392, 694 (1998)], while capacity fade was significant during cycling at room temperature. However, both the absolute capacity and cycleability were improved at 55°C. © 1999 The Electrochemical Society. All rights reserved.
Read moreAbstract Sn‐BEA zeolite is known to catalyze the aldose‐to‐ketose isomerization of xylose and glucose; however, the selectivity to pentose and hexose isomers is not stoichiometric, suggesting the formation of other products. In the present study, we have observed near‐complete conversion of all pentose and hexose isomers when xylose and glucose were reacted in the presence of Sn‐BEA at 140 °C and 200 °C, respectively. The previously unidentified products were identified by nuclear magnetic resonance and mass spectrometry to be hydroxyalkanoic acids and their derivatives. The hydroxyl‐rich acids comprise a significant fraction of the converted sugars and are potential monomers for the synthesis of hyper‐crosslinked, biodegradable polymers.
Read moreThe next generation of alternative fuels is being investigated through advanced chemical and biological production techniques for the purpose of finding suitable replacements to diesel and gasoline while lowering production costs and increasing process yields. Chemical conversion of biomass to fuels provides a plethora of pathways with a variety of fuel molecules, both novel and traditional, which may be targeted. In the search for new fuels, an initial, intuition-driven evaluation of fuel compounds with desired properties is required. Due to the high cost and significant production time needed to synthesize these materials at a scale sufficient for exhaustive testing, a predictive model would allow chemists to preemptively screen fuel properties of potentially desirable fuel candidates. Recent work has shown that predictive models, in this case artificial neural networks (ANN’s) analyzing quantitative structure property relationships (QSPR’s), can predict the cetane number (CN) of a proposed fuel molecule with relatively small error. A fuel’s CN is a measure of its ignition quality, typically defined using prescribed ASTM standards and a cetane testing engine. Alternatively, the analogous derived cetane number (DCN), obtained using an Ignition Quality Tester (IQT), is a direct measurement alternative to the CN that uses an empirical inverse relationship to the ignition delay found in the constant volume combustion chamber apparatus. DCN data points acquired using an IQT were utilized for model validation and expansion of the experimental database used in this study. The present work improves on an existing model by optimizing the model architecture along with the key learning variables of the ANN and by making the model more generalizable to a wider variety of fuel candidate types, specifically the class of furans and furan derivatives, by including specific molecules for the model to incorporate. The new molecules considered include tetrahydrofuran, 2-methylfuran, 2-methyltetrahydrofuran, 5,5'-(furan-2-ylmethylene)bis(2-methylfuran), 5,5'-((tetrahydrofuran-2-yl)methylene)bis(2-methyltetrahydrofuran), tris(5-methylfuran-2-yl)methane, and tris(5-methyltetrahydrofuran-2-yl)methane. Model architecture adjustments improved the overall root-mean-squared error (RMSE) of the base database predictions by 5.54%. Additionally, through the targeted database expansion, it is shown that the predicted cetane number of the furan-based molecules improves on average by 49.21% (3.74 CN units) and significantly for a few of the individual molecules. This indicates that a selected subset of representative molecules can be used to extend the model’s predictive accuracy to new molecular classes. The approach, bolstered by the improvements presented in this paper, enables chemists to focus on promising molecules by eliminating less favorable candidates in relation to their ignition quality.
Read moreThe ability to efficiently locate transition states is critically important to the widespread adoption of theoretical chemistry techniques for their ability to accurately predict kinetic constants. Existing surface walking techniques to locate such transition states typically require an extremely good initial guess that is often beyond human intuition to estimate. To alleviate this problem, automated techniques to locate transition state guesses have been created that take the known reactant and product endpoint structures as inputs. In this work, we present a simple method to build an approximate reaction path through a combination of interpolation and optimization. Starting from the known reactant and product structures, new nodes are interpolated inwards towards the transition state, partially optimized orthogonally to the reaction path, and then frozen before a new pair of nodes is added. The algorithm is stopped once the string ends connect. For the practical user, this method provides a quick and convenient way to generate transition state structure guesses. Tests on three reactions (cyclization of cis,cis-2,4-hexadiene, alanine dipeptide conformation transition, and ethylene dimerization in a Ni-exchanged zeolite) show that this “freezing string” method is an efficient way to identify complex transition states with significant cost savings over existing methods, particularly when high quality linear synchronous transit interpolation is employed.
Read more