skip to content

Department of Chemistry


Computational physical chemistry

The broad area of the research activity in the group is computational physical chemistry. The main computational tool is Density Functional Theory based Molecular Dynamics (DFTMD or "Car-Parrinello"). Reactive species and the condensed phase environment are treated at strictly the same level of theory ("all-atom"). Using this fundamental approach we study reactivity in solutions, solids and at liquid/solid interfaces computing equilibrium constants, activation energies, electronic and vibrational spectra and other physical quantities such as dielectric properties. The DFTMD code we use is CP2K( and occasionally also Quantum Espresso(

Interfacial electrochemistry

Computational electrochemistry is an application where all-atom DFTMD methods are particularly instructive and has become a major research topic for us. Our favourite model systems are the transition metal oxide electrodes used in energy conversion and storage. Proton coupled electron transfer (PCET) is a key process in these reactions. Unified treatment of oxidation/reduction and (de)protonation reactions is therefore crucial. Combining DFTMD and free energy perturbation (FTP) methods we have developed such a scheme. The central computational tool is reversible insertion of electrons and protons in periodic DFTMD model systems. Applying this method we can refer the computed redox free energies and electronic energy levels to the Standard Hydrogen Electrode (SHE) without having to introduce an additional interface with vacuum. We are using this approach to study the level alignment at electrochemical interfaces of typical transition metal oxides such as TiO2 and MnO2. The aim is to correlate the catalytic properties of the electrochemical interface to the electronic structure.

Figure: DFTMD model system of a TiO2 rutile (110) water interface at the point of zero net proton charge (PZC) consisting of 48 TiO2 units (Ti ions in yellow) and 71 H2O molecules in a periodically replicated orthorhombic MD cell (size 11.9 × 13.2 × 24.2 Å). The system has been reduced by inserting one electron and one proton. The excess electron visualized by the green spin density contours occupies a semilocalized orbital in the TiO2 slab. The proton is attached to the highlighted water molecule. The free energy of the double insertion gives the electron affinity of TiO2 (effectively the conduction band minimum) in contact with water relative to SHE.


Surface acidity and complexation

The surface of metal oxides exchanges protons with aqueous electrolytic solutions and binds also other ions. The excess surface charge is compensated by counterions in an electrical double layer. Surface acidity is therefore an important factor for the understanding of the reactivity. Since the solvation free energy of the proton is also the reference for acidity constants, the proton insertion method can be equally applied to compute the pKa of aqueous species in homogeneous solution as well as surfaces. We have used this method to compute surface pKa's of TiO2 and MnO2 and of main group oxides such as SiO2 (quartz) and Al(OH)3 (Gibbsite). Computational investigation of the properties of electric double layers and the effect on surface reactivity is a major challenge which is adressed in current research.


Electronic Energy Levels and Band Alignment for Aqueous Phenol and Phenolate from First Principles
D Opalka, TA Pham, M Sprik, G Galli
– J Phys Chem B
The temperature dependence of the symmetry factor for a model Fe<sup>3+</sup>(aq)/Fe<sup>2+</sup>(aq) redox half reaction
C Drechsel-Grau, M Sprik
– Molecular Physics
Temperature dependence of interfacial structures and acidity of clay edge surfaces
X Liu, X Lu, J Cheng, M Sprik, R Wang
– Geochimica et Cosmochimica Acta
Interfacial structures and acidity of edge surfaces of ferruginous smectites
X Liu, J Cheng, M Sprik, X Lu, R Wang
– Geochimica et Cosmochimica Acta
Molecular Simulation Study of Hydrated Na-Rectorite
J Zhou, ES Boek, J Zhu, X Lu, M Sprik, H He
– Langmuir
Reductive Hydrogenation of the Aqueous Rutile TiO<inf>2</inf>(110) Surface
J Cheng, X Liu, J VandeVondele, M Sprik
– Electrochimica Acta
Aqueous transition-metal cations as impurities in a wide gap oxide: the Cu(2+)/Cu(+) and Ag(2+)/Ag(+) redox couples revisited.
X Liu, J Cheng, M Sprik
– J Phys Chem B
Frontispiece: Aligning Electronic and Protonic Energy Levels of Proton-Coupled Electron Transfer in Water Oxidation on Aqueous TiO 2
J Cheng, X Liu, JA Kattirtzi, J VandeVondele, M Sprik
– Angew Chem Int Ed Engl
Redox potentials and acidity constants from density functional theory based molecular dynamics
J Cheng, X Liu, J VandeVondele, M Sulpizi, M Sprik
– Acc Chem Res
Aligning electronic and protonic energy levels of proton-coupled electron transfer in water oxidation on aqueous TiO2
J Cheng, X Liu, X Liu, JA Kattirtzi, J VandeVondele, M Sprik
– Angewandte Chemie - International Edition
  • <
  • 2 of 22
  • >

Research Group

Research Interest Groups

Telephone number

01223 336422 (fax)

Email address