1
|
Cheong O, Bornhake T, Zhu X, Eikerling MH. Stay Hydrated! Impact of Solvation Phenomena on the CO 2 Reduction Reaction at Pb(100) and Ag(100) surfaces. CHEMSUSCHEM 2023; 16:e202300885. [PMID: 37539768 DOI: 10.1002/cssc.202300885] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Revised: 08/03/2023] [Accepted: 08/04/2023] [Indexed: 08/05/2023]
Abstract
Herein, a comprehensive computational study of the impact of solvation on the reduction reaction of CO2 to formic acid (HCOOH) and carbon monoxide on Pb(100) and Ag(100) surfaces is presented. Results further the understanding of how solvation phenomena influence the adsorption energies of reaction intermediates. We applied an explicit solvation scheme harnessing a combined density functional theory (DFT)/microkinetic modeling approach for the CO2 reduction reaction. This approach reveals high selectivities for CO formation at Ag and HCOOH formation on Pb, resolving the prior disparity between ab initio calculations and experimental observations. Furthermore, the detailed analysis of adsorption energies of relevant reaction intermediates shows that the total number of hydrogen bonds formed by HCOO plays a primary role for the adsorption strength of intermediates and the electrocatalytic activity. Results emphasize the importance of explicit solvation for adsorption and electrochemical reaction phenomena on metal surfaces.
Collapse
Affiliation(s)
- Oskar Cheong
- Institute of Energy and Climate Research (IEK-13), Forschungszentrum Jülich, Wilhelm-Johnen-Straße, 52425, Jülich, Germany
- Chair of Theory and Computation of Energy Materials, Faculty of Georesources and Materials Enginering, RWTH Aachen University, Intzestrasse 5, 52072, Aachen, Germany
- JARA Energy & Center for Simulation and Data Science (CSD), 52425, Jülich, Germany
| | - Thomas Bornhake
- Institute of Energy and Climate Research (IEK-13), Forschungszentrum Jülich, Wilhelm-Johnen-Straße, 52425, Jülich, Germany
- JARA Energy & Center for Simulation and Data Science (CSD), 52425, Jülich, Germany
| | - Xinwei Zhu
- Institute of Energy and Climate Research (IEK-13), Forschungszentrum Jülich, Wilhelm-Johnen-Straße, 52425, Jülich, Germany
- Chair of Theory and Computation of Energy Materials, Faculty of Georesources and Materials Enginering, RWTH Aachen University, Intzestrasse 5, 52072, Aachen, Germany
- JARA Energy & Center for Simulation and Data Science (CSD), 52425, Jülich, Germany
| | - Michael H Eikerling
- Institute of Energy and Climate Research (IEK-13), Forschungszentrum Jülich, Wilhelm-Johnen-Straße, 52425, Jülich, Germany
- Chair of Theory and Computation of Energy Materials, Faculty of Georesources and Materials Enginering, RWTH Aachen University, Intzestrasse 5, 52072, Aachen, Germany
- JARA Energy & Center for Simulation and Data Science (CSD), 52425, Jülich, Germany
| |
Collapse
|
2
|
Chen BWJ, Zhang X, Zhang J. Accelerating explicit solvent models of heterogeneous catalysts with machine learning interatomic potentials. Chem Sci 2023; 14:8338-8354. [PMID: 37564405 PMCID: PMC10411631 DOI: 10.1039/d3sc02482b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2023] [Accepted: 07/11/2023] [Indexed: 08/12/2023] Open
Abstract
Realistically modelling how solvents affect catalytic reactions is a longstanding challenge due to its prohibitive computational cost. Typically, an explicit atomistic treatment of the solvent molecules is needed together with molecular dynamics (MD) simulations and enhanced sampling methods. Here, we demonstrate the utility of machine learning interatomic potentials (MLIPs), coupled with active learning, to enable fast and accurate explicit solvent modelling of adsorption and reactions on heterogeneous catalysts. MLIPs trained on-the-fly were able to accelerate ab initio MD simulations by up to 4 orders of magnitude while reproducing with high fidelity the geometrical features of water in the bulk and at metal-water interfaces. Using these ML-accelerated simulations, we accurately predicted key catalytic quantities such as the adsorption energies of CO*, OH*, COH*, HCO*, and OCCHO* on Cu surfaces and the free energy barriers of C-H scission of ethylene glycol over Cu and Pd surfaces, as validated with ab initio calculations. We envision that such simulations will pave the way towards detailed and realistic studies of solvated catalysts at large time- and length-scales.
Collapse
Affiliation(s)
- Benjamin W J Chen
- Institute of High Performance Computing (IHPC), Agency for Science, Technology and Research (A*STAR) 1 Fusionopolis Way, #16-16 Connexis Singapore 138632 Singapore
| | - Xinglong Zhang
- Institute of High Performance Computing (IHPC), Agency for Science, Technology and Research (A*STAR) 1 Fusionopolis Way, #16-16 Connexis Singapore 138632 Singapore
| | - Jia Zhang
- Institute of High Performance Computing (IHPC), Agency for Science, Technology and Research (A*STAR) 1 Fusionopolis Way, #16-16 Connexis Singapore 138632 Singapore
| |
Collapse
|
3
|
Bramley GA, Beynon OT, Stishenko PV, Logsdail AJ. The application of QM/MM simulations in heterogeneous catalysis. Phys Chem Chem Phys 2023; 25:6562-6585. [PMID: 36810655 DOI: 10.1039/d2cp04537k] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
The QM/MM simulation method is provenly efficient for the simulation of biological systems, where an interplay of extensive environment and delicate local interactions drives a process of interest through a funnel on a complex energy landscape. Recent advances in quantum chemistry and force-field methods present opportunities for the adoption of QM/MM to simulate heterogeneous catalytic processes, and their related systems, where similar intricacies exist on the energy landscape. Herein, the fundamental theoretical considerations for performing QM/MM simulations, and the practical considerations for setting up QM/MM simulations of catalytic systems, are introduced; then, areas of heterogeneous catalysis are explored where QM/MM methods have been most fruitfully applied. The discussion includes simulations performed for adsorption processes in solvent at metallic interfaces, reaction mechanisms within zeolitic systems, nanoparticles, and defect chemistry within ionic solids. We conclude with a perspective on the current state of the field and areas where future opportunities for development and application exist.
Collapse
Affiliation(s)
- Gabriel Adrian Bramley
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, CF10 3AT, UK.
| | - Owain Tomos Beynon
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, CF10 3AT, UK.
| | | | - Andrew James Logsdail
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, CF10 3AT, UK.
| |
Collapse
|
4
|
Experimental and theoretical study unveiling the role of solvents on CO activation and hydrogenation to methanol in three-phase reactor system. J Mol Struct 2023. [DOI: 10.1016/j.molstruc.2022.134392] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
5
|
Abstract
Adsorption energy (AE) of reactive intermediate is currently the most important descriptor for electrochemical reactions (e.g., water electrolysis, hydrogen fuel cell, electrochemical nitrogen fixation, electrochemical carbon dioxide reduction, etc.), which can bridge the gap between catalyst's structure and activity. Tracing the history and evolution of AE can help to understand electrocatalysis and design optimal electrocatalysts. Focusing on oxygen electrocatalysis, this review aims to provide a comprehensive introduction on how AE is selected as the activity descriptor, the intrinsic and empirical relationships related to AE, how AE links the structure and electrocatalytic performance, the approaches to obtain AE, the strategies to improve catalytic activity by modulating AE, the extrinsic influences on AE from the environment, and the methods in circumventing linear scaling relations of AE. An outlook is provided at the end with emphasis on possible future investigation related to the obstacles existing between adsorption energy and electrocatalytic performance.
Collapse
Affiliation(s)
- Junming Zhang
- School of Chemical and Biomedical Engineering, Nanyang Technological University, Singapore 637459, Singapore
| | - Hong Bin Yang
- School of Chemical and Biomedical Engineering, Nanyang Technological University, Singapore 637459, Singapore
| | - Daojin Zhou
- State Key Laboratory of Chemical Resource Engineering, Beijing Advanced Innovation Center for Soft Matter Science and Engineering, College of Chemistry, Beijing University of Chemical Technology, Beijing 100029, P. R. China.,Department of Electrical and Computer Engineering, University of Toronto, 35 St. George Street, Toronto, Ontario M5S 1A4, Canada
| | - Bin Liu
- School of Chemical and Biomedical Engineering, Nanyang Technological University, Singapore 637459, Singapore
| |
Collapse
|
6
|
Zare M, Saleheen MS, Singh N, Uline MJ, Faheem M, Heyden A. Liquid-Phase Effects on Adsorption Processes in Heterogeneous Catalysis. JACS AU 2022; 2:2119-2134. [PMID: 36186571 PMCID: PMC9516566 DOI: 10.1021/jacsau.2c00389] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/08/2022] [Revised: 07/21/2022] [Accepted: 07/26/2022] [Indexed: 06/16/2023]
Abstract
Aqueous solvation free energies of adsorption have recently been measured for phenol adsorption on Pt(111). Endergonic solvent effects of ∼1 eV suggest solvents dramatically influence a metal catalyst's activity with significant implications for the catalyst design. However, measurements are indirect and involve adsorption isotherm models, which potentially reduces the reliability of the extracted energy values. Computational, implicit solvation models predict exergonic solvation effects for phenol adsorption, failing to agree with measurements even qualitatively. In this study, an explicit, hybrid quantum mechanical/molecular mechanical approach for computing solvation free energies of adsorption is developed, solvation free energies of phenol adsorption are computed, and experimental data for solvation free energies of phenol adsorption are reanalyzed using multiple adsorption isotherm models. Explicit solvation calculations predict an endergonic solvation free energy for phenol adsorption that agrees well with measurements to within the experimental and force field uncertainties. Computed adsorption free energies of solvation of carbon monoxide, ethylene glycol, benzene, and phenol over the (111) facet of Pt and Cu suggest that liquid water destabilizes all adsorbed species, with the largest impact on the largest adsorbates.
Collapse
Affiliation(s)
- Mehdi Zare
- Department
of Chemical Engineering, University of South
Carolina, 301 Main Street, Columbia, South Carolina 29208, United States
| | - Mohammad S. Saleheen
- Department
of Chemical Engineering, University of South
Carolina, 301 Main Street, Columbia, South Carolina 29208, United States
| | - Nirala Singh
- Department
of Chemical Engineering and Catalysis Science and Technology Institute, University of Michigan, Ann Arbor, Michigan 48109-2136, United States
| | - Mark J. Uline
- Department
of Chemical Engineering, University of South
Carolina, 301 Main Street, Columbia, South Carolina 29208, United States
| | - Muhammad Faheem
- Department
of Chemical Engineering, University of South
Carolina, 301 Main Street, Columbia, South Carolina 29208, United States
- Department
of Chemical Engineering, University of Engineering
& Technology, Lahore 54890, Pakistan
| | - Andreas Heyden
- Department
of Chemical Engineering, University of South
Carolina, 301 Main Street, Columbia, South Carolina 29208, United States
| |
Collapse
|
7
|
Muramoto E, Patel DA, Chen W, Sautet P, Sykes ECH, Madix RJ. Direct Observation of Solvent–Reaction Intermediate Interactions in Heterogeneously Catalyzed Alcohol Coupling. J Am Chem Soc 2022; 144:17387-17398. [DOI: 10.1021/jacs.2c02199] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Eri Muramoto
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Dipna A. Patel
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, United States
| | - Wei Chen
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
- Center for Functional Nanomaterials, Brookhaven National Laboratory, Upton, New York 11973, United States
| | - Philippe Sautet
- Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, California 90095, United States
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
| | - E. Charles H. Sykes
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, United States
| | - Robert J. Madix
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| |
Collapse
|
8
|
Geng Y, Li H. Hydrogen Spillover-Enhanced Heterogeneously Catalyzed Hydrodeoxygenation for Biomass Upgrading. CHEMSUSCHEM 2022; 15:e202102495. [PMID: 35230748 DOI: 10.1002/cssc.202102495] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2021] [Revised: 02/25/2022] [Indexed: 06/14/2023]
Abstract
Hydrodeoxygenation (HDO) is regarded as a promising technology for biomass upgrading to obtain sustainable and competitive chemicals and fuels. In fact, biomass HDO over heterogeneous solid catalysts is often accompanied by the phenomenon of hydrogen spillover, which further affects the catalytic performance. Thus, it is necessary to gain in-depth understand the promoting effect of hydrogen spillover in the biomass HDO process to obtain desired conversion and selectivity. This Review summarized the extensive research on hydrogen spillover in biomass refining and discussed in detail the regulation mechanism of hydrogen spillover in biomass HDO process, mainly by regulating different active center sites on catalyst supports, such as metal sites, acid sites, surface functional groups, and defective sites, which exhibit independent and synergistic characteristics promoting catalyst activity, selectivity, and stability. Finally, the prospective of hydrogen spillover in biomass HDO applications was critically evaluated, and the key technical challenges in developing "hydrogen-free" HDO and upgrading biofuels were highlighted. The presentation of hydrogen spillover-enhanced catalytic biomass HDO in this Review will hopefully provide insight and guidance for further development of efficient catalysts and preparation of high-value chemicals in the future.
Collapse
Affiliation(s)
- Yanyan Geng
- National-Local Joint Engineering Laboratory for Energy Conservation in Chemical Process Integration and Resources Utilization, Tianjin Key Laboratory of Chemical Process Safety, School of Chemical Engineering and Technology, Hebei University of Technology, 8 Guangrong Road, Tianjin, 300130, P. R. China
| | - Hao Li
- National-Local Joint Engineering Laboratory for Energy Conservation in Chemical Process Integration and Resources Utilization, Tianjin Key Laboratory of Chemical Process Safety, School of Chemical Engineering and Technology, Hebei University of Technology, 8 Guangrong Road, Tianjin, 300130, P. R. China
| |
Collapse
|
9
|
Korpelin V, Kiljunen T, Melander MM, Caro MA, Kristoffersen HH, Mammen N, Apaja V, Honkala K. Addressing Dynamics at Catalytic Heterogeneous Interfaces with DFT-MD: Anomalous Temperature Distributions from Commonly Used Thermostats. J Phys Chem Lett 2022; 13:2644-2652. [PMID: 35297635 PMCID: PMC8959310 DOI: 10.1021/acs.jpclett.2c00230] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2022] [Accepted: 03/11/2022] [Indexed: 05/28/2023]
Abstract
Density functional theory-based molecular dynamics (DFT-MD) has been widely used for studying the chemistry of heterogeneous interfacial systems under operational conditions. We report frequently overlooked errors in thermostated or constant-temperature DFT-MD simulations applied to study (electro)catalytic chemistry. Our results demonstrate that commonly used thermostats such as Nosé-Hoover, Berendsen, and simple velocity-rescaling methods fail to provide a reliable temperature description for systems considered. Instead, nonconstant temperatures and large temperature gradients within the different parts of the system are observed. The errors are not a "feature" of any particular code but are present in several ab initio molecular dynamics implementations. This uneven temperature distribution, due to inadequate thermostatting, is well-known in the classical MD community, where it is ascribed to the failure in kinetic energy equipartition among different degrees of freedom in heterogeneous systems (Harvey et al. J. Comput. Chem. 1998, 726-740) and termed the flying ice cube effect. We provide tantamount evidence that interfacial systems are susceptible to substantial flying ice cube effects and demonstrate that the traditional Nosé-Hoover and Berendsen thermostats should be applied with care when simulating, for example, catalytic properties or structures of solvated interfaces and supported clusters. We conclude that the flying ice cube effect in these systems can be conveniently avoided using Langevin dynamics.
Collapse
Affiliation(s)
- Ville Korpelin
- Department
of Chemistry, Nanoscience Center, University
of Jyväskylä, P.O. Box 35 (YN), FI-40014 Jyväskylä, Finland
| | - Toni Kiljunen
- Department
of Chemistry, Nanoscience Center, University
of Jyväskylä, P.O. Box 35 (YN), FI-40014 Jyväskylä, Finland
| | - Marko M. Melander
- Department
of Chemistry, Nanoscience Center, University
of Jyväskylä, P.O. Box 35 (YN), FI-40014 Jyväskylä, Finland
| | - Miguel A. Caro
- Department
of Electrical Engineering and Automation, Aalto University, FIN-02150 Espoo, Finland
| | | | - Nisha Mammen
- Department
of Physics,Nanoscience Center, University
of Jyväskylä, P.O. Box
35 (YN), FI-40014 Jyväskylä, Finland
| | - Vesa Apaja
- Department
of Physics,Nanoscience Center, University
of Jyväskylä, P.O. Box
35 (YN), FI-40014 Jyväskylä, Finland
| | - Karoliina Honkala
- Department
of Chemistry, Nanoscience Center, University
of Jyväskylä, P.O. Box 35 (YN), FI-40014 Jyväskylä, Finland
| |
Collapse
|