1
|
Zhai Z, Chen Q, Wang Y, Ren W, Guo P. Orientational dynamics of the water layer adjacent to Au surface accelerated by polarization effect. J Chem Phys 2024; 160:234704. [PMID: 38884408 DOI: 10.1063/5.0198777] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2024] [Accepted: 05/28/2024] [Indexed: 06/18/2024] Open
Abstract
The orientation and rearrangement of water on a gold electrode significantly influences its physicochemical heterogeneous performance. Despite numerous experimental and theoretical studies aimed at uncovering the structural characteristics of interfacial water, the orientational behavior resulting from electrode-induced rearrangements remains a subject of ongoing debate. Here, we employed molecular dynamics simulations to investigate the adaptive structure and dynamics properties of interfacial water on Au(111) and Au(100) surfaces by considering a polarizable model for Au atoms in comparison with the non-polarizable model. Compared to the nonpolarizable systems, the polarization effect can enhance the interaction between water molecules and the gold surface. Unexpectedly, the rotational dynamics directly associated with the orientational behavior of water adjacent to the gold surface is accelerated, thereby reducing the hydrogen bond lifetime. The underlying mechanism for this anomalous phenomenon originates from the polarization effect, which induces the attraction of the positive hydrogen atoms to the surface by the negative image charge. This leads to a change in orientation that disrupts the hydrogen bonds in the first water layer and subsequently accelerates reorientation dynamics of water molecules adjacent to the gold surface. These results shed light on the intricate interplay between polarization effects and water molecule dynamics on metal surfaces, establishing the foundation for the rational regulation of the orientation of interfacial water.
Collapse
Affiliation(s)
- Zhidong Zhai
- Department of Physics, Shanghai Key Laboratory of High Temperature Superconductors, International Centre of Quantum and Molecular Structures, Shanghai University, Shanghai 200444, China
| | - Qun Chen
- Department of Physics, Shanghai Key Laboratory of High Temperature Superconductors, International Centre of Quantum and Molecular Structures, Shanghai University, Shanghai 200444, China
| | - Yin Wang
- Department of Physics, Shanghai Key Laboratory of High Temperature Superconductors, International Centre of Quantum and Molecular Structures, Shanghai University, Shanghai 200444, China
| | - Wei Ren
- Department of Physics, Shanghai Key Laboratory of High Temperature Superconductors, International Centre of Quantum and Molecular Structures, Shanghai University, Shanghai 200444, China
| | - Pan Guo
- Department of Physics, Shanghai Key Laboratory of High Temperature Superconductors, International Centre of Quantum and Molecular Structures, Shanghai University, Shanghai 200444, China
| |
Collapse
|
2
|
Sriram A, Choi S, Yu X, Brabson LM, Das A, Ulissi Z, Uyttendaele M, Medford AJ, Sholl DS. The Open DAC 2023 Dataset and Challenges for Sorbent Discovery in Direct Air Capture. ACS CENTRAL SCIENCE 2024; 10:923-941. [PMID: 38799660 PMCID: PMC11117325 DOI: 10.1021/acscentsci.3c01629] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
Direct air capture (DAC) of CO2 with porous adsorbents such as metal-organic frameworks (MOFs) has the potential to aid large-scale decarbonization. Previous screening of MOFs for DAC relied on empirical force fields and ignored adsorbed H2O and MOF deformation. We performed quantum chemistry calculations overcoming these restrictions for thousands of MOFs. The resulting data enable efficient descriptions using machine learning.
Collapse
Affiliation(s)
- Anuroop Sriram
- Fundamental AI Research,
Meta AI, Meta, Menlo Park, California 94025, United States
| | - Sihoon Choi
- Fundamental AI Research,
Meta AI, Meta, Menlo Park, California 94025, United States
- School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Xiaohan Yu
- School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Logan M. Brabson
- School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Abhishek Das
- Fundamental AI Research,
Meta AI, Meta, Menlo Park, California 94025, United States
| | - Zachary Ulissi
- Fundamental AI Research,
Meta AI, Meta, Menlo Park, California 94025, United States
| | - Matt Uyttendaele
- Fundamental AI Research,
Meta AI, Meta, Menlo Park, California 94025, United States
| | - Andrew J. Medford
- School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - David S. Sholl
- School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
- Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-2008, United States
| |
Collapse
|
3
|
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
|
4
|
Wu X, Steinmann SN, Michel C. Gaussian attractive potential for carboxylate/cobalt surface interactions. J Chem Phys 2023; 159:164115. [PMID: 37902224 DOI: 10.1063/5.0173351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2023] [Accepted: 10/09/2023] [Indexed: 10/31/2023] Open
Abstract
Ligand-decorated metal surfaces play a pivotal role in various areas of chemistry, particularly in selective catalysis. Molecular dynamics simulations at the molecular mechanics level of theory are best adapted to gain complementary insights to experiments regarding the structure and dynamics of such organic films. However, standard force fields tend to capture only weak physisorption interactions. This is inadequate for ligands that are strongly adsorbed such as carboxylates on metal surfaces. To address this limitation, we employ the Gaussian Lennard-Jones (GLJ) potential, which incorporates an attractive Gaussian potential between the surface and ligand atoms. Here, we develop this approach for the interaction between cobalt surfaces and carboxylate ligands. The accuracy of the GLJ approach is validated through the analysis of the interaction of oxygen with two distinct cobalt surfaces. The accuracy of this method reaches a root mean square deviation (RMSD) of about 3 kcal/mol across all probed configurations, which corresponds to a percentage error of roughly 4%. Application of the GLJ force field to the dynamics of the organic layer on these surfaces reveals how the ligand concentration influences the film order, and highlights differing mobility in the x and y directions, attributable to surface corrugation on Co(112̄0). GLJ is versatile, suitable for a broad range of metal/ligand systems, and can, subsequently, be utilized to study the organic film on the adsorption/desorption of reactants and products during a catalytic process.
Collapse
Affiliation(s)
- Xiaojing Wu
- École Normale Supérieure de Lyon, CNRS, Laboratoire de Chimie UMR 5182, 46 allée d'Italie, F-69364 Lyon, France
| | - Stephan N Steinmann
- École Normale Supérieure de Lyon, CNRS, Laboratoire de Chimie UMR 5182, 46 allée d'Italie, F-69364 Lyon, France
| | - Carine Michel
- École Normale Supérieure de Lyon, CNRS, Laboratoire de Chimie UMR 5182, 46 allée d'Italie, F-69364 Lyon, France
| |
Collapse
|
5
|
Abidi N, Steinmann SN. An Electrostatically Embedded QM/MM Scheme for Electrified Interfaces. ACS APPLIED MATERIALS & INTERFACES 2023; 15:25009-25017. [PMID: 37163568 DOI: 10.1021/acsami.3c01430] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
Atomistic modeling of electrified interfaces remains a major issue for detailed insights in electrocatalysis, corrosion, electrodeposition, batteries, and related devices such as pseudocapacitors. In these domains, the use of grand-canonical density functional theory (GC-DFT) in combination with implicit solvation models has become popular. GC-DFT can be conveniently applied not only to metallic surfaces but also to semiconducting oxides and sulfides and is, furthermore, sufficiently robust to achieve a consistent description of reaction pathways. However, the accuracy of implicit solvation models for solvation effects at interfaces is in general unknown. One promising way to overcome the limitations of implicit solvents is going toward hybrid quantum mechanical (QM)/molecular mechanics (MM) models. For capturing the electrochemical potential dependence, the key quantity is the capacitance, i.e., the relation between the surface charge and the electrochemical potential. In order to retrieve the electrochemical potential from a QM/MM hybrid scheme, an electrostatic embedding is required. Furthermore, the charge of the surface and of the solvent regions has to be strictly opposite in order to consistently simulate charge-neutral unit cells in MM and in QM. To achieve such a QM/MM scheme, we present the implementation of electrostatic embedding in the VASP code. This scheme is broadly applicable to any neutral or charged solid/liquid interface. Here, we demonstrate its use in the context of GC-DFT for the hydrogen evolution reaction (HER) over a noble-metal-free electrocatalyst, MoS2. We investigate the effect of electrostatic embedding compared to the implicit solvent model for three contrasting active sites on MoS2: (i) the sulfur vacancy defect, which is rather apolar; (ii) a Mo antisite defect, where the active site is a surface bound highly polar OH group; and (iii) a reconstructed edge site, which is generally believed to be responsible for most of the catalytic activity. According to our results, the electrostatic embedding leads to almost indistinguishable results compared to the implicit solvent for the apolar system but has a significant effect on polar sites. This demonstrates the reliability of the hybrid QM/MM, electrostatically embedded solvation model for electrified interfaces.
Collapse
Affiliation(s)
- Nawras Abidi
- Ecole Normale Supérieure de Lyon, CNRS, Laboratoire de Chimie UMR 5182, 46 allée d'Italie, F-69364 Lyon, France
| | - Stephan N Steinmann
- Ecole Normale Supérieure de Lyon, CNRS, Laboratoire de Chimie UMR 5182, 46 allée d'Italie, F-69364 Lyon, France
| |
Collapse
|
6
|
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
|
7
|
Steinmann SN, Wang Q, Seh ZW. How machine learning can accelerate electrocatalysis discovery and optimization. MATERIALS HORIZONS 2023; 10:393-406. [PMID: 36541226 DOI: 10.1039/d2mh01279k] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
Advances in machine learning (ML) provide the means to bypass bottlenecks in the discovery of new electrocatalysts using traditional approaches. In this review, we highlight the currently achieved work in ML-accelerated discovery and optimization of electrocatalysts via a tight collaboration between computational models and experiments. First, the applicability of available methods for constructing machine-learned potentials (MLPs), which provide accurate energies and forces for atomistic simulations, are discussed. Meanwhile, the current challenges for MLPs in the context of electrocatalysis are highlighted. Then, we review the recent progress in predicting catalytic activities using surrogate models, including microkinetic simulations and more global proxies thereof. Several typical applications of using ML to rationalize thermodynamic proxies and predict the adsorption and activation energies are also discussed. Next, recent developments of ML-assisted experiments for catalyst characterization, synthesis optimization and reaction condition optimization are illustrated. In particular, the applications in ML-enhanced spectra analysis and the use of ML to interpret experimental kinetic data are highlighted. Additionally, we also show how robotics are applied to high-throughput synthesis, characterization and testing of electrocatalysts to accelerate the materials exploration process and how this equipment can be assembled into self-driven laboratories.
Collapse
Affiliation(s)
| | - Qing Wang
- Univ Lyon, ENS de Lyon, CNRS, Laboratoire de Chimie UMR 5182, Lyon, France.
| | - Zhi Wei Seh
- Institute of Materials Research and Engineering, Agency for Science, Technology and Research (A*STAR), 2 Fusionopolis Way, Innovis, 138634, Singapore.
| |
Collapse
|
8
|
Clabaut P, Beisert M, Michel C, Steinmann SN. Beyond single-crystal surfaces: The GAL21 water/metal force field. J Chem Phys 2022; 157:194705. [DOI: 10.1063/5.0130368] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
Solvent effects are notoriously difficult to describe for metallic nanoparticles (NPs). Here, we introduce GAL21 which is the first pairwise additive force field that is specifically designed to modulate the near chemisorption energy of water as a function of the coordination numbers of the metallic atoms. We find a quadratic dependence to be most suitable for capturing the dependence of the adsorption energy of water on the generalized coordination number (GCN) of the metal atoms. GAL21 has been fitted against DFT adsorption energies for Cu, Ag, Au, Ni, Pd, Pt, and Co on 500 configurations and validated on about 3000 configurations for each metal, constructed on five surfaces with GCNs varying from 2.5 to 11.25. Depending on the metals, the root mean square deviation is found between 0.7 kcal mol−1 (Au) to 1.6 kcal mol−1 (Ni). Using GAL21, as implemented in the open-source code CP2K, we then evaluate the solvation energy of Au55 and Pt55 NPs in water using thermodynamic integration. The solvation free energy is found to be larger for Pt than for Au and systematically larger than 200 kcal mol−1, demonstrating the large impact of solvent on the surface energetics of NPs. Still, given that the amorphous NPs are both, the most stable and the most solvated ones, we do not predict a change in the preferred morphology between the gas-phase and in water. Finally, based on a linear regression on three sizes of NPs (from 38 to 147), the solvation energy for Au and Pt surface atoms is found to be −5.2 and −9.9 kcal mol−1, respectively.
Collapse
Affiliation(s)
- Paul Clabaut
- Ecole Normale Supérieure de Lyon, CNRS, Laboratoire de Chimie UMR 5182, 46 allée d’Italie, F-69364 Lyon, France
| | - Matthieu Beisert
- Ecole Normale Supérieure de Lyon, CNRS, Laboratoire de Chimie UMR 5182, 46 allée d’Italie, F-69364 Lyon, France
| | - Carine Michel
- Ecole Normale Supérieure de Lyon, CNRS, Laboratoire de Chimie UMR 5182, 46 allée d’Italie, F-69364 Lyon, France
| | - Stephan N. Steinmann
- Ecole Normale Supérieure de Lyon, CNRS, Laboratoire de Chimie UMR 5182, 46 allée d’Italie, F-69364 Lyon, France
| |
Collapse
|
9
|
Steinmann SN, Michel C. How to Gain Atomistic Insights on Reactions at the Water/Solid Interface? ACS Catal 2022. [DOI: 10.1021/acscatal.2c00594] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Stephan N. Steinmann
- Ecole Normale Supérieure de Lyon, CNRS, Laboratoire de Chimie
UMR 5182, 46 allée d’Italie, F-69364 Lyon, France
| | - Carine Michel
- Ecole Normale Supérieure de Lyon, CNRS, Laboratoire de Chimie
UMR 5182, 46 allée d’Italie, F-69364 Lyon, France
| |
Collapse
|
10
|
Mikkelsen AEG, Kristoffersen HH, Schiøtz J, Vegge T, Hansen HA, Jacobsen KW. Structure and energetics of liquid water-hydroxyl layers on Pt(111). Phys Chem Chem Phys 2022; 24:9885-9890. [PMID: 35416202 DOI: 10.1039/d2cp00190j] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The interactions between liquid water and hydroxyl species on Pt(111) surfaces have been intensely investigated due to their importance to fuel cell electrocatalysis. Here we present a molecular dynamics study of their structure and energetics using an ensemble of neural network potentials, which allow us to obtain unprecedented statistical sampling. We first study the energetics of hydroxyl formation, where we find a near-linear adsorption energy profile, which exhibits a soft and gradual increase in the differential adsorption energy at high hydroxyl coverages. This is strikingly different from the predictions of the conventional bilayer model, which displays a kink at 1/3ML OH coverage indicating a sizeable jump in differential adsorption energy, but within the statistical uncertainty of previously reported ab initio molecular dynamics studies. We then analyze the structure of the interface, where we provide evidence for the water-OH/Pt(111) interface being hydrophobic at high hydroxyl coverages. We furthermore explain the observed adsorption energetics by analyzing the hydrogen bonding in the water-hydroxyl adlayers, where we argue that the increase in differential adsorption energy at high OH coverage can be explained by a reduction in the number of hydrogen bonds from the adsorbed water molecules to the hydroxyls.
Collapse
Affiliation(s)
- August E G Mikkelsen
- Department of Energy Conversion and Storage, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark.
| | | | - Jakob Schiøtz
- CAMD, Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Tejs Vegge
- Department of Energy Conversion and Storage, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark.
| | - Heine A Hansen
- Department of Energy Conversion and Storage, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark.
| | - Karsten W Jacobsen
- CAMD, Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| |
Collapse
|
11
|
Lee JD, Miller JB, Shneidman AV, Sun L, Weaver JF, Aizenberg J, Biener J, Boscoboinik JA, Foucher AC, Frenkel AI, van der Hoeven JES, Kozinsky B, Marcella N, Montemore MM, Ngan HT, O'Connor CR, Owen CJ, Stacchiola DJ, Stach EA, Madix RJ, Sautet P, Friend CM. Dilute Alloys Based on Au, Ag, or Cu for Efficient Catalysis: From Synthesis to Active Sites. Chem Rev 2022; 122:8758-8808. [PMID: 35254051 DOI: 10.1021/acs.chemrev.1c00967] [Citation(s) in RCA: 28] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
The development of new catalyst materials for energy-efficient chemical synthesis is critical as over 80% of industrial processes rely on catalysts, with many of the most energy-intensive processes specifically using heterogeneous catalysis. Catalytic performance is a complex interplay of phenomena involving temperature, pressure, gas composition, surface composition, and structure over multiple length and time scales. In response to this complexity, the integrated approach to heterogeneous dilute alloy catalysis reviewed here brings together materials synthesis, mechanistic surface chemistry, reaction kinetics, in situ and operando characterization, and theoretical calculations in a coordinated effort to develop design principles to predict and improve catalytic selectivity. Dilute alloy catalysts─in which isolated atoms or small ensembles of the minority metal on the host metal lead to enhanced reactivity while retaining selectivity─are particularly promising as selective catalysts. Several dilute alloy materials using Au, Ag, and Cu as the majority host element, including more recently introduced support-free nanoporous metals and oxide-supported nanoparticle "raspberry colloid templated (RCT)" materials, are reviewed for selective oxidation and hydrogenation reactions. Progress in understanding how such dilute alloy catalysts can be used to enhance selectivity of key synthetic reactions is reviewed, including quantitative scaling from model studies to catalytic conditions. The dynamic evolution of catalyst structure and composition studied in surface science and catalytic conditions and their relationship to catalytic function are also discussed, followed by advanced characterization and theoretical modeling that have been developed to determine the distribution of minority metal atoms at or near the surface. The integrated approach demonstrates the success of bridging the divide between fundamental knowledge and design of catalytic processes in complex catalytic systems, which can accelerate the development of new and efficient catalytic processes.
Collapse
Affiliation(s)
- Jennifer D Lee
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Jeffrey B Miller
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Anna V Shneidman
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Lixin Sun
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Jason F Weaver
- Department of Chemical Engineering, University of Florida, Gainesville, Florida 32611, United States
| | - Joanna Aizenberg
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States.,John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States.,Wyss Institute for Biologically Inspired Engineering, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Juergen Biener
- Nanoscale Synthesis and Characterization Laboratory, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - J Anibal Boscoboinik
- Center for Functional Nanomaterials, Brookhaven National Laboratory, Upton, New York 11973, United States
| | - Alexandre C Foucher
- Department of Materials Science and Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Anatoly I Frenkel
- Department of Materials Science and Chemical Engineering, Stony Brook University, Stony Brook, New York 11794, United States.,Division of Chemistry, Brookhaven National Laboratory, Upton, New York 11973, United States
| | - Jessi E S van der Hoeven
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States.,John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Boris Kozinsky
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Nicholas Marcella
- Department of Materials Science and Chemical Engineering, Stony Brook University, Stony Brook, New York 11794, United States
| | - Matthew M Montemore
- Department of Chemical and Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, United States
| | - Hio Tong Ngan
- Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Christopher R O'Connor
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Cameron J Owen
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States.,John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Dario J Stacchiola
- Center for Functional Nanomaterials, Brookhaven National Laboratory, Upton, New York 11973, United States
| | - Eric A Stach
- Department of Materials Science and Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Robert J Madix
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Philippe Sautet
- Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, California 90095, United States.,Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Cynthia M Friend
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States.,John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| |
Collapse
|
12
|
Ringe S, Hörmann NG, Oberhofer H, Reuter K. Implicit Solvation Methods for Catalysis at Electrified Interfaces. Chem Rev 2021; 122:10777-10820. [PMID: 34928131 PMCID: PMC9227731 DOI: 10.1021/acs.chemrev.1c00675] [Citation(s) in RCA: 49] [Impact Index Per Article: 16.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
Implicit solvation
is an effective, highly coarse-grained approach
in atomic-scale simulations to account for a surrounding liquid electrolyte
on the level of a continuous polarizable medium. Originating in molecular
chemistry with finite solutes, implicit solvation techniques are now
increasingly used in the context of first-principles modeling of electrochemistry
and electrocatalysis at extended (often metallic) electrodes. The
prevalent ansatz to model the latter electrodes and the reactive surface
chemistry at them through slabs in periodic boundary condition supercells
brings its specific challenges. Foremost this concerns the difficulty
of describing the entire double layer forming at the electrified solid–liquid
interface (SLI) within supercell sizes tractable by commonly employed
density functional theory (DFT). We review liquid solvation methodology
from this specific application angle, highlighting in particular its
use in the widespread ab initio thermodynamics approach
to surface catalysis. Notably, implicit solvation can be employed
to mimic a polarization of the electrode’s electronic density
under the applied potential and the concomitant capacitive charging
of the entire double layer beyond the limitations of the employed
DFT supercell. Most critical for continuing advances of this effective
methodology for the SLI context is the lack of pertinent (experimental
or high-level theoretical) reference data needed for parametrization.
Collapse
Affiliation(s)
- Stefan Ringe
- Department of Energy Science and Engineering, Daegu Institute of Science and Technology (DGIST), Daegu 42988, Republic of Korea.,Energy Science & Engineering Research Center, Daegu Institute of Science and Technology (DGIST), Daegu 42988, Republic of Korea
| | - Nicolas G Hörmann
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany.,Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstraße 4, D-85747 Garching, Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstraße 4, D-85747 Garching, Germany.,Chair for Theoretical Physics VII and Bavarian Center for Battery Technology (BayBatt), University of Bayreuth, Universitätsstraße 30, 95447 Bayreuth, Germany
| | - Karsten Reuter
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| |
Collapse
|
13
|
Mikkelsen AEG, Schiøtz J, Vegge T, Jacobsen KW. Is the water/Pt(111) interface ordered at room temperature? J Chem Phys 2021; 155:224701. [PMID: 34911304 DOI: 10.1063/5.0077580] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
The structure of the water/Pt(111) interface has been a subject of debate over the past decades. Here, we report the results of a room temperature molecular dynamics study based on neural network potentials, which allow us to access long time scale simulations while retaining ab initio accuracy. We find that the water/Pt(111) interface is characterized by a double layer composed of a primary, strongly bound adsorption layer with a coverage of ∼0.15 ML, which is coupled to a secondary, weakly bound adsorption layer with a coverage of ∼0.58 ML. By studying the order of the primary adsorption layer, we find that there is an effective repulsion between the adsorbed water molecules, which gives rise to a dynamically changing, semi-ordered interfacial structure, where the water molecules in the primary adsorption layer are distributed homogeneously across the interface, forming frequent hydrogen bonds to water molecules in the secondary adsorption layer. We further show that these conclusions are beyond the time scales accessible to ab initio molecular dynamics.
Collapse
Affiliation(s)
- August E G Mikkelsen
- Department of Energy Conversion and Storage, Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Jakob Schiøtz
- CAMD, Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Tejs Vegge
- Department of Energy Conversion and Storage, Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Karsten W Jacobsen
- CAMD, Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| |
Collapse
|
14
|
Parvulescu VI, Epron F, Garcia H, Granger P. Recent Progress and Prospects in Catalytic Water Treatment. Chem Rev 2021; 122:2981-3121. [PMID: 34874709 DOI: 10.1021/acs.chemrev.1c00527] [Citation(s) in RCA: 67] [Impact Index Per Article: 22.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Presently, conventional technologies in water treatment are not efficient enough to completely mineralize refractory water contaminants. In this context, the implementation of catalytic processes could be an alternative. Despite the advantages provided in terms of kinetics of transformation, selectivity, and energy saving, numerous attempts have not yet led to implementation at an industrial scale. This review examines investigations at different scales for which controversies and limitations must be solved to bridge the gap between fundamentals and practical developments. Particular attention has been paid to the development of solar-driven catalytic technologies and some other emerging processes, such as microwave assisted catalysis, plasma-catalytic processes, or biocatalytic remediation, taking into account their specific advantages and the drawbacks. Challenges for which a better understanding related to the complexity of the systems and the coexistence of various solid-liquid-gas interfaces have been identified.
Collapse
Affiliation(s)
- Vasile I Parvulescu
- Department of Organic Chemistry, Biochemistry and Catalysis, University of Bucharest, B-dul Regina Elisabeta 4-12, Bucharest 030016, Romania
| | - Florence Epron
- Université de Poitiers, CNRS UMR 7285, Institut de Chimie des Milieux et Matériaux de Poitiers (IC2MP), 4 rue Michel Brunet, TSA 51106, 86073 Poitiers Cedex 9, France
| | - Hermenegildo Garcia
- Instituto Universitario de Tecnología Química, Universitat Politecnica de Valencia-Consejo Superior de Investigaciones Científicas, Universitat Politencia de Valencia, Av. de los Naranjos s/n, 46022 Valencia, Spain
| | - Pascal Granger
- CNRS, Centrale Lille, Univ. Artois, UMR 8181 - UCCS - Unité de Catalyse et Chimie du Solide, Univ. Lille, F-59000 Lille, France
| |
Collapse
|
15
|
Dávila López AC, Eggert T, Reuter K, Hörmann NG. Static and dynamic water structures at interfaces: A case study with focus on Pt(111). J Chem Phys 2021; 155:194702. [PMID: 34800953 DOI: 10.1063/5.0067106] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An accurate atomistic treatment of aqueous solid-liquid interfaces necessitates the explicit description of interfacial water ideally via ab initio molecular dynamics simulations. Many applications, however, still rely on static interfacial water models, e.g., for the computation of (electro)chemical reaction barriers and focus on a single, prototypical structure. In this work, we systematically study the relation between density functional theory-derived static and dynamic interfacial water models with specific focus on the water-Pt(111) interface. We first introduce a general construction protocol for static 2D water layers on any substrate, which we apply to the low index surfaces of Pt. Subsequently, we compare these with structures from a broad selection of reference works based on the Smooth Overlap of Atomic Positions descriptor. The analysis reveals some structural overlap between static and dynamic water ensembles; however, static structures tend to overemphasize the in-plane hydrogen bonding network. This feature is especially pronounced for the widely used low-temperature hexagonal ice-like structure. In addition, a complex relation between structure, work function, and adsorption energy is observed, which suggests that the concentration on single, static water models might introduce systematic biases that are likely reduced by averaging over consistently created structural ensembles, as introduced here.
Collapse
Affiliation(s)
| | - Thorben Eggert
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Karsten Reuter
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Nicolas G Hörmann
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| |
Collapse
|
16
|
Rey J, Blanck S, Clabaut P, Loehlé S, Steinmann SN, Michel C. Transferable Gaussian Attractive Potentials for Organic/Oxide Interfaces. J Phys Chem B 2021; 125:10843-10853. [PMID: 34533310 DOI: 10.1021/acs.jpcb.1c05156] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Organic/oxide interfaces play an important role in many areas of chemistry and in particular for lubrication and corrosion. Molecular dynamics simulations are the method of choice for providing complementary insight to experiments. However, the force fields used to simulate the interaction between molecules and oxide surfaces tend to capture only weak physisorption interactions, discarding the stabilizing Lewis acid/base interactions. We here propose a simple complement to the straightforward molecular mechanics description based on "out-of-the-box" Lennard-Jones potentials and electrostatic interactions: the addition of an attractive Gaussian potential between reactive sites of the surface and heteroatoms of adsorbed organic molecules, leading to the Gaussian Lennard-Jones (GLJ) potential. The interactions of four oxygenated and four amine molecules with the typical and widespread hematite and γ-alumina surfaces are investigated. The root mean square deviation (RMSD) for all probed molecules is only 5.7 kcal/mol, which corresponds to an error of 23% over hematite. On γ-alumina, the RMSD is 11.2 kcal/mol using a single parameter for all five chemically inequivalent surface aluminum atoms. Applying GLJ to the simulation of organic films on oxide surfaces demonstrates that the mobility of the surfactants is overestimated by the simplistic LJ potential, while GLJ and other qualitatively correct potentials show a strong structuration and slow dynamics of the surface films, as could be expected from the first-principles adsorption energies for model head groups.
Collapse
Affiliation(s)
- Jérôme Rey
- Université de Lyon, École Normale Supérieure de Lyon, CNRS UMR 5182, Laboratoire de Chimie, 46 allée d'Italie, Lyon F69364, France
| | - Sarah Blanck
- Université de Lyon, École Normale Supérieure de Lyon, CNRS UMR 5182, Laboratoire de Chimie, 46 allée d'Italie, Lyon F69364, France.,Total Marketing & Services, Chemin du Canal-BP 22, Solaize 69360, France
| | - Paul Clabaut
- Université de Lyon, École Normale Supérieure de Lyon, CNRS UMR 5182, Laboratoire de Chimie, 46 allée d'Italie, Lyon F69364, France
| | - Sophie Loehlé
- Total Marketing & Services, Chemin du Canal-BP 22, Solaize 69360, France
| | - Stephan N Steinmann
- Université de Lyon, École Normale Supérieure de Lyon, CNRS UMR 5182, Laboratoire de Chimie, 46 allée d'Italie, Lyon F69364, France
| | - Carine Michel
- Université de Lyon, École Normale Supérieure de Lyon, CNRS UMR 5182, Laboratoire de Chimie, 46 allée d'Italie, Lyon F69364, France
| |
Collapse
|
17
|
Granda-Marulanda LP, McCrum IT, Koper MTM. A simple method to calculate solution-phase free energies of charged species in computational electrocatalysis. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:204001. [PMID: 33761487 DOI: 10.1088/1361-648x/abf19d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/13/2021] [Accepted: 03/24/2021] [Indexed: 06/12/2023]
Abstract
Determining the adsorption potential of adsorbed ions in the field of computational electrocatalysis is of great interest to study their interaction with the electrode material and the solvent, and to map out surface phase diagrams and reaction pathways. Calculating the adsorption potentials of ions with density functional theory and comparing across various ions requires an accurate reference energy of the ion in solution and electrons at the same electrochemical scale. Here we highlight a previously used method for determining the reference free energy of solution phase ions using a simple electrochemical thermodynamic cycle, which allows this free energy to be calculated from that of a neutral gas-phase or solid species and an experimentally measured equilibrium potential, avoiding the need to model solvent around the solution phase ion in the electronic structure calculations. While this method is not new, we describe its use and utility in detail and show that this same method can be used to find the free energy of any ion from any reaction, as long as the half-cell equilibrium potential is known, even for reactions that do not transfer the same number of protons and electrons. To illustrate its usability, we compare the adsorption potentials obtained with DFT of I*, Br*, Cl*, and SO4*on Pt(111) and Au(111) and OH*and Ag*on Pt(111) with those measured experimentally and find that this simple and computationally affordable method reproduces the experimental trends.
Collapse
Affiliation(s)
| | - Ian T McCrum
- Leiden Institute of Chemistry, Leiden University, PO Box 9502, 2300 RA Leiden, The Netherlands
- Department of Chemical & Biomolecular Engineering, Clarkson University, 8 Clarkson Ave., Potsdam, NY 13699, United States of America
| | - Marc T M Koper
- Leiden Institute of Chemistry, Leiden University, PO Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
18
|
Ripani G, Flachmüller A, Peter C, Palleschi A. Coarse-Grained Simulation of the Adsorption of Water on Au(111) Surfaces Using a Modified Stillinger-Weber Potential. ACS OMEGA 2020; 5:31055-31059. [PMID: 33324813 PMCID: PMC7726742 DOI: 10.1021/acsomega.0c04071] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/23/2020] [Accepted: 10/22/2020] [Indexed: 06/12/2023]
Abstract
For reproducing the behavior of water molecules adsorbed on gold surfaces in terms of density of both bulk and interfacial water and in terms of structuring of water on top of gold atoms, the implementation of a multibody potential is necessary, thus the Stillinger-Weber potential was tested. The goal is using a single nonbonded potential for coarse-grained models, without the usage of explicit charges. In order to modify the angular part of the Stillinger-Weber potential from a single cosine to a piecewise function accounting for multiple equilibrium angles, employed for Au-Au-Au and Au-Au-water triplets, it is necessary to create a version of the simulation package LAMMPS that supports the assignment of multiple favored angles. This novel approach is able to reproduce the data obtained using quantum mechanical calculations and density profiles of both bulk and adsorbed water molecules obtained using classical polarizable force fields.
Collapse
Affiliation(s)
- Giorgio Ripani
- Department
of Chemical Science and Technologies, University
of Rome “Tor Vergata”, Via della Ricerca Scientifica, Rome 00133, Italy
| | - Alexander Flachmüller
- Theoretical
Chemistry, University of Konstanz, Konstanz 78547, Baden-Württemberg, Germany
| | - Christine Peter
- Theoretical
Chemistry, University of Konstanz, Konstanz 78547, Baden-Württemberg, Germany
| | - Antonio Palleschi
- Department
of Chemical Science and Technologies, University
of Rome “Tor Vergata”, Via della Ricerca Scientifica, Rome 00133, Italy
| |
Collapse
|
19
|
Clabaut P, Schweitzer B, Götz AW, Michel C, Steinmann SN. Solvation Free Energies and Adsorption Energies at the Metal/Water Interface from Hybrid Quantum-Mechanical/Molecular Mechanics Simulations. J Chem Theory Comput 2020; 16:6539-6549. [PMID: 32931268 DOI: 10.1021/acs.jctc.0c00632] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Modeling adsorption at metal/water interfaces is a cornerstone toward an improved understanding in a variety of fields from heterogeneous catalysis to corrosion. We propose and validate a hybrid scheme that combines the adsorption free energies obtained in the gas phase at the density functional theory level with the variation in solvation from the bulk phase to the interface evaluated using a MM-based alchemical transformation, denoted MMsolv. Using the GAL17 force field for the platinum/water interaction, we retrieve a qualitatively correct interaction energy of the water solvent at the interface. This interaction is of near chemisorption character and thus challenging, both for the alchemical transformation and also for the fixed point-charge electrostatics. Our scheme passes through a state characterized by a well-behaved physisorption potential for the Pt(111)/H2O interaction to converge the free energy difference. The workflow is implemented in the freely available SolvHybrid package. We first assess the adsorption of a water molecule at the Pt/water interface, which turns out to be a stringent test. The intrinsic error of our quantum-mechanical/molecular mechanics (QM/MM) hybrid scheme is limited to 6 kcal mol-1 through the introduction of a correction term to attenuate the electrostatic interaction between near-chemisorbed water molecules and the underlying Pt atoms. Next, we show that the MMsolv solvation free energy of Pt (-0.46 J m-2) is in good agreement with the experimental estimate (-0.32 J m-2). Furthermore, we show that the entropy contribution at room temperature is roughly of equal magnitude as the free energy but with an opposite sign. Finally, we compute the adsorption energy of benzene and phenol at the Pt(111)/water interface, one of the rare systems for which experimental data are available. In qualitative agreement with the experiment, but in stark contrast with a standard implicit solvent model, the adsorption of these aromatic molecules is strongly reduced (i.e., less exothermic by ∼30 and 40 kcal mol-1 for our QM/MM hybrid scheme and experiment, respectively, but ∼0 with the implicit solvent) at the solid/liquid interface compared to the solid/gas interface. This reduction occurs mainly because of the competition between the organic adsorbate and the solvent for adsorption on the metallic surface. The semiquantitative agreement with experimental estimates for the adsorption energy of aromatic molecules thus validates the soundness of our hybrid QM/MM scheme.
Collapse
Affiliation(s)
- Paul Clabaut
- Univ. Lyon, ENS de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratorie de Chimie, F-69342 Lyon, France
| | - Benjamin Schweitzer
- Univ. Lyon, ENS de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratorie de Chimie, F-69342 Lyon, France
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| | - Carine Michel
- Univ. Lyon, ENS de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratorie de Chimie, F-69342 Lyon, France
| | - Stephan N Steinmann
- Univ. Lyon, ENS de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratorie de Chimie, F-69342 Lyon, France
| |
Collapse
|
20
|
Abidi N, Lim KRG, Seh ZW, Steinmann SN. Atomistic modeling of electrocatalysis: Are we there yet? WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1499] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Nawras Abidi
- Univ Lyon, Ens de Lyon, CNRS UMR 5182 Université Claude Bernard Lyon 1, Laboratoire de Chimie, F69342, Lyon France
| | - Kang Rui Garrick Lim
- Institute of Materials Research and Engineering, Agency for Science, Technology and Research (A*STAR) Singapore
| | - Zhi Wei Seh
- Institute of Materials Research and Engineering, Agency for Science, Technology and Research (A*STAR) Singapore
| | - Stephan N. Steinmann
- Univ Lyon, Ens de Lyon, CNRS UMR 5182 Université Claude Bernard Lyon 1, Laboratoire de Chimie, F69342, Lyon France
| |
Collapse
|
21
|
Clabaut P, Staub R, Galiana J, Antonetti E, Steinmann SN. Water adlayers on noble metal surfaces: Insights from energy decomposition analysis. J Chem Phys 2020; 153:054703. [DOI: 10.1063/5.0013040] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Affiliation(s)
- Paul Clabaut
- Univ Lyon, Ecole Normale Supérieure de Lyon, CNRS Université Lyon 1, Laboratoire de Chimie UMR 5182, 46 allée d’Italie, F-69364 Lyon, France
| | - Ruben Staub
- Univ Lyon, Ecole Normale Supérieure de Lyon, CNRS Université Lyon 1, Laboratoire de Chimie UMR 5182, 46 allée d’Italie, F-69364 Lyon, France
| | - Joachim Galiana
- Univ Lyon, Ecole Normale Supérieure de Lyon, CNRS Université Lyon 1, Laboratoire de Chimie UMR 5182, 46 allée d’Italie, F-69364 Lyon, France
| | - Elise Antonetti
- Univ Lyon, Ecole Normale Supérieure de Lyon, CNRS Université Lyon 1, Laboratoire de Chimie UMR 5182, 46 allée d’Italie, F-69364 Lyon, France
| | - Stephan N. Steinmann
- Univ Lyon, Ecole Normale Supérieure de Lyon, CNRS Université Lyon 1, Laboratoire de Chimie UMR 5182, 46 allée d’Italie, F-69364 Lyon, France
| |
Collapse
|
22
|
Clabaut P, Fleurat-Lessard P, Michel C, Steinmann SN. Ten Facets, One Force Field: The GAL19 Force Field for Water-Noble Metal Interfaces. J Chem Theory Comput 2020; 16:4565-4578. [PMID: 32413265 DOI: 10.1021/acs.jctc.0c00091] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Understanding the structure of the water/metal interfaces plays an important role in many areas ranging from surface chemistry to environmental processes. The size, required phase-space sampling, and the slow diffusion of molecules at the water/metal interfaces motivate the development of accurate force fields. We develop and parametrize GAL19, a novel force field, to describe the interaction of water with two facets (111 and 100) of five metals (Pt, Pd, Au, Ag, Cu). To increase transferability compared to its predecessor GAL17, the water-metal interaction is described as a sum of pairwise terms. The interaction energy has three contributions: (i) physisorption is described via a Tang and Toennies potential, (ii) chemisorption and surface corrugation rely on an attractive Gaussian term, and (iii) the angular dependence is explicitly included as a truncated Fourier series. Thirteen parameters are used for each metal surface and were fitted on 250 water adsorption energies computed at the PBE+dDsC level. The performance of GAL19 was evaluated on a set of more than 600 DFT adsorption energies for each surface, leading to an average root-mean-square deviation of only 1 kcal/mol, correctly reproducing the adsorption trends: strong on Pt and Pd but weaker on Ag, Au, and Cu. This force field was then used to simulate the water/metal interface for all ten surfaces for 1 ns. Structural analyses reveal similar tendencies for all surfaces: a first, dense water layer that is mostly adsorbed on the metal top sites and a second layer up to around 6 Å, which is less structured. On Pt and Pd, the first layer is strongly organized with water lying flat on the surface. The pairwise additive functional form allows one to simulate the water adsorption on alloys, which is demonstrated at the example of Ag/Cu and Au/Pt alloys. The water/Ag-Cu interface is predicted to be disordered with water mostly adsorbed on Cu which should exacerbate the Ag reactivity. On the contrary, incorporating Pt into Au materials leads to a structuring of the water interface. Our promising results make GAL19 an ideal candidate to get representative sampling of complex metal/water interfaces as a first step toward accurate estimation of free energies of reactions in solution at the metal interface.
Collapse
Affiliation(s)
- Paul Clabaut
- Univ. Lyon, ENS de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratorie de Chimie, 46 allée d'Italie, F-69364 Lyon, France
| | - Paul Fleurat-Lessard
- Institut de Chimie Moléculaire de l'Université de Bourgogne (ICMUB), UMR CNRS 6302, Université de Bourgogne Franche-Comté (UBFC), 9 avenue Alain Savary 21078 Dijon, France
| | - Carine Michel
- Univ. Lyon, ENS de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratorie de Chimie, 46 allée d'Italie, F-69364 Lyon, France
| | - Stephan N Steinmann
- Univ. Lyon, ENS de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratorie de Chimie, 46 allée d'Italie, F-69364 Lyon, France
| |
Collapse
|
23
|
Kühne TD, Iannuzzi M, Del Ben M, Rybkin VV, Seewald P, Stein F, Laino T, Khaliullin RZ, Schütt O, Schiffmann F, Golze D, Wilhelm J, Chulkov S, Bani-Hashemian MH, Weber V, Borštnik U, Taillefumier M, Jakobovits AS, Lazzaro A, Pabst H, Müller T, Schade R, Guidon M, Andermatt S, Holmberg N, Schenter GK, Hehn A, Bussy A, Belleflamme F, Tabacchi G, Glöß A, Lass M, Bethune I, Mundy CJ, Plessl C, Watkins M, VandeVondele J, Krack M, Hutter J. CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations. J Chem Phys 2020; 152:194103. [PMID: 33687235 DOI: 10.1063/5.0007045] [Citation(s) in RCA: 856] [Impact Index Per Article: 214.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.
Collapse
Affiliation(s)
- Thomas D Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Marcella Iannuzzi
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Mauro Del Ben
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Vladimir V Rybkin
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Patrick Seewald
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Frederick Stein
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Teodoro Laino
- IBM Research Europe, CH-8803 Rüschlikon, Switzerland
| | - Rustam Z Khaliullin
- Department of Chemistry, McGill University, CH-801 Sherbrooke St. West, Montreal, Quebec H3A 0B8, Canada
| | - Ole Schütt
- Department of Materials, ETH Zürich, CH-8092 Zürich, Switzerland
| | | | - Dorothea Golze
- Department of Applied Physics, Aalto University, Otakaari 1, FI-02150 Espoo, Finland
| | - Jan Wilhelm
- Institute of Theoretical Physics, University of Regensburg, Universitätsstraße 31, D-93053 Regensburg, Germany
| | - Sergey Chulkov
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln, United Kingdom
| | | | - Valéry Weber
- IBM Research Europe, CH-8803 Rüschlikon, Switzerland
| | | | | | | | | | - Hans Pabst
- Intel Extreme Computing, Software and Systems, Zürich, Switzerland
| | - Tiziano Müller
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Robert Schade
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Manuel Guidon
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Samuel Andermatt
- Integrated Systems Laboratory, ETH Zürich, CH-8092 Zürich, Switzerland
| | - Nico Holmberg
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076 Aalto, Finland
| | - Gregory K Schenter
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Anna Hehn
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Augustin Bussy
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Fabian Belleflamme
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Gloria Tabacchi
- Department of Science and High Technology, University of Insubria and INSTM, via Valleggio 9, I-22100 Como, Italy
| | - Andreas Glöß
- BASF SE, Carl-Bosch-Straße 38, D-67056 Ludwigshafen am Rhein, Germany
| | - Michael Lass
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Iain Bethune
- Hartree Centre, Science and Technology Facilities Council, Sci-Tech Daresbury, Warrington WA4 4AD, United Kingdom
| | - Christopher J Mundy
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Christian Plessl
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Matt Watkins
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln, United Kingdom
| | - Joost VandeVondele
- Swiss National Supercomputing Centre (CSCS), ETH Zürich, Zürich, Switzerland
| | - Matthias Krack
- Laboratory for Scientific Computing and Modelling, Paul Scherrer Institute, CH-5232 Villigen PSI, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
24
|
Zhang X, Savara A, Getman RB. A Method for Obtaining Liquid-Solid Adsorption Rates from Molecular Dynamics Simulations: Applied to Methanol on Pt(111) in H 2O. J Chem Theory Comput 2020; 16:2680-2691. [PMID: 32134649 DOI: 10.1021/acs.jctc.9b01249] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Adsorption is an important step in heterogeneous catalysis as it predetermines how many reactant molecules can participate in a surface reaction per unit time. While the rate of adsorption processes is well studied in gas-solid adsorption in both theory and experiment, such rates are still not well studied for liquid-solid adsorption. This is partly because the ever-changing configurations of liquid-phase solvent molecules impede the ability to study a molecule approaching a surface from a liquid phase by either experiment or theory. In this work, we develop a method using molecular dynamics (MD) simulations to study the rate of adsorption in liquid-solid adsorption processes. Specifically, we use MD to model the diffusion of a methanol molecule in aqueous solvent and its adsorption to a Pt(111) surface. We find that by approximating the solute motion as following the same displacement rates as a random walk model, the adsorbed and non-adsorbed states of the methanol molecule near the Pt(111) surface can be discerned and quantified. In particular, this methodology enables extracting a sticking coefficient and a macroscopically relatable adsorption rate. This method can be applied to arbitrary types of reactants and surfaces, as well as different liquid environments, thus providing a general tool for predicting quantitative adsorption rates of liquid-solid adsorption systems.
Collapse
Affiliation(s)
- Xiaohong Zhang
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634-0909, United States
| | - Aditya Savara
- Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Rachel B Getman
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634-0909, United States
| |
Collapse
|
25
|
Hardikar RP, Mondal U, Thakkar FM, Roy S, Ghosh P. Theoretical investigations of a platinum-water interface using quantum-mechanics-molecular-mechanics based molecular dynamics simulations. Phys Chem Chem Phys 2019; 21:24345-24353. [PMID: 31663549 DOI: 10.1039/c9cp03558c] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Pt-water interfaces have been of immense interest in the field of energy storage and conversion. Studying this interface using both experimental and theoretical tools is challenging. On the theoretical front, typically one uses classical molecular dynamics (MD) simulations to handle large system sizes or time scales while for a more accurate quantum mechanical description Born Oppenheimer MD (BOMD) is typically used. The latter is limited to smaller system sizes and time-scales. In this study using quantum-mechanics-molecular-mechanics (QMMM), we have performed atomistic MD simulations to have a microscopic understanding of the structure of the Pt-water interface using a system size that is much larger than that accessible when using BOMD simulations. In contrast to recent reports using BOMD simulations, our study reveals that the water molecules typically form two distinct layers above the Pt-surface before they form bulk like structures. Further, we also find that a significant fraction of the water molecules at the interface are pointed towards the surface thereby disrupting the H-bond network. Consistent with this observation, the layer resolved oxygen-oxygen radial distribution function for the water molecules belonging to the solvating water layer shows a high density liquid like behaviour even though the overall water behaves like a low density liquid. A charge transfer analysis reveals that this solvating water layer donates electrons to the Pt atoms in contact with it thereby resulting in the formation of an interface dipole that is pointing towards the surface. Our results suggest that, using QMMM-MD, on one hand it is possible to study more realistic models of solid-liquid interfaces that are inaccessible with BOMD, while on the other hand one also has access to information about such systems that are not obtained from conventional classical MD simulations.
Collapse
Affiliation(s)
- R P Hardikar
- Department of Physics, Indian Institute of Science Education and Research (IISER), Pune 411008, Maharashtra, India.
| | | | | | | | | |
Collapse
|
26
|
Theory and experiments join forces to characterize the electrocatalytic interface. Proc Natl Acad Sci U S A 2019; 116:7611-7613. [PMID: 30923115 DOI: 10.1073/pnas.1903412116] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
|
27
|
|
28
|
NAKANOA H, SATO H. Classical Molecular Dynamics Simulation of Metal Electrodes-Electrolyte Interface. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2019. [DOI: 10.2477/jccj.2018-0040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Hiroshi NAKANOA
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
- Elements Strategy Initiative for Catalysts and Batteries, Kyoto University
| | - Hirofumi SATO
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
- Elements Strategy Initiative for Catalysts and Batteries, Kyoto University
| |
Collapse
|
29
|
Schweitzer B, Steinmann SN, Michel C. Can microsolvation effects be estimated from vacuum computations? A case-study of alcohol decomposition at the H2O/Pt(111) interface. Phys Chem Chem Phys 2019; 21:5368-5377. [DOI: 10.1039/c8cp06331a] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Activation and reaction energies of alcohol decomposition at Pt(111) are barely modified by a PCM, in contrast to adding a single water molecule, whose effect can be predicted based on vacuum computations.
Collapse
Affiliation(s)
- Benjamin Schweitzer
- Univ Lyon
- Ens de Lyon
- CNRS UMR 5182
- Université Claude Bernard Lyon 1
- Laboratoire de Chimie
| | - Stephan N. Steinmann
- Univ Lyon
- Ens de Lyon
- CNRS UMR 5182
- Université Claude Bernard Lyon 1
- Laboratoire de Chimie
| | - Carine Michel
- Univ Lyon
- Ens de Lyon
- CNRS UMR 5182
- Université Claude Bernard Lyon 1
- Laboratoire de Chimie
| |
Collapse
|
30
|
Staub R, Iannuzzi M, Khaliullin RZ, Steinmann SN. Energy Decomposition Analysis for Metal Surface-Adsorbate Interactions by Block Localized Wave Functions. J Chem Theory Comput 2018; 15:265-275. [PMID: 30462497 DOI: 10.1021/acs.jctc.8b00957] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The energy decomposition analysis based on block localized wave functions (BLW-EDA) allows one to gain physical insight into the nature of chemical bonding, decomposing the interaction energy in (1) a "frozen" term, accounting for the attraction due to electrostatic and dispersion interactions, modulated by Pauli repulsion, (2) the variationally assessed polarization energy, and (3) the charge transfer. This method has so far been applied to gas- and condensed-phase molecular systems. However, its standard version is not compatible with fractionally occupied orbitals (i.e., electronic smearing) and, as a consequence, cannot be applied to metallic surfaces. In this work, we propose a simple and practical extension of BLW-EDA to fractionally occupied orbitals, termed Ensemble BLW-EDA. As illustrative examples, we have applied the developed method to analyze the nature of the interaction of various adsorbates on Pt(111), ranging from physisorbed water to strongly chemisorbed ethylene. Our results show that polarization and charge transfer both contribute significantly at the adsorption minimum for all studied systems. The energy decomposition analysis provides details with respect to competing adsorption sites (e.g., CO on atop vs hollow sites) and elucidates the respective importance of polarization and charge transfer for the increased adsorption energy of H2S compared to H2O. Our development will enable a deeper understanding of the impact of charge transfer on catalytic processes in general.
Collapse
Affiliation(s)
- Ruben Staub
- Univ Lyon, Ecole Normale Supérieure de Lyon, CNRS Université Lyon 1 , Laboratoire de Chimie UMR 5182 , 46 allée d'Italie , Lyon , F-69364 , France
| | - Marcella Iannuzzi
- Institut für Chemie , University of Zurich , Winterthurerstrasse 190 , Zurich , CH-8057 , Switzerland
| | - Rustam Z Khaliullin
- Department of Chemistry , McGill University , 801 Sherbrooke Street West , Montreal , Québec H3A 0B8 , Canada
| | - Stephan N Steinmann
- Univ Lyon, Ecole Normale Supérieure de Lyon, CNRS Université Lyon 1 , Laboratoire de Chimie UMR 5182 , 46 allée d'Italie , Lyon , F-69364 , France
| |
Collapse
|