In silico catalyst design is a grand challenge of chemistry. Traditional computational approaches have been limited by the need to compute properties for an intractably large number of possible ...catalysts. Recently, inverse design methods have emerged, starting from a desired property and optimizing a corresponding chemical structure. Techniques used for exploring chemical space include gradient-based optimization, alchemical transformations, and machine learning. Though the application of these methods to catalysis is in its early stages, further development will allow for robust computational catalyst design. This review provides an overview of the evolution of inverse design approaches and their relevance to catalysis. The strengths and limitations of existing techniques are highlighted, and suggestions for future research are provided.
Efficient photoelectrochemical water oxidation may open a way to produce energy from renewable solar power. In biology, generation of fuel due to water oxidation happens efficiently on an immense ...scale during the light reactions of photosynthesis. To oxidize water, photosynthetic organisms have evolved a highly conserved protein complex, Photosystem II. Within that complex, water oxidation happens at the CaMn4O5 inorganic catalytic cluster, the so-called oxygen-evolving complex (OEC), which cycles through storage “S” states as it accumulates oxidizing equivalents and produces molecular oxygen. In recent years, there has been significant progress in understanding the OEC as it evolves through the catalytic cycle. Studies have combined conventional and femtosecond X-ray crystallography with extended X-ray absorption fine structure (EXAFS) and quantum mechanics/molecular mechanics (QM/MM) methods and have addressed changes in protonation states of μ-oxo bridges and the coordination of substrate water through the analysis of ammonia binding as a chemical analog of water. These advances are thought to be critical to understanding the catalytic cycle since protonation states regulate the relative stability of different redox states and the geometry of the OEC. Therefore, establishing the mechanism for substrate water binding and the nature of protonation/redox state transitions in the OEC is essential for understanding the catalytic cycle of O2 evolution. The structure of the dark-stable S1 state has been a target for X-ray crystallography for the past 15 years. However, traditional X-ray crystallography has been hampered by radiation-induced reduction of the OEC. Very recently, a revolutionary X-ray free electron laser (XFEL) technique was applied to PSII to reveal atomic positions at 1.95 Å without radiation damage, which brought us closer than ever to establishing the ultimate structure of the OEC in the S1 state. However, the atom positions in this crystal structure are still not consistent with high-resolution EXAFS spectroscopy, partially due to the poorly resolved oxygen positions next to Mn centers and partial reduction due to extended dark adaptation of the sample. These inconsistencies led to the new models of the OEC with an alternative low oxidation state and raised questions on the protonation state of the cluster, especially the O5 μ-oxo bridge. This Account summarizes the most recent models of the OEC that emerged from QM/MM, EXAFS and femtosecond X-ray crystallography methods. When PSII in the S1 state is exposed to light, the S1 state is advanced to the higher oxidation states and eventually binds substrate water molecules. Identifying the substrate waters is of paramount importance for establishing the water-oxidation mechanism but is complicated by a large number of spectroscopically similar waters. Water analogues can, therefore, be helpful because they serve as spectroscopic markers that help to track the motion of the substrate waters. Due to a close structural and electronic similarity to water, ammonia has been of particular interest. We review three competing hypotheses on substrate water/ammonia binding and compile theoretical and experimental evidence to support them. Binding of ammonia as a sixth ligand to Mn4 during the S1 → S2 transition seems to satisfy most of the criteria, especially the most compelling recent EPR data on D1-D61A mutated PSII. Such a binding mode suggests delivery of water from the “narrow” channel through a “carousel” rearrangement of waters around Mn4 upon the S2 → S3 transition. An alternative hypothesis suggests water delivery through the “large” channel on the Ca side. However, both water delivery paths lead to a similar S3 structure, seemingly reaching consensus on the nature of the last detectable S-state intermediate in the Kok cycle before O2 evolution.
Tungsten carbide is one of the most promising electrocatalysts for the hydrogen evolution reaction, although it exhibits sluggish kinetics due to a strong tungsten-hydrogen bond. In addition, ...tungsten carbide's catalytic activity toward the oxygen evolution reaction has yet to be reported. Here, we introduce a superaerophobic nitrogen-doped tungsten carbide nanoarray electrode exhibiting high stability and activity toward hydrogen evolution reaction as well as driving oxygen evolution efficiently in acid. Nitrogen-doping and nanoarray structure accelerate hydrogen gas release from the electrode, realizing a current density of -200 mA cm
at the potential of -190 mV vs. reversible hydrogen electrode, which manifest one of the best non-noble metal catalysts for hydrogen evolution reaction. Under acidic conditions (0.5 M sulfuric acid), water splitting catalyzed by nitrogen-doped tungsten carbide nanoarray starts from about 1.4 V, and outperforms most other water splitting catalysts.
Water oxidation is the key kinetic bottleneck of photoelectrochemical devices for fuel synthesis. Despite advances in the identification of intermediates, elucidating the catalytic mechanism of this ...multi-redox reaction on metal-oxide photoanodes remains a significant experimental and theoretical challenge. Here, we report an experimental analysis of water oxidation kinetics on four widely studied metal oxides, focusing particularly on haematite. We observe that haematite is able to access a reaction mechanism that is third order in surface-hole density, which is assigned to equilibration between three surface holes and M(OH)-O-M(OH) sites. This reaction exhibits low activation energy (E
≈ 60 meV). Density functional theory is used to determine the energetics of charge accumulation and O-O bond formation on a model haematite (110) surface. The proposed mechanism shows parallels with the function of the oxygen evolving complex of photosystem II, and provides new insights into the mechanism of heterogeneous water oxidation on a metal oxide surface.
Efficient photoelectrochemical water oxidation may open a way to produce energy from renewable solar power. In biology, generation of fuel due to water oxidation happens efficiently on an immense ...scale during the light reactions of photosynthesis. To oxidize water, photosynthetic organisms have evolved a highly conserved protein complex, Photosystem II. Within that complex, water oxidation happens at the CaMn
O
inorganic catalytic cluster, the so-called oxygen-evolving complex (OEC), which cycles through storage "S" states as it accumulates oxidizing equivalents and produces molecular oxygen. In recent years, there has been significant progress in understanding the OEC as it evolves through the catalytic cycle. Studies have combined conventional and femtosecond X-ray crystallography with extended X-ray absorption fine structure (EXAFS) and quantum mechanics/molecular mechanics (QM/MM) methods and have addressed changes in protonation states of μ-oxo bridges and the coordination of substrate water through the analysis of ammonia binding as a chemical analog of water. These advances are thought to be critical to understanding the catalytic cycle since protonation states regulate the relative stability of different redox states and the geometry of the OEC. Therefore, establishing the mechanism for substrate water binding and the nature of protonation/redox state transitions in the OEC is essential for understanding the catalytic cycle of O
evolution. The structure of the dark-stable S
state has been a target for X-ray crystallography for the past 15 years. However, traditional X-ray crystallography has been hampered by radiation-induced reduction of the OEC. Very recently, a revolutionary X-ray free electron laser (XFEL) technique was applied to PSII to reveal atomic positions at 1.95 Å without radiation damage, which brought us closer than ever to establishing the ultimate structure of the OEC in the S
state. However, the atom positions in this crystal structure are still not consistent with high-resolution EXAFS spectroscopy, partially due to the poorly resolved oxygen positions next to Mn centers and partial reduction due to extended dark adaptation of the sample. These inconsistencies led to the new models of the OEC with an alternative low oxidation state and raised questions on the protonation state of the cluster, especially the O5 μ-oxo bridge. This Account summarizes the most recent models of the OEC that emerged from QM/MM, EXAFS and femtosecond X-ray crystallography methods. When PSII in the S
state is exposed to light, the S
state is advanced to the higher oxidation states and eventually binds substrate water molecules. Identifying the substrate waters is of paramount importance for establishing the water-oxidation mechanism but is complicated by a large number of spectroscopically similar waters. Water analogues can, therefore, be helpful because they serve as spectroscopic markers that help to track the motion of the substrate waters. Due to a close structural and electronic similarity to water, ammonia has been of particular interest. We review three competing hypotheses on substrate water/ammonia binding and compile theoretical and experimental evidence to support them. Binding of ammonia as a sixth ligand to Mn4 during the S
→ S
transition seems to satisfy most of the criteria, especially the most compelling recent EPR data on D1-D61A mutated PSII. Such a binding mode suggests delivery of water from the "narrow" channel through a "carousel" rearrangement of waters around Mn4 upon the S
→ S
transition. An alternative hypothesis suggests water delivery through the "large" channel on the Ca side. However, both water delivery paths lead to a similar S
structure, seemingly reaching consensus on the nature of the last detectable S-state intermediate in the Kok cycle before O
evolution.
Restructuring-induced catalytic activity is an intriguing phenomenon of fundamental importance to rational design of high-performance catalyst materials. We study three copper-complex materials for ...electrocatalytic carbon dioxide reduction. Among them, the copper(II) phthalocyanine exhibits by far the highest activity for yielding methane with a Faradaic efficiency of 66% and a partial current density of 13 mA cm
at the potential of - 1.06 V versus the reversible hydrogen electrode. Utilizing in-situ and operando X-ray absorption spectroscopy, we find that under the working conditions copper(II) phthalocyanine undergoes reversible structural and oxidation state changes to form ~ 2 nm metallic copper clusters, which catalyzes the carbon dioxide-to-methane conversion. Density functional calculations rationalize the restructuring behavior and attribute the reversibility to the strong divalent metal ion-ligand coordination in the copper(II) phthalocyanine molecular structure and the small size of the generated copper clusters under the reaction conditions.
Exploration of heterogeneous molecular catalysts combining the atomic-level tunability of molecular structures and the practical handling advantages of heterogeneous catalysts represents an ...attractive approach to developing high-performance catalysts for important and challenging chemical reactions such as electrochemical carbon dioxide reduction which holds the promise for converting emissions back to fuels utilizing renewable energy. Thus, far, efficient and selective electroreduction of CO2 to deeply reduced products such as hydrocarbons remains a big challenge. Here, we report a molecular copper-porphyrin complex (copper(II)-5,10,15,20-tetrakis(2,6-dihydroxyphenyl)porphyrin) that can be used as a heterogeneous electrocatalyst with high activity and selectivity for reducing CO2 to hydrocarbons in aqueous media. At −0.976 V vs the reversible hydrogen electrode, the catalyst is able to drive partial current densities of 13.2 and 8.4 mA cm–2 for methane and ethylene production from CO2 reduction, corresponding to turnover frequencies of 4.3 and 1.8 molecules·site–1·s–1 for methane and ethylene, respectively. This represents the highest catalytic activity to date for hydrocarbon production over a molecular CO2 reduction electrocatalyst. The unprecedented catalytic performance is attributed to the built-in hydroxyl groups in the porphyrin structure and the reactivity of the copper(I) metal center.
We introduce the “tensor-train split-operator Fourier transform” (TT-SOFT) method for simulations of multidimensional nonadiabatic quantum dynamics. TT-SOFT is essentially the grid-based SOFT method ...implemented in dynamically adaptive tensor-train representations. In the same spirit of all matrix product states, the tensor-train format enables the representation, propagation, and computation of observables of multidimensional wave functions in terms of the grid-based wavepacket tensor components, bypassing the need of actually computing the wave function in its full-rank tensor product grid space. We demonstrate the accuracy and efficiency of the TT-SOFT method as applied to propagation of 24-dimensional wave packets, describing the S 1/S 2 interconversion dynamics of pyrazine after UV photoexcitation to the S 2 state. Our results show that the TT-SOFT method is a powerful computational approach for simulations of quantum dynamics of polyatomic systems since it avoids the exponential scaling problem of full-rank grid-based representations.
We introduce a new family of complexes with the general formula Ru n (tda)(py)2 m+ (n = 2, m = 0, 1; n = 3, m = 1, 2 +; n = 4, m = 2, 3 2+), with tda2– being ...2,2′:6′,2″-terpyridine-6,6″-dicarboxylate, including complex RuIV(OH)(tda-κ-N3O)(py)2+, 4H +, which we find to be an impressive water oxidation catalyst, formed by hydroxo coordination to 3 2+ under basic conditions. The complexes are synthesized, isolated, and thoroughly characterized by analytical, spectroscopic (UV–vis, nuclear magnetic resonance, electron paramagnetic resonance), computational, and electrochemical techniques (cyclic voltammetry, differential pulse voltammetry, coulometry), including solid-state monocrystal X-ray diffraction analysis. In oxidation state IV, the Ru center is seven-coordinated and diamagnetic, whereas in oxidation state II, the complex has an unbonded dangling carboxylate and is six-coordinated while still diamagnetic. With oxidation state III, the coordination number is halfway between the coordination of oxidation states II and IV. Species generated in situ have also been characterized by spectroscopic, computational, and electrochemical techniques, together with the related species derived from a different degree of protonation and oxidation states. 4H + can be generated potentiometrically, or voltammetrically, from 3 2+, and both coexist in solution. While complex 3 2+ is not catalytically active, the catalytic performance of complex 4H + is characterized by the foot of the wave analysis, giving an impressive turnover frequency record of 8000 s–1 at pH 7.0 and 50 000 s–1 at pH 10.0. Density functional theory calculations provide a complete description of the water oxidation catalytic cycle of 4H +, manifesting the key functional role of the dangling carboxylate in lowering the activation free energies that lead to O–O bond formation.
Copper oxides have been of considerable interest as electrocatalysts for CO2 reduction (CO2R) in aqueous electrolytes. However, their role as an active catalyst in reducing the required overpotential ...and improving the selectivity of reaction compared with that of polycrystalline copper remains controversial. Here, we introduce the use of selected-ion flow tube mass spectrometry, in concert with chronopotentiometry, in situ Raman spectroscopy, and computational modeling, to investigate CO2R on Cu2O nanoneedles, Cu2O nanocrystals, and Cu2O nanoparticles. We show experimentally that the selective formation of gaseous C2 products (i.e., ethylene) in CO2R is preceded by the reduction of the copper oxide (Cu2OR) surface to metallic copper. On the basis of density functional theory modeling, CO2R products are not formed as long as Cu2O is present at the surface because Cu2OR is kinetically and energetically more favorable than CO2R.