1
|
Zare M, Sahsah D, Saleheen M, Behler J, Heyden A. Hybrid Quantum Mechanical, Molecular Mechanical, and Machine Learning Potential for Computing Aqueous-Phase Adsorption Free Energies on Metal Surfaces. J Chem Theory Comput 2024. [PMID: 39254514 DOI: 10.1021/acs.jctc.4c00869] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/11/2024]
Abstract
Performing reliable computer simulations of elementary processes occurring at metal-water interfaces is pivotal for novel catalyst design in sustainable energy applications. Computational catalyst design hinges on the ability to reliably and efficiently compute the potential energy surface (PES) of the system. Due to the large system sizes needed for studying processes at liquid water-metal interfaces, these systems can currently not be described using density functional theory (DFT). In this work, we used a hybrid quantum mechanical, molecular mechanical, and machine learning potential for studying the adsorption behavior of phenol, atomic hydrogen, 2-butanol, and 2-butanone on the (0001) facet of Ru under reducing conditions when Ru is not oxidized. Specifically, we describe the adsorbate and the surrounding metal atoms at the DFT level of theory. Here, we also considered the electrostatic field effect of the water molecules on adsorbate-metal interactions. Next, for the water-water and water-adsorbate interactions, we used established classical force fields. Finally, for the water-Ru surface interaction, for which no reliable force fields have been published, we used Behler-Parrinello high-dimensional neural network potentials (HDNNPs). Employing this setup, we used our explicit solvation for metal surface (eSMS) approach to compute the aqueous-phase effect on the low-coverage adsorption of selected molecules and atoms on the (0001) facet of Ru. In agreement with previous experimental and computational studies of oxygenated molecules over transition metal facets, we found that liquid water destabilizes the tested adsorbates on Ru(0001). Interestingly, our findings indicate that adsorbates on Ru are less affected by the presence of an aqueous phase than on other transition metals (e.g., Pt), highlighting the necessity of experimental investigations of Ru-based catalytic systems in liquid water.
Collapse
Affiliation(s)
- Mehdi Zare
- Department of Chemical Engineering, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Dia Sahsah
- Department of Chemical Engineering, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Mohammad Saleheen
- Department of Chemical Engineering, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Jörg Behler
- Lehrstuhl für Theoretische Chemie II, Ruhr-Universität Bochum, Bochum 44780, Germany
- Research Center Chemical Sciences and Sustainability, Research Alliance Ruhr, Bochum 44780, Germany
| | - Andreas Heyden
- Department of Chemical Engineering, University of South Carolina, Columbia, South Carolina 29208, United States
| |
Collapse
|
2
|
Wittmann L, Gordiy I, Friede M, Helmich-Paris B, Grimme S, Hansen A, Bursch M. Extension of the D3 and D4 London dispersion corrections to the full actinides series. Phys Chem Chem Phys 2024; 26:21379-21394. [PMID: 39092890 DOI: 10.1039/d4cp01514b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/04/2024]
Abstract
Efficient dispersion corrections are an indispensable component of modern density functional theory, semi-empirical quantum mechanical, and even force field methods. In this work, we extend the well established D3 and D4 London dispersion corrections to the full actinides series, francium, and radium. To keep consistency with the existing versions, the original parameterization strategy of the D4 model was only slightly modified. This includes improved reference Hirshfeld atomic partial charges at the ωB97M-V/ma-def-TZVP level to fit the required electronegativity equilibration charge (EEQ) model. In this context, we developed a new actinide data set called AcQM, which covers the most common molecular actinide compound space. Furthermore, the efficient calculation of dynamic polarizabilities that are needed to construct CAB6 dispersion coefficients was implemented into the ORCA program package. The extended models are assessed for the computation of dissociation curves of actinide atoms and ions, geometry optimizations of crystal structure cutouts, gas-phase structures of small uranium compounds, and an example extracted from a small actinide complex protein assembly. We found that the novel parameterizations perform on par with the computationally more demanding density-dependent VV10 dispersion correction. With the presented extension, the excellent cost-accuracy ratio of the D3 and D4 models can now be utilized in various fields of computational actinide chemistry and, e.g., in efficient composite DFT methods such as r2SCAN-3c. They are implemented in our freely available standalone codes (dftd4, s-dftd3) and the D4 version will be also available in the upcoming ORCA 6.0 program package.
Collapse
Affiliation(s)
- Lukas Wittmann
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Igor Gordiy
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Marvin Friede
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Markus Bursch
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
- FACCTs GmbH, 50677, Köln, Germany
| |
Collapse
|
3
|
Domínguez-Flores F, Kiljunen T, Groß A, Sakong S, Melander MM. Metal-water interface formation: Thermodynamics from ab initio molecular dynamics simulations. J Chem Phys 2024; 161:044705. [PMID: 39056392 DOI: 10.1063/5.0220576] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2024] [Accepted: 07/08/2024] [Indexed: 07/28/2024] Open
Abstract
Metal-water interfaces are central to many electrochemical, (electro)catalytic, and materials science processes and systems. However, our current understanding of their thermodynamic properties is limited by the scarcity of accurate experimental and computational data and procedures. In this work, thermodynamic quantities for metal-water interface formation are computed for a range of FCC(111) surfaces (Pd, Pt, Au, Ag, Rh, and PdAu) through extensive density functional theory based molecular dynamics and the two-phase entropy model. We show that metal-water interface formation is thermodynamically favorable and that most metal surfaces studied in this work are completely wettable, i.e., have contact angles of zero. Interfacial water has higher entropy than bulk water due to the increased population of low-frequency translational modes. The entropic contributions also correlate with the orientational water density, and the highest solvation entropies are observed for interfaces with a moderately ordered first water layer; the entropic contributions account for up to ∼25% of the formation free energy. Water adsorption energy correlates with the water orientation and structure and is found to be a good descriptor of the internal energy part of the interface formation free energy, but it alone cannot satisfactorily explain the interfacial thermodynamics; the interface formation is driven by the competition between energetic and entropic contributions. The obtained results and insight can be used to develop, parameterize, and benchmark theoretical and computational methods for studying metal-water interfaces. Overall, our study yields benchmark-quality data and fundamental insight into the thermodynamic forces driving metal-water interface formation.
Collapse
Affiliation(s)
- Fabiola Domínguez-Flores
- Institute of Theoretical Chemistry, Ulm University, 89081 Ulm, Germany
- Department of Chemistry, Nanoscience Center, University of Jyväskylä, P.O. Box 35 (YN), FI-40014 Jyväskylä, Finland
| | - Toni Kiljunen
- Department of Chemistry, Nanoscience Center, University of Jyväskylä, P.O. Box 35 (YN), FI-40014 Jyväskylä, Finland
| | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, 89081 Ulm, Germany
| | - Sung Sakong
- Institute of Theoretical Chemistry, Ulm University, 89081 Ulm, Germany
| | - Marko M Melander
- Department of Chemistry, Nanoscience Center, University of Jyväskylä, P.O. Box 35 (YN), FI-40014 Jyväskylä, Finland
| |
Collapse
|
4
|
Xu M, Liu S, Vijay S, Bligaard T, Kastlunger G. Benchmarking water adsorption on metal surfaces with ab initio molecular dynamics. J Chem Phys 2024; 160:244707. [PMID: 38920400 DOI: 10.1063/5.0205552] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2024] [Accepted: 06/10/2024] [Indexed: 06/27/2024] Open
Abstract
Solid-water interfaces are ubiquitous in nature and technology. In particular, technologies evolving in the green transition, such as electrocatalysis, heavily rely on the junction of an electrolyte and an electrode as a central part of the device. For the understanding of atomic-scale processes taking place at the electrolyte-electrode interface, density functional theory (DFT) has become the de facto standard. The validation of DFT's ability to simulate the interfacial solid/water interaction is crucial, and ideal simulation setups need to be identified in order to prevent avoidable systematic errors. Here, we develop a rigorous sampling protocol for benchmarking the adsorption/desorption energetics of water on metallic surfaces against experimental temperature programmed desorption, single crystal adsorption calorimetry, and thermal energy atom scattering. We screened DFT's quality on a series of transition metal surfaces, applying three of the most common exchange-correlation approximations: PBE-D3, RPBE-D3, and BEEF-vdW. We find that all three xc-functionals reflect the pseudo-zeroth order desorption of water rooted in the combination of attractive adsorbate-adsorbate interactions and their saturation at low and intermediate coverages, respectively. However, both RPBE-D3 and BEEF-vdW lead to more accurate water adsorption strengths, while PBE-D3 clearly overbinds near-surface water. We relate the variations in binding strength to specific variations in water-metal and water-water interactions, highlighting the structural consequences inherent in an uninformed choice of simulation parameters. Our study gives atomistic insight into water's complex adsorption equilibrium. Furthermore, it represents a guideline for future DFT-based simulations of solvated solid interfaces by providing an assessment of systematic errors in specific setups.
Collapse
Affiliation(s)
- Mianle Xu
- Catalysis Theory Center, Department of Physics, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Sihang Liu
- Catalysis Theory Center, Department of Physics, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Sudarshan Vijay
- Catalysis Theory Center, Department of Physics, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Thomas Bligaard
- Catalysis Theory Center, Department of Physics, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Georg Kastlunger
- Catalysis Theory Center, Department of Physics, Technical University of Denmark, Kongens Lyngby, Denmark
| |
Collapse
|
5
|
Deng H, Huang Y, Li J. Orientational Water Bonding on Pt(111): Beyond the Frontier Orbital Principle. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2023. [PMID: 37494475 DOI: 10.1021/acs.langmuir.3c01545] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
For decades, our understanding of water-metal bonding has been dominated by the frontier orbital principle in which globally stable water-metal interactions are ruled by HOMO interacting with metal surfaces. Using density functional theory calculations, herein, we have revealed that the frontier orbital principle cannot be applied to metastable water bonding on Pt(111), where the decisive role of HOMO is replaced by HOMO-1 in terms of the greatest orbital shifts and depopulations as the two different bonding indicators. Unlike the stable water configuration in which both HOMO-1 and HOMO prefer to overlap with metal states through σ-like orbital interactions, metastable configurations exhibit delicate competition or balance between σ-like and π-like orbital interactions exerted by HOMO-1 and HOMO, respectively. These findings have significantly deepened our understanding of orbital roles in water-metal bonding interactions and bridged the gap between theoretical understanding of electrified waters at electrochemical interfaces and water science on metal surfaces.
Collapse
Affiliation(s)
- Haochang Deng
- School of Materials Science and Engineering, Xiangtan University, Xiangtan 411105, China
| | - Yongli Huang
- School of Materials Science and Engineering, Xiangtan University, Xiangtan 411105, China
| | - Jibiao Li
- School of Materials Science and Engineering, Yangtze Normal University, Chongqing 408100, China
| |
Collapse
|
6
|
Xu D, Pei Z, Yang X, Li Q, Zhang F, Zhu R, Cheng X, Ma L. A Review of Trends in Corrosion-Resistant Structural Steels Research-From Theoretical Simulation to Data-Driven Directions. MATERIALS (BASEL, SWITZERLAND) 2023; 16:ma16093396. [PMID: 37176277 PMCID: PMC10179958 DOI: 10.3390/ma16093396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/23/2023] [Revised: 04/18/2023] [Accepted: 04/20/2023] [Indexed: 05/15/2023]
Abstract
This paper provides a review of models commonly used over the years in the study of microscopic models of material corrosion mechanisms, data mining methods and the corrosion-resistant performance control of structural steels. The virtual process of material corrosion is combined with experimental data to reflect the microscopic mechanism of material corrosion from a nano-scale to macro-scale, respectively. Data mining methods focus on predicting and modeling the corrosion rate and corrosion life of materials. Data-driven control of the corrosion resistance of structural steels is achieved through micro-alloying and organization structure control technology. Corrosion modeling has been used to assess the effects of alloying elements, grain size and organization purity on corrosion resistance, and to determine the contents of alloying elements.
Collapse
Affiliation(s)
- Di Xu
- Institute for Advanced Materials and Technology, University of Science and Technology Beijing, Beijing 100083, China
| | - Zibo Pei
- Institute for Advanced Materials and Technology, University of Science and Technology Beijing, Beijing 100083, China
| | - Xiaojia Yang
- Institute for Advanced Materials and Technology, University of Science and Technology Beijing, Beijing 100083, China
- Shunde Graduate School, University of Science and Technology Beijing, Foshan 528399, China
| | - Qing Li
- Institute for Advanced Materials and Technology, University of Science and Technology Beijing, Beijing 100083, China
| | - Fan Zhang
- Key State Laboratories, Wuhan Research Institute of Materials Protection, Wuhan 430030, China
| | - Renzheng Zhu
- Institute for Advanced Materials and Technology, University of Science and Technology Beijing, Beijing 100083, China
| | - Xuequn Cheng
- Institute for Advanced Materials and Technology, University of Science and Technology Beijing, Beijing 100083, China
- National Materials Corrosion and Protection Scientific Data Center, University of Science and Technology Beijing, Beijing 100083, China
| | - Lingwei Ma
- Institute for Advanced Materials and Technology, University of Science and Technology Beijing, Beijing 100083, China
- National Materials Corrosion and Protection Scientific Data Center, University of Science and Technology Beijing, Beijing 100083, China
| |
Collapse
|
7
|
Zhang KX, Liu ZP. Electrochemical hydrogen evolution on Pt-based catalysts from a theoretical perspective. J Chem Phys 2023; 158:141002. [PMID: 37061480 DOI: 10.1063/5.0142540] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
Hydrogen evolution reaction (HER) by splitting water is a key technology toward a clean energy society, where Pt-based catalysts were long known to have the highest activity under acidic electrochemical conditions but suffer from high cost and poor stability. Here, we overview the current status of Pt-catalyzed HER from a theoretical perspective, focusing on the methodology development of electrochemistry simulation, catalytic mechanism, and catalyst stability. Recent developments in theoretical methods for studying electrochemistry are introduced, elaborating on how they describe solid-liquid interface reactions under electrochemical potentials. The HER mechanism, the reaction kinetics, and the reaction sites on Pt are then summarized, which provides an atomic-level picture of Pt catalyst surface dynamics under reaction conditions. Finally, state-of-the-art experimental solutions to improve catalyst stability are also introduced, which illustrates the significance of fundamental understandings in the new catalyst design.
Collapse
Affiliation(s)
- Ke-Xiang Zhang
- Collaborative Innovation Center of Chemistry for Energy Material, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Key Laboratory of Computational Physical Science, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Zhi-Pan Liu
- Collaborative Innovation Center of Chemistry for Energy Material, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Key Laboratory of Computational Physical Science, Department of Chemistry, Fudan University, Shanghai 200433, China
| |
Collapse
|
8
|
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
|
9
|
Affiliation(s)
- Andrew J. Medford
- School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | | | - Karsten Wedel Jacobsen
- CAMD, Department of Physics, Technical University of Denmark, Kongens Lyngby DK-2800, Denmark
| | - Andrew A. Peterson
- School of Engineering, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
10
|
Muñoz-Santiburcio D. Accurate diffusion coefficients of the excess proton and hydroxide in water via extensive ab initio simulations with different schemes. J Chem Phys 2022; 157:024504. [PMID: 35840376 DOI: 10.1063/5.0093958] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Despite its simple molecular formula, obtaining an accurate in silico description of water is far from straightforward. Many of its very peculiar properties are quite elusive, and in particular, obtaining good estimations of the diffusion coefficients of the solvated proton and hydroxide at a reasonable computational cost has been an unsolved challenge until now. Here, I present extensive results of several unusually long ab initio molecular dynamics (MD) simulations employing different combinations of the Born-Oppenheimer and second-generation Car-Parrinello MD propagation methods with different ensembles (NVE and NVT) and thermostats, which show that these methods together with the RPBE-D3 functional provide a very accurate estimation of the diffusion coefficients of the solvated H3O+ and OH- ions, together with an extremely accurate description of several properties of neutral water (such as the structure of the liquid and its diffusion and shear viscosity coefficients). In addition, I show that the estimations of DH3O+ and DOH- depend dramatically on the simulation length, being necessary to reach timescales in the order of hundreds of picoseconds to obtain reliable results.
Collapse
Affiliation(s)
- Daniel Muñoz-Santiburcio
- CIC nanoGUNE BRTA, Tolosa Hiribidea 76, 20018 San Sebastián, Spain and Instituto de Fusión Nuclear "Guillermo Velarde," Universidad Politécnica de Madrid, C/ José Gutiérrez Abascal 2, 28006 Madrid, Spain
| |
Collapse
|
11
|
Schienbein P, Blumberger J. Nanosecond solvation dynamics of the hematite/liquid water interface at hybrid DFT accuracy using committee neural network potentials. Phys Chem Chem Phys 2022; 24:15365-15375. [PMID: 35703465 DOI: 10.1039/d2cp01708c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Metal oxide/water interfaces play an important role in biology, catalysis, energy storage and photocatalytic water splitting. The atomistic structure at these interfaces is often difficult to characterize by experimental techniques, whilst results from ab initio molecular dynamics simulations tend to be uncertain due to the limited length and time scales accessible. In this work, we train a committee neural network potential to simulate the hematite/water interface at the hybrid DFT level of theory to reach the nanosecond timescale and systems containing more than 3000 atoms. The NNP enables us to converge dynamical properties, not possible with brute-force ab initio molecular dynamics. Our simulations uncover a rich solvation dynamics at the hematite/water interface spanning three different time scales: picosecond H-bond dynamics between surface hydroxyls and the first water layer, in-plane/out-of-plane tilt motion of surface hydroxyls on the 10 ps time scale, and diffusion of water molecules from the oxide surface characterized by a mean residence lifetime of about 60 ps. Calculation of vibrational spectra confirm that H-bonds between surface hydroxyls and first layer water molecules are stronger than H-bonds in bulk water. Our study showcases how state of the art machine learning approaches can routinely be utilized to explore the structural dynamics at transition metal oxide interfaces with complex electronic structure. It foreshadows that c-NNPs are a promising tool to tackle the sampling problem in ab initio electrochemistry with explicit solvent molecules.
Collapse
Affiliation(s)
- Philipp Schienbein
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London, WC1E 6BT, UK.
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London, WC1E 6BT, UK.
| |
Collapse
|
12
|
Thomsen B, Shiga M. Structures of liquid and aqueous water isotopologues at ambient temperature from ab initio path integral simulations. Phys Chem Chem Phys 2022; 24:10851-10859. [PMID: 35504275 DOI: 10.1039/d2cp00499b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The heavy hydrogen isotopes D and T are found in trace amounts in water, and when their concentration increases they can play an intricate role in modulating the physical properties of the liquid. We present an analysis of the microscopic structures of ambient light water (H2O(l)), heavy water (D2O(l)), T2O(l), HDO(aq) and HTO(aq) studied by ab initio path integral molecular dynamics (PIMD). Unlike previous ab initio PIMD investigations of H2O(l) and D2O(l) [Chen et al., Phys. Rev. Lett., 2003, 91, 215503] [Machida et al., J. Chem. Phys., 2017, 148, 102324] we find that D2O(l) is more structured than H2O(l), as is predicted by the experiment. The agreement between the experiment and our simulation for H2O(l) and D2O(l) allows us to accurately predict the intra- and intermolecular structures of T2O(l) HDO(aq) and HTO(aq). T2O(l) is found to have a similar intermolecular structure to that of D2O(l), while the intramolecular structure is more compact, giving rise to a smaller dipole moment than those of H2O(l) and D2O(l). For the mixed isotope species, HDO(aq) and HTO(aq), we find smaller dipole moments and fewer hydrogen bonds when compared with the pure species H2O and D2O. We can attribute this effect to the relative compactness of the mixed isotope species, which results in a lower dipole moment than that of the pure species.
Collapse
Affiliation(s)
- Bo Thomsen
- CCSE, Japan Atomic Energy Agency, 178-4-4, Wakashiba, Kashiwa, Chiba, 277-0871, Japan.
| | - Motoyuki Shiga
- CCSE, Japan Atomic Energy Agency, 178-4-4, Wakashiba, Kashiwa, Chiba, 277-0871, Japan.
| |
Collapse
|
13
|
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]
|
14
|
Mattsson S, Paulus B. First principle calculations including ab initio molecular dynamics studies for the activation of hydrogen fluoride on Ni(111). Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111469] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
15
|
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
|
16
|
Lin X, Shao A, Hua M, Tian X. A first principle study of water adsorbed on flat and stepped silver surfaces. Phys Chem Chem Phys 2022; 24:6803-6810. [DOI: 10.1039/d1cp04618g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The structural, electronic and vibrational properties of a water layer on Ag(100) and Ag(511) have been studied by first principles calculations and ab initio molecular dynamics simulations. The most stable...
Collapse
|
17
|
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: 51] [Impact Index Per Article: 17.0] [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
|
18
|
In situ Raman spectroscopy reveals the structure and dissociation of interfacial water. Nature 2021; 600:81-85. [PMID: 34853456 DOI: 10.1038/s41586-021-04068-z] [Citation(s) in RCA: 283] [Impact Index Per Article: 94.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2020] [Accepted: 09/29/2021] [Indexed: 11/08/2022]
Abstract
Understanding the structure and dynamic process of water at the solid-liquid interface is an extremely important topic in surface science, energy science and catalysis1-3. As model catalysts, atomically flat single-crystal electrodes exhibit well-defined surface and electric field properties, and therefore may be used to elucidate the relationship between structure and electrocatalytic activity at the atomic level4,5. Hence, studying interfacial water behaviour on single-crystal surfaces provides a framework for understanding electrocatalysis6,7. However, interfacial water is notoriously difficult to probe owing to interference from bulk water and the complexity of interfacial environments8. Here, we use electrochemical, in situ Raman spectroscopic and computational techniques to investigate the interfacial water on atomically flat Pd single-crystal surfaces. Direct spectral evidence reveals that interfacial water consists of hydrogen-bonded and hydrated Na+ ion water. At hydrogen evolution reaction (HER) potentials, dynamic changes in the structure of interfacial water were observed from a random distribution to an ordered structure due to bias potential and Na+ ion cooperation. Structurally ordered interfacial water facilitated high-efficiency electron transfer across the interface, resulting in higher HER rates. The electrolytes and electrode surface effects on interfacial water were also probed and found to affect water structure. Therefore, through local cation tuning strategies, we anticipate that these results may be generalized to enable ordered interfacial water to improve electrocatalytic reaction rates.
Collapse
|
19
|
Hydroalcoholic mixtures: the effect of the solvent on the electronic spectrum of the indicator methyl orange. MONATSHEFTE FUR CHEMIE 2021. [DOI: 10.1007/s00706-021-02875-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
20
|
Thomsen B, Shiga M. Ab initio study of nuclear quantum effects on sub- and supercritical water. J Chem Phys 2021; 155:194107. [PMID: 34800944 DOI: 10.1063/5.0071857] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The structures of water in the ambient, subcritical, and supercritical conditions at various densities were studied systematically by ab initio path integral molecular dynamics simulations. It was found that the nuclear quantum effects (NQEs) have a significant impact on the structure of hydrogen bonds in close contact, not only in the ambient condition but also in the sub- and supercritical conditions. The NQEs on the structure beyond the hydrogen bond contact are important in ambient water, but not much for water in the sub- and supercritical conditions. The NQEs are furthermore important for determining the number of hydrogen bonds in the ambient conditions, and this role is, however, diminished in the sub- and supercritical conditions. The NQEs do, nevertheless, show their importance in determining the intramolecular structure of water and the close contact structures of the hydrogen bonds, even at sub- and supercritical conditions. Using the RPBE-D3 functional, the computed radial distribution functions for ambient water are in excellent agreement with experimental data, upgrading our previous results using the BLYP-D2 functional [Machida et al., J. Chem. Phys. 148, 102324 (2018)]. The computed radial distribution functions for water in the sub- and supercritical conditions were carefully compared with experiment. In particular, we found that the first peak in hydrogen pair distribution functions matches only when the NQEs are taken into account.
Collapse
Affiliation(s)
- Bo Thomsen
- CCSE, Japan Atomic Energy Agency, 178-4-4, Wakashiba, Kashiwa, Chiba 277-0871, Japan
| | - Motoyuki Shiga
- CCSE, Japan Atomic Energy Agency, 178-4-4, Wakashiba, Kashiwa, Chiba 277-0871, Japan
| |
Collapse
|
21
|
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
|
22
|
Tesch R, Kowalski PM, Eikerling MH. Properties of the Pt(111)/electrolyte electrochemical interface studied with a hybrid DFT-solvation approach. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:444004. [PMID: 34348250 DOI: 10.1088/1361-648x/ac1aa2] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/26/2021] [Accepted: 08/04/2021] [Indexed: 06/13/2023]
Abstract
Self-consistent modeling of the interface between solid metal electrode and liquid electrolyte is a crucial challenge in computational electrochemistry. In this contribution, we adopt the effective screening medium reference interaction site method (ESM-RISM) to study the charged interface between a Pt(111) surface that is partially covered with chemisorbed oxygen and an aqueous acidic electrolyte. This method proves to be well suited to describe the chemisorption and charging state of the interface at controlled electrode potential. We present an in-depth assessment of the ESM-RISM parameterization and of the importance of computing near-surface water molecules explicitly at the quantum mechanical level. We found that ESM-RISM is able to reproduce some key interface properties, including the peculiar, non-monotonic charging relation of the Pt(111)/electrolyte interface. The comparison with independent theoretical models and explicit simulations of the interface reveals strengths and limitations of ESM-RISM for modeling electrochemical interfaces.
Collapse
Affiliation(s)
- Rebekka Tesch
- Institute of Energy and Climate Research, Theory and Computation of Energy Materials (IEK-13), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
- Chair of Theory and Computation of Energy Materials, Faculty of Georesources and Materials Engineering, RWTH Aachen University, 52062 Aachen, Germany
- Jülich Aachen Research Alliance, JARA-CSD and JARA-ENERGY, 52425 Jülich, Germany
| | - Piotr M Kowalski
- Institute of Energy and Climate Research, Theory and Computation of Energy Materials (IEK-13), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
- Jülich Aachen Research Alliance, JARA-CSD and JARA-ENERGY, 52425 Jülich, Germany
| | - Michael H Eikerling
- Institute of Energy and Climate Research, Theory and Computation of Energy Materials (IEK-13), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
- Chair of Theory and Computation of Energy Materials, Faculty of Georesources and Materials Engineering, RWTH Aachen University, 52062 Aachen, Germany
- Jülich Aachen Research Alliance, JARA-CSD and JARA-ENERGY, 52425 Jülich, Germany
| |
Collapse
|
23
|
Segalina A, Lebègue S, Rocca D, Piccinin S, Pastore M. Structure and Energetics of Dye-Sensitized NiO Interfaces in Water from Ab Initio MD and Large-Scale GW Calculations. J Chem Theory Comput 2021; 17:5225-5238. [PMID: 34324810 DOI: 10.1021/acs.jctc.1c00354] [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/29/2022]
Abstract
The energy-level alignment across solvated molecule/semiconductor interfaces is a crucial property for the correct functioning of dye-sensitized photoelectrodes, where, following the absorption of solar light, a cascade of interfacial hole/electron transfer processes has to efficiently take place. In light of the difficulty of performing X-ray photoelectron spectroscopy measurements at the molecule/solvent/metal-oxide interface, being able to accurately predict the level alignment by first-principles calculations on realistic structural models would represent an important step toward the optimization of the device. In this respect, dye/NiO surfaces, employed in p-type dye-sensitized solar cells, are undoubtedly challenging for ab initio methods and, also for this reason, much less investigated than the n-type dye/TiO2 counterpart. Here, we consider the C343-sensitized NiO surface in water and combine ab initio molecular dynamics (AIMD) simulations with GW (G0W0) calculations, performed along the MD trajectory to reliably describe the structure and energetics of the interface when explicit solvation and finite temperature effects are accounted for. We show that the differential perturbative correction on the NiO and molecule states obtained at the GW level is mandatory to recover the correct (physical) interfacial energetics, allowing hole transfer from the semiconductor valence band to the highest occupied molecular orbital (HOMO) of the dye. Moreover, the calculated average driving force quantitatively agrees with the experimental estimate.
Collapse
Affiliation(s)
- Alekos Segalina
- Université de Lorraine & CNRS, LPCT, UMR 7019, F-54000 Nancy, France
| | - Sébastien Lebègue
- Université de Lorraine & CNRS, LPCT, UMR 7019, F-54000 Nancy, France
| | - Dario Rocca
- Université de Lorraine & CNRS, LPCT, UMR 7019, F-54000 Nancy, France
| | - Simone Piccinin
- Consiglio Nazionale delle Ricerche, Istituto Officina dei Materiali, 34136 Trieste, Italy
| | | |
Collapse
|
24
|
Li H, Brémond E, Sancho-García JC, Adamo C. Pairing double hybrid functionals with a tailored basis set for an accurate thermochemistry of hydrocarbons. RSC Adv 2021; 11:26073-26082. [PMID: 35479441 PMCID: PMC9037073 DOI: 10.1039/d1ra04108h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2021] [Accepted: 07/19/2021] [Indexed: 11/21/2022] Open
Abstract
A collection of five challenging datasets, including noncovalent interactions, reaction barriers and electronic rearrangements of medium-sized hydrocarbons, has been selected to verify the robustness of double-hybrid functionals used in conjunction with the small DH-SVPD basis set, especially developed for noncovalent interactions. The analysis is completed by other, more standard functionals, for a total of 17 models, including also empirical corrections for dispersion. The obtained results show that the chemical accuracy threshold, that is an error lower than 1.0 kcal mol−1, can be obtained by pairing the nonempirical PBE-QIDH functional with the DH-SVPD basis set, as well as by other semi-empirical functionals, such as DSD-PBEP86, using larger basis sets and empirical corrections. More in general, a significant improvement can be obtained using the DH-SVPD basis set with DHs, without resorting to any empirical corrections. This choice leads to a fast computational protocol that, avoiding any empirical potential, remains on a fully quantum ground. The pairing of the PBE-QIDH double-hybrid functional with a tailored split-valence basis set leads to a fast computational protocol for the accurate evaluation of hydrocarbon thermochemistry, without resorting to any empirical correction.![]()
Collapse
Affiliation(s)
- Hanwei Li
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Health and Life Sciences F-75005 Paris France
| | - Eric Brémond
- Université de Paris, ITODYS, CNRS F-75006 Paris France
| | | | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Health and Life Sciences F-75005 Paris France .,Institut Universitaire de France 103 Boulevard Saint Michel F-75005 Paris France
| |
Collapse
|
25
|
Al‐Mahayni H, Wang X, Harvey J, Patience GS, Seifitokaldani A. Experimental methods in chemical engineering: Density functional theory. CAN J CHEM ENG 2021. [DOI: 10.1002/cjce.24127] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
| | - Xiao Wang
- Chemical Engineering McGill University Montréal Québec Canada
| | | | | | | |
Collapse
|
26
|
Kolling I, Hölzl C, Imoto S, Alfarano SR, Vondracek H, Knake L, Sebastiani F, Novelli F, Hoberg C, Brubach JB, Roy P, Forbert H, Schwaab G, Marx D, Havenith M. Aqueous TMAO solution under high hydrostatic pressure. Phys Chem Chem Phys 2021; 23:11355-11365. [PMID: 33972970 DOI: 10.1039/d1cp00703c] [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/21/2022]
Abstract
Trimethylamine N-oxide (TMAO) is a well known osmolyte in nature, which is used by deep sea fish to stabilize proteins against High Hydrostatic Pressure (HHP). We present a combined ab initio molecular dynamics, force field molecular dynamics, and THz absorption study of TMAO in water up to 12 kbar to decipher its solvation properties upon extreme compression. On the hydrophilic oxygen side of TMAO, AIMD simulations at 1 bar and 10 kbar predict a change of the coordination number from a dominating TMAO·(H2O)3 complex at ambient conditions towards an increased population of a TMAO·(H2O)4 complex at HHP conditions. This increase of the TMAO-oxygen coordination number goes in line with a weakening of the local hydrogen bond network, spectroscopic shifts and intensity changes of the corresponding intermolecular THz bands. Using a pressure-dependent HHP force field, FFMD simulations predict a significant increase of hydrophobic hydration from 1 bar up to 4-5 kbar, which levels off at higher pressures up to 10 kbar. THz spectroscopic data reveal two important pressure regimes with spectroscopic inflection points of the dominant intermolecular modes: The first regime (1.5-2 kbar) is barely recognizable in the simulation data. However, it relates well with the observation that the apparent molar volume of solvated TMAO is nearly constant in the biologically relevant pressure range up to 1 kbar as found in the deepest habitats on Earth in the ocean. The second inflection point around 4-5 kbar is related to the amount of hydrophobic hydration as predicted by the FFMD simulations. In particular, the blueshift of the intramolecular CNC bending mode of TMAO at about 390 cm-1 is the spectroscopic signature of increasingly pronounced pressure-induced changes in the solvation shell of TMAO. Thus, the CNC bend can serve as local pressure sensor in the multi-kbar pressure regime.
Collapse
Affiliation(s)
- Inga Kolling
- Lehrstuhl für Physikalische Chemie II, Ruhr-Universität Bochum, 44780 Bochum, Germany.
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
27
|
Lansford JL, Vlachos DG. Spectroscopic Probe Molecule Selection Using Quantum Theory, First-Principles Calculations, and Machine Learning. ACS NANO 2020; 14:17295-17307. [PMID: 33196162 DOI: 10.1021/acsnano.0c07408] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Probe molecule vibrational spectra have a long history of being used to characterize materials including metals, oxides, metal-organic frameworks, and even human proteins. Furthermore, recent advances in machine learning have enabled computationally generated spectra to aid in detailed characterization of complex surfaces with probe molecules. Despite widespread use of probe molecules, the science of probe molecule selection is underdeveloped. Here, we develop physical concepts, including orbital interaction energy and the energy overlap integral, to explain and predict the ability of probe molecules to discriminate structural descriptors. We resolve the crystal orbital overlap population (COOP) to specific molecular orbitals and quantify their bonding character, which directly influences vibrational frequencies. Using only a single adsorbate calculation from density function theory (DFT), we compute the interaction energy of individual adsorbate molecular orbitals with adsorption site atomic orbitals across many different sites. Combining the molecular orbital resolved COOP and changes in orbital interaction energy enables probe molecule selection for improved discrimination of various sites. We demonstrate these concepts by comparing the predicted effectiveness of carbon monoxide (CO), nitric oxide (NO), and ethylene (C2H4) to probe Pt adsorption sites. Finally, using a previously developed machine learning framework, we show that models trained on hundreds of thousands of C2H4 spectra, computed from DFT, which regress surface binding-type and generalized coordination number, outperform those trained using CO and NO spectra. A python package, pDOS_overlap, for implementing the electron density-based analysis on any combination of adsorbates and materials, is also made available.
Collapse
Affiliation(s)
- Joshua L Lansford
- Department of Chemical and Biomolecular Engineering, University of Delaware, 150 Academy Street, Newark, Delaware 19716, United States
| | - Dionisios G Vlachos
- Department of Chemical and Biomolecular Engineering, University of Delaware, 150 Academy Street, Newark, Delaware 19716, United States
- Catalysis Center for Energy Innovation, University of Delaware, 221 Academy Street, Newark, Delaware 19716, United States
| |
Collapse
|
28
|
Abstract
Thinking about water is inextricably linked to hydrogen bonds, which are highly directional in character and determine the unique structure of water, in particular its tetrahedral H-bond network. Here, we assess if this common connotation also holds for supercritical water. We employ extensive ab initio molecular dynamics simulations to systematically monitor the evolution of the H-bond network mode of water from room temperature, where it is the hallmark of its fluctuating three-dimensional network structure, to supercritical conditions. Our simulations reveal that the oscillation period required for H-bond vibrations to occur exceeds the lifetime of H-bonds in supercritical water by far. Instead, the corresponding low-frequency intermolecular vibrations of water pairs as seen in supercritical water are found to be well represented by isotropic van-der-Waals interactions only. Based on these findings, we conclude that water in its supercritical phase is not a H-bonded fluid.
Collapse
Affiliation(s)
- Philipp Schienbein
- Lehrstuhl für Theoretische ChemieRuhr-Universität Bochum44780BochumGermany
| | - Dominik Marx
- Lehrstuhl für Theoretische ChemieRuhr-Universität Bochum44780BochumGermany
| |
Collapse
|
29
|
Affiliation(s)
- Philipp Schienbein
- Lehrstuhl für Theoretische Chemie Ruhr-Universität Bochum 44780 Bochum Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie Ruhr-Universität Bochum 44780 Bochum Germany
| |
Collapse
|
30
|
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
|
31
|
Zhou G, Huang L. A review of recent advances in computational and experimental analysis of first adsorbed water layer on solid substrate. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1786086] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Guobing Zhou
- School of Chemical, Biological and Materials Engineering, University of Oklahoma, Norman, OK, USA
| | - Liangliang Huang
- School of Chemical, Biological and Materials Engineering, University of Oklahoma, Norman, OK, USA
| |
Collapse
|
32
|
Heenen HH, Gauthier JA, Kristoffersen HH, Ludwig T, Chan K. Solvation at metal/water interfaces: An ab initio molecular dynamics benchmark of common computational approaches. J Chem Phys 2020; 152:144703. [DOI: 10.1063/1.5144912] [Citation(s) in RCA: 63] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Hendrik H. Heenen
- Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Joseph A. Gauthier
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, USA
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | | | - Thomas Ludwig
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, USA
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - Karen Chan
- Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| |
Collapse
|
33
|
Abstract
The diffusion of water molecules and clusters across the surfaces of materials is important to a wide range of processes. Interestingly, experiments have shown that on certain substrates, water dimers can diffuse more rapidly than water monomers. Whilst explanations for anomalously fast diffusion have been presented for specific systems, the general underlying physical principles are not yet established. We investigate this through a systematic ab initio study of water monomer and dimer diffusion on a range of surfaces. Calculations reveal different mechanisms for fast water dimer diffusion, which is found to be more widespread than previously anticipated. The key factors affecting diffusion are the balance of water-water versus water-surface bonding and the ease with which hydrogen-bond exchange can occur (either through a classical over-the-barrier process or through quantum-mechanical tunnelling). We anticipate that the insights gained will be useful for understanding future experiments on the diffusion and clustering of hydrogen-bonded adsorbates. The experimental observation that water dimers diffuse more rapidly than monomers across materials’ surfaces is yet to be clarified. Here the authors show by ab initio calculations classical and quantum mechanical mechanisms for faster water dimer diffusion on a broad range of metal and non-metal surfaces.
Collapse
|
34
|
Lansford JL, Vlachos DG. Infrared spectroscopy data- and physics-driven machine learning for characterizing surface microstructure of complex materials. Nat Commun 2020; 11:1513. [PMID: 32251293 PMCID: PMC7089992 DOI: 10.1038/s41467-020-15340-7] [Citation(s) in RCA: 52] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Accepted: 02/20/2020] [Indexed: 11/17/2022] Open
Abstract
There is a need to characterize complex materials and their dynamics under reaction conditions to accelerate materials design. Adsorbate vibrational excitations are selective to adsorbate/surface interactions and infrared (IR) spectra associated with activating adsorbate vibrational modes are accurate, capture details of most modes, and can be obtained operando. Current interpretation depends on heuristic peak assignments for simple spectra, precluding the possibility of obtaining detailed structural information. Here, we combine data-based approaches with chemistry-dependent problem formulation to develop physics-driven surrogate models that generate synthetic IR spectra from first-principles calculations. Using synthetic IR spectra of carbon monoxide on platinum, we implement multinomial regression via neural network ensembles to learn probability distributions functions (pdfs) that describe adsorption sites and quantify uncertainty. We use these pdfs to infer detailed surface microstructure from experimental spectra and extend this methodology to other systems as a first step towards characterizing complex interfaces and closing the materials gap.
Collapse
Affiliation(s)
- Joshua L Lansford
- Department of Chemical Biomolecular Engineering, University of Delaware, 150 Academy Street, Newark, DE, 19716, USA
| | - Dionisios G Vlachos
- Department of Chemical Biomolecular Engineering, University of Delaware, 150 Academy Street, Newark, DE, 19716, USA.
- Catalysis Center for Energy Innovation, University of Delaware, 221 Academy Street, Newark, DE, 19716, USA.
| |
Collapse
|
35
|
Mattsson S, Paulus B. A theoretical study on the binding and electrolytic splitting of hydrogen fluoride on Ni(111) and Ni(211). Phys Chem Chem Phys 2020; 22:4407-4415. [DOI: 10.1039/c9cp06348j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Anhydrous hydrogen fluoride (HF) is electrochemically activated on Ni electrodes in the Simons process. We study the electrode/electrolyte interface using DFT methods and find the splitting of HF to be exothermic even at low cell potentials.
Collapse
Affiliation(s)
- Stefan Mattsson
- Institut für Chemie und Biochemie
- Freie Universität Berlin
- 14195 Berlin
- Germany
| | - Beate Paulus
- Institut für Chemie und Biochemie
- Freie Universität Berlin
- 14195 Berlin
- Germany
| |
Collapse
|
36
|
Schienbein P, Marx D. Assessing the properties of supercritical water in terms of structural dynamics and electronic polarization effects. Phys Chem Chem Phys 2020; 22:10462-10479. [DOI: 10.1039/c9cp05610f] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Evolution of water's structural dynamics from ambient liquid to supercritical dense liquid-like and dilute gas-like conditions.
Collapse
Affiliation(s)
- Philipp Schienbein
- Lehrstuhl für Theoretische Chemie
- Ruhr-Universität Bochum
- 44780 Bochum
- Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie
- Ruhr-Universität Bochum
- 44780 Bochum
- Germany
| |
Collapse
|
37
|
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
|
38
|
Lindgren P, Kastlunger G, Peterson AA. A Challenge to the G ∼ 0 Interpretation of Hydrogen Evolution. ACS Catal 2019. [DOI: 10.1021/acscatal.9b02799] [Citation(s) in RCA: 112] [Impact Index Per Article: 22.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Per Lindgren
- School of Engineering, Brown University, Providence, Rhode Island 02912, United States
| | - Georg Kastlunger
- School of Engineering, Brown University, Providence, Rhode Island 02912, United States
| | - Andrew A. Peterson
- School of Engineering, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
39
|
Hölzl C, Kibies P, Imoto S, Noetzel J, Knierbein M, Salmen P, Paulus M, Nase J, Held C, Sadowski G, Marx D, Kast SM, Horinek D. Structure and thermodynamics of aqueous urea solutions from ambient to kilobar pressures: From thermodynamic modeling, experiments, and first principles simulations to an accurate force field description. Biophys Chem 2019; 254:106260. [PMID: 31522071 DOI: 10.1016/j.bpc.2019.106260] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2019] [Revised: 08/29/2019] [Accepted: 08/29/2019] [Indexed: 12/17/2022]
Abstract
Molecular simulations based on classical force fields are a powerful method for shedding light on the complex behavior of biomolecules in solution. When cosolutes are present in addition to water and biomolecules, subtle balances of weak intermolecular forces have to be accounted for. This imposes high demands on the quality of the underlying force fields, and therefore force field development for small cosolutes is still an active field. Here, we present the development of a new urea force field from studies of urea solutions at ambient and elevated hydrostatic pressures based on a combination of experimental and theoretical approaches. Experimental densities and solvation shell properties from ab initio molecular dynamics simulations at ambient conditions served as the target properties for the force field optimization. Since urea is present in many marine life forms, elevated hydrostatic pressure was rigorously addressed: densities at high pressure were measured by vibrating tube densitometry up to 500 bar and by X-ray absorption up to 5 kbar. Densities were determined by the perturbed-chain statistical associating fluid theory equation of state. Solvation properties were determined by embedded cluster integral equation theory and ab initio molecular dynamics. Our new force field is able to capture the properties of urea solutions at high pressures without further high-pressure adaption, unlike trimethylamine-N-oxide, for which a high-pressure adaption is necessary.
Collapse
Affiliation(s)
- Christoph Hölzl
- Institut für Physikalische und Theoretische Chemie, Universität Regensburg, 93040 Regensburg, Germany
| | - Patrick Kibies
- Physikalische Chemie III, Technische Universität Dortmund, 44227 Dortmund, Germany
| | - Sho Imoto
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Jan Noetzel
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Michael Knierbein
- Laboratory of Thermodynamics, Technische Universität Dortmund, 44227 Dortmund, Germany
| | - Paul Salmen
- Fakultät Physik/DELTA, Technische Universität Dortmund, 44221 Dortmund, Germany
| | - Michael Paulus
- Fakultät Physik/DELTA, Technische Universität Dortmund, 44221 Dortmund, Germany
| | - Julia Nase
- Fakultät Physik/DELTA, Technische Universität Dortmund, 44221 Dortmund, Germany
| | - Christoph Held
- Laboratory of Thermodynamics, Technische Universität Dortmund, 44227 Dortmund, Germany
| | - Gabriele Sadowski
- Laboratory of Thermodynamics, Technische Universität Dortmund, 44227 Dortmund, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, 44227 Dortmund, Germany.
| | - Dominik Horinek
- Institut für Physikalische und Theoretische Chemie, Universität Regensburg, 93040 Regensburg, Germany.
| |
Collapse
|
40
|
Qian J, Ye Y, Yang H, Yano J, Crumlin EJ, Goddard WA. Initial Steps in Forming the Electrode-Electrolyte Interface: H 2O Adsorption and Complex Formation on the Ag(111) Surface from Combining Quantum Mechanics Calculations and Ambient Pressure X-ray Photoelectron Spectroscopy. J Am Chem Soc 2019; 141:6946-6954. [PMID: 30945541 DOI: 10.1021/jacs.8b13672] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The interaction of water with metal surfaces is at the heart of electrocatalysis. But there remain enormous uncertainties about the atomistic interactions at the electrode-electrolyte interface (EEI). As the first step toward an understanding of the EEI, we report here the details of the initial steps of H2O adsorption and complex formation on a Ag(111) surface, based on coupling quantum mechanics (QM) and ambient-pressure X-ray photoelectron spectroscopy (APXPS) experiments. We find a close and direct comparison between simulation and experiment, validated under various isotherm and isobar conditions. We identify five observable oxygen-containing species whose concentrations depend sensitively on temperature and pressure: chemisorbed O* and OH*, H2O* stabilized by hydrogen bond interactions with OH* or O*, and multilayer H2O*. We identify the species experimentally by their O 1s core-level shift that we calculate with QM along with the structures and free energies as a function of temperature and pressure. This leads to a chemical reaction network (CRN) that we use to predict the time evolution of their concentrations over a wide range of temperature (298-798 K) and pressure conditions (10-6-1 Torr), which agree well with the populations determined from APXPS. This multistep simulation CRN protocol should be useful for other heterogeneous catalytic systems.
Collapse
Affiliation(s)
- Jin Qian
- Joint Center for Artificial Photosynthesis , California Institute of Technology , Pasadena , California 91125 , United States.,Materials and Process Simulation Center , California Institute of Technology , Pasadena , California 91125 , United States
| | - Yifan Ye
- Joint Center for Artificial Photosynthesis , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States.,Advanced Light Source , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States.,Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| | - Hao Yang
- Materials and Process Simulation Center , California Institute of Technology , Pasadena , California 91125 , United States
| | - Junko Yano
- Joint Center for Artificial Photosynthesis , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| | - Ethan J Crumlin
- Advanced Light Source , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States.,Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| | - William A Goddard
- Joint Center for Artificial Photosynthesis , California Institute of Technology , Pasadena , California 91125 , United States.,Materials and Process Simulation Center , California Institute of Technology , Pasadena , California 91125 , United States
| |
Collapse
|
41
|
Mahlberg D, Sakong S, Forster-Tonigold K, Groß A. Improved DFT Adsorption Energies with Semiempirical Dispersion Corrections. J Chem Theory Comput 2019; 15:3250-3259. [DOI: 10.1021/acs.jctc.9b00035] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- David Mahlberg
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany
| | - Sung Sakong
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany
| | - Katrin Forster-Tonigold
- Helmholtz Institute Ulm (HIU) for Electrochemical Energy Storage, 89069 Ulm, Germany
- Karlsruhe Institute of Technology (KIT), P.O. Box
3640, 76021 Karlsruhe, Germany
| | - Axel Groß
- Institute of Theoretical Chemistry, Ulm University, 89069 Ulm, Germany
- Helmholtz Institute Ulm (HIU) for Electrochemical Energy Storage, 89069 Ulm, Germany
| |
Collapse
|
42
|
Inagaki T, Takenaka N, Nagaoka M. The crucial role of electron transfer from interfacial molecules in the negative potential shift of Au electrode immersed in ionic liquids. Phys Chem Chem Phys 2018; 20:29362-29373. [PMID: 30451258 DOI: 10.1039/c8cp04594a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Potential of zero charge (PZC) is essential in electrochemistry to understand physical and chemical phenomena at the interface between an electrode and a solution. A negative potential shift from the work function to the PZC has been experimentally observed in a metal/ionic liquid (IL) system, but the mechanism remains unclear and controversial. In this paper we provide valuable insight into the mechanism on the potential shift in the Au/IL (1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)amide: [BMIM][TFSA]) system using a computational approach combining classical molecular dynamics simulations and first-principles calculations. By separately estimating some contributions to the potential shift, the shift is calculated in an easy-to-understand manner. The resultant PZC is shown to be in good agreement with the experimental one. Among the contributions, the electron redistribution at the Au/IL interface is found to provide the largest negative potential change. This indicates that the redistribution plays a crucial role in determining the potential shift of the Au electrode immersed in the IL. Detailed analyses suggest that the redistribution corresponds to the electron transfer not only from the anionic TFSA but also from the cationic BMIM molecules to the Au electrode surface. This unique observation is understood to originate from the interfacial structure where the IL molecules are in very close proximity to the electrode surface via dispersion interaction.
Collapse
Affiliation(s)
- Taichi Inagaki
- Graduate School of Informatics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan.
| | | | | |
Collapse
|
43
|
Gauthier JA, Ringe S, Dickens CF, Garza AJ, Bell AT, Head-Gordon M, Nørskov JK, Chan K. Challenges in Modeling Electrochemical Reaction Energetics with Polarizable Continuum Models. ACS Catal 2018. [DOI: 10.1021/acscatal.8b02793] [Citation(s) in RCA: 114] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Joseph A. Gauthier
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Stefan Ringe
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Colin F. Dickens
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Alejandro J. Garza
- The Joint Center for Artificial Photosynthesis, Lawrence Berkeley National Laboratory, Berkeley California 94720, United States
| | - Alexis T. Bell
- The Joint Center for Artificial Photosynthesis, Lawrence Berkeley National Laboratory, Berkeley California 94720, United States
- Department of Chemical and Biomolecular Engineering, University of California at Berkeley, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- The Joint Center for Artificial Photosynthesis, Lawrence Berkeley National Laboratory, Berkeley California 94720, United States
- Department of Chemistry, University of California at Berkeley, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Jens K. Nørskov
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Karen Chan
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
44
|
García‐Muelas R, Rellán‐Piñeiro M, Li Q, López N. Developments in the Atomistic Modelling of Catalytic Processes for the Production of Platform Chemicals from Biomass. ChemCatChem 2018. [DOI: 10.1002/cctc.201801271] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Rodrigo García‐Muelas
- Institute of Chemical Research of Catalonia, ICIQThe Barcelona Institute of Science and Technology Av. Països Catalans 16 Tarragona 43007 Spain
| | - Marcos Rellán‐Piñeiro
- Institute of Chemical Research of Catalonia, ICIQThe Barcelona Institute of Science and Technology Av. Països Catalans 16 Tarragona 43007 Spain
| | - Qiang Li
- Institute of Chemical Research of Catalonia, ICIQThe Barcelona Institute of Science and Technology Av. Països Catalans 16 Tarragona 43007 Spain
| | - Núria López
- Institute of Chemical Research of Catalonia, ICIQThe Barcelona Institute of Science and Technology Av. Països Catalans 16 Tarragona 43007 Spain
| |
Collapse
|
45
|
Buchner F, Uhl B, Forster-Tonigold K, Bansmann J, Groß A, Behm RJ. Structure formation and surface chemistry of ionic liquids on model electrode surfaces-Model studies for the electrode | electrolyte interface in Li-ion batteries. J Chem Phys 2018; 148:193821. [PMID: 30307189 DOI: 10.1063/1.5012878] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
Ionic liquids (ILs) are considered as attractive electrolyte solvents in modern battery concepts such as Li-ion batteries. Here we present a comprehensive review of the results of previous model studies on the interaction of the battery relevant IL 1-butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide ([BMP]+[TFSI]-) with a series of structurally and chemically well-defined model electrode surfaces, which are increasingly complex and relevant for battery applications [Ag(111), Au(111), Cu(111), pristine and lithiated highly oriented pyrolytic graphite (HOPG), and rutile TiO2(110)]. Combining surface science techniques such as high resolution scanning tunneling microscopy and X-ray photoelectron spectroscopy for characterizing surface structure and chemical composition in deposited (sub-)monolayer adlayers with dispersion corrected density functional theory based calculations, this work aims at a molecular scale understanding of the fundamental processes at the electrode | electrolyte interface, which are crucial for the development of the so-called solid electrolyte interphase (SEI) layer in batteries. Performed under idealized conditions, in an ultrahigh vacuum environment, these model studies provide detailed insights on the structure formation in the adlayer, the substrate-adsorbate and adsorbate-adsorbate interactions responsible for this, and the tendency for chemically induced decomposition of the IL. To mimic the situation in an electrolyte, we also investigated the interaction of adsorbed IL (sub-)monolayers with coadsorbed lithium. Even at 80 K, postdeposited Li is found to react with the IL, leading to decomposition products such as LiF, Li3N, Li2S, LixSOy, and Li2O. In the absence of a [BMP]+[TFSI]- adlayer, it tends to adsorb, dissolve, or intercalate into the substrate (metals, HOPG) or to react with the substrate (TiO2) above a critical temperature, forming LiOx and Ti3+ species in the latter case. Finally, the formation of stable decomposition products was found to sensitively change the equilibrium between surface Li and Li+ intercalated in the bulk, leading to a deintercalation from lithiated HOPG in the presence of an adsorbed IL adlayer at >230 K. Overall, these results provide detailed insights into the surface chemistry at the solid | electrolyte interface and the initial stages of SEI formation at electrode surfaces in the absence of an applied potential, which is essential for the further improvement of future Li-ion batteries.
Collapse
Affiliation(s)
- Florian Buchner
- Helmholtz Institute Ulm Electrochemical Energy Storage (HIU), Helmholtzstraße 11, D-89081 Ulm, Germany
| | - Benedikt Uhl
- Institute of Surface Chemistry and Catalysis, Ulm University, Albert-Einstein-Allee 47, D-89081 Ulm, Germany
| | - Katrin Forster-Tonigold
- Helmholtz Institute Ulm Electrochemical Energy Storage (HIU), Helmholtzstraße 11, D-89081 Ulm, Germany
| | - Joachim Bansmann
- Institute of Surface Chemistry and Catalysis, Ulm University, Albert-Einstein-Allee 47, D-89081 Ulm, Germany
| | - Axel Groß
- Helmholtz Institute Ulm Electrochemical Energy Storage (HIU), Helmholtzstraße 11, D-89081 Ulm, Germany
| | - R Jürgen Behm
- Helmholtz Institute Ulm Electrochemical Energy Storage (HIU), Helmholtzstraße 11, D-89081 Ulm, Germany
| |
Collapse
|
46
|
Sakong S, Groß A. The electric double layer at metal-water interfaces revisited based on a charge polarization scheme. J Chem Phys 2018; 149:084705. [DOI: 10.1063/1.5040056] [Citation(s) in RCA: 94] [Impact Index Per Article: 15.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sung Sakong
- Institute of Theoretical Chemistry, 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
|
47
|
The structure of metal-water interface at the potential of zero charge from density functional theory-based molecular dynamics. J Electroanal Chem (Lausanne) 2018. [DOI: 10.1016/j.jelechem.2017.09.002] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
48
|
Imoto S, Forbert H, Marx D. Aqueous TMAO solutions as seen by theoretical THz spectroscopy: hydrophilic versus hydrophobic water. Phys Chem Chem Phys 2018; 20:6146-6158. [DOI: 10.1039/c7cp07003a] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
All THz resonances of aqueous TMAO solutions are computed and assigned based on ab initio molecular dynamics simulations.
Collapse
Affiliation(s)
- Sho Imoto
- Lehrstuhl für Theoretische Chemie
- Ruhr-Universität Bochum
- 44780 Bochum
- Germany
| | - Harald Forbert
- Center for Solvation Science ZEMOS
- Ruhr-Universität Bochum
- 44780 Bochum
- Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie
- Ruhr-Universität Bochum
- 44780 Bochum
- Germany
| |
Collapse
|
49
|
Cao L, Raciti D, Li C, Livi KJT, Rottmann PF, Hemker KJ, Mueller T, Wang C. Mechanistic Insights for Low-Overpotential Electroreduction of CO2 to CO on Copper Nanowires. ACS Catal 2017. [DOI: 10.1021/acscatal.7b03107] [Citation(s) in RCA: 86] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Liang Cao
- Department
of Chemical and Biomolecular Engineering, ‡Department of Materials Science
and Engineering §Department of Physics and Astronomy, and ∥Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - David Raciti
- Department
of Chemical and Biomolecular Engineering, ‡Department of Materials Science
and Engineering §Department of Physics and Astronomy, and ∥Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Chenyang Li
- Department
of Chemical and Biomolecular Engineering, ‡Department of Materials Science
and Engineering §Department of Physics and Astronomy, and ∥Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Kenneth J. T. Livi
- Department
of Chemical and Biomolecular Engineering, ‡Department of Materials Science
and Engineering §Department of Physics and Astronomy, and ∥Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Paul F. Rottmann
- Department
of Chemical and Biomolecular Engineering, ‡Department of Materials Science
and Engineering §Department of Physics and Astronomy, and ∥Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Kevin J. Hemker
- Department
of Chemical and Biomolecular Engineering, ‡Department of Materials Science
and Engineering §Department of Physics and Astronomy, and ∥Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Tim Mueller
- Department
of Chemical and Biomolecular Engineering, ‡Department of Materials Science
and Engineering §Department of Physics and Astronomy, and ∥Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Chao Wang
- Department
of Chemical and Biomolecular Engineering, ‡Department of Materials Science
and Engineering §Department of Physics and Astronomy, and ∥Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
50
|
Schienbein P, Marx D. Liquid–Vapor Phase Diagram of RPBE-D3 Water: Electronic Properties along the Coexistence Curve and in the Supercritical Phase. J Phys Chem B 2017; 122:3318-3329. [DOI: 10.1021/acs.jpcb.7b09761] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Philipp Schienbein
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|