1
|
van Bree RAB, Gerrits N, Kroes GJ. Dissociative chemisorption of O 2 on Al(111): dynamics on a potential energy surface computed with a non-self-consistent screened hybrid density functional approach. Faraday Discuss 2024; 251:361-381. [PMID: 38787655 DOI: 10.1039/d3fd00165b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2024]
Abstract
Density functional theory (DFT) at the generalized gradient approximation (GGA) level is often considered the best compromise between feasibility and accuracy for reactions of molecules on metal surfaces. Recent work, however, strongly suggests that density functionals (DFs) based on GGA exchange are not able to describe molecule-metal surface reactions for which the work function of the metal surface minus the electron affinity of the molecule is less than 7 eV. Systems for which this is true exhibit an increased charge transfer from the metal to the molecule at the transition state, increasing the delocalisation of the electron density. This enlarged delocalisation can cause GGA-DFT to underestimate energy values relative to the gas-phase and thus underestimate the barrier height, similar to what has been observed for several gas-phase reactions. An example of such a molecule-metal surface system is O2 + Al(111). Following a similar strategy as for gas-phase reactions, previous work showed results of increased accuracy when using a screened hybrid DF for O2 + Al(111). However, even screened hybrid DFs are computationally expensive to use for metal surfaces. To resolve this, we test a non-self-consistent field (NSCF) screened hybrid DF approach. This approach computes screened hybrid DFT energies based on self-consistent-field (SCF) GGA electronic densities. Here, we explore the accuracy of the NSCF screened hybrid DF approach by implementing the NSCF HSE03-1/3x@RPBE DF for O2 + Al(111). We compute and analyse molecular beam sticking probabilities as well as a set of sticking probabilities for rotationally aligned O2. Our results show that the NSCF approach results in reaction probability curves that reproduce SCF results with near-chemical accuracy, suggesting that the NSCF approach can be used advantageously for exploratory purposes. An analysis of the potential energy surface and the barriers gives insight into the cause of the disagreement between the SCF and NSCF reaction probabilities and into the changes needed in theoretical modelling to further improve the description of the O2 + Al(111) system. Finally, the hole model yields fair agreement with dynamics results for the reaction probability curve, but results in an increased slope of the reaction probability curve compared to the molecular dynamics, with a shift to lower or higher energies depending on whether the vibrational energy of the molecule is included in the initial energy of the molecule or not.
Collapse
Affiliation(s)
- Robert A B van Bree
- Leiden Institute of Chemisty, Leiden University, Gorlaeus Labaratories, P.O. Box 9502 2300 RA, Leiden, The Netherlands.
| | - Nick Gerrits
- Leiden Institute of Chemisty, Leiden University, Gorlaeus Labaratories, P.O. Box 9502 2300 RA, Leiden, The Netherlands.
| | - Geert-Jan Kroes
- Leiden Institute of Chemisty, Leiden University, Gorlaeus Labaratories, P.O. Box 9502 2300 RA, Leiden, The Netherlands.
| |
Collapse
|
2
|
Gallagher C, Siddiqui W, Arnold T, Cheng C, Su E, Zhao Q. Benchmarking a Molecular Flake Model on the Road to Programmable Graphene-Based Single-Atom Catalysts. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2024; 128:2876-2883. [PMID: 38414836 PMCID: PMC10895666 DOI: 10.1021/acs.jpcc.3c07681] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/21/2023] [Revised: 01/19/2024] [Accepted: 01/30/2024] [Indexed: 02/29/2024]
Abstract
Single-atom catalysts (SACs) of embedding an active metal in nitrogen-doped graphene are emergent catalytic materials in various applications. The rational design of efficient SACs necessitates an electronic and mechanistic understanding of those materials with reliable quantum mechanical simulations. Conventional computational methods of modeling SACs involve using an infinite slab model with periodic boundary condition, limiting to the selection of generalized gradient approximations as the exchange correlation (XC) functional within density functional theory (DFT). However, these DFT approximations suffer from electron self-interaction error and delocalization error, leading to errors in predicted charge-transfer energetics. An alternative strategy is using a molecular flake model, which carved out the important catalytic center by cleaving C-C bonds and employing a hydrogen capping scheme to saturate the innocent dangling bonds at the molecular boundary. By doing so, we can afford more accurate hybrid XC functionals, or even high-level correlated wavefunction theory, to study those materials. In this work, we compared the structural, electronic, and catalytic properties of SACs simulated using molecular flake models and periodic slab models with first-row transition metals as the active sites. Molecular flake models successfully reproduced structural properties, including both global distortion and local metal-coordination environment, as well as electronic properties, including spin magnetic moments and metal partial charges, for all transition metals studied. In addition, we calculated CO binding strength as a descriptor for electrochemical CO2 reduction reactivity and noted qualitatively similar trends between two models. Using the computationally efficient molecular flake models, we investigated the effect of tuning Hartree-Fock exchange in a global hybrid functional on the CO binding strength and observed system-dependent sensitivities. Overall, our calculations provide valuable insights into the development of accurate and efficient computational tools to simulate SACs.
Collapse
Affiliation(s)
- Colin Gallagher
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Wali Siddiqui
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Tyler Arnold
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Carmen Cheng
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Eric Su
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Qing Zhao
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| |
Collapse
|
3
|
Powell A, Gerrits N, Tchakoua T, Somers MF, Busnengo HF, Meyer J, Kroes GJ, Doblhoff-Dier K. Best-of-Both-Worlds Predictive Approach to Dissociative Chemisorption on Metals. J Phys Chem Lett 2024; 15:307-315. [PMID: 38169287 PMCID: PMC10788952 DOI: 10.1021/acs.jpclett.3c02972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2023] [Revised: 12/21/2023] [Accepted: 12/28/2023] [Indexed: 01/05/2024]
Abstract
Predictive capability, accuracy, and affordability are essential features of a theory that is capable of describing dissociative chemisorption on a metal surface. This type of reaction is important for heterogeneous catalysis. Here we present an approach in which we use diffusion Monte Carlo (DMC) to pin the minimum barrier height and construct a density functional that reproduces this value. This predictive approach allows the construction of a potential energy surface at the cost of density functional theory while retaining near DMC accuracy. Scrutinizing effects of energy dissipation and quantum tunneling, dynamics calculations suggest the approach to be of near chemical accuracy, reproducing molecular beam sticking experiments for the showcase H2 + Al(110) system to ∼1.4 kcal/mol.
Collapse
Affiliation(s)
- Andrew
D. Powell
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Nick Gerrits
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Theophile Tchakoua
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Mark F. Somers
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Heriberto F. Busnengo
- Instituto
de Física Rosario (IFIR), CONICET-UNR, 2000 Rosario, Argentina
- Facultad
de Ciencias Exatas, Ingeniería y
Agrimensura, UNR, 2000 Rosario, Argentina
| | - Jörg Meyer
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Katharina Doblhoff-Dier
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| |
Collapse
|
4
|
Tchakoua T, Jansen T, van Nies Y, van den Elshout RFA, van Boxmeer BAB, Poort SP, Ackermans MG, Beltrão GS, Hildebrand SA, Beekman SEJ, van der Drift T, Kaart S, Šantić A, Spuijbroek EE, Gerrits N, Somers MF, Kroes GJ. Constructing Mixed Density Functionals for Describing Dissociative Chemisorption on Metal Surfaces: Basic Principles. J Phys Chem A 2023; 127:10481-10498. [PMID: 38051300 PMCID: PMC10726370 DOI: 10.1021/acs.jpca.3c01932] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- Théophile Tchakoua
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Tim Jansen
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Youri van Nies
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | | | - Bart A B van Boxmeer
- Faculty of Applied Engineering, University of Antwerp, Groenenborgerlaan 171, 2020 Antwerpen, Belgium
| | - Saskia P Poort
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Michelle G Ackermans
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Gabriel Spiller Beltrão
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Stefan A Hildebrand
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Steijn E J Beekman
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Thijs van der Drift
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Sam Kaart
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Anthonie Šantić
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Esmee E Spuijbroek
- Department of Chemical Engineering, Delft University of Technology, 2629 HZ Delft, The Netherlands
| | - Nick Gerrits
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Mark F Somers
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
5
|
Wei Z, Martirez JMP, Carter EA. Introducing the embedded random phase approximation: H2 dissociative adsorption on Cu(111) as an exemplar. J Chem Phys 2023; 159:194108. [PMID: 37971031 DOI: 10.1063/5.0181229] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2023] [Accepted: 10/16/2023] [Indexed: 11/19/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Ziyang Wei
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, USA
| | - John Mark P Martirez
- Applied Materials and Sustainability Sciences, Princeton Plasma Physics Laboratory, Princeton, New Jersey 08540-6655, USA
| | - Emily A Carter
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, USA
- Applied Materials and Sustainability Sciences, Princeton Plasma Physics Laboratory, Princeton, New Jersey 08540-6655, USA
- Andlinger Center for Energy and the Environment and Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
6
|
Kim SY, Park JW. Approximate Excited-State Geometry Optimization with the State-Averaged Adaptive Sampling Configuration Interaction Algorithm with Large Active Spaces. J Chem Theory Comput 2023; 19:7260-7272. [PMID: 37800852 DOI: 10.1021/acs.jctc.3c00808] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/07/2023]
Abstract
The selected configuration interaction (SCI) wave function is a useful approximation to the full configuration interaction (FCI) one. The adaptive sampling CI (ASCI) method is a deterministic SCI method. By combining ASCI and orbital optimization, the ASCI self-consistent field (ASCI-SCF) method, which is an approximation of the complete active space self-consistent field (CASSCF) method, can be formulated as well. However, their applicability has been tested mainly on the systems in their electronically ground states. In this work, we implement the state-average (SA) ansatz in ASCI-SCF calculations to calculate excited states. We also derive expressions for the approximate analytical gradient and implement them as a computer program. We demonstrate the applicability of the current method for calculating vertical and adiabatic excitation energies and optimizing the molecular geometries of thermally activated delayed fluorescence (TADF) molecules.
Collapse
Affiliation(s)
- So Yeon Kim
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| | - Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|
7
|
He N, Huang M, Evangelista FA. CO Inversion on a NaCl(100) Surface: A Multireference Quantum Embedding Study. J Phys Chem A 2023; 127:1975-1987. [PMID: 36799901 PMCID: PMC9986868 DOI: 10.1021/acs.jpca.2c05844] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
Abstract
We develop a multireference quantum embedding model to investigate a recent experimental observation of the isomerization of vibrationally excited CO molecules on a NaCl(100) surface [Science 2020, 367, 175-178]. To explore this mechanism, we built a reduced potential energy surface of CO interacting with NaCl(100) using a second-order multireference perturbation theory, modeling the adsorbate-surface interaction with our previously developed active space embedding theory (ASET). We considered an isolated CO molecule on NaCl(100) and a high-coverage CO monolayer (1/1), and for both we generated potential energy surfaces parametrized by the CO stretching, adsorption, and inversion coordinates. These surfaces are used to determine stationary points and adsorption energies and to perform a vibrational analysis of the states relevant to the inversion mechanism. We found that for near-equilibrium bond lengths, CO adsorbed in the C-down configuration is lower in energy than in the O-down configuration. Stretching of the C-O bond reverses the energetic order of these configurations, supporting the accepted isomerization mechanism. The vibrational constants obtained from these potential energy surfaces show a small (< 10 cm-1) blue- and red-shift for the C-down and O-down configurations, respectively, in agreement with experimental assignments and previous theoretical studies. Our vibrational analysis of the monolayer case suggests that the O-down configuration is energetically more stable than the C-down one beyond the 16th vibrational excited state of CO, a value slightly smaller than the one from quasi-classical trajectory simulations (22nd) and consistent with the experiment. Our analysis suggests that CO-CO interactions in the monolayer play an important role in stabilizing highly vibrationally excited states in the O-down configuration and reducing the barrier between the C-down and O-down geometries, therefore playing a crucial role in the inversion mechanism.
Collapse
Affiliation(s)
- Nan He
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Meng Huang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
8
|
Chen Z, Liu Z, Xu X. Accurate descriptions of molecule-surface interactions in electrocatalytic CO 2 reduction on the copper surfaces. Nat Commun 2023; 14:936. [PMID: 36807556 PMCID: PMC9941474 DOI: 10.1038/s41467-023-36695-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2022] [Accepted: 02/14/2023] [Indexed: 02/22/2023] Open
Abstract
Copper-based catalysts play a pivotal role in many industrial processes and hold a great promise for electrocatalytic CO2 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 CO2 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.
Collapse
Affiliation(s)
- Zheng Chen
- grid.8547.e0000 0001 0125 2443Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, 200433 Shanghai, People’s Republic of China
| | - Zhangyun Liu
- grid.8547.e0000 0001 0125 2443Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, 200433 Shanghai, People’s Republic of China
| | - Xin Xu
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, 200433, Shanghai, People's Republic of China. .,Hefei National Laboratory, 230088, Hefei, P. R. China.
| |
Collapse
|
9
|
Zhang X, Zhu L, Hou Q, Guan J, Lu Y, Keal TW, Buckeridge J, Catlow CRA, Sokol AA. Toward a Consistent Prediction of Defect Chemistry in CeO 2. CHEMISTRY OF MATERIALS : A PUBLICATION OF THE AMERICAN CHEMICAL SOCIETY 2023; 35:207-227. [PMID: 36644213 PMCID: PMC9835833 DOI: 10.1021/acs.chemmater.2c03019] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Revised: 11/25/2022] [Indexed: 05/10/2023]
Abstract
Polarizable shell-model potentials are widely used for atomic-scale modeling of charged defects in solids using the Mott-Littleton approach and hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) embedded-cluster techniques. However, at the pure MM level of theory, the calculated defect energetics may not satisfy the requirement of quantitative predictions and are limited to only certain charged states. Here, we proposed a novel interatomic potential development scheme that unifies the predictions of all relevant charged defects in CeO2 based on the Mott-Littleton approach and QM/MM electronic-structure calculations. The predicted formation energies of oxygen vacancies accompanied by different excess electron localization patterns at the MM level of theory reach the accuracy of density functional theory (DFT) calculations using hybrid functionals. The new potential also accurately reproduces a wide range of physical properties of CeO2, showing excellent agreement with experimental and other computational studies. These findings provide opportunities for accurate large-scale modeling of the partial reduction and nonstoichiometry in CeO2, as well as a prototype for developing robust interatomic potentials for other defective crystals.
Collapse
Affiliation(s)
- Xingfan Zhang
- Kathleen
Lonsdale Materials Chemistry, Department of Chemistry, University College London, LondonWC1H 0AJ, United Kingdom
| | - Lei Zhu
- Kathleen
Lonsdale Materials Chemistry, Department of Chemistry, University College London, LondonWC1H 0AJ, United Kingdom
| | - Qing Hou
- Kathleen
Lonsdale Materials Chemistry, Department of Chemistry, University College London, LondonWC1H 0AJ, United Kingdom
- Institute
of Photonic Chips, University of Shanghai
for Science and Technology, Shanghai200093, China
| | - Jingcheng Guan
- Kathleen
Lonsdale Materials Chemistry, Department of Chemistry, University College London, LondonWC1H 0AJ, United Kingdom
| | - You Lu
- Scientific
Computing Department, STFC Daresbury Laboratory, Warrington, CheshireWA4 4AD, United Kingdom
| | - Thomas W. Keal
- Scientific
Computing Department, STFC Daresbury Laboratory, Warrington, CheshireWA4 4AD, United Kingdom
| | - John Buckeridge
- School
of Engineering, London South Bank University, LondonSE1 OAA, United Kingdom
| | - C. Richard A. Catlow
- Kathleen
Lonsdale Materials Chemistry, Department of Chemistry, University College London, LondonWC1H 0AJ, United Kingdom
- School
of Chemistry, Cardiff University, Park Place, CardiffCF10 1AT, United
Kingdom
| | - Alexey A. Sokol
- Kathleen
Lonsdale Materials Chemistry, Department of Chemistry, University College London, LondonWC1H 0AJ, United Kingdom
| |
Collapse
|
10
|
Tchakoua T, Gerrits N, Smeets EWF, Kroes GJ. SBH17: Benchmark Database of Barrier Heights for Dissociative Chemisorption on Transition Metal Surfaces. J Chem Theory Comput 2022; 19:245-270. [PMID: 36529979 PMCID: PMC9835835 DOI: 10.1021/acs.jctc.2c00824] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Accurate barriers for rate controlling elementary reactions on metal surfaces are key to understanding, controlling, and predicting the rate of heterogeneously catalyzed processes. While barrier heights for gas phase reactions have been extensively benchmarked, dissociative chemisorption barriers for the reactions of molecules on metal surfaces have received much less attention. The first database called SBH10 and containing 10 entries was recently constructed based on the specific reaction parameter approach to density functional theory (SRP-DFT) and experimental results. We have now constructed a new and improved database (SBH17) containing 17 entries based on SRP-DFT and experiments. For this new SBH17 benchmark study, we have tested three algorithms (high, medium, and light) for calculating barrier heights for dissociative chemisorption on metals, which we have named for the amount of computational effort involved in their use. We test the performance of 14 density functionals at the GGA, GGA+vdW-DF, and meta-GGA rungs. Our results show that, in contrast with the previous SBH10 study where the BEEF-vdW-DF2 functional seemed to be most accurate, the workhorse functional PBE and the MS2 density functional are the most accurate of the GGA and meta-GGA functionals tested. Of the GGA+vdW functionals tested, the SRP32-vdW-DF1 functional is the most accurate. Additionally, we found that the medium algorithm is accurate enough for assessing the performance of the density functionals tested, while it avoids geometry optimizations of minimum barrier geometries for each density functional tested. The medium algorithm does require metal lattice constants and interlayer distances that are optimized separately for each functional. While these are avoided in the light algorithm, this algorithm is found not to give a reliable description of functional performance. The combination of relative ease of use and demonstrated reliability of the medium algorithm will likely pave the way for incorporation of the SBH17 database in larger databases used for testing new density functionals and electronic structure methods.
Collapse
Affiliation(s)
- T. Tchakoua
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands
| | - N. Gerrits
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands,PLASMANT,
Department of Chemistry, University of Antwerp, BE-2610Antwerp, Belgium
| | - E. W. F. Smeets
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands,ALTEN
Nederland, Technology, Fascinatio Boulevard 582, 2909 VACapelle a/d IJssel, The Netherlands
| | - G.-J. Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RALeiden, The Netherlands,E-mail: . Phone: +31 71 527 4396
| |
Collapse
|
11
|
Charting C–C coupling pathways in electrochemical CO
2
reduction on Cu(111) using embedded correlated wavefunction theory. Proc Natl Acad Sci U S A 2022; 119:e2202931119. [PMID: 36306330 PMCID: PMC9636923 DOI: 10.1073/pnas.2202931119] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The electrochemical CO
2
reduction reaction (CO
2
RR) powered by excess zero-carbon-emission electricity to produce especially multicarbon (C
2+
) products could contribute to a carbon-neutral to carbon-negative economy. Foundational to the rational design of efficient, selective CO
2
RR electrocatalysts is mechanistic analysis of the best metal catalyst thus far identified, namely, copper (Cu), via quantum mechanical computations to complement experiments. Here, we apply embedded correlated wavefunction (ECW) theory, which regionally corrects the electron exchange-correlation error in density functional theory (DFT) approximations, to examine multiple C–C coupling steps involving adsorbed CO (*CO) and its hydrogenated derivatives on the most ubiquitous facet, Cu(111). We predict that two adsorbed hydrogenated CO species, either *COH or *CHO, are necessary precursors for C–C bond formation. The three kinetically feasible pathways involving these species yield all three possible products: *COH–CHO, *COH–*COH, and *OCH–*OCH. The most kinetically favorable path forms *COH–CHO. In contrast, standard DFT approximations arrive at qualitatively different conclusions, namely, that only *CO and *COH will prevail on the surface and their C–C coupling paths produce only *COH–*COH and *CO–*CO, with a preference for the first product. This work demonstrates the importance of applying qualitatively and quantitatively accurate quantum mechanical method to simulate electrochemistry in order ultimately to shed light on ways to enhance selectivity toward C
2+
product formation via CO
2
RR electrocatalysts.
Collapse
|
12
|
Vitillo JG, Cramer CJ, Gagliardi L. Multireference Methods are Realistic and Useful Tools for Modeling Catalysis. Isr J Chem 2022. [DOI: 10.1002/ijch.202100136] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Jenny G. Vitillo
- Department of Science and High Technology and INSTM Università degli Studi dell'Insubria Via Valleggio 9 I-22100 Como Italy
| | - Christopher J. Cramer
- Underwriters Laboratories Inc. 333 Pfingsten Road Northbrook Illinois 60602 United States
| | - Laura Gagliardi
- Department of Chemistry Pritzker School of Molecular Engineering James Franck Institute University of Chicago Chicago Illinois 60637 United States
| |
Collapse
|
13
|
Abstract
Quantum embedding schemes are a promising way to extend multireference computations to large molecules with strong correlation effects localized on a small number of atoms. This work introduces a second-order active-space embedding theory [ASET(2)] which improves upon mean-field frozen embedding by treating fragment-environment interactions via an approximate canonical transformation. The canonical transformation employed in ASET(2) is formulated using the driven similarity renormalization group. The ASET(2) scheme is benchmarked on the N═N bond dissociation in pentyldiazene, the S0 to S1 excitation in 1-octene, and the interaction energy of the O2-benzene complex. The ASET(2) explicit treatment of fragment-environment interactions beyond the mean-field level generally improves the accuracy of embedded computations, and it becomes necessary to achieve an accurate description of excitation energies of 1-octene and the singlet-triplet gap of the O2-benzene complex.
Collapse
Affiliation(s)
- Nan He
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chenyang Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
14
|
Gerhards L, Klüner T. Theoretical investigation of CH-bond activation by photocatalytic excited SO 2 and the effects of C-, N-, S-, and Se-doped TiO 2. Phys Chem Chem Phys 2022; 24:2051-2069. [PMID: 35014643 DOI: 10.1039/d1cp04335h] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The photocatalytic sulfoxidation on TiO2 discovered by Parrino et al. represents a new, interesting and lower energy route for the synthesis of sulfonic acids. Sulfonic acids are important precursors for dyes, detergents and drugs. In the commonly known industrial process, SO2 and a specific hydrocarbon are converted into sulfonic acids using high-energy UV light. In this reaction, SO2 is excited into a metastable triplet state (3SO2), which has the potential to activate a CH-bond of hydrocarbons and start a radical reaction cycle. By introducing TiO2 as a photocatalyst, it has been shown that visible light can be used for the synthesis. This offers the potential to be a cost-effective reaction approach for industrial use. However, experimental studies indicate that the initial excitation mechanism of SO2 on TiO2 is significantly different from the catalyst-free mechanism. Parrino et al. were able to reveal first evidence for the existence of a charge-transfer process from SO2 to the TiO2 surface by means of electrochemical experiments. First theoretical investigations from first principles were able to further substantiate the existence of a charge-transfer. However, to fully understand this mechanism, more accurate methods such as Time Dependent Density Functional Theory (TD-DFT) or ab initio multireference methods such as the Complete Active Space Self Consistent Field (CASSCF) method are required. Furthermore, after understanding the charge-transfer mechanism, the introduction of dopants into TiO2 can be investigated in order to possibly redshift the excitation energy. This might open the route to using lower energy light for the sulfoxidation of hydrocarbons on TiO2 as a new potential industrial reaction for the synthesis of sulfonic acids. In this work, we will study the initial step of the photocatalytic sulfoxidation of hydrocarbons using the TD-DFT and CASSCF methods by using a combined approach consisting of calculations with periodic boundary conditions and a newly constructed embedded cluster model. Furthermore, we will explore the effects of doping by introducing four heteroatoms (C, N, S, and Se) into the TiO2 surfaces anatase[101] and rutile[110] to find a possible enhancement of the photocatalytic reactivity by lowering the electronic excitation energy.
Collapse
Affiliation(s)
- Luca Gerhards
- School of Mathematics and Science Chemistry Department Carl von Ossietzky University of Oldenburg, 26111 Oldenburg, Germany.
| | - Thorsten Klüner
- School of Mathematics and Science Chemistry Department Carl von Ossietzky University of Oldenburg, 26111 Oldenburg, Germany.
| |
Collapse
|
15
|
Mejuto-Zaera C, Tzeli D, Williams-Young D, Tubman NM, Matoušek M, Brabec J, Veis L, Xantheas SS, de Jong WA. The Effect of Geometry, Spin, and Orbital Optimization in Achieving Accurate, Correlated Results for Iron-Sulfur Cubanes. J Chem Theory Comput 2022; 18:687-702. [PMID: 35034448 DOI: 10.1021/acs.jctc.1c00830] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Iron-sulfur clusters comprise an important functional motif in the catalytic centers of biological systems, capable of enabling important chemical transformations at ambient conditions. This remarkable capability derives from a notoriously complex electronic structure that is characterized by a high density of states that is sensitive to geometric changes. The spectral sensitivity to subtle geometric changes has received little attention from correlated, large active space calculations, owing partly to the exceptional computational complexity for treating these large and correlated systems accurately. To provide insight into this aspect, we report the first Complete Active Space Self Consistent Field (CASSCF) calculations for different geometries of the [Fe(II/III)4S4(SMe)4]-2 clusters using two complementary, correlated solvers: spin-pure Adaptive Sampling Configuration Interaction (ASCI) and Density Matrix Renormalization Group (DMRG). We find that the previously established picture of a double-exchange driven magnetic structure, with minute energy gaps (<1 mHa) between consecutive spin states, has a weak dependence on the underlying geometry. However, the spin gap between the singlet and the spin state 2S + 1 = 19, corresponding to a maximal number of Fe-d electrons being unpaired and of parallel spin, is strongly geometry dependent, changing by a factor of 3 upon slight deformations that are still within biologically relevant parameters. The CASSCF orbital optimization procedure, using active spaces as large as 86 electrons in 52 orbitals, was found to reduce this gap compared to typical mean-field orbital approaches. Our results show the need for performing large active space calculations to unveil the challenging electronic structure of these complex catalytic centers and should serve as accurate starting points for fully correlated treatments upon inclusion of dynamical correlation outside the active space.
Collapse
Affiliation(s)
- Carlos Mejuto-Zaera
- University of California, Berkeley, California 94720, United States.,Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Demeter Tzeli
- Laboratory of Physical Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, Panepistimiopolis Zografou, Athens 15784, Greece.,Theoretical and Physical Chemistry Institute, National Hellenic Research Foundation, Vas. Constantinou 48, Athens 11635, Greece
| | - David Williams-Young
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Norm M Tubman
- Quantum Artificial Intelligence Lab. (QuAIL), Exploration Technology Directorate, NASA Ames Research Center, Moffett Field, California 94035, United States
| | - Mikuláš Matoušek
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jiri Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Sotiris S Xantheas
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, United States.,Department of Chemistry, University of Washington, Seattle, Washington 98185, United States
| | - Wibe A de Jong
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
16
|
Geng X, Liu X, Mawella-Vithanage L, Hewa-Rahinduwage CC, Zhang L, Brock SL, Tan T, Luo L. Photoexcited NO 2 Enables Accelerated Response and Recovery Kinetics in Light-Activated NO 2 Gas Sensing. ACS Sens 2021; 6:4389-4397. [PMID: 34784175 DOI: 10.1021/acssensors.1c01694] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Slow response and recovery kinetics is a major challenge for practical room-temperature NO2 gas sensing. Here, we report the use of visible light illumination to significantly shorten the response and recovery times of a PbSe quantum dot (QD) gel sensor by 21% (to 27 s) and 63% (to 102 s), respectively. When combined with its high response (0.04%/ppb) and ultralow limit of detection (3 ppb), the reduction in response and recovery time makes the PbSe QD gel sensor among the best p-type room-temperature NO2 sensors reported to date. A combined experimental and theoretical investigation reveals that the accelerated response and recovery time is primarily due to photoexcitation of NO2 gaseous molecules and adsorbed NO2 on the gel surface, rather than the excitation of the semiconductor sensing material, as suggested by the currently prevailing light-activated gas-sensing theory. Furthermore, we find that the extent of improvement attained in the recovery speed also depends on the distribution of excited electrons in the adsorbed NO2/QD gel system. This work suggests that the design of light-activated sensor platforms may benefit from a careful assessment of the photophysics of the analyte in the gas phase and when adsorbed onto the semiconductor surface.
Collapse
Affiliation(s)
- Xin Geng
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - Xiaolong Liu
- Laboratory of Theoretical and Computational Nanoscience, CAS Center for Excellence in Nanoscience, National Center for Nanoscience and Technology, Chinese Academy of Sciences, Beijing 100190, China
| | | | | | - Liang Zhang
- School of Vehicle and Mobility, Tsinghua University, Beijing 100084, China
- Center for Combustion Energy, Tsinghua University, Beijing 100084, China
| | - Stephanie L. Brock
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - Ting Tan
- Laboratory of Theoretical and Computational Nanoscience, CAS Center for Excellence in Nanoscience, National Center for Nanoscience and Technology, Chinese Academy of Sciences, Beijing 100190, China
| | - Long Luo
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| |
Collapse
|
17
|
Rybkin VV. Formulation and Implementation of Density Functional Embedding Theory Using Products of Basis Functions. J Chem Theory Comput 2021; 17:3995-4005. [PMID: 34048247 DOI: 10.1021/acs.jctc.1c00175] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The representation of embedding potential using products of atomic orbital basis functions has been developed in the context of density functional embedding theory. The formalism allows to treat pseudopotential and all-electron calculations on the same footing and enables simple transfer of the embedding potential in a compact matrix form. In addition, a cost-reduction procedure for the basis set and potential reduction based on population analysis has been proposed. Implemented for the condensed-phase and molecular systems within Gaussian and plane-waves and Gaussian and augmented plane-waves formalisms, the scheme has been tested for proton-transfer reactions in the cluster and the condensed phase and projected density of states of carbon monoxide adsorbed on platinum surface. With the computational scaling of the embedding potential optimization similar to that of hybrid density functional theory with a significantly reduced prefactor, the method allows for large-scale applications to extended systems.
Collapse
Affiliation(s)
- Vladimir V Rybkin
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, Zürich 8057, Switzerland
| |
Collapse
|
18
|
Martirez JMP, Carter EA. Projector-Free Capped-Fragment Scheme within Density Functional Embedding Theory for Covalent and Ionic Compounds. J Chem Theory Comput 2021; 17:4105-4121. [DOI: 10.1021/acs.jctc.1c00285] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- John Mark P. Martirez
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Emily A. Carter
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
19
|
Martirez JMP, Carter EA. Metal-to-Ligand Charge-Transfer Spectrum of a Ru-Bipyridine-Sensitized TiO 2 Cluster from Embedded Multiconfigurational Excited-State Theory. J Phys Chem A 2021; 125:4998-5013. [PMID: 34077662 DOI: 10.1021/acs.jpca.1c02628] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Understanding optical properties of the dye molecule in dye-sensitized solar cells (DSSCs) from first-principles quantum mechanics can contribute to improving the efficiency of such devices. While density functional theory (DFT) and time-dependent DFT have been pivotal in simulating optoelectronic properties of photoanodes used in DSSCs at the atomic scale, questions remain regarding DFT's adequacy and accuracy to furnish critical information needed to understand the various excited-state processes involved. Here, we simulate the absorption spectra of a dye-sensitized solar cell analogue, comprised of a Ru-bipyridine (Ru-bpy) dye molecule and a small TiO2 cluster via DFT and via an accurate embedded correlated wavefunction (CW) theory. We generated CW spectra for the adsorbed Ru-bpy dye via a recently introduced capped density functional embedding theory or capped-DFET (to generate the embedding potential that accounts for the interaction of the molecule and the TiO2 cluster). We then combined capped-DFET with the accurate but expensive multiconfigurational complete active space second-order perturbation theory (CASPT2)-embedded CASPT2. Because the CW theory is conducted on only a portion of the total system in the presence of an embedding potential that describes that portion's interaction with its environment, we efficiently obtain CW-quality predictions that reflect local properties of the entire system. Specifically, for example, with capped-DFET and embedded CW theory, we can simulate accurately a plethora of metal-to-ligand charge-transfer excited properties at a manageable computational cost. Here, we predict detailed electronic spectra within the visible region, featuring the lowest three singlet and triplet excited states, along with predictions of the singlets' lifetimes. We illustrated these results using a Jablonski diagram that show the relative energy position of the singlet and longer-lived triplet excited states and analyzed and proposed relaxation paths for the excited state corresponding to the most intense but short-lived absorption (interconversion, intersystem crossing, fluorescence, and phosphorescence) that may lead to longer-lived excited states necessary for efficient charge separation required to generate current in solar cells.
Collapse
Affiliation(s)
- John Mark P Martirez
- Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, California 90095-1592, United States
| | - Emily A Carter
- Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, California 90095-1592, United States.,Office of the Chancellor, University of California, Los Angeles, Box 951405, Los Angeles, California 90095-1405, United States.,Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, United States
| |
Collapse
|
20
|
Chilkuri VG, Neese F. Comparison of Many-Particle Representations for Selected Configuration Interaction: II. Numerical Benchmark Calculations. J Chem Theory Comput 2021; 17:2868-2885. [PMID: 33886300 PMCID: PMC8279407 DOI: 10.1021/acs.jctc.1c00081] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
The present work
is the second part in our three-part series on
the comparison of many-particle representations for the selected configuration
interaction (CI) method. In this work, we present benchmark calculations
based on our selected CI program called the iterative configuration
expansion (ICE) that is inspired by the CIPSI method of Malrieu and
co-workers (Malrieu1973, 58, ( (12), ), 5745−5759). We describe the main parameters
that enter in this algorithm and perform benchmark calculations on
a set of 21 small molecules and compare ground state energies with
full configuration interaction (FCI) results (FCI21 test set). The
focus is the comparison of the performance of three different types
of many-particle basis functions (MPBFs): (1) individual Slater determinants
(DETS), (2) individual spin-adapted configuration state functions
(CSFs), and (3) all CSFs of a given total spin that can be generated
from spatial configurations (CFGs). An analysis of the cost of the
calculation in terms of the number of wavefunction parameters and
the energy error is evaluated for the DET-, CFG-, and CSF-based ICE.
The main differences for the three many-particle basis representations
show up in the number of wavefunction parameters and the rate of convergence
toward the FCI limit with the thresholds of the ICE. Next, we analyze
the best way to extrapolate the ICE energies toward the FCI results
as a function of the thresholds. The efficiency of the extrapolation
is investigated relative to the FCI21 test set as well as near FCI
calculations on three moderately sized hydrocarbon molecules CH4, C2H4, and C4H6. Finally, we comment on the size-inconsistency error for the three
many-particle representations and compare it with the error in the
total energy. The implication for selected CI implementations with
any of the three many-particle representations is discussed.
Collapse
Affiliation(s)
- Vijay Gopal Chilkuri
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, Mülheim an der Ruhr 45470, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, Mülheim an der Ruhr 45470, Germany
| |
Collapse
|
21
|
Zhao Q, Martirez JMP, Carter EA. Revisiting Understanding of Electrochemical CO 2 Reduction on Cu(111): Competing Proton-Coupled Electron Transfer Reaction Mechanisms Revealed by Embedded Correlated Wavefunction Theory. J Am Chem Soc 2021; 143:6152-6164. [PMID: 33851840 DOI: 10.1021/jacs.1c00880] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Copper (Cu) electrodes, as the most efficacious of CO2 reduction reaction (CO2RR) electrocatalysts, serve as prototypes for determining and validating reaction mechanisms associated with electrochemical CO2 reduction to hydrocarbons. As in situ electrochemical mechanism determination by experiments is still out of reach, such mechanistic analysis typically is conducted using density functional theory (DFT). The semilocal exchange-correlation (XC) approximations most often used to model such catalysis unfortunately engender a basic error: predicting the wrong adsorption site for CO (a key CO2RR intermediate) on the most ubiquitous facet of Cu, namely, Cu(111). This longstanding inconsistency casts lingering doubt on previous DFT predictions of the attendant CO2RR kinetics. Here, we apply embedded correlated wavefunction (ECW) theory, which corrects XC functional error, to study the CO2RR on Cu(111) via both surface hydride (*H) transfer and proton-coupled electron transfer (PCET). We predict that adsorbed CO (*CO) reduces almost equally to two intermediates, namely, hydroxymethylidyne (*COH) and formyl (*CHO) at -0.9 V vs the RHE. In contrast, semilocal DFT approximations predict a strong preference for *COH. With increasing applied potential, the dominance of *COH (formed via potential-independent surface *H transfer) diminishes, switching to the competitive formation of both *CHO and *COH (both formed via potential-dependent PCET). Our results also demonstrate the importance of including explicitly modeled solvent molecules in predicting electron-transfer barriers and reveal the pitfalls of overreliance on simple surface *H transfer models of reduction reactions.
Collapse
Affiliation(s)
- Qing Zhao
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, United States
| | - John Mark P Martirez
- Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, California 90095-1592, United States
| | - Emily A Carter
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, United States.,Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, California 90095-1592, United States.,Office of the Chancellor, Box 951405, University of California, Los Angeles, Los Angeles, California 90095-1405, United States
| |
Collapse
|
22
|
Morales‐García Á, Viñes F, Gomes JRB, Illas F. Concepts, models, and methods in computational heterogeneous catalysis illustrated through
CO
2
conversion. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1530] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Ángel Morales‐García
- Departament de Ciència de Materials i Química Física & Institut de Química Teòrica i Computacional (IQTCUB) Universitat de Barcelona Barcelona Spain
| | - Francesc Viñes
- Departament de Ciència de Materials i Química Física & Institut de Química Teòrica i Computacional (IQTCUB) Universitat de Barcelona Barcelona Spain
| | - José R. B. Gomes
- CICECO—Aveiro Institute of Materials, Department of Chemistry University of Aveiro Aveiro Portugal
| | - Francesc Illas
- Departament de Ciència de Materials i Química Física & Institut de Química Teòrica i Computacional (IQTCUB) Universitat de Barcelona Barcelona Spain
| |
Collapse
|
23
|
Park JW. Second-Order Orbital Optimization with Large Active Spaces Using Adaptive Sampling Configuration Interaction (ASCI) and Its Application to Molecular Geometry Optimization. J Chem Theory Comput 2021; 17:1522-1534. [DOI: 10.1021/acs.jctc.0c01292] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|