Computational Investigation of Polar Surface/Electrolyte Interfaces

Lead Research Organisation: University of Cambridge
Department Name: Chemistry

Abstract

Polar surfaces of solids are surfaces with a net effective charge density. Polar surfaces can be created by terminating an ionic solid by a plane with a net charge. Also surfaces of ferroelectric compounds are polar when the surface is not parallel to the bulk polarization. Gas phase polar surfaces are unstable against non-stoichiometric reconstruction or depolarization in case of ferroelectric compounds. A first objective of the project is to verify whether polar surfaces can be stabilized by ions in an electrolyte screening the surface charge (electric double layer formation). The second aim is to investigate the surface chemistry of such polar surface-electrolyte interfaces by computing chemisorption energies, surface acidities and redox potentials of electro-active species. These (free) energies will be compared to the corresponding energies of non-polar surfaces of the same compound.

The computational method to be used is Density Functional Theory based Molecular Dynamics (DFTMD) extended with a recently developed finite field technique for the treatment of electric double layers in periodic model systems (arXiv preprint arXiv:1609.07754 (2016)). This study is a first step in a computational exploration of the potential of polar surfaces in heterogeneous catalysis. Model systems that will be considered are main group oxides, such as Al2O3, d0 transition metal oxides such as ZrO2 and ferroelectric solids such as PbTiO3 used in heterogeneous acid base chemistry.

Publications

10 25 50

Studentship Projects

Project Reference Relationship Related To Start End Student Name
EP/N509620/1 01/10/2016 30/09/2022
1800796 Studentship EP/N509620/1 01/10/2016 31/08/2020 Thomas Sayer
 
Title Research data supporting "Finite electric displacement simulations of polar ionic solid-electrolyte interfaces: Application to NaCl(111)/aqueous NaCl solution" 
Description GROMACS files and CP2K files, both input and output, required to reproduce the figures of the publication to which these data belong. That is to say classical and ab initio MD trajectories of a polar NaCl (111) slab in contact with NaCl(aq) electrolyte under a plethora of applied field. See readme.txt for detailed breakdown. GROMACS files were run on a local machine, CP2K was run on Archer with 384 cores and/or Beskow with 256 cores, again all such .pbs input files are provided. 
Type Of Material Database/Collection of data 
Year Produced 2018 
Provided To Others? Yes  
 
Title Stabilization of AgI's polar surfaces by the aqueous environment, and its implications for ice formation (data set) 
Description See also the README file. This dataset contains three sub-directories: (1) D0o0 -- contains input files for performing simulations at D = 0 (2) DCNC -- contains input files for performing simulations at DCNC (2) ECNC -- contains input files for performing simulations at ECNC All simulations can be run with the LAMMPS code (16 Mar 2018). See also https://github.com/uccasco/FiniteFields for additional source code required to apply the E and D fields. Each directory contains the followings files: (a) init.data -- initial structure for pure water in contact with AgI. (b) in.tip4p2005.equil -- input file for performing the initial equilibration of the system at 252K. (c) in.tip4p2005.cool -- input file for performing the cooling ramp simulation between 252K and 242K. (d) in.tip4p2005.constT -- input file for performing a constant T simulation at 242K. The above files perform simulations with an immobile AgI crystal. The DCNC additionally contains a file "in.tip4p2005.constT.mob" which demonstrates the changes needed to perform a simulation with a mobile AgI crystal. (The other equilibration and cooling input files can be similarly adapted.) The ECNC and DCNC directories also contain a file "init.data.electrolyte" which contains an initial structure for NaCl electrolyte in contact with AgI. Please see Table S1 of the article for values of D and E fields used. AgI.table -- tabulated interatomic potential for AgI crystal. 
Type Of Material Database/Collection of data 
Year Produced 2019 
Provided To Others? Yes