Diffusion Barriers for Carbon Monoxide on the Cu(001) Surface Using Many-Body Perturbation Theory and Various Density Functionals

Facet-Defined Dilute Metal Alloy Nanorods for Efficient Electroreduction of CO2 to n-Propanol

Electroreduction of CO2 into liquid fuels is a compelling strategy for storing intermittent renewable energy. Here, we introduce a family of facet-defined dilute copper alloy nanocrystals as catalysts to improve the electrosynthesis of n-propanol from CO2 and H2O. We show that substituting a dilute amount of weak-CO-binding metals into the Cu(100) surface improves CO2-to-n-propanol activity and selectivity by modifying the electronic structure of catalysts to facilitate C1-C2 coupling while preserving the (100)-like 4-fold Cu ensembles which favor C1-C1 coupling. With the Au0.02Cu0.98 champion catalyst, we achieve an n-propanol Faradaic efficiency of 18.2 ± 0.3% at a low potential of -0.41 V versus the reversible hydrogen electrode and a peak production rate of 16.6 mA·cm-2. This study demonstrates that shape-controlled dilute-metal-alloy nanocrystals represent a new frontier in electrocatalyst design, and precise control of the host and minority metal distributions is crucial for elucidating structure-composition-property relationships and attaining superior catalytic performance.

Constructing Mixed Density Functionals for Describing Dissociative Chemisorption on Metal Surfaces: Basic Principles

The production of a majority of chemicals involves heterogeneous catalysis at some stage, and the rates of many heterogeneously catalyzed processes are governed by transition states for dissociative chemisorption on metals. Accurate values of barrier heights for dissociative chemisorption on metals are therefore important to benchmarking electronic structure theory in general and density functionals in particular. Such accurate barriers can be obtained using the semiempirical specific reaction parameter (SRP) approach to density functional theory. However, this approach has thus far been rather ad hoc in its choice of the generic expression of the SRP functional to be used, and there is a need for better heuristic approaches to determining the mixing parameters contained in such expressions. Here we address these two issues. We investigate the ability of several mixed, parametrized density functional expressions combining exchange at the generalized gradient approximation (GGA) level with either GGA or nonlocal correlation to reproduce barrier heights for dissociative chemisorption on metal surfaces. For this, seven expressions of such mixed density functionals are tested on a database consisting of results for 16 systems taken from a recently published slightly larger database called SBH17. Three expressions are derived that exhibit high tunability and use correlation functionals that are either of the PBE GGA form or of one of two limiting nonlocal forms also describing the attractive van der Waals interaction in an approximate way. We also find that, for mixed density functionals incorporating GGA correlation, the optimum fraction of repulsive RPBE GGA exchange obtained with a specific GGA density functional is correlated with the charge-transfer parameter, which is equal to the difference in the work function of the metal surface and the electron affinity of the molecule. However, the correlation is generally not large and not large enough to obtain accurate guesses of the mixing parameter for the systems considered, suggesting that it does not give rise to a very effective search strategy.

Introducing the embedded random phase approximation: H2 dissociative adsorption on Cu(111) as an exemplar

The random phase approximation (RPA) as a means of treating electron correlation recently has been shown to outperform standard density functional theory (DFT) approximations in a variety of cases. However, the computational cost of the RPA is substantially more than DFT, especially when aiming to study extended surfaces. Properly accounting for sufficient surface ensemble size, Brillouin zone sampling, and vacuum separation of periodic images in standard periodic-planewave-based DFT code raises the cost to achieve converged results. Here, we show that sub-system embedding schemes enable use of the RPA for modeling heterogeneous reactions at reduced computational cost. We explore two different embedded RPA (emb-RPA) approaches, periodic emb-RPA and cluster emb-RPA. We use the (experimentally and theoretically) well-studied H2 dissociative adsorption on Cu(111) as our exemplar, and first perform full periodic RPA calculations as a benchmark. The full RPA results match well the semi-empirical barrier fit to experimental observables and others derived from high-level computations, e.g., from recent embedded n-electron valence second order perturbation theory [Zhao et al., J. Chem. Theory Comput. 16(11), 7078-7088 (2020)] and quantum Monte Carlo [Doblhoff-Dier et al., J. Chem. Theory Comput. 13(7), 3208-3219 (2017)] simulations. Among the two emb-RPA approaches tested, the cluster emb-RPA accurately reproduces the energy profile (maximum error of 50 meV along the reaction pathway) while reducing the computational cost by approximately two orders of magnitude. We therefore expect that the embedded cluster approach will enable wider RPA implementation in heterogeneous catalysis.

Accurate descriptions of molecule-surface interactions in electrocatalytic CO2 reduction on the copper surfaces

Copper-based catalysts play a pivotal role in many industrial processes and hold a great promise for electrocatalytic CO₂ reduction reaction into valuable chemicals and fuels. Towards the rational design of catalysts, the growing demand on theoretical study is seriously at odds with the low accuracy of the most widely used functionals of generalized gradient approximation. Here, we present results using a hybrid scheme that combines the doubly hybrid XYG3 functional and the periodic generalized gradient approximation, whose accuracy is validated against an experimental set on copper surfaces. A near chemical accuracy is established for this set, which, in turn, leads to a substantial improvement for the calculated equilibrium and onset potentials as against the experimental values for CO₂ reduction to CO on Cu(111) and Cu(100) electrodes. We anticipate that the easy use of the hybrid scheme will boost the predictive power for accurate descriptions of molecule-surface interactions in heterogeneous catalysis.

PtOx Cly (OH)z (H2 O)n Complexes under Oxidative and Reductive Conditions: Impact of the Level of Theory on Thermodynamic Stabilities

Platinum-based catalysts with Cl^- , OH^- , O^2- and H₂ O ligands, are involved in many industrial processes. Their final chemical properties are impacted by calcination and reduction applied during the preparation and activation steps. We investigate their stability under these reactive conditions with density functional theory (DFT). We benchmark various functionals (PBE-dDsC, optPBE, B3LYP, HSE06, PBE0, TPSS, RTPSS and SCAN) against ACFDT-RPA. PBE-dDsC is well adapted, although hybrid functionals are more accurate for redox reactions. Thermodynamic phase diagrams are determined by computing the chemical potential of the species as a function of temperature and partial pressures of H₂ O, HCl, O₂ and H₂ . The stability and nature of the Pt species are highly sensitive to the activation conditions. Under O₂ , high temperatures favour PtO₂ while under H₂ , platinum is easily reduced to Pt(0). Chlorine modifies the coordination sphere of platinum during calcination by stabilizing PtCl₄ and shifts the reduction of platinum to higher temperatures under H₂ .

Modeling Electrochemical Processes with Grand Canonical Treatment of Many-Body Perturbation Theory

Electrocatalysis plays a key role in sustainable energy conversion and storage. It is critical to model the grand canonical treatment of electrons, which accounts for the electrochemical potential explicitly, at the atomic scale and understand these reactions at electrified interfaces. However, such a grand canonical treatment for electrocatalytic modeling is in practice restricted to a treatment of electronic structure with density functional theory, and more accurate methods are in many cases desirable. Here, we develop an original workflow combining the grand canonical treatment of electrons with many-body perturbation theory in its random phase approximation (RPA). Using the potential dependent adsorption of carbon monoxide on the copper (100) facet, we show that the grand canonical RPA energetics provide the correct on-top Cu geometry for CO at reducing potential. The match with experimental results is significantly improved compared to the functionals at the generalized gradient approximation level, which is the most commonly used approximation for electrochemical applications. We expect this development to pave the way to further electrochemical applications using RPA.

Correlation-Consistent Gaussian Basis Sets for Solids Made Simple

The rapidly growing interest in simulating condensed-phase materials using quantum chemistry methods calls for a library of high-quality Gaussian basis sets suitable for periodic calculations. Unfortunately, most standard Gaussian basis sets commonly used in molecular simulation show significant linear dependencies when used in close-packed solids, leading to severe numerical issues that hamper the convergence to the complete basis set (CBS) limit, especially in correlated calculations. In this work, we revisit Dunning's strategy for construction of correlation-consistent basis sets and examine the relationship between accuracy and numerical stability in periodic settings. We find that limiting the number of primitive functions avoids the appearance of problematic small exponents while still providing smooth convergence to the CBS limit. As an example, we generate double-, triple-, and quadruple-ζ correlation-consistent Gaussian basis sets for periodic calculations with Goedecker-Teter-Hutter (GTH) pseudopotentials. Our basis sets cover the main-group elements from the first three rows of the periodic table. Especially for atoms on the left side of the periodic table, our basis sets are less diffuse than those used in molecular calculations. We verify the fast and reliable convergence to the CBS limit in both Hartree-Fock and post-Hartree-Fock (MP2) calculations, using a diverse test set of 19 semiconductors and insulators.