Development of a Transferable Density-Functional Tight-Binding Model for Organic Molecules at the Water/Platinum Interface.

Journal: Journal of chemical theory and computation
Published Date:

Abstract

A computationally efficient and transferable approach for modeling reactions at metal/water interfaces could significantly accelerate our understanding and ultimately the development of new catalytic transformations, particularly in the context of the emerging field of biomass conversion. Here, we present a parametrization of Pt- ( = H, O, C) density-functional tight-binding (DFTB) for addressing this need. We first constructed Pt-H, Pt-O, and Pt-C repulsive potential splines. These pairwise parameters were then augmented to include many-body interactions using the Chebyshev Interaction Model for Efficient Simulation (ChIMES). We compare the geometrical and energetic performances of both DFTB and DFTB/ChIMES methods with DFT reference data across a variety of organic molecules at the platinum surface from nanoparticles to single-crystal surfaces. DFTB shows limited transferability between extended crystal surfaces and small nanoparticles. This transferability is significantly improved through the introduction of three-body interactions with Pt in DFTB/ChIMES, which provides consistent results across various systems, with reductions in the RMSD from around 30 kcal/mol in DFTB to around 10 kcal/mol. We demonstrate the stability and reliability of the obtained parameters by performing metadynamic simulations for the adsorption of phenol on Pt(111). We observe that DFTB itself is undersolvating the surface, leading to only one or two chemisorbed water molecules in a c(4 × 6) unit cell. In contrast, DFTB/ChIMES leads to a coverage of about 0.5 ML and successfully captures the chemisorbed mode of phenol at both the solid/liquid and the solid/gas interfaces. Furthermore, in agreement with experimental measurements, the adsorption at the solid/liquid interface is significantly weaker than that at the solid/gas interface. Furthermore, we highlight that even with DFTB, where we can accumulate dynamics for more than 1 ns for a given system, the simulations are not fully converged.

Authors

  • Qing Wang
    School of Chemistry and Chemical Engineering, Southwest Petroleum University, Chengdu 610500, China. qwang@163.com.
  • Mingjun Gu
    ENS de Lyon, CNRS, LCH, UMR 5182, 69342 Lyon cedex 07, France.
  • Carine Michel
    CNRS, ENS de Lyon, LCH, UMR 5182, 69342 Lyon cedex 07, France.
  • Nir Goldman
    Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550, United States.
  • Thomas Niehaus
    Institut Lumière Matière, UMR 5306 CNRS, Université Claude Bernard Lyon 1, 10 rue Ada Byron, 69622 Villeurbanne, France.
  • Stephan N Steinmann
    CNRS, ENS de Lyon, LCH, UMR 5182, 69342 Lyon cedex 07, France.

Keywords

No keywords available for this article.