1
|
Sokolov M, Doblhoff-Dier K, Exner KS. Best practices of modeling complex materials in electrocatalysis, exemplified by oxygen evolution reaction on pentlandites. Phys Chem Chem Phys 2024; 26:22359-22370. [PMID: 39158931 DOI: 10.1039/d4cp01792g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/20/2024]
Abstract
Pentlandites are natural ores with structural properties comparable to that of [FeNi] hydrogenases. While this class of transition-metal sulfide materials - (Fe,Ni)9S8 - with a variable Fe : Ni ratio has been proven to be an active electrode material for the hydrogen evolution reaction, it is also discussed as electrocatalyst for the alkaline oxygen evolution reaction (OER), corresponding to the bottleneck of anion exchange membrane electrolyzers for green hydrogen production. Despite the experimental evidence for the use of (Fe,Ni)9S8 as an OER catalyst, a detailed investigation of the elementary reaction steps, including consideration of adsorbate coverages and limiting steps under anodic polarizing conditions, is still missing. We address this gap in the present manuscript by gaining atomistic insights into the OER on an Fe4.5Ni4.5S8(111) surface through density functional theory calculations combined with a descriptor-based analysis. We use this system to introduce best practices for modeling this rather complex material by pointing out hidden pitfalls that can arise when using the popular computational hydrogen electrode approach to describe electrocatalytic processes at the electrified solid/liquid interface for energy conversion and storage.
Collapse
Affiliation(s)
- Maksim Sokolov
- Faculty of Chemistry, Theoretical Inorganic Chemistry, University Duisburg-Essen, Universitätsstraße 5, 45141 Essen, Germany.
- Cluster of Excellence RESOLV, 44801 Bochum, Germany
| | - Katharina Doblhoff-Dier
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, Leiden 2300, RA, The Netherlands
| | - Kai S Exner
- Faculty of Chemistry, Theoretical Inorganic Chemistry, University Duisburg-Essen, Universitätsstraße 5, 45141 Essen, Germany.
- Cluster of Excellence RESOLV, 44801 Bochum, Germany
- Center for Nanointegration (CENIDE) Duisburg-Essen, 47057 Duisburg, Germany
| |
Collapse
|
2
|
Illner H, Sakong S, Groß A, Wintterlin J. Solution of the structure of the high-coverage CO layer on the Ru(0001) surface-A combined study by density functional theory and scanning tunneling microscopy. J Chem Phys 2024; 161:014703. [PMID: 38949590 DOI: 10.1063/5.0215872] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2024] [Accepted: 06/12/2024] [Indexed: 07/02/2024] Open
Abstract
Structures formed by dense CO adsorption layers can provide information about the balance between molecule-surface and molecule-molecule interactions. However, in many cases, the structure models are not clear. Using density functional theory (DFT) and scanning tunneling microscopy (STM), we have investigated the high-coverage CO layer on the Ru(0001) surface. Previous investigations by low-energy electron diffraction (LEED) and vibrational spectroscopy led to conflicting results about the structure. In the present study, 88 models with coverages between 0.58 and 0.77 monolayers have been analyzed by DFT. The most stable structures consist of small, compact CO clusters with an internal pseudo 1×1 structure. The CO molecules in the cluster centers occupy on-top sites in an upright position, whereas the molecules farther outside are slightly shifted from these sites and tilted outward. STM data of the CO-saturated surface at low temperatures, corresponding to a coverage of 0.66 monolayers, show a quasi-hexagonal pattern of features with an internal hexagonal fine structure. Simulated images based on the cluster model agree with the experimental data. It is concluded that the high-coverage CO layer consists of the close-packed clusters predicted by DFT as the most stable structure elements. In the experiment, the sizes and shapes of the clusters vary. However, the arrangement is not random but follows defined tiling rules. The structure remains ordered, almost up to room temperature. The LEED data are re-interpreted on the basis of the Fourier transforms of the STM data, solving the long-standing conflict about the structure.
Collapse
Affiliation(s)
- Hannah Illner
- Department of Chemistry, Ludwig-Maximilians-Universität München, 81377 Munich, Germany
| | - Sung Sakong
- Institute of Theoretical Chemistry, University of Ulm, 89081 Ulm, Germany
| | - Axel Groß
- Institute of Theoretical Chemistry, University of Ulm, 89081 Ulm, Germany
| | - Joost Wintterlin
- Department of Chemistry, Ludwig-Maximilians-Universität München, 81377 Munich, Germany
- Center for NanoScience, Schellingstr. 4, 80799 Munich, Germany
| |
Collapse
|
3
|
Mancera LA, Groß A, Behm RJ. Stability, electronic properties and CO adsorption properties of bimetallic PtAg/Pt(111) surfaces. Phys Chem Chem Phys 2024; 26:18435-18448. [PMID: 38916054 DOI: 10.1039/d4cp01640h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/26/2024]
Abstract
Aiming at a better fundamental understanding of the chemistry of bimetallic PtAg/Pt(111) surfaces, we have investigated the stability, electronic properties and CO adsorption properties of bimetallic PtAg surfaces, including pseudomorphic Ag film covered Pt(111) surfaces and PtxAg1-x/Pt(111) monolayer surface alloys, using periodic density functional theory calculations. The data provide detailed insights into the relative stabilities of different surface configurations, as indicated by their formation enthalpies and surface energies, and changes in their electronic properties, i.e., in the projected local densities of states and shifts in the d-band center. The adsorption properties of different Ptn ensembles were systematically tested using CO as a probe molecule. In addition to electronic ligand and strain effects, we were particularly interested in the role of different adsorption sites and of the local COad coverage, given by the number of CO molecules per Pt surface atom in the Ptn ensemble. Different from PdAg surfaces, variations in the adsorption energy with adsorption sites and with increasing local coverage are small up to one COad per Pt surface atom. Finally, formation of multicarbonyl species with more than one COad per Pt surface atom was tested for separated Pt1 monomers and can be excluded at finite temperatures. General trends and aspects are derived by comparison with comparable data for PdAg bimetallic surfaces. Fundamental insights relevant for applications of bimetallic Pt catalysts, specifically PtAg catalysts, are briefly discussed.
Collapse
Affiliation(s)
- Luis A Mancera
- Institute of Theoretical Chemistry, Ulm University, Oberberghof 7, D-89081 Ulm, Germany.
| | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, Oberberghof 7, D-89081 Ulm, Germany.
| | - R Jürgen Behm
- Institute of Theoretical Chemistry, Ulm University, Oberberghof 7, D-89081 Ulm, Germany.
| |
Collapse
|
4
|
Stottmeister D, Wildersinn L, Maibach J, Hofmann A, Jeschull F, Groß A. Unraveling Propylene Oxide Formation in Alkali Metal Batteries. CHEMSUSCHEM 2024; 17:e202300995. [PMID: 37820026 DOI: 10.1002/cssc.202300995] [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/10/2023] [Revised: 10/11/2023] [Accepted: 10/11/2023] [Indexed: 10/13/2023]
Abstract
The increasing need for electrochemical energy storage drives the development of post-lithium battery systems. Among the most promising new battery types are sodium-based battery systems. However, like its lithium predecessor, sodium batteries suffer from various issues like parasitic side reactions, which lead to a loss of active sodium inventory, thus reducing the capacity over time. Some problems in sodium batteries arise from an unstable solid electrolyte interphase (SEI) reducing its protective power e. g., due to increased solubility of SEI components in sodium battery systems. While it is known that the electrolyte affects the SEI structure, the exact formation mechanism of the SEI is not yet fully understood. In this study, we follow the initial SEI formation on a piece of sodium metal submerged in propylene carbonate with and without the electrolyte salt sodium perchlorate. We combine X-ray photoelectron spectroscopy, gas chromatography, and density functional theory to unravel the sudden emergence of propylene oxide after adding sodium perchlorate to the electrolyte solvent. We identify the formation of a sodium chloride layer as a crucial step in forming propylene oxide by enabling precursors formed from propylene carbonate on the sodium metal surface to undergo a ring-closing reaction. Based on our combined theoretical and experimental approach, we identify changes in the electrolyte decomposition process, propose a reaction mechanism to form propylene oxide and discuss alternatives based on known synthesis routes.
Collapse
Affiliation(s)
| | - Leonie Wildersinn
- Karlsruher Institut für Technologie, Institut für Angewandte Materialien (IAM), Herrmann-von-Helmholtz Platz 1, 76344, Eggenstein - Leopoldshafen, Germany
| | - Julia Maibach
- Karlsruher Institut für Technologie, Institut für Angewandte Materialien (IAM), Herrmann-von-Helmholtz Platz 1, 76344, Eggenstein - Leopoldshafen, Germany
- Department of Physics, Chalmers University of Technology, SE - 412 96, Gothenburg, Sweden
| | - Andreas Hofmann
- Karlsruher Institut für Technologie, Institut für Angewandte Materialien (IAM), Herrmann-von-Helmholtz Platz 1, 76344, Eggenstein - Leopoldshafen, Germany
| | - Fabian Jeschull
- Karlsruher Institut für Technologie, Institut für Angewandte Materialien (IAM), Herrmann-von-Helmholtz Platz 1, 76344, Eggenstein - Leopoldshafen, Germany
| | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, 89069, Ulm, Germany
- Helmholtz Institute Ulm (HIU) Electrochemical Energy Storage, Helmholtzstr. 11, 89069, Ulm, Germany
| |
Collapse
|
5
|
Putra MH, Bagemihl B, Rau S, Groß A. Prediction of Strong Solvatochromism in a Molecular Photocatalyst. Chemistry 2024; 30:e202302643. [PMID: 37754665 DOI: 10.1002/chem.202302643] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2023] [Revised: 09/25/2023] [Accepted: 09/25/2023] [Indexed: 09/28/2023]
Abstract
Based on quantum chemical calculations, we predict strong solvatochromism in a light-driven molecular photocatalyst for hydrogen generation, that is we show that the electronic and optical properties of the photocatalyst strongly depend on the solvent it is dissolved in. Our calculations in particular indicate a solvent-dependent relocation of the highest occupied molecular orbital (HOMO). Ground-state density functional theory and linear response time-dependent density functional theory calculations were applied in order to investigate the influence of implicit solvents on the structural, electronic and optical properties of a molecular photocatalyst. Only at high dielectric constants of the solvent, is the HOMO located at the metal center of the photosensitizer, whereas at low dielectric constants the HOMO is centered at the metal atom of the catalytically active complex. We elucidate the electronic origins of this strong solvatochromic effect and sketch the consequences of these insights for the use of photocatalysts in different environments.
Collapse
Affiliation(s)
| | - Benedikt Bagemihl
- Institute of Inorganic Chemistry I, Materials and Catalysis, Ulm University, 89069, Ulm, Germany
| | - Sven Rau
- Institute of Inorganic Chemistry I, Materials and Catalysis, Ulm University, 89069, Ulm, Germany
| | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, 89069, Ulm, Germany
- Helmholtz Institute Ulm (HIU), Electrochemical Energy Storage, 89069, Ulm, Germany
| |
Collapse
|
6
|
Ye P, Zhang Y, Tong T, Ao L, Chen Z, Huang H, Hussain A, Ramiere A, Cai X, Liu D, Shen J. 3D Lithiophilic CuZrAg Metallic Glass Based-Current Collector for High-Performance Lithium Metal Anode. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2023; 19:e2304373. [PMID: 37649179 DOI: 10.1002/smll.202304373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/25/2023] [Revised: 08/20/2023] [Indexed: 09/01/2023]
Abstract
Lithium metal anodes face several challenges in practical applications, such as dendrite growth, poor cycle efficiency, and volume variation. 3D hosts with lithiophilic surfaces have emerged as a promising design strategy for anodes. In this study, inspiration from the intrinsic isotropy, chemical heterogeneity, and wide tunability of metallic glass (MG) is drew to develop a 3D mesoporous host with a lithiophilic surface. The CuZrAg MG is prepared using the scalable melt-spinning technique and subsequently treated with a simple one-step chemical dealloying method. This resultes in the creation of a host with a homogeneously distributed abundance of lithium affinity sites on the surface. The excellent lithiophilic property and capability for uniform lithium deposition of the 3D CuZrAg electrode have been confirmed through theoretical calculations. Therefore, the 3D CuZrAg electrode displays excellent cyclic stability for over 400 cycles with 96% coulomb efficiency, and ultra-low overpotentials of 5 mV for over 2000 h at 1.0 mA cm-2 and 1.0 mAh cm-2 . Additionally, the full cells partied with either LiFePO4 or LiNi0.8 Co0.1 Mn0.1 O2 cathode deliver exceptional long-term cyclability and rate capability. This work demonstrates the great potential of metallic glass in lithium metal anode application.
Collapse
Affiliation(s)
- Pengfei Ye
- College of Mechatronics and Control Engineering, Shenzhen University, Shenzhen, 518060, P. R. China
| | - Yanhui Zhang
- State Key Lab of Metastable Materials Science and Technology, and, College of Materials Science and Engineering, Yanshan University, Qinhuangdao, Hebei, 066004, P. R. China
| | - Tong Tong
- College of Mechatronics and Control Engineering, Shenzhen University, Shenzhen, 518060, P. R. China
| | - Lihong Ao
- College of Physics and Optoelectronic Engineering, Shenzhen University, Shenzhen, 518060, P. R. China
| | - Zhe Chen
- College of Mechatronics and Control Engineering, Shenzhen University, Shenzhen, 518060, P. R. China
| | - Huayu Huang
- College of Mechatronics and Control Engineering, Shenzhen University, Shenzhen, 518060, P. R. China
| | - Arshad Hussain
- Institute for Advanced Study, Shenzhen University, Shenzhen, 518060, P. R. China
| | - Aymeric Ramiere
- College of Physics and Optoelectronic Engineering, Shenzhen University, Shenzhen, 518060, P. R. China
| | - Xingke Cai
- Institute for Advanced Study, Shenzhen University, Shenzhen, 518060, P. R. China
| | - Dongqing Liu
- College of Mechatronics and Control Engineering, Shenzhen University, Shenzhen, 518060, P. R. China
| | - Jun Shen
- College of Mechatronics and Control Engineering, Shenzhen University, Shenzhen, 518060, P. R. China
| |
Collapse
|
7
|
Richard A, Corà F. Influence of Dispersion Interactions on the Polymorphic Stability of Crystalline Oxides. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:10766-10776. [PMID: 37313119 PMCID: PMC10259254 DOI: 10.1021/acs.jpcc.3c01013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Revised: 05/11/2023] [Indexed: 06/15/2023]
Abstract
The accurate determination of relative phase stabilities using DFT methods is a significant challenge when some of these can vary by only a few kJ/mol. Here, we demonstrate that for a selection of oxides (TiO2, MnO2, and ZnO) the inclusion of dispersion interactions, accomplished using the DFT-D3 correction scheme, allows for the correct ordering and an improved calculation of the energy differences between polymorphic phases. The energetic correction provided is of the same order of magnitude as the energy difference between phases. D3-corrected hybrid functionals systematically yield results closest to experiment. We propose that the inclusion of dispersion interactions makes a significant contribution to the relative energetics of polymorphic phases, especially those with different densities, and should therefore be included for calculations of relative energies using DFT methods.
Collapse
|
8
|
Zhou Y, Ouyang Y, Zhang Y, Li Q, Wang J. Machine Learning Assisted Simulations of Electrochemical Interfaces: Recent Progress and Challenges. J Phys Chem Lett 2023; 14:2308-2316. [PMID: 36847421 DOI: 10.1021/acs.jpclett.2c03288] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
The electrochemical interface, where the adsorption of reactants and electrocatalytic reactions take place, has long been a focus of attention. Some of the important processes on it tend to possess relatively slow kinetic characteristics, which are usually beyond the scope of ab initio molecular dynamics. The newly emerging technique, machine learning methods, provides an alternative approach to achieve thousands of atoms and nanosecond time scale while ensuring precision and efficiency. In this Perspective, we summarize in detail the recent progress and achievements made by the introduction of machine learning to simulate electrochemical interfaces, and focus on the limitations of current machine learning models, such as accurate descriptions of long-range electrostatic interactions and the kinetics of the electrochemical reactions occurring at the interface. Finally, we further point out the future directions for machine learning to expand in the field of electrochemical interfaces.
Collapse
Affiliation(s)
- Yipeng Zhou
- School of Physics, Southeast University, Nanjing 211189, China
| | - Yixin Ouyang
- School of Physics, Southeast University, Nanjing 211189, China
| | - Yehui Zhang
- School of Physics, Southeast University, Nanjing 211189, China
| | - Qiang Li
- School of Physics, Southeast University, Nanjing 211189, China
| | - Jinlan Wang
- School of Physics, Southeast University, Nanjing 211189, China
| |
Collapse
|
9
|
Liu P, Wang J, Avargues N, Verdi C, Singraber A, Karsai F, Chen XQ, Kresse G. Combining Machine Learning and Many-Body Calculations: Coverage-Dependent Adsorption of CO on Rh(111). PHYSICAL REVIEW LETTERS 2023; 130:078001. [PMID: 36867825 DOI: 10.1103/physrevlett.130.078001] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/22/2022] [Revised: 11/27/2022] [Accepted: 01/10/2023] [Indexed: 06/18/2023]
Abstract
Adsorption of carbon monoxide (CO) on transition-metal surfaces is a prototypical process in surface sciences and catalysis. Despite its simplicity, it has posed great challenges to theoretical modeling. Pretty much all existing density functionals fail to accurately describe surface energies and CO adsorption site preference as well as adsorption energies simultaneously. Although the random phase approximation (RPA) cures these density functional theory failures, its large computational cost makes it prohibitive to study the CO adsorption for any but the simplest ordered cases. Here, we address these challenges by developing a machine-learned force field (MLFF) with near RPA accuracy for the prediction of coverage-dependent adsorption of CO on the Rh(111) surface through an efficient on-the-fly active learning procedure and a Δ-machine learning approach. We show that the RPA-derived MLFF is capable to accurately predict the Rh(111) surface energy and CO adsorption site preference as well as adsorption energies at different coverages that are all in good agreement with experiments. Moreover, the coverage-dependent ground-state adsorption patterns and adsorption saturation coverage are identified.
Collapse
Affiliation(s)
- Peitao Liu
- University of Vienna, Faculty of Physics and Center for Computational Materials Science, Kolingasse 14-16, A-1090 Vienna, Austria
- Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, 110016 Shenyang, China
| | - Jiantao Wang
- Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, 110016 Shenyang, China
| | - Noah Avargues
- University of Vienna, Faculty of Physics and Center for Computational Materials Science, Kolingasse 14-16, A-1090 Vienna, Austria
| | - Carla Verdi
- University of Vienna, Faculty of Physics and Center for Computational Materials Science, Kolingasse 14-16, A-1090 Vienna, Austria
| | | | - Ferenc Karsai
- VASP Software GmbH, Sensengasse 8, A-1090 Vienna, Austria
| | - Xing-Qiu Chen
- Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, 110016 Shenyang, China
| | - Georg Kresse
- University of Vienna, Faculty of Physics and Center for Computational Materials Science, Kolingasse 14-16, A-1090 Vienna, Austria
- VASP Software GmbH, Sensengasse 8, A-1090 Vienna, Austria
| |
Collapse
|
10
|
Delarmelina M, Dlamini MW, Pattisson S, Davies PR, Hutchings GJ, Catlow CRA. The effect of dissolved chlorides on the photocatalytic degradation properties of titania in wastewater treatment. Phys Chem Chem Phys 2023; 25:4161-4176. [PMID: 36655703 DOI: 10.1039/d2cp03140j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
We investigate the effect of chlorides on the photocatalytic degradation of phenol by titania polymorphs (anatase and rutile). We demonstrate how solubilised chlorides can affect the hydroxyl radical formation on both polymorphs with an overall effect on their photodegradative activity. Initially, the photocatalytic activity of anatase and rutile for phenol degradation is investigated in both standard water and brines. With anatase, a significant reduction of the phenol conversion rate is observed (from a pseudo-first-order rate constant k = 5.3 × 10-3 min-1 to k = 3.5 × 10-3 min-1). In contrast, the presence of solubilised chlorides results in enhancement of rutile activity under the same reaction conditions (from 2.3 × 10-3 min-1 to 4.8 × 10-3 min-1). Periodic DFT methods are extensively employed and we show that after the generation of charge separation in the modelled titania systems, adsorbed chlorides are the preferential site for partial hole localisation, although small energy differences are computed between partially localised hole densities over adsorbed chloride or hydroxyl. Moreover, chlorides can reduce or inhibit the ability of r-TiO2 (110) and a-TiO2 (101) systems to localise polarons in the slab structure. These results indicate that both mechanisms - hole scavenging and the inhibition of hole localisation - can be the origin of the effect of chlorides on photocatalytic activity of both titania polymorphs. These results provide fundamental insight into the photocatalytic properties of titania polymorphs and elucidate the effect of adsorbed anions over radical formation and oxidative decomposition of organic pollutants.
Collapse
Affiliation(s)
- Maicon Delarmelina
- School of Chemistry, Cardiff University, Cardiff, CF10 3AT, UK. .,UK Catalysis Hub, Research Complex at Harwell, STFC Rutherford Appleton Laboratory, Didcot, Oxfordshire, OX11 0FA, UK
| | - Mbongiseni W Dlamini
- UK Catalysis Hub, Research Complex at Harwell, STFC Rutherford Appleton Laboratory, Didcot, Oxfordshire, OX11 0FA, UK.,Max Planck-Cardiff Centre on the Fundamentals of Heterogeneous Catalysis FUNCAT, Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff, CF10 3AT, UK
| | - Samuel Pattisson
- Max Planck-Cardiff Centre on the Fundamentals of Heterogeneous Catalysis FUNCAT, Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff, CF10 3AT, UK
| | - Philip R Davies
- UK Catalysis Hub, Research Complex at Harwell, STFC Rutherford Appleton Laboratory, Didcot, Oxfordshire, OX11 0FA, UK.,Max Planck-Cardiff Centre on the Fundamentals of Heterogeneous Catalysis FUNCAT, Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff, CF10 3AT, UK
| | - Graham J Hutchings
- UK Catalysis Hub, Research Complex at Harwell, STFC Rutherford Appleton Laboratory, Didcot, Oxfordshire, OX11 0FA, UK.,Max Planck-Cardiff Centre on the Fundamentals of Heterogeneous Catalysis FUNCAT, Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff, CF10 3AT, UK
| | - C Richard A Catlow
- School of Chemistry, Cardiff University, Cardiff, CF10 3AT, UK. .,UK Catalysis Hub, Research Complex at Harwell, STFC Rutherford Appleton Laboratory, Didcot, Oxfordshire, OX11 0FA, UK.,Department of Chemistry, University College London, 20 Gordon St., London WC1 HOAJ, UK
| |
Collapse
|
11
|
Fingerhut J, Lecroart L, Borodin D, Schwarzer M, Hörandl S, Kandratsenka A, Auerbach DJ, Wodtke AM, Kitsopoulos TN. Binding Energy and Diffusion Barrier of Formic Acid on Pd(111). J Phys Chem A 2022; 127:142-152. [PMID: 36583672 PMCID: PMC9841570 DOI: 10.1021/acs.jpca.2c07414] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Velocity-resolved kinetics is used to measure the thermal rate of formic acid desorption from Pd(111) between 228 and 273 K for four isotopologues: HCOOH, HCOOD, DCOOH, DCOOD. Upon molecular adsorption, formic acid undergoes decomposition to CO2 and H2 and thermal desorption. To disentangle the contributions of individual processes, we implement a mass-balance-based calibration procedure from which the branching ratio between desorption and decomposition for formic acid is determined. From experimentally derived elementary desorption rate constants, we obtain the binding energy 639 ± 8 meV and the diffusion barrier 370 ± 130 meV using the detailed balance rate model (DBRM). The DBRM explains the observed kinetic isotope effects.
Collapse
Affiliation(s)
- Jan Fingerhut
- Institute
for Physical Chemistry, Georg-August University
of Goettingen, Goettingen 37077, Germany
| | - Loïc Lecroart
- Department
of Dynamics at Surfaces, Max Planck Institute
for Multidisciplinary Sciences, Goettingen 37077, Germany
| | - Dmitriy Borodin
- Institute
for Physical Chemistry, Georg-August University
of Goettingen, Goettingen 37077, Germany,Department
of Dynamics at Surfaces, Max Planck Institute
for Multidisciplinary Sciences, Goettingen 37077, Germany,Email
| | - Michael Schwarzer
- Institute
for Physical Chemistry, Georg-August University
of Goettingen, Goettingen 37077, Germany
| | - Stefan Hörandl
- Institute
for Physical Chemistry, Georg-August University
of Goettingen, Goettingen 37077, Germany
| | - Alexander Kandratsenka
- Department
of Dynamics at Surfaces, Max Planck Institute
for Multidisciplinary Sciences, Goettingen 37077, Germany
| | - Daniel J. Auerbach
- Department
of Dynamics at Surfaces, Max Planck Institute
for Multidisciplinary Sciences, Goettingen 37077, Germany
| | - Alec M. Wodtke
- Institute
for Physical Chemistry, Georg-August University
of Goettingen, Goettingen 37077, Germany,Department
of Dynamics at Surfaces, Max Planck Institute
for Multidisciplinary Sciences, Goettingen 37077, Germany,International
Center for Advanced Studies of Energy Conversion, Georg-August University of Goettingen, Goettingen 37077, Germany
| | - Theofanis N. Kitsopoulos
- Institute
for Physical Chemistry, Georg-August University
of Goettingen, Goettingen 37077, Germany,Department
of Dynamics at Surfaces, Max Planck Institute
for Multidisciplinary Sciences, Goettingen 37077, Germany,Department
of Chemistry, University of Crete, Heraklion 715 00, Greece,Institute
of Electronic Structure and Laser − FORTH, Heraklion 70013, Greece,Email
| |
Collapse
|
12
|
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
|
13
|
Borodin D, Galparsoro O, Rahinov I, Fingerhut J, Schwarzer M, Hörandl S, Auerbach DJ, Kandratsenka A, Schwarzer D, Kitsopoulos TN, Wodtke AM. Steric Hindrance of NH 3 Diffusion on Pt(111) by Co-Adsorbed O-Atoms. J Am Chem Soc 2022; 144:21791-21799. [DOI: 10.1021/jacs.2c10458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Affiliation(s)
- Dmitriy Borodin
- Institute for Physical Chemistry, Georg-August University of Goettingen, Tammannstraße 6, Goettingen37077, Germany
- Department of Dynamics at Surfaces, Max Planck Institute for Multidisciplinary Sciences, Am Fassberg 11, Goettingen37077, Germany
| | - Oihana Galparsoro
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, Donostia-San Sebastián20018, Spain
- Kimika Fakultatea, Euskal Herriko Unibertsitatea UPV/EHU, P.K. 1072Donostia-San Sebastián20018, Spain
| | - Igor Rahinov
- Department of Natural Sciences, The Open University of Israel, Raanana4353701, Israel
| | - Jan Fingerhut
- Institute for Physical Chemistry, Georg-August University of Goettingen, Tammannstraße 6, Goettingen37077, Germany
| | - Michael Schwarzer
- Institute for Physical Chemistry, Georg-August University of Goettingen, Tammannstraße 6, Goettingen37077, Germany
| | - Stefan Hörandl
- Institute for Physical Chemistry, Georg-August University of Goettingen, Tammannstraße 6, Goettingen37077, Germany
| | - Daniel J. Auerbach
- Department of Dynamics at Surfaces, Max Planck Institute for Multidisciplinary Sciences, Am Fassberg 11, Goettingen37077, Germany
| | - Alexander Kandratsenka
- Department of Dynamics at Surfaces, Max Planck Institute for Multidisciplinary Sciences, Am Fassberg 11, Goettingen37077, Germany
| | - Dirk Schwarzer
- Department of Dynamics at Surfaces, Max Planck Institute for Multidisciplinary Sciences, Am Fassberg 11, Goettingen37077, Germany
| | - Theofanis N. Kitsopoulos
- Institute for Physical Chemistry, Georg-August University of Goettingen, Tammannstraße 6, Goettingen37077, Germany
- Department of Dynamics at Surfaces, Max Planck Institute for Multidisciplinary Sciences, Am Fassberg 11, Goettingen37077, Germany
- Department of Chemistry, University of Crete, Heraklion71500, Greece
- Institute of Electronic Structure and Laser − FORTH, Heraklion70013, Greece
| | - Alec M. Wodtke
- Institute for Physical Chemistry, Georg-August University of Goettingen, Tammannstraße 6, Goettingen37077, Germany
- Department of Dynamics at Surfaces, Max Planck Institute for Multidisciplinary Sciences, Am Fassberg 11, Goettingen37077, Germany
- International Center for Advanced Studies of Energy Conversion, Georg-August University of Goettingen, Tammannstraße 6, Goettingen37077, Germany
| |
Collapse
|
14
|
Adsorption energies on transition metal surfaces: towards an accurate and balanced description. Nat Commun 2022; 13:6853. [PMID: 36369277 PMCID: PMC9652424 DOI: 10.1038/s41467-022-34507-y] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2021] [Accepted: 10/27/2022] [Indexed: 11/13/2022] Open
Abstract
Density functional theory predictions of binding energies and reaction barriers provide invaluable data for analyzing chemical transformations in heterogeneous catalysis. For high accuracy, effects of band structure and coverage, as well as the local bond strength in both covalent and non-covalent interactions, must be reliably described and much focus has been put on improving functionals to this end. Here, we show that a correction from higher-level calculations on small metal clusters can be applied to improve periodic band structure adsorption energies and barriers. We benchmark against 38 reliable experimental covalent and non-covalent adsorption energies and five activation barriers with mean absolute errors of 2.2 kcal mol-1, 2.7 kcal mol-1, and 1.1 kcal mol-1, respectively, which are lower than for functionals widely used and tested for surface science evaluations, such as BEEF-vdW and RPBE.
Collapse
|
15
|
Didar BR, Groß A. Solvation structure and dynamics of Li and LiO2 and their transformation in non-aqueous organic electrolyte solvents from first-principles simulations. CHINESE JOURNAL OF CATALYSIS 2022. [DOI: 10.1016/s1872-2067(22)64098-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/31/2022]
|
16
|
Younis U, Qayyum F, Muhammad I, Yaseen M, Sun Q. A Stable Three‐Dimensional Porous Carbon as a High‐Performance Anode Material for Lithium, Sodium, and Potassium Ion Batteries. ADVANCED THEORY AND SIMULATIONS 2022. [DOI: 10.1002/adts.202200230] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Umer Younis
- School of Material Science and Engineering Center for Applied Physics and Technology (CAPT) Peking University Beijing 10087 China
| | - Fizzah Qayyum
- National Centre for Nanoscience and Technology (NCNST), Beiyitiao Zhongguancun Beijing 100190 China
- University of Chinese Academy of Science no. 19A Yuquan Road Beijing 100049 China
| | - Imran Muhammad
- School of Material Science and Engineering Center for Applied Physics and Technology (CAPT) Peking University Beijing 10087 China
| | - Muhammad Yaseen
- Department of Physics University of Agriculture Faisalabad Faisalabad 38000 Pakistan
| | - Qiang Sun
- School of Material Science and Engineering Center for Applied Physics and Technology (CAPT) Peking University Beijing 10087 China
| |
Collapse
|
17
|
Sundararaman R, Vigil-Fowler D, Schwarz K. Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface. Chem Rev 2022; 122:10651-10674. [PMID: 35522135 PMCID: PMC10127457 DOI: 10.1021/acs.chemrev.1c00800] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. We discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.
Collapse
Affiliation(s)
- Ravishankar Sundararaman
- Department of Materials Science and Engineering, Rensselaer Polytechnic Institute, 110 8th Street, Troy, New York 12180, United States
| | - Derek Vigil-Fowler
- Materials, Chemical, and Computational Science Directorate, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Kathleen Schwarz
- Material Measurement Laboratory, National Institute of Standards and Technology, 100 Bureau Drive, Gaithersburg, Maryland 20899, United States
| |
Collapse
|
18
|
Structure of PtRu/Ru(0 0 0 1) and AgPd/Pd(1 1 1) surface alloys: A kinetic Monte Carlo study. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2021.111428] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
19
|
Abstract
Structures and processes at water/metal interfaces play an important technological role in electrochemical energy conversion and storage, photoconversion, sensors, and corrosion, just to name a few. However, they are also of fundamental significance as a model system for the study of solid-liquid interfaces, which requires combining concepts from the chemistry and physics of crystalline materials and liquids. Particularly interesting is the fact that the water-water and water-metal interactions are of similar strength so that the structures at water/metal interfaces result from a competition between these comparable interactions. Because water is a polar molecule and water and metal surfaces are both polarizable, explicit consideration of the electronic degrees of freedom at water/metal interfaces is mandatory. In principle, ab initio molecular dynamics simulations are thus the method of choice to model water/metal interfaces, but they are computationally still rather demanding. Here, ab initio simulations of water/metal interfaces will be reviewed, starting from static systems such as the adsorption of single water molecules, water clusters, and icelike layers, followed by the properties of liquid water layers at metal surfaces. Technical issues such as the appropriate first-principles description of the water-water and water-metal interactions will be discussed, and electrochemical aspects will be addressed. Finally, more approximate but numerically less demanding approaches to treat water at metal surfaces from first-principles will be briefly discussed.
Collapse
Affiliation(s)
- Axel Groß
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany.,Electrochemical Energy Storage, Helmholtz Institute Ulm (HIU), 89069 Ulm, Germany
| | - Sung Sakong
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany
| |
Collapse
|
20
|
Hinsch JJ, Liu J, Wang Y. Reinvestigating oxygen adsorption on Ag(111) by using strongly constrained and appropriately normed semi-local density functional with the revised Vydrov van Voorhis van der Waals force correction. J Chem Phys 2021; 155:234704. [PMID: 34937376 DOI: 10.1063/5.0073407] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
While density functional theory (DFT) at the generalized gradient approximation (GGA) level has made great success in catalysis, it fails in some important systems such as the adsorption of the oxygen molecule on the Ag(111) surface. Previous DFT studies at the GGA level revealed theoretical inconsistencies on the adsorption energies and dissociation barriers of O2 on Ag(111) in comparison with the experimental conclusion. In this study, the strongly constrained and appropriately normed-revised Vydrov van Voorhis van der Waals correction functional (SCAN-rVV10) method at the meta-GGA level with the nonlocal van der Waals (vdW) force correction was used to reinvestigate the adsorption properties of O2 on the Ag(111) surface. The SCAN-rVV10 results successfully confirm the experimental observation that both molecular and dissociative adsorptions can exist for oxygen on Ag(111). The calculated adsorption energy for the physisorption state and the relevant dissociation energy barrier are close to the experimental data. It demonstrates that SCAN-rVV10 can outperform functionals at the GGA level for O2/Ag(111). Therefore, our findings suggest that SCAN-rVV10 can be the desired method for systems where the correct description of intermediate-ranged vdW forces is essential, such as the physisorption of small molecules on the solid surface.
Collapse
Affiliation(s)
- Jack J Hinsch
- School of Environment and Science, Centre for Catalysis and Clean Energy, Griffith University, Gold Coast QLD 4222, Australia
| | - Junxian Liu
- School of Environment and Science, Centre for Catalysis and Clean Energy, Griffith University, Gold Coast QLD 4222, Australia
| | - Yun Wang
- School of Environment and Science, Centre for Catalysis and Clean Energy, Griffith University, Gold Coast QLD 4222, Australia
| |
Collapse
|
21
|
Putra MH, Seidenath S, Kupfer S, Gräfe S, Groß A. Coupling of photoactive transition metal complexes to a functional polymer matrix*. Chemistry 2021; 27:17104-17114. [PMID: 34761834 PMCID: PMC9299502 DOI: 10.1002/chem.202102776] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2021] [Indexed: 11/10/2022]
Abstract
Conductive polymers represent a promising alternative to semiconducting oxide electrodes typically used in dye-sensitized cathodes as they more easily allow a tuning of the physicochemical properties. This can then also be very beneficial for using them in light-driven catalysis. In this computational study, we address the coupling of Ru-based photosensitizers to a polymer matrix by combining two different first-principles electronic structure approaches. We use a periodic density functional theory code to properly account for the delocalized nature of the electronic states in the polymer. These ground state investigations are complemented by time-dependent density functional theory simulations to assess the Franck-Condon photophysics of the present photoactive hybrid material based on a molecular model system. Our results are consistent with recent experimental observations and allow to elucidate the light-driven redox chemical processes - eventually leading to charge separation - in the present functional hybrid systems with potential application as photocathode materials.
Collapse
Affiliation(s)
| | - Sebastian Seidenath
- Institute for Physical Chemistry (IPC) and Abbe Center of Photonics, Friedrich Schiller University Jena, Helmholtzweg 4, 07743, Jena
| | - Stephan Kupfer
- Institute for Physical Chemistry (IPC) and Abbe Center of Photonics, Friedrich Schiller University Jena, Helmholtzweg 4, 07743, Jena
| | - Stefanie Gräfe
- Institute for Physical Chemistry (IPC) and Abbe Center of Photonics, Friedrich Schiller University Jena, Helmholtzweg 4, 07743, Jena
| | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, 89069, Ulm, Germany.,Helmholtz Institute Ulm (HIU), Electrochemical Energy Storage, 89069, Ulm, Germany
| |
Collapse
|
22
|
Didar BR, Yashina L, Groß A. First-Principles Study of the Surfaces and Equilibrium Shape of Discharge Products in Li-Air Batteries. ACS APPLIED MATERIALS & INTERFACES 2021; 13:24984-24994. [PMID: 34009936 DOI: 10.1021/acsami.1c05863] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Li-air batteries are a promising alternative to Li-ion batteries as they theoretically provide the highest possible specific energy density. Mainly, Li2O2 (lithium peroxide) and to a lesser extent, Li2O (lithium oxide) are assumed to be the discharge products of these batteries formed with the soluble LiO2 (lithium superoxide) considered to be an intermediate product. Bulk Li2O2 is an electronic insulator, and the precipitation of this compound on the cathode is thought to be the main limiting factor in achieving high capacities in lithium-oxygen cells. For the most promising electrolytes including solvents with high donor numbers, microscopy observations frequently reveal crystallite morphologies of Li2O2 compounds, rather than uniform layers covering the electrode surface. The precise morphologies of Li2O and Li2O2 particles, and their effect and their extent of contact with the electrode, which may all affect the capacity and rechargeability, however, remain largely undetermined. Here, we address the stability of various Li2O and Li2O2 surfaces and consequently, their crystallite morphologies using density functional theory calculations and ab initio thermodynamics. In contrast to previous studies, we also consider high-index surface terminations, which exhibit surprisingly low surface energies. We carefully analyze the reasons for the stability of these high-index surfaces, which also prominently influence the equilibrium shape of the particles, at least for Li2O2, and discuss the consequences for the observed morphology of the reaction products.
Collapse
Affiliation(s)
| | - Lada Yashina
- Lomonosov Moscow State University, 119991 Moscow, Russia
| | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany
- Helmholtz Institute Ulm (HIU), Electrochemical Energy Storage, 89069 Ulm, Germany
| |
Collapse
|
23
|
Kroes GJ. Computational approaches to dissociative chemisorption on metals: towards chemical accuracy. Phys Chem Chem Phys 2021; 23:8962-9048. [PMID: 33885053 DOI: 10.1039/d1cp00044f] [Citation(s) in RCA: 42] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (Eb, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate Eb can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for Eb for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH4. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of Eb for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate Eb for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O2, H2O, NH3, CO2, and CH3OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands.
| |
Collapse
|
24
|
Forster-Tonigold K, Kim J, Bansmann J, Groß A, Buchner F. Model Studies on the Formation of the Solid Electrolyte Interphase: Reaction of Li with Ultrathin Adsorbed Ionic-Liquid Films and Co 3 O 4 (111) Thin Films. Chemphyschem 2021; 22:441-454. [PMID: 33373085 PMCID: PMC7986933 DOI: 10.1002/cphc.202001033] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2020] [Indexed: 11/15/2022]
Abstract
In this work we aim towards the molecular understanding of the solid electrolyte interphase (SEI) formation at the electrode electrolyte interface (EEI). Herein, we investigated the interaction between the battery‐relevant ionic liquid (IL) 1‐butyl‐1‐methylpyrrolidinium bis(trifluoromethylsulfonyl)imide (BMP‐TFSI), Li and a Co3O4(111) thin film model anode grown on Ir(100) as a model study of the SEI formation in Li‐ion batteries (LIBs). We employed mostly X‐ray photoelectron spectroscopy (XPS) in combination with dispersion‐corrected density functional theory calculations (DFT‐D3). If the surface is pre‐covered by BMP‐TFSI species (model electrolyte), post‐deposition of Li (Li+ ion shuttle) reveals thermodynamically favorable TFSI decomposition products such as LiCN, Li2NSO2CF3, LiF, Li2S, Li2O2, Li2O, but also kinetic products like Li2NCH3C4H9 or LiNCH3C4H9 of BMP. Simultaneously, Li adsorption and/or lithiation of Co3O4(111) to LinCo3O4 takes place due to insertion via step edges or defects; a partial transformation to CoO cannot be excluded. Formation of Co0 could not be observed in the experiment indicating that surface reaction products and inserted/adsorbed Li at the step edges may inhibit or slow down further Li diffusion into the bulk. This study provides detailed insights of the SEI formation at the EEI, which might be crucial for the improvement of future batteries.
Collapse
Affiliation(s)
- Katrin Forster-Tonigold
- Helmholtz Institute Ulm Electrochemical Energy Storage (HIU), Helmholtzstraße 11, 89081, Ulm, Germany.,Karlsruhe Institute of Technology (KIT), P.O. Box 3640, 76021, Karlsruhe, Germany
| | - Jihyun Kim
- Institute of Surface Chemistry and Catalysis, Ulm University, Albert-Einstein-Allee 47, 89081, Ulm, Germany
| | - Joachim Bansmann
- Institute of Surface Chemistry and Catalysis, Ulm University, Albert-Einstein-Allee 47, 89081, Ulm, Germany
| | - Axel Groß
- Helmholtz Institute Ulm Electrochemical Energy Storage (HIU), Helmholtzstraße 11, 89081, Ulm, Germany.,Institute of Theoretical Chemistry, Ulm University, Albert-Einstein-Allee 11, 89081, Ulm, Germany
| | - Florian Buchner
- Institute of Surface Chemistry and Catalysis, Ulm University, Albert-Einstein-Allee 47, 89081, Ulm, Germany
| |
Collapse
|
25
|
Gossenberger F, Juarez F, Groß A. Sulfate, Bisulfate, and Hydrogen Co-adsorption on Pt(111) and Au(111) in an Electrochemical Environment. Front Chem 2020; 8:634. [PMID: 32850652 PMCID: PMC7411137 DOI: 10.3389/fchem.2020.00634] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2020] [Accepted: 06/18/2020] [Indexed: 11/13/2022] Open
Abstract
The co-adsorption of sulfate, bisulfate and hydrogen on Pt(111) and Au(111) electrodes was studied based on periodic density functional calculations with the aqueous electrolyte represented by both explicit and implicit solvent models. The influence of the electrochemical control parameters such as the electrode potential and pH was taken into account in a grand-canonical approach. Thus, phase diagrams of the stable coadsorption phases as a function of the electrochemical potential and Pourbaix diagrams have been derived which well reproduce experimental findings. We demonstrate that it is necessary to include explicit water molecules in order to determine the stable adsorbate phases as the (bi)sulfate adsorbates rows become significantly stabilized by bridging water molecules.
Collapse
Affiliation(s)
| | - Fernanda Juarez
- Institute of Theoretical Chemistry, Ulm University, Ulm, Germany
| | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, Ulm, Germany
| |
Collapse
|
26
|
Wei J, Chen YX, Magnussen OM. Vacancy dynamics on CO-covered Pt(111) electrodes. Chem Commun (Camb) 2020; 56:8309-8312. [PMID: 32573589 DOI: 10.1039/d0cc04046k] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
In situ video-STM studies of Pt(111) electrodes in CO-saturated 0.1 M H2SO4 solution are presented, which reveal the presence of defined point defects in the CO pre-oxidation regime, where the surface is covered by a highly dynamic apparent (1 × 1)-CO adlayer. These defects are generated at the Pt steps and can switch between a mobile and an immobile state. They are assigned to diffusing vacancies in the Pt surface layer, induced by interaction of the highly mobile CO adsorbates with Pt surface atoms at potentials as low as 0.30 VAg/AgCl.
Collapse
Affiliation(s)
- Jie Wei
- Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei 230026, China
| | | | | |
Collapse
|
27
|
Delarmelina M, Quesne MG, Catlow CRA. Modelling the bulk properties of ambient pressure polymorphs of zirconia. Phys Chem Chem Phys 2020; 22:6660-6676. [PMID: 32159203 DOI: 10.1039/d0cp00032a] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report a detailed survey of the calculated bulk properties of zirconia using GGA and meta-GGA functionals (PBE, PBEsol, RPBE, and TPSS), dispersion (Grimme's D2 and D3 approach), and on-site Coulomb repulsion correction (U = 2-8 eV). Structural, elastic, mechanical, and dielectric properties, as well as energetics, electronic structure, and phonon dispersion curves were computed and compared to previous investigations to identify the best DFT approach for a consistent in silico description of zirconia polymorphs. In general, inclusion of dispersion corrections led to only small changes in the calculated properties, whereas DFT+U (U = 2 or 4 eV) reduced the deviations of calculated properties from the experimental results, although deterioration of the structure and relative stabilities may be observed in some cases. Standard PBEsol, RPBE+U, and PBE+U were the best methodologies for a simultaneous description of the three polymorphs of ZrO2. RPBE+U, however, was the only functional to conserve the distinct structures and stabilities of c-, t-, and m-ZrO2 when U = 4 eV was used, resulting in the best in silico replication of the band gaps of ZrO2, whilst outperforming the other methodologies in the description of elastic, mechanical, and dielectric properties of this material. Overall, these results provide insight into the most appropriate DFT methodology for in silico investigations of ZrO2, and show that simultaneous description of all three ambient pressure zirconia polymorphs by DFT techniques with acceptable levels of accuracy can be achieved only when the correct choice of methodology is applied.
Collapse
Affiliation(s)
- Maicon Delarmelina
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, UK. and UK Catalysis Hub, Research Complex at Harwell, STFC Rutherford Appleton Laboratory, Didcot, Oxfordshire OX11 0FA, UK
| | - Matthew G Quesne
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, UK. and UK Catalysis Hub, Research Complex at Harwell, STFC Rutherford Appleton Laboratory, Didcot, Oxfordshire OX11 0FA, UK
| | - C Richard A Catlow
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, UK. and UK Catalysis Hub, Research Complex at Harwell, STFC Rutherford Appleton Laboratory, Didcot, Oxfordshire OX11 0FA, UK
| |
Collapse
|
28
|
Sakong S, Groß A. Water structures on a Pt(111) electrode from ab initio molecular dynamic simulations for a variety of electrochemical conditions. Phys Chem Chem Phys 2020; 22:10431-10437. [PMID: 31976502 DOI: 10.1039/c9cp06584a] [Citation(s) in RCA: 44] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
A structural analysis of solvating water layers on a Pt(111) electrode has been performed based on extensive ab initio molecular dynamics simulations. We have emulated different electrochemical conditions by varying the concentration of hydrogen ions in the water layers, which effectively corresponds to a variation in the electrode potential. We present a detailed analysis of the arrangement and orientation of the water molecules and also address their mobility in the solvation layer.
Collapse
Affiliation(s)
- Sung Sakong
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany.
| | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, and Helmholtz Institute Ulm (HIU) Electrochemical Energy Storage, 89069 Ulm, Germany.
| |
Collapse
|
29
|
Zhao Q, Kulik HJ. Stable Surfaces That Bind Too Tightly: Can Range-Separated Hybrids or DFT+U Improve Paradoxical Descriptions of Surface Chemistry? J Phys Chem Lett 2019; 10:5090-5098. [PMID: 31411023 PMCID: PMC6748670 DOI: 10.1021/acs.jpclett.9b01650] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Accepted: 08/14/2019] [Indexed: 05/25/2023]
Abstract
Approximate, semilocal density functional theory (DFT) suffers from delocalization error that can lead to a paradoxical model of catalytic surfaces that both overbind adsorbates yet are also too stable. We investigate the effect of two widely applied approaches for delocalization error correction, (i) affordable DFT+U (i.e., semilocal DFT augmented with a Hubbard U) and (ii) hybrid functionals with an admixture of Hartree-Fock (HF) exchange, on surface and adsorbate energies across a range of rutile transition metal oxides widely studied for their promise as water-splitting catalysts. We observe strongly row- and period-dependent trends with DFT+U, which increases surface formation energies only in early transition metals (e.g., Ti and V) and decreases adsorbate energies only in later transition metals (e.g., Ir and Pt). Both global and local hybrids destabilize surfaces and reduce adsorbate binding across the periodic table, in agreement with higher-level reference calculations. Density analysis reveals why hybrid functionals correct both quantities, whereas DFT+U does not. We recommend local, range-separated hybrids for the accurate modeling of catalysis in transition metal oxides at only a modest increase in computational cost over semilocal DFT.
Collapse
Affiliation(s)
- Qing Zhao
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, Cambridge, Massachusetts 02139, United States
- Department
of Mechanical Engineering, Massachusetts
Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|