DFT-MD of the (110)-Co3O4 cobalt oxide semiconductor in contact with liquid water, preliminary chemical and physical insights into the electrochemical environment

Water-air interface revisited by means of path-integral ab initio molecular dynamics

Although nuclear quantum effects (NQEs) have been considered on bulk liquid water, the impact of these latter on the air-water interface has not yet been reported. Herein, by performing and comparing ab initio molecular dynamics (AIMD) and path integral AIMD (PI-AIMD) simulations, we reveal the impact of NQEs on structural, dynamical and electronic properties as well as IR spectra of the air-water interface at room temperature. NQEs, being able to describe a more accurate proton delocalization in H-bonded system than AIMD, reveal a different structural arrangement and dynamical behaviour of both bulk and interfacial water molecules in comparison to AIMD results. A more de-structured and de-bound water arrangement and coordination are identified when the quantum nature of nuclei are considered for both bulk and interfacial water molecules. Structural properties, such as inter-/intra-molecular bond lengths, coordination numbers and H-bonding angles of bulk and interfacial water molecules here calculated, are affected by NQEs mitigating the overstructured description given by AIMD. Further evidences of an AIMD overstructured description of bulk water are in the computed IR spectra, where an increased absorption peak intensity and an increased strength of the hydrogen-bond network are alleviated by NQEs. In addition, NQEs show a valuable impact on the electronic structure of the air-water interface, reducing the total valence bandwidth and the electronic energy band-gap when passing from bulk to interfacial water. This work proves how NQEs significantly affect properties and features of the air-water interface, that are essential to accurately describe H-bonded interfacial systems.

Co3O4 (111) surfaces in contact with water: molecular dynamics study of the surface chemistry and structure at room temperature

In this work, we have used ab initio molecular dynamics at room temperature to study the adsorption and dissociation of a thin water film on Co3O4 (111) surfaces, considering the O-rich and Co-rich terminations named as the A-type and B-type surface terminations, respectively. We investigate the occupation of active sites, the hydrogen bond network at the interface and the structural response of the surfaces to water adsorption. On both terminations, water adsorbs via a partial dissociative mode. The contact layer is populated by molecular water as well as OH groups and surface OH resulting from proton transfer to the surface. The B-termination is more reactive, with a higher degree of dissociation in the contact layer with water (46%). On the B-terminated surface, water barely adsorbs on the Co2+ sites and almost exclusively binds and dissociates on the Co3+ sites. The interaction with the surface consists mostly of Co3+-Ow bonds and proton transfer exclusively to the 3-fold unsaturated surface Os1. Hydrogen bonds between water molecules in the aqueous film dominate the hydrogen bond network and no hydrogen bonds between water and the surface are observed. The A-terminated surface is less reactive. Water binds covalently on Co2+ sites, with a dissociation degree of 13%. Proton transfer occurs mostly on the 3-fold unsaturated surface oxygens Os1. Besides, short-lived surface OH arising from proton transfer to 3-fold unsaturated surface oxygens Os2 is observed. H-bonding to surface Os1 and Os2 constitutes 12.7% and 19.8% of the H-bond network, respectively, and the largest contribution is found among the water molecules (67.4%). On both surfaces, the coordination number of the active sites drives the relaxations of the outermost atom positions to the their bulk counterparts. The occupation of active sites on B-termination could reach up to 3 adsorbates per Co3+ leading to a binding motif in which the Co is octahedrally coordinated and which was observed experimentally.

Liquid-Solid Slip on Charged Walls: The Dramatic Impact of Charge Distribution

Nanofluidic systems show great promise for applications in energy conversion, where their performance can be enhanced by nanoscale liquid-solid slip. However, efficiency is also controlled by surface charge, which is known to reduce slip. Combining molecular dynamics simulations and analytical developments, we show the dramatic impact of surface charge distribution on the slip-charge coupling. Homogeneously charged graphene exhibits a very favorable slip-charge relation (rationalized with a new theoretical model correcting some weaknesses of the existing ones), leading to giant electrokinetic energy conversion. In contrast, slip is strongly affected on heterogeneously charged surfaces, due to the viscous drag induced by counterions trapped on the surface. In that case slip should depend on the detailed physical chemistry of the interface controlling the fraction of bound ions. Our numerical results and theoretical models provide new fundamental insight into the molecular mechanisms of liquid-solid slip, and practical guidelines for searching new functional interfaces with optimal energy conversion properties, e.g., for blue energy or waste heat harvesting.

Implications of the fractional charge of hydroxide at the electrochemical interface

Rational design of materials that efficiently convert electrical energy into chemical bonds will ultimately depend on a thorough understanding of the electrochemical interface at the atomic level. Towards this goal, the use of density functional theory (DFT) at the generalized gradient approximation (GGA) level has been applied widely in the past 15 years. In the calculation of electrochemical reaction energetics using GGA-DFT, it is frequently implicitly assumed that ions in the Helmholtz plane have unit charge. However, the ion charge is observed to be fractional near the interface through both a capacitor model and through Bader charge partitioning. In this work, we show that this spurious charge transfer can be effectively mitigated by continuum charging of the electrolyte. We then show that, similar to hydronium, the observed fractional charge of hydroxide is not due to a GGA level self-interaction error, as the partial charge is observed even when using hybrid level exchange-correlation functionals.

Structure and sum-frequency generation spectra of water on uncharged Q4 silica surfaces: a molecular dynamics study

The structural characteristics and sum-frequency generation (SFG) spectra of water near neutral Q₄ silica surfaces were investigated by molecular dynamics simulations. The interactions of water molecules with atoms of the solid were described by different potential models, in particular by the CLAYFF [Cygan et al., J. Phys. Chem. B, 2004, 108, 1255] and INTERFACE [Heinz et al. Langmuir, 2013, 29, 1754] force fields. The calculations of the contact angle of water have shown that the silica surface modeled with CLAYFF behaves as macroscopically hydrophilic, in contrast to the surface described with the INTERFACE model. The hydrophilicity of CLAYFF stems from too attractive electrostatic surface-water interactions. Regardless of the surface's affinity for water, the aqueous phase has a layered structure in the direction perpendicular to the surface with density fluctuations decaying within a distance of 10 Å from the surface. The orientational ordering of H₂O molecules was found to be more short-range than the density fluctuations, especially for the hydrophobic surfaces. Modeling the SFG spectra has shown that the spectra of all studied hydrophobic silica-water interfaces are similar and have features in common with the spectrum of the water-vapor interface. The spectra fairly agree with experimental results obtained for the silica-water interface at low pH conditions [Myalitsin et al., J. Phys. Chem. C, 2016, 120, 9357]. The spectral response for the hydrophobic interface was computed to primarily arise from the topmost molecules of the first layer of interfacial water. In contrast, the SFG signal from the hydrophilic silica-water interface is accumulated over a greater distance extending for several water layers due to more long-range perturbation of the structure by the surface.

Enhanced conductivity of water at the electrified air–water interface: a DFT-MD characterization

DFT-based molecular dynamics simulations of the electrified air–liquid water interface are presented, where a homogeneous field is applied parallel to the surface plane (i.e. parallel to the 2D-HBonded-Network/2DN).

Coupling of Surface Chemistry and Electric Double Layer at TiO2 Electrochemical Interfaces

Surfaces of metal oxides at working conditions are usually electrified because of the acid-base chemistry. The charged interface compensated with counterions forms the so-called electric double layer. The coupling of surface chemistry and the electric double layer is considered to be crucial but is poorly understood because of the lack of information at the atomistic scale. Here, we used the latest development in density functional theory-based finite-field molecular dynamics simulation to investigate the pH dependence of the Helmholtz capacitance at electrified rutile TiO₂(110)-NaCl electrolyte interfaces. It is found that, because of competing forces from surface adsorption and from the electric double layer, water molecules have a stronger structural fluctuation at high pH, and this leads to a much larger capacitance. It is also seen that interfacial proton transfers at low pH increase significantly the capacitance value. These findings elucidate the microscopic origin of the same trend observed in titration experiments.

Special Topic on Interfacial Electrochemistry and Photo(electro)catalysis

Interfacial electrochemistry and photo(electro)catalysis are key processes that convert the energy of photons or electrons to chemical bonds in many energy conversion and storage technologies. Achieving a molecular level understanding of the fundamental interfacial structure, energetics, dynamics, and reaction mechanisms that govern these processes represents a broad frontier for chemical physics and physical chemistry. This Special Topic contains a collection of articles that range from the development of new experimental and computational techniques to the novel application of those techniques for mechanistic studies, as the principal investigators seek a fundamental molecular understanding of both electrode/electrolyte interfaces and the relevant electrocatalytic, photocatalytic, and photoelectrochemical reactions taking place thereabout. Altogether, this collection of articles captures the current state of this very active, frontier research field and highlights the current and remaining key scientific challenges and opportunities.

Graphene-covered transition metal halide molecules as efficient and durable electrocatalysts for oxygen reduction and evolution reactions

Graphene-covered halides are designed as durable and efficient electrocatalysts in acid media. A design principle has been established through the DFT calculations, from which the best catalysts could be predicted for fuel cells.