Display options
Share it on

J Chem Theory Comput. 2013 Nov 12;9(11):5086-97. doi: 10.1021/ct400698y. Epub 2013 Oct 11.

Simulation of Adsorption Processes at Metallic Interfaces: An Image Charge Augmented QM/MM Approach.

Journal of chemical theory and computation

Dorothea Golze, Marcella Iannuzzi, Manh-Thuong Nguyen, Daniele Passerone, Jürg Hutter

Affiliations

  1. Institute of Physical Chemistry, University of Zürich , Winterthurerstrasse 190, CH-8057 Zürich, Switzerland.
  2. The Abdus Salam International Centre for Theoretical Physics , Strada Costiera 11, I-34151 Trieste, Italy.
  3. Empa, Swiss Federal Laboratories for Materials Science and Technology , nanotech@surfaces Laboratory, Überlandstrasse 129, CH-8600 Dübendorf, Switzerland.

PMID: 26583423 DOI: 10.1021/ct400698y

Abstract

A novel method for including polarization effects within hybrid quantum mechanics/molecular mechanics (QM/MM) simulations of adsorbate-metal systems is presented. The interactions between adsorbate (QM) and metallic substrate (MM) are described at the MM level of theory. Induction effects are additionally accounted for by applying the image charge formulation. The charge distribution induced within the metallic substrate is modeled by a set of Gaussian charges (image charges) centered at the metal atoms. The image charges and the electrostatic response of the QM potential are determined self-consistently by imposing the constant-potential condition within the metal. The implementation is embedded in a highly efficient Gaussian and plane wave framework and is naturally suited for periodic systems. Even though the electronic properties of the metallic substrate are not taken into account explicitly, the augmented QM/MM scheme can reproduce characteristic polarization effects of the adsorbate. The method is assessed through the investigation of structural and electronic properties of benzene, nitrobenzene, thymine, and guanine on Au(111). The study of small water clusters adsorbed on Pt(111) is also reported in order to demonstrate that the approach provides a sizable correction of the MM-based interactions between adsorbate and substrate. Large-scale molecular dynamics (MD) simulations of a water film in contact with a Pt(111) surface show that the method is suitable for simulations of liquid/metal interfaces at reduced computational cost.

Publication Types