1
|
Ye HZ, Berkelbach TC. Adsorption and vibrational spectroscopy of CO on the surface of MgO from periodic local coupled-cluster theory. Faraday Discuss 2024. [PMID: 39049598 DOI: 10.1039/d4fd00041b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/27/2024]
Abstract
The adsorption of CO on the surface of MgO has long been a model problem in surface chemistry. Here, we report periodic Gaussian-based calculations for this problem using second-order perturbation theory (MP2) and coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)], with the latter two performed using a recently developed extension of the local natural orbital approximation to problems with periodic boundary conditions. The low cost of periodic local correlation calculations allows us to calculate the full CCSD(T) binding curve of CO approaching the surface of MgO (and thus the adsorption energy) and the two-dimensional potential energy surface (PES) as a function of the distance from the surface and the CO stretching coordinate. From the PES, we obtain the fundamental vibrational frequency of CO on MgO, whose shift from the gas phase value is a common experimental probe of surface adsorption. We find that CCSD(T) correctly predicts a positive frequency shift upon adsorption of +14.7 cm-1, in excellent agreement with the experimental shift of +14.3 cm-1. We use our CCSD(T) results to assess the accuracy of MP2, CCSD, and several density functional theory (DFT) approximations, including exchange correlation functionals and dispersion corrections. We find that MP2 and CCSD yield reasonable binding energies and frequency shifts, whereas many DFT calculations overestimate the magnitude of the adsorption energy by 5-15 kJ mol-1 and predict a negative frequency shift of about -20 cm-1, which we attribute to self-interaction-induced delocalization errors that are mildly ameliorated with hybrid functionals. Our findings highlight the accuracy and computational efficiency of the periodic local correlation for the simulation of surface chemistry with accurate wavefunction methods.
Collapse
Affiliation(s)
- Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, NY 10027, USA.
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, NY 10027, USA.
- Initiative for Computational Catalysis, Flatiron Institute, New York, NY 10010, USA
| |
Collapse
|
2
|
Sun Y. Accelerating Density Matrix Embedding with Stochastic Density Fitting Theory: An Application to Hydrogen Bonded Clusters. J Chem Theory Comput 2024. [PMID: 39029091 DOI: 10.1021/acs.jctc.4c00489] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/21/2024]
Abstract
In this work, we demonstrate how using semistochastic density fitting (ss-DF) can accelerate self-consistent density matrix embedding theory (DMET) calculations by reducing the number of auxiliary orbitals in the three-indexed DF integrals. This reduction results in significant time savings when building the Hartree-Fock (HF) Coulomb and Exchange Matrices and in transforming integrals from the atomic orbital (AO) basis to the embedding orbital (EO) basis. We apply ss-DF to a range of hydrogen-bonded clusters to showcase its effectiveness. First, we examine how the amount of deterministic space impacts the quality of the calculation in a (H2O)10 cluster. Next, we test the computational efficiency of ss-DF compared to deterministic DF (d-DF) in water clusters containing 6-30 water molecules using a triple-ζ basis set. Finally, we perform numerical structural optimizations on water and hydrogen fluoride clusters, revealing that DMET can recover weak interactions using a back-transformed energy formula. This work demonstrates the potential of using stochastic resolution of identity in quantum embedding theories and highlights its capability to recover weak interactions effectively.
Collapse
Affiliation(s)
- Yi Sun
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
3
|
Yang MY, Wu XP. Level-Shifted Embedded Cluster Method for Modeling the Chemistry of Metal Oxides. J Chem Theory Comput 2024. [PMID: 38300767 DOI: 10.1021/acs.jctc.3c01123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2024]
Abstract
The embedded cluster method has been used extensively in the study of the chemical and physical properties of metal oxides. This method has been a popular tool due to its relatively high accuracy and low computational cost. An even more promising option may entail integrating the embedded cluster method with the combined quantum mechanical and molecular mechanical (QM/MM) approach, thereby enabling further consideration of interactions within the entire system for superior results. We aim to accurately model the chemistry of metal oxides using this combined scheme. Here, using the prototypical MgO(100) surface as a test system, with Mg9O14 as the cluster in the quantum mechanical region, we show that the embedded cluster with untailored boundary effective core potentials (ECPs) can have frontier orbital energy levels that substantially deviate from the quantum mechanical reference results. This occurs even when Mg9O9, which retains the stoichiometry of MgO, is used as the cluster in the quantum mechanical region. As a result, the chemical properties of the embedded cluster models differ from those of the quantum mechanical reference model. To address this issue, we propose a new variant of the embedded cluster method called the level-shifted embedded cluster (LSEC) method, which allows the energy levels to be shifted to match the reference levels by tuning the boundary ECPs. Our validation calculations on the adsorption of various adsorbates with different properties on the MgO(100) surface show that the overall performance of QM/MM with the LSEC method is excellent for the adsorption energies, geometries, and charge properties. The excellent performance holds for both the nonstoichiometric and stoichiometric clusters (i.e., Mg9O14 and Mg9O9, respectively), demonstrating the robustness of the LSEC method. We expect that the LSEC method can be combined with QM/MM or used separately for future chemical studies of metal oxides and other ionically bonded systems.
Collapse
Affiliation(s)
- Ming-Yu Yang
- State Key Laboratory of Green Chemical Engineering and Industrial Catalysis, Centre for Computational Chemistry and Research Institute of Industrial Catalysis, School of Chemistry and Molecular Engineering, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, P.R. China
| | - Xin-Ping Wu
- State Key Laboratory of Green Chemical Engineering and Industrial Catalysis, Centre for Computational Chemistry and Research Institute of Industrial Catalysis, School of Chemistry and Molecular Engineering, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, P.R. China
| |
Collapse
|
4
|
Shi B, Zen A, Kapil V, Nagy PR, Grüneis A, Michaelides A. Many-Body Methods for Surface Chemistry Come of Age: Achieving Consensus with Experiments. J Am Chem Soc 2023; 145:25372-25381. [PMID: 37948071 PMCID: PMC10683001 DOI: 10.1021/jacs.3c09616] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2023] [Revised: 10/15/2023] [Accepted: 10/17/2023] [Indexed: 11/12/2023]
Abstract
The adsorption energy of a molecule onto the surface of a material underpins a wide array of applications, spanning heterogeneous catalysis, gas storage, and many more. It is the key quantity where experimental measurements and theoretical calculations meet, with agreement being necessary for reliable predictions of chemical reaction rates and mechanisms. The prototypical molecule-surface system is CO adsorbed on MgO, but despite intense scrutiny from theory and experiment, there is still no consensus on its adsorption energy. In particular, the large cost of accurate many-body methods makes reaching converged theoretical estimates difficult, generating a wide range of values. In this work, we address this challenge, leveraging the latest advances in diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] to obtain accurate predictions for CO on MgO. These reliable theoretical estimates allow us to evaluate the inconsistencies in published temperature-programed desorption experiments, revealing that they arise from variations in employed pre-exponential factors. Utilizing this insight, we derive new experimental estimates of the (electronic) adsorption energy with a (more) precise pre-exponential factor. As a culmination of all of this effort, we are able to reach a consensus between multiple theoretical calculations and multiple experiments for the first time. In addition, we show that our recently developed cluster-based CCSD(T) approach provides a low-cost route toward achieving accurate adsorption energies. This sets the stage for affordable and reliable theoretical predictions of chemical reactions on surfaces to guide the realization of new catalysts and gas storage materials.
Collapse
Affiliation(s)
- Benjamin
X. Shi
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CB2 1EW Cambridge, U.K.
| | - Andrea Zen
- Dipartimento
di Fisica Ettore Pancini, Università
di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
- Department
of Earth Sciences, University College London, Gower Street, WC1E 6BT London, U.K.
| | - Venkat Kapil
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CB2 1EW Cambridge, U.K.
| | - Péter R. Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Müegyetem rkp. 3, H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111 Budapest, Hungary
| | - Andreas Grüneis
- Institute
for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, 1040 Vienna, Austria
| | - Angelos Michaelides
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CB2 1EW Cambridge, U.K.
| |
Collapse
|
5
|
Di Liberto G, Pacchioni G. Modeling Single-Atom Catalysis. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2023; 35:e2307150. [PMID: 37749881 DOI: 10.1002/adma.202307150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Revised: 08/17/2023] [Indexed: 09/27/2023]
Abstract
Electronic structure calculations represent an essential complement of experiments to characterize single-atom catalysts (SACs), consisting of isolated metal atoms stabilized on a support, but also to predict new catalysts. However, simulating SACs with quantum chemistry approaches is not as simple as often assumed. In this work, the essential factors that characterize a reliable simulation of SACs activity are examined. The Perspective focuses on the importance of precise atomistic characterization of the active site, since even small changes in the metal atom's surroundings can result in large changes in reactivity. The dynamical behavior and stability of SACs under working conditions, as well as the importance of adopting appropriate methods to solve the Schrödinger equation for a quantitative evaluation of reaction energies are addressed. The Perspective also focuses on the relevance of the model adopted. For electrocatalysis this must include the effects of the solvent, the presence of electrolytes, the pH, and the external potential. Finally, it is discussed how the similarities between SACs and coordination compounds may result in reaction intermediates that usually are not observed on metal electrodes. When these aspects are not adequately considered, the predictive power of electronic structure calculations is quite limited.
Collapse
Affiliation(s)
- Giovanni Di Liberto
- Dipartimento di Scienza dei Materiali, Università degli studi di Milano Bicocca, Via R. Cozzi 55, Milano, 20125, Italy
| | - Gianfranco Pacchioni
- Dipartimento di Scienza dei Materiali, Università degli studi di Milano Bicocca, Via R. Cozzi 55, Milano, 20125, Italy
| |
Collapse
|
6
|
Verma S, Mitra A, Jin Y, Haldar S, Vorwerk C, Hermes MR, Galli G, Gagliardi L. Optical Properties of Neutral F Centers in Bulk MgO with Density Matrix Embedding. J Phys Chem Lett 2023; 14:7703-7710. [PMID: 37606586 DOI: 10.1021/acs.jpclett.3c01875] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
The optical spectra of neutral oxygen vacancies (F0 centers) in the bulk MgO lattice are investigated using density matrix embedding theory. The impurity Hamiltonian is solved with the complete active space self-consistent field and second-order n-electron valence state perturbation theory (NEVPT2-DMET) multireference methods. To estimate defect-localized vertical excitation energies at the nonembedding and thermodynamic limits, a double extrapolation scheme is employed. The extrapolated NEVPT2-DMET vertical excitation energy value of 5.24 eV agrees well with the experimental absorption maxima at 5.03 eV, whereas the excitation energy value of 2.89 eV at the relaxed triplet defect-localized state geometry overestimates the experimental emission at 2.4 eV by only nearly 0.5 eV, indicating the involvement of the triplet-singlet decay pathway.
Collapse
Affiliation(s)
- Shreya Verma
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Abhishek Mitra
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Yu Jin
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Soumi Haldar
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Christian Vorwerk
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew R Hermes
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Giulia Galli
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Materials Science Division and Center for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Laura Gagliardi
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Argonne National Laboratory, 9700 S. Cass Avenue, Lemont, Illinois 60439, United States
| |
Collapse
|
7
|
Nusspickel M, Ibrahim B, Booth GH. Effective Reconstruction of Expectation Values from Ab Initio Quantum Embedding. J Chem Theory Comput 2023; 19:2769-2791. [PMID: 37155201 DOI: 10.1021/acs.jctc.2c01063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
Quantum embedding is an appealing route to fragment a large interacting quantum system into several smaller auxiliary "cluster" problems to exploit the locality of the correlated physics. In this work, we critically review approaches to recombine these fragmented solutions in order to compute nonlocal expectation values, including the total energy. Starting from the democratic partitioning of expectation values used in density matrix embedding theory, we motivate and develop a number of alternative approaches, numerically demonstrating their efficiency and improved accuracy as a function of increasing cluster size for both energetics and nonlocal two-body observables in molecular and solid state systems. These approaches consider the N-representability of the resulting expectation values via an implicit global wave function across the clusters, as well as the importance of including contributions to expectation values spanning multiple fragments simultaneously, thereby alleviating the fundamental locality approximation of the embedding. We clearly demonstrate the value of these introduced functionals for reliable extraction of observables and robust and systematic convergence as the cluster size increases, allowing for significantly smaller clusters to be used for a desired accuracy compared to traditional approaches in ab initio wave function quantum embedding.
Collapse
Affiliation(s)
- Max Nusspickel
- Department of Physics, King's College London, Strand, London WC2R 2LS, United Kingdom
| | - Basil Ibrahim
- Department of Physics, King's College London, Strand, London WC2R 2LS, United Kingdom
| | - George H Booth
- Department of Physics, King's College London, Strand, London WC2R 2LS, United Kingdom
| |
Collapse
|
8
|
Baranov O, Bazaka K, Belmonte T, Riccardi C, Roman HE, Mohandas M, Xu S, Cvelbar U, Levchenko I. Recent innovations in the technology and applications of low-dimensional CuO nanostructures for sensing, energy and catalysis. NANOSCALE HORIZONS 2023; 8:568-602. [PMID: 36928662 DOI: 10.1039/d2nh00546h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
Low-dimensional copper oxide nanostructures are very promising building blocks for various functional materials targeting high-demanded applications, including energy harvesting and transformation systems, sensing and catalysis. Featuring a very high surface-to-volume ratio and high chemical reactivity, these materials have attracted wide interest from researchers. Currently, extensive research on the fabrication and applications of copper oxide nanostructures ensures the fast progression of this technology. In this article we briefly outline some of the most recent, mostly within the past two years, innovations in well-established fabrication technologies, including oxygen plasma-based methods, self-assembly and electric-field assisted growth, electrospinning and thermal oxidation approaches. Recent progress in several key types of leading-edge applications of CuO nanostructures, mostly for energy, sensing and catalysis, is also reviewed. Besides, we briefly outline and stress novel insights into the effect of various process parameters on the growth of low-dimensional copper oxide nanostructures, such as the heating rate, oxygen flow, and roughness of the substrates. These insights play a key role in establishing links between the structure, properties and performance of the nanomaterials, as well as finding the cost-and-benefit balance for techniques that are capable of fabricating low-dimensional CuO with the desired properties and facilitating their integration into more intricate material architectures and devices without the loss of original properties and function.
Collapse
Affiliation(s)
- Oleg Baranov
- Department of Theoretical Mechanics, Engineering and Robomechanical Systems, National Aerospace University, Kharkiv 61070, Ukraine.
- Department of Gaseous Electronics, Jozef Stefan Institute, Ljubljana 1000, Slovenia
| | - Kateryna Bazaka
- School of Engineering, The Australian National University, Canberra, ACT 2601, Australia
| | | | - Claudia Riccardi
- Dipartimento di Fisica "Giuseppe Occhialini", Università degli Studi di Milano-Bicocca, Piazza della Scienza 3, I20126 Milan, Italy
| | - H Eduardo Roman
- Dipartimento di Fisica "Giuseppe Occhialini", Università degli Studi di Milano-Bicocca, Piazza della Scienza 3, I20126 Milan, Italy
| | - Mandhakini Mohandas
- Center for Nanoscience and Technology, Anna University, Chennai, 600 025, India
| | - Shuyan Xu
- Plasma Sources and Application Centre, NIE, Nanyang Technological University, 637616, Singapore.
| | - Uroš Cvelbar
- Department of Gaseous Electronics, Jozef Stefan Institute, Ljubljana 1000, Slovenia
| | - Igor Levchenko
- Plasma Sources and Application Centre, NIE, Nanyang Technological University, 637616, Singapore.
| |
Collapse
|
9
|
Haldar S, Mitra A, Hermes MR, Gagliardi L. Local Excitations of a Charged Nitrogen Vacancy in Diamond with Multireference Density Matrix Embedding Theory. J Phys Chem Lett 2023; 14:4273-4280. [PMID: 37126760 DOI: 10.1021/acs.jpclett.3c00551] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
We investigate the negatively charged nitrogen-vacancy center in diamond using periodic density matrix embedding theory (pDMET). To describe the strongly correlated excited states of this system, the complete active space self-consistent field (CASSCF) followed by n-electron valence state second-order perturbation theory (NEVPT2) was used as the impurity solver. Since the NEVPT2-DMET energies show a linear dependence on the inverse of the size of the embedding subspace, we performed an extrapolation of the excitation energies to the nonembedding limit using a linear regression. The extrapolated NEVPT2-DMET first triplet-triplet excitation energy is 2.31 eV and that for the optically inactive singlet-singlet transition is 1.02 eV, both in agreement with the experimentally observed vertical excitation energies of ∼2.18 eV and ∼1.26 eV, respectively. This is the first application of pDMET to a charged periodic system and the first investigation of the NV- defect using NEVPT2 for periodic supercell models.
Collapse
Affiliation(s)
- Soumi Haldar
- Department of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Abhishek Mitra
- Department of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew R Hermes
- Department of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
- Argonne National Laboratory, 9700 S. Cass Avenue, Lemont, Illinois 60439, United States
| |
Collapse
|
10
|
Ai Y, Sun Q, Jiang H. Efficient Multiconfigurational Quantum Chemistry Approach to Single-Ion Magnets Based on Density Matrix Embedding Theory. J Phys Chem Lett 2022; 13:10627-10634. [PMID: 36350106 DOI: 10.1021/acs.jpclett.2c02890] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
Density matrix embedding theory (DMET) provides a systematic framework to combine low-level (e.g., Hartree-Fock approximation) and high-level correlated quantum chemistry methods to treat strongly correlated systems. In this work, we propose an efficient quantum embedding approach that combines DMET with the complete active space self-consistent field and subsequent state interaction treatment of spin-orbit coupling (CASSI-SO) and apply it to a theoretical description of single-ion magnets (SIMs). We have developed a novel regularized direct inversion of iterative subspace (R-DIIS) technique that ensures restricted open-shell Hartree-Fock converging to a physically correct ground state, which is found to be crucial for the efficacy of subsequent CASSI-SO calculation. The DMET+CASSI-SO approach can produce reliable zero-field splitting parameters in typical 3d-SIMs with dramatically reduced computational cost compared to its all-electron counterpart. This work therefore demonstrates the great potential of DMET-based multiconfigurational approaches for efficient ab initio study of magneto-structural correlations in complex molecular magnetic systems.
Collapse
Affiliation(s)
- Yuhang Ai
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Qiming Sun
- Axiomquant Investment Management LLC, Beijing 100871, China
| | - Hong Jiang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|