<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn> <publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">604799</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2020.604799</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Impact of Solvation on the Structure and Reactivity of the Co<sub>3</sub>O<sub>4</sub> (001)/H<sub>2</sub>O Interface: Insights From Molecular Dynamics Simulations</article-title>
<alt-title alt-title-type="left-running-head">Kox et al.</alt-title>
<alt-title alt-title-type="right-running-head">Co3O4 (001)/H2O Interface</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Kox</surname>
<given-names>T.</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">
<sup>&#x2a;</sup>
</xref>
<uri xlink:href="http://loop.frontiersin.org/people/1112416/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Spohr</surname>
<given-names>E.</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">
<sup>&#x2a;</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Kenmoe</surname>
<given-names>S.</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">
<sup>&#x2a;</sup>
</xref>
<uri xlink:href="http://loop.frontiersin.org/people/1084003/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<label>
<sup>1</sup>
</label>Department of Theoretical Chemistry, University of Duisburg-Essen, <addr-line>Essen</addr-line>, <country>Germany</country>
</aff>
<aff id="aff2">
<label>
<sup>2</sup>
</label>Center of Computational Science and Simulation, University of Duisburg-Essen, <addr-line>Essen</addr-line>, <country>Germany</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/872307/overview">Jun Huang</ext-link>, Central South University, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1095255/overview">Mohammad J. Eslamibidgoli</ext-link>, Forschungszentrum J&#xfc;lich, Germany</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1096546/overview">Zhengda He</ext-link>, Helmholtz Association of German Research Centers (HZ), Germany</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: T. Kox, <email>tim.kox@uni-due.de</email>E. Spohr, <email>eckhard.spohr@uni-due.de</email>S. Kenmoe, <email>stephane.kenmoe@uni-due.de</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Electrochemical Energy Conversion and Storage, a section of the journal Frontiers in Energy Research</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>26</day>
<month>11</month>
<year>2020</year>
</pub-date> <pub-date pub-type="collection">
<year>2020</year>
</pub-date>
<volume>8</volume>
<elocation-id>604799</elocation-id>
<history>
<date date-type="received">
<day>10</day>
<month>09</month>
<year>2020</year>
</date>
<date date-type="accepted">
<day>21</day>
<month>10</month>
<year>2020</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2020 Kox, Spohr and Kenmoe</copyright-statement>
<copyright-holder>Kox, Spohr and Kenmoe</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>The spinel Co<sub>3</sub>O<sub>4</sub> has many beneficial properties for potential use in catalysis. In operando, water is always present and alters the properties of the catalyst. We have used ab initio molecular dynamics to understand the effect of water and solvation on the structure and reactivity of the Co<sub>3</sub>O<sub>4</sub> (001) A-type and B-type surface terminations. Water adsorbs on both terminations via a partial dissociative mode, and the A-termination is seen to be more reactive. On this surface, a higher degree of dissociation is observed in the topmost layers of the crystal in contact with water. Water dissociates more frequently on the Co<sup>2&#x2b;</sup> sites (about 75%) than on the adjacent Co<sup>3&#x2b;</sup> sites, where the degree of dissociation is about 50%. Increasing water coverage does not change the degree of water dissociation significantly. OH<sup>&#x2212;</sup> adsorption on the Co<sup>2&#x2b;</sup> sites leads to a reduction of the amount of reconstruction and relaxation of the surface relative to the clean surface at room temperature. Proton transfer within the water films and between water molecules and surface has localized character. The B-terminated interface is less dynamic, and water forms epitactic layers on top of the Co<sup>3&#x2b;</sup> sites, with a dissociation degree of about 25% in the contact layer.</p>
</abstract>
<kwd-group>
<kwd>solvation</kwd>
<kwd>reactivity</kwd>
<kwd>ab initio molecular dynamics</kwd>
<kwd>surface termination</kwd>
<kwd>coordination number, Cobalt oxide, Spinel</kwd>
</kwd-group>
<counts>
<page-count count="0"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>Recent years have witnessed a rise of interest for heterogeneous catalysis in the liquid phase. Its paramount importance in industrial chemical process technology has prompted intense efforts to gain an atomistic-level understanding of the fundamental process steps. Due to their abundance, the use of transition metal oxides as heterogeneous catalysts is seen as an important low-cost alternative to costly noble metal-based ones. Unlike in the latter case, deep insight into the interplay of factors defining the activity and selectivity of metal oxide catalysts is still missing. More knowledge is still required about the different factors governing their efficiency in order to allow their rational design and application. Several metal oxides have been found relevant for catalysis, and their interactions with water have captured particular attention as water is always present in operando. Co<sub>3</sub>O<sub>4</sub> is one of them, as it shows electronic, magnetic and redox properties that make it a potential candidate for many catalytic reactions. For example, Co<sub>3</sub>O<sub>4</sub> nanoparticles have been successfully used for water and carbon monoxide oxidation, selective 2-propanol oxidation, or steam reforming of ethanol (<xref ref-type="bibr" rid="B10">Llorca et al., 2004</xref>; <xref ref-type="bibr" rid="B6">Deng and T&#xfc;ys&#xfc;z, 2014</xref>; <xref ref-type="bibr" rid="B1">Anke et al., 2019</xref>; <xref ref-type="bibr" rid="B2">Cai et al., 2019</xref>).</p>
<p>Zasada and co-workers combined DFT calculations and high resolution scanning transmission electron microscopy (HRSTEM) to investigate the structure and shape of cobalt spinel nanostructures (<xref ref-type="bibr" rid="B20">Zasada et al., 2011</xref>). Their calculations of surface energies predicted the following order of stability for the low index facets: (001) &#x3e;(111) &#x3e;(110), in agreement with <xref ref-type="bibr" rid="B12">Montoya et al., (2011)</xref> who came to the same conclusion in a separate DFT study. Based on these surfaces energies, a Wulff construction was applied to predict the shape of nanoparticles. They showed rhombicuboctahedral grains exposing mainly the (001) and (111) facets, with a percentage of exposure of 48% and 41%, respectively, and the (101) facet with a minor percentage of 11%. These predictions compared well with HRSTEM images of synthesized spinel nanocrystals. The impact of the vapor pressure of water on the equilibrium shape of the nanoparticles was studied by <xref ref-type="bibr" rid="B19">Zasada et al. (2010)</xref>. Their study reports no faceting transformation upon water adsorption, and the percentage of exposure of the least stable (110) facet is slightly increased to 14%. The presence of water (in the submonolayer and monolayer regime) slightly changes the relative abundance of exposure between the (001) and (111) facets to 39% and 47%, respectively. However, as mentioned in the same study, the (001) facet remains the most attractive to catalyse liquid phase processes. The reason is that on the (111) facet, because of the highly unsaturated nature of the surface Co ions, water molecules do not desorb even at elevated temperatures of more than 450&#xb0;C, which implies that water may block active sites and inhibit catalytic processes taking place at lower temperatures.</p>
<p>To enhance the catalytic response of Co<sub>3</sub>O<sub>4</sub> nano-catalysts, several recipes have been proposed over the years to allow a selective exposure of the more reactive crystal facets for particular catalytic reactions(<xref ref-type="bibr" rid="B15">Sun et al., 2013</xref>). Depending on the recipe, Co<sub>3</sub>O<sub>4</sub> nanoparticles assume the shapes of nanocubes, octahedrons, or nanorods, exposing, respectively, the (001), the (111), or the (101) facets only. These nanoparticle shapes were proven to be beneficial for many catalytical reactions such as hydrodesulfurization of carbonyl sulfide (<xref ref-type="bibr" rid="B17">Wang et al., 2011</xref>) or ethylene oxidation(<xref ref-type="bibr" rid="B11">Ma et al., 2010</xref>). For these reactions, the least-exposed (110) facet was found to be the most reactive one, and, therefore, the reaction was found to be selectively promoted on nanorods. Observations of this nature have motivated studies on the interaction between this particular facet and water (<xref ref-type="bibr" rid="B18">Xu and Li, 2011</xref>; <xref ref-type="bibr" rid="B3">Chen and Selloni, 2012</xref>; <xref ref-type="bibr" rid="B5">Creazzo et al., 2019</xref>), which have provided some atomistic understanding of the water monolayer regime (<xref ref-type="bibr" rid="B3">Chen and Selloni, 2012</xref>) and, very recently, also the inferface with liquid water (<xref ref-type="bibr" rid="B5">Creazzo et al., 2019</xref>), under electrochemical conditions.</p>
<p>Unlike for the (110) surface, atomic-level information for the (001)/H<sub>2</sub>O interface of Co<sub>3</sub>O<sub>4</sub> is scarce, although this facet is reported as being the most stable one and the most relevant for many catalytical reactions. To our knowlegde, only an early attempt by Zasada et al. has been made to understand the water sub- and monolayer adsorption regime on an oxygen-poor Co<sub>3</sub>O<sub>4</sub> surface(<xref ref-type="bibr" rid="B19">Zasada et al., 2010</xref>). Thus, several open issues need to be addressed for the (001) surface. The current paper aims at providing insight into the structure and reactivity of the (001) surface upon water adsorption. We use ab initio molecular dynamics simulations to study the following properties of interfacial water: 1) the active binding sites, 2) the adsorption mode, 3) the degree of dissociation, 4) the fundamental interactions that stabilize the adsorption layer, 5) the structural response of the substrate and 6) the hydrogen bond network. Specifically, we have studied the impact of increasing water content between 16 and 32 water molecules on the A-terminated surface, which contains 16 formal Co<sup>3&#x2b;</sup> and eight formal Co<sup>3&#x2b;</sup> ions, and on the B-terminated surface with 16 formal Co<sup>2&#x2b;</sup> ions, and we compared these simulations with room temperature simulations of the water-free systems.</p>
</sec>
<sec id="s2">
<title>Computational Details</title>
<p>Co<sub>3</sub>O<sub>4</sub> crystallises in a spinel structure with a face-centered cubic unit cell. The primitive unit cell contains 14 atoms (2 Co<sup>2&#x2b;</sup> , 4 Co<sup>3&#x2b;</sup> , and 8 O<sup>2&#x2212;</sup> ). Four such primitive cells form the conventional face-centered cubic unit cell of the antiferromagnetic spinel structure (<xref ref-type="bibr" rid="B4">Chen et al., 2011</xref>), consisting of 56 atoms (8 Co<sup>2&#x2b;</sup> , 16 Co<sup>3&#x2b;</sup> , and 32 O<sup>2&#x2212;</sup> ions). When cutting the bulk in the (001) direction, two (001) surfaces are possible for different termination layers. A surface layer of 16 Co<sup>3&#x2b;</sup> ions forms the B-terminated structure. Slightly above this layer, an additional layer containing 8 Co<sup>2&#x2b;</sup> ions forms the A-termination. <xref ref-type="fig" rid="F1">Figure 1</xref> shows side-views of the two terminations. Note the exposed (green) Co<sup>2&#x2b;</sup> ions on the A-terminated surface (top figure), whereas the Co<sup>3&#x2b;</sup> ions (purple) are, on both ideal surfaces, embedded in the top-most oxygen layer.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption>
<p>A- (top) and B-terminations (bottom) of the Co<sub>3</sub>O<sub>4</sub> (001) surface. Co<sup>2&#x2b;</sup> ions are shown in green, Co<sup>3&#x2b;</sup> in purple, oxygen in red. On the left, various layers are defined, and atom density profiles along the z axis are shown on the right. Layers L6-L10 and L8-L13 were kept rigid in all studies of A- and B-terminated surfaces, respectively.</p>
</caption>
<graphic xlink:href="fenrg-08-604799-g001.tif"/>
</fig>
<p>The Co<sub>3</sub>O<sub>4</sub> (001) surfaces were modeled as slabs of 11 and 13 atomic layers for A- and B-termination, respectively. The choice of non-stoichiometric slabs with an odd number of layers was motivated by the requirement that two equivalent surfaces on top and bottom minimize the dipole moment within the supercell, at least in the water-free case. The atoms in the bottom five or six slab layers were kept immobile at their bulk positions for A- and B-termination (see <xref ref-type="fig" rid="F1">Figure 1</xref>). All other atoms were free to move. Orthorhombic supercells with (2 &#xd7; 2) periodicity in the lateral directions (x,y) and dimensions <inline-formula id="inf1">
<mml:math>
<mml:mrow>
<mml:mn>16.18</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> <inline-formula id="inf2">
<mml:math>
<mml:mo>&#xd7;</mml:mo>
</mml:math>
</inline-formula> <inline-formula id="inf3">
<mml:math>
<mml:mrow>
<mml:mn>16.18</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> <inline-formula id="inf4">
<mml:math>
<mml:mo>&#xd7;</mml:mo>
</mml:math>
</inline-formula> <inline-formula id="inf5">
<mml:math>
<mml:mrow>
<mml:mn>24.9</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf6">
<mml:math>
<mml:mrow>
<mml:mn>16.18</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> <inline-formula id="inf7">
<mml:math>
<mml:mo>&#xd7;</mml:mo>
</mml:math>
</inline-formula> <inline-formula id="inf8">
<mml:math>
<mml:mrow>
<mml:mn>16.18</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> <inline-formula id="inf9">
<mml:math>
<mml:mo>&#xd7;</mml:mo>
</mml:math>
</inline-formula> <inline-formula id="inf10">
<mml:math>
<mml:mrow>
<mml:mn>27.2</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> for the A- and B-terminated surfaces, respectively, were used to study the interaction of water films with the slabs. Water molecules were adsorbed asymmetrically only on the (mobile) top surface and an additional vacuum region of at least 10&#xa0;&#xc5; thickness was introduced, resulting in the box dimensions as given above. All calculations were performed at the <inline-formula id="inf11">
<mml:math>
<mml:mtext>&#x393;</mml:mtext>
</mml:math>
</inline-formula> point. Because of the asymmetry imposed by adsorption and the constrained atoms in the bottommost layers, a dipole correction was applied to cancel the electric field gradient in the vacuum.</p>
<p>Spin polarized Born-Oppenheimer molecular dynamics (MD) simulations were performed using the CP2K/Quickstep package (<xref ref-type="bibr" rid="B16">The CP2K developers group, 2016</xref>). The Generalized Gradient Approximation (GGA) within its PBE formulation (<xref ref-type="bibr" rid="B13">Perdew et al., 1996</xref>) was used to treat the exchange correlation effects. For a more efficient description of dispersion interactions, particularly in the water films, the Grimme D3 correction was added to the PBE exchange-correlation functional (<xref ref-type="bibr" rid="B7">Grimme et al., 2010</xref>). A Hubbard correction (<xref ref-type="bibr" rid="B8">Hubbard and Flowers, 1963</xref>) term (<italic>U</italic>) was added for a correct description of Co 3d states. In a previous study, different <italic>U</italic> values have been used for Co<sub>3</sub>O<sub>4</sub>. These values strongly depend on the computational procedure (<xref ref-type="bibr" rid="B4">Chen et al., 2011</xref>; <xref ref-type="bibr" rid="B12">Montoya et al., 2011</xref>; <xref ref-type="bibr" rid="B3">Chen and Selloni, 2012</xref>; <xref ref-type="bibr" rid="B5">Creazzo et al., 2019</xref>). For our system, the Hubbard correction was set to <inline-formula id="inf12">
<mml:math>
<mml:mrow>
<mml:mi>U</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>&#xa0;eV. With this value, the calculated band gap (1.6&#xa0;eV) agrees well with experiment (<xref ref-type="bibr" rid="B9">Kim and Park, 2013</xref>; <xref ref-type="bibr" rid="B14">Shinde et al., 2006</xref>). The 3&#xa0;s, 3&#xa0;p, 3&#xa0;d, and 4&#xa0;s electrons of Co atoms and the 2&#xa0;s and 2&#xa0;p electrons of O atoms were considered as valence electrons, and the Goedecker-Teter-Hutter (GTH) pseudopotentials were utilized to treat the core electrons. The mixed-repesentation basis sets consist of double-&#x3b6; quality local basis functions with one set of polarization functions (DZVP), together with plane waves with a cutoff of <inline-formula id="inf13">
<mml:math>
<mml:mrow>
<mml:mn>500</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>Ry</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>A Nos&#xe9;-Hoover thermostat with a target temperature of <inline-formula id="inf14">
<mml:math>
<mml:mrow>
<mml:mn>300</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>K</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> and a time constant of <inline-formula id="inf15">
<mml:math>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>ps</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> was used to impose NVT conditions on the system. A time step of 0.5&#xa0;fs was used to propagate the Born-Oppenheimer molecular dynamics trajectories. Total simulation time was 20&#xa0;ps for the water-covered films (after which the system reaches thermal equilibrium, see <bold>Supplementary</bold> <xref ref-type="sec" rid="s8">Figure S1</xref>) and <inline-formula id="inf16">
<mml:math>
<mml:mrow>
<mml:mn>10</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>ps</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> for the water-free systems.</p>
</sec>
<sec sec-type="results|discussion" id="s3">
<title>Results and Discussions</title>
<sec id="s3-1">
<title>Water-Induced Relaxation</title>
<p>We first studied the structural response of the A- and B-terminated substrates to water adsorption. The clean surfaces were taken as references. <xref ref-type="fig" rid="F2">Figure 2</xref> shows each clean surface termination before and after a <inline-formula id="inf17">
<mml:math>
<mml:mrow>
<mml:mn>10</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>ps</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> MD run at room temperature.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption>
<p>Top views of the initial (top row) and final (after 10&#xa0;ps; bottom row) surface structures for A- (left) and B-terminated water-free Co<sub>3</sub>O<sub>4</sub> (001) surfaces, with atom colouring (Co<sup>2&#x2b;</sup>: green, Co<sup>3&#x2b;</sup>: purple, O: red) as in <xref ref-type="fig" rid="F1">Figure 1</xref>. Co<sup>2&#x2b;</sup> ions below the top surface layers are partially grayed out.</p>
</caption>
<graphic xlink:href="fenrg-08-604799-g002.tif"/>
</fig>
</sec>
<sec id="s3-2">
<title>A-Termination</title>
<p>
<xref ref-type="fig" rid="F2">Figure 2</xref> shows that the A-terminated water-free surface undergoes reconstruction at room temperature. The relaxation to the final structure is almost instantaneous at the begin of the simulation. The reconstructed A-terminated surface evolves from the unreconstructed one by (green) Co<sup>2&#x2b;</sup> ions moving along the diagonal rows in the direction to one of the empty sites. The occupation of these four-fold hollow sites has already been reported by <xref ref-type="bibr" rid="B12">Montoya et al. (2011)</xref>. <xref ref-type="fig" rid="F3">Figure 3A</xref> shows the radial distribution function (RDF) of the surface oxygen atoms (those in the L1 layer) to the outermost Co<sup>2&#x2b;</sup> ions in layer L0 (see <xref ref-type="fig" rid="F1">Figure 1</xref>). At 300&#xa0;K, the first maximum of the RDF shifts towards larger distances from the corresponding bulk crystal value of <inline-formula id="inf18">
<mml:math>
<mml:mrow>
<mml:mn>1.95</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> (black curve). The shift is largest for the clean surface (grey), and is smaller, when water is present. In the presence of water, the reconstruction of the A-terminated surface (<xref ref-type="fig" rid="F2">Figure 2B</xref>) is lifted and the final surface positions look similar to the ones in <xref ref-type="fig" rid="F2">Figure 2A</xref> again (not shown), with some minor relaxations. On all A-terminated surfaces, the RDF in <xref ref-type="fig" rid="F3">Figure 3A</xref> shows a diffuse distance spectrum in the range up to about 3&#xa0;&#xc5;, irrespective of actual degree of relaxation and reconstruction, and independent of water coverage. The RDF between surface oxygens and Co<sup>3&#x2b;</sup> ions in layer L1 (<xref ref-type="fig" rid="F3">Figure 3B</xref>) shows a similar behavior of the first peak on the surface relative to the bulk, namely a broadening and large shift in the water-free case and smaller shifts and broadening in the hydrated case, with no significant dependence on water content.</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption>
<p>Radial distribution functions of surface oxygens to outermost (A) Co<sup>2&#x2b;</sup> of the A-termination (layer L0), (B) Co<sup>3&#x2b;</sup> of the A-termination (L1) and (C) the Co<sup>3&#x2b;</sup> of the B-termination (L1). Data are obtained by averaging over the last 10 ps.</p>
</caption>
<graphic xlink:href="fenrg-08-604799-g003.tif"/>
</fig>
<p>The calculated average coordination numbers (defining the cut-off distance as 2.3&#xa0;&#xc5; and averaging over the last half of the simulations, see <bold>Supplementary</bold> <xref ref-type="sec" rid="s8">Figure S2</xref> for illustration) of Co<sup>2&#x2b;</sup> ions in L0 and Co<sup>3&#x2b;</sup> ions in L1 reveal the reason for the lifting of the reconstruction upon water adsorption. Relative to the ideal bulk lattice, coordination in the L0 layer of the water-free surface increases from 2 to 3.5, while the coordination in L1 (Co<sup>3&#x2b;</sup>) is reduced slightly from 5 to 4.9. With the adsorption of 16 water molecules, which mostly (see below) solvate the Co<sup>3&#x2b;</sup> ions in L1, the coordination number in L1 increases to 5.6, at the expense of a slight reduction of the average coordination in L0 from 3.5 to 3.0, indicative also of the lifting of the reconstruction addressed above. Increasing water content to 24 molecules solvates the L0 layer, by increasing the average coordination number of Co<sup>2&#x2b;</sup> to 3.7, while hardly changing the average coordination number in the L1 layer from 5.6 to 5.5. Increasing the water content to 32 water molecules leads to a very minor increase in the coordination number of Co<sup>2&#x2b;</sup> ions in the L0 layer to 3.8, while the coordination number in the L1 layer is stable at 5.5. Taken together with the snapshots in <xref ref-type="fig" rid="F2">Figure 2</xref>, the results indicate that water adsorption leads to an (expected) increase in overall cobalt coordination number, which alleviates the need for reconstruction (driven by an increase in average coordination number), thereby reducing the strain energy of the crystal.</p>
</sec>
<sec id="s3-3">
<title>B-Termination</title>
<p>The free B-terminated surface shows no reconstruction during the 10 ps molecular dynamics runs, as can be observed by comparing <xref ref-type="fig" rid="F2">Figures 2C, D</xref>. Relaxation effects are less prounounced than on the A-termination. The maximum of the RDF for Co<sup>3&#x2b;</sup> ions (layer L1) in <xref ref-type="fig" rid="F3">Figure 3C</xref> shifts towards smaller distance on the free surface and is sharper than for the A-termination, where the shift is in the opposite direction (<xref ref-type="fig" rid="F3">Figure 3B</xref>). Water adsorption reduces this shift relative to the bulk, which can be interpreted as the additional solvation of surface Co<sup>3&#x2b;</sup> ions, which increases their coordination number from 5 towards 6, thereby creating a more bulk-like local oxygen environment for the surface ions.</p>
</sec>
<sec id="s3-4">
<title>The Structure of Interfacial Water</title>
<p>
<xref ref-type="fig" rid="F4">Figure 4</xref> shows top, front and side views of the final snapshots of the simulated water films with 16, 24, and 32 molecules on the A- and B-terminated Co<sub>3</sub>O<sub>4</sub> (001) surfaces. On both terminations, it is seen that, for all coverages, water adsorbs via a partial dissociative mode; the dissociation degree within the contact layer is by and large independent of water content. To gain some statistical insights into the varying properties of the water layers, their interfacial bonds and hydrogen bond network, we analyse density profiles, characteristic radial distribution functions and OH orientations of adsorbed species on the different surface terminations.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption>
<p>Snapshots of final configurations of each simulation. Co<sup>2&#x2b;</sup>: green, Co<sup>3&#x2b;</sup>: purple, slab oxygens: red, water/hydroxide oxygen: blue, hydrogen: white. Greyed-out atoms are located in lower layers.</p>
</caption>
<graphic xlink:href="fenrg-08-604799-g004.tif"/>
</fig>
</sec>
<sec id="s3-5">
<title>A-Termination</title>
<p>On the A-termination, water can bind either to the Co<sup>2&#x2b;</sup> ions in layer L0 or to the Co<sup>3&#x2b;</sup> surface ions in the adjacent L1 layer. Altogether, 32 possible adsorption sites are, thus, present in the 2 &#xd7; 2 cell (2 per Co<sup>2&#x2b;</sup> and one per Co<sup>3&#x2b;</sup> ion). Hence, 16, 24 and 32 adsorbed water molecules correspond, respectively, to 0.5, 0.75 and 1&#xa0;ML coverage. <xref ref-type="fig" rid="F5">Figure 5A</xref> shows the out-of-plane density profiles of water oxygen (thick lines) and hydrogen atoms (thin lines) at different coverages. All oxygen densities have two maxima, one around 2&#xa0;&#xc5; and another one around 2.5&#xa0;&#xc5; from the position of the L1 layer, indicating a buckled water overlayer. At the coverage of 0.5&#xa0;ML the majority of the water molecules are found in the lower sublayer. As the first 16 water molecules have been found to preferentially solvate the Co<sup>3&#x2b;</sup> ions in the L1 layer, one possible explanation for the bimodality might be that the closer sublayer solvates Co<sup>3&#x2b;</sup> and the more distant sublayer solvates Co<sup>2&#x2b;</sup> ions. Integrating the density peak from 0 to 2.3&#xa0;&#xc5; and from 2.3 to 3.0&#xa0;&#xc5; yields on average approximately 12.8 water molecules or hydroxide species in the closer and 3.2 molecules in the more distant sublayer, which does not agree well with the visual inspection of the interface (see also <xref ref-type="fig" rid="F6">Figure 6</xref>), which leads to an estimate of about six molecules interacting with the Co<sup>2&#x2b;</sup> ions in the L0 layer. Thus, it is quite likely that, in addition to the difference between Co<sup>3&#x2b;</sup> and Co<sup>2&#x2b;</sup> solvation, the intrinsic (static and dynamic) inhomogeneity of the L0 layer itself may contribute to the observed oxygen density (see discussion of <xref ref-type="fig" rid="F7">Figure 7</xref>). Increasing water coverage to 0.75&#xa0;ML increases the contribution of the closer subpeak of the oxygen density only slightly to 12.9, whereas most of the additional water molecules are found in the more distant sublayer, in agreement with the progressive solvation of Co<sup>2&#x2b;</sup> ions in L0. Even further increase of water content to 1&#xa0;ML coverage yields 12.6 and 15.7 molecules solvating the closer and more distant layers, respectively. Thus, in summary, these data and coordination numbers suggest that at 1&#xa0;ML coverage most Co<sup>3&#x2b;</sup> ions are solvated by one water molecule or hydroxy group, whereas on average more than 1.5 water molecules or hydroxy species interact with each Co<sup>2&#x2b;</sup> ion directly. The remaining water molecules at the higher water coverages are located in the diffuse range beyond 3&#xa0;&#xc5;.</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption>
<p>Density profiles of oxygen atoms of water and hydroxide species (thick lines) and all hydrogen atoms (thin lines), averaged over the last 10 ps. The center of mass of layer L1 (see <xref ref-type="fig" rid="F1">Figure 1</xref>) is taken as reference.</p>
</caption>
<graphic xlink:href="fenrg-08-604799-g005.tif"/>
</fig>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption>
<p>Trajectory traces of interfacial Co ions (Co<sup>2&#x2b;</sup>: green, Co<sup>3&#x2b;</sup>: purple) and oxygen atoms of hydroxyl groups (red: oxygen originally from water molecule, orange: oxygen original from lattice) and water oxygen atoms (grey). Water oxygen atoms have been shifted upwards (along the y-direction) by 0.5&#xa0;&#xc5; for better discernibility.</p>
</caption>
<graphic xlink:href="fenrg-08-604799-g006.tif"/>
</fig>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption>
<p>Density profiles of the Co<sup>2&#x2b;</sup> atoms for A-terminated surfaces. The center of mass of layer L1 (see <xref ref-type="fig" rid="F1">Figure 1</xref>) is taken as reference position <inline-formula id="inf19">
<mml:math>
<mml:mrow>
<mml:mi>z</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
</caption>
<graphic xlink:href="fenrg-08-604799-g007.tif"/>
</fig>
<p>The hydrogen density profiles show a peak centered around 1 to 1.1&#xa0;&#xc5;, which corresponds to protons transferred from water molecules to the oxide ions. By integrating this peak one can conclude that on average 10.9, 12.0, and 14.0 water molecules transfer a proton to the oxygen atoms in the crystal lattice for 0.5, 0.75, and 1&#xa0;ML coverage, respectively. The second hydrogen peak (clearly visible around 2.2-2.3&#xa0;&#xc5; for the two lower coverages) indicates that the hydroxyl species formed from the water oxygen atoms are oriented on average at a substantial angle to the surface normal with the oxygen atom pointing down, which can be seen from the fact that this peak is further away than the oxygen peak (at around 2&#xa0;&#xc5;). At the highest coverage (red line), the H peak starts rising at approximately the same distance from the reference surface, but the maximum is shifted further away from the surface, probably because the distribution of hydrogen atoms from OH overlaps with the distribution of water hydrogen atoms. Integrating this second hydrogen peak over the range from 1.7 to 2.6&#xa0;&#xc5; yields 15.4, 18.9, and 25.0 hydrogen atoms. This also supports the conclusion that water and OH hydrogens are part of this second peak.</p>
<p>
<xref ref-type="fig" rid="F6">Figure 6</xref> shows trajectory traces of surface cobalt ions (Co<sup>2&#x2b;</sup>: green, Co<sup>3&#x2b;</sup>: purple) and oxygen atoms from hydroxyl groups (red: O originally from water, orange: O originally from oxide) and from undissociated water molecules (grey, shifted upwards by 0.5&#xa0;&#xc5; for better discernibility). First one notes, in agreement with the conclusions drawn before, that there is no visible reconstruction of the A-sites. Most oxygen traces for the A-terminated surface (left column) are localized, indicating a relatively immobile water/hydroxy overlayer. Already at 0.5&#xa0;ML coverage (<xref ref-type="fig" rid="F6">Figure 6</xref>) one notes a substantial amount of dissociated water molecules that transfer a proton to an (orange-colored) oxygen atom in the lattice and produce a second (red) hydroxy group. Undissociated water molecules are visible in grey. In line with the coordination numbers discussed above, many water molecules or hydroxy groups are attached to the Co<sup>3&#x2b;</sup> ions. However, not all of the Co<sup>3&#x2b;</sup> and not all of the Co<sup>2&#x2b;</sup> ions interact with these oxygen species. With increasing water content (<xref ref-type="fig" rid="F6">Figure 6</xref>), all cobalt ions become solvated at 1&#xa0;ML (32 water molecules), with the exception of a single Co<sup>3&#x2b;</sup> ion (at the upper right corner). Also note that there are several traces where red hydroxy and grey water traces are located roughly at the same spot, indicating proton equilibria between OH and H<sub>2</sub>O species. At 1&#xa0;ML coverage, Co<sup>2&#x2b;</sup> ions are not fully solvated from the top but tend to have one or two patches of OH in the off-normal direction, indicating that also these ions are well solvated, often by two OH groups and/or water molecules.</p>
<p>As shown in <xref ref-type="fig" rid="F7">Figure 7</xref>, the structure of the topmost Co<sup>2&#x2b;</sup> ion layer (L0) is not uniform. The density profiles are bimodal with two maxima separated by about 0.6&#xa0;&#xc5;. The bistable nature of the Co<sup>2&#x2b;</sup> density has no obvious origin: While there is a larger tendency for ions solvated by a water molecule to be closer to the surface and a tendency for ions solvated by two hydroxide ions to be further from the surface, all other configurations can also be observed, at least at one of the three coverages. In addition, several ions dynamically switch from an oscillatory motion around one of the configurations to one around the other for extended periods of time. Certainly, the Co<sup>2&#x2b;</sup> ions and their environment form the most dynamic part of the A-terminated interface.</p>
<p>Proton transfer takes place both from water films to the surface and within the water films. In the former case, the dissociated H atoms are transferred to surface oxygens across the neighboring trench which are bound to Co<sup>2&#x2b;</sup> ions underneath. In the latter case, proton transfer is made via Grothuss mechanism between water molecules binding to the surface and the ones bridging two such molecules. This mechanism is more likely to occur with increasing coverage, as more bridging molecules are present on the surface. However, the high dissociation degree in the contact Co<sup>2&#x2b;</sup> layer remains almost the same (see <xref ref-type="fig" rid="F4">Figures 4B,C</xref>).</p>
<p>
<xref ref-type="fig" rid="F8">Figure 8A</xref> shows the radial distribution functions of Co<sup>2&#x2b;</sup> ions in layer L0 at the surface and oxygen atoms of water and hydroxide molecules. At 0.5 ML, the Co<sup>2&#x2b;</sup>-O<sub>w</sub> radial distribution function shows a sharp peak centered at <inline-formula id="inf20">
<mml:math>
<mml:mrow>
<mml:mn>1.85</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> (<xref ref-type="fig" rid="F8">Figure 8A</xref>), as only hydroxide molecules coordinate the Co<sup>2&#x2b;</sup> sites, with short bond lengths. With increasing water coverage, the peak broadens, and its center shifts to <inline-formula id="inf21">
<mml:math>
<mml:mrow>
<mml:mn>2.1</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> due to the coexistence of molecular water in the solvation shell (<xref ref-type="fig" rid="F8">Figure 8A</xref>). On the Co<sup>3&#x2b;</sup> sites below (in layer L1; <xref ref-type="fig" rid="F8">Figure 8B</xref>), both intact water molecules and OH groups are present at 0.5&#xa0;ML in the peak at <inline-formula id="inf22">
<mml:math>
<mml:mrow>
<mml:mn>1.9</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>, which extends up to <inline-formula id="inf23">
<mml:math>
<mml:mrow>
<mml:mn>2.2</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>. As water coverage increases, a second peak, centered at <inline-formula id="inf24">
<mml:math>
<mml:mrow>
<mml:mn>2.4</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>, arises. This can be attributed to additional water molecules that physisorb close to Co<sup>3&#x2b;</sup> sites, since the shape and height of the first Co<sup>3&#x2b;</sup>-O<sub>w</sub> peak barely changes.</p>
<fig id="F8" position="float">
<label>Figure 8</label>
<caption>
<p>Radial distribution function of topmost Co<sup>2&#x2b;</sup>/Co<sup>3&#x2b;</sup> to oxygen atoms of water and hydroxide ions.</p>
</caption>
<graphic xlink:href="fenrg-08-604799-g008.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F9">Figure 9</xref> shows the distribution of the cosine of the angle &#x3b8; between the surface normal and the OH vector of intact water molecules, (top), OH groups formed from water molecules (center) and OH groups formed from the proton transferred from water to an oxygen atom of the surface (bottom). The plane of molecular water which adsorbs at 0.5&#xa0;ML on Co<sup>3&#x2b;</sup> lies mostly parallel to the surface. Increased water coverage increases the number of configurations with protons pointing slightly towards the surface. In addition, some OH groups point almost straight up or straight down towards the surface. The predominantly flat orientation allows for optimum hydrogen bonding at the limited water coverage, even at 1&#xa0;ML coverage.</p>
<fig id="F9" position="float">
<label>Figure 9</label>
<caption>
<p>Orientation of OH vectors of water molecules (top), hydroxide species with oxygen from water molecule (<inline-formula id="inf25">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mtext>O</mml:mtext>
<mml:mtext>W</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>; middle) and hydroxide species formed from the proton transferred from a water molecule to a surface oxygen atom <inline-formula id="inf26">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mtext>O</mml:mtext>
<mml:mtext>S</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> (bottom).</p>
</caption>
<graphic xlink:href="fenrg-08-604799-g009.tif"/>
</fig>
<p>Hydroxide species, which are formed by proton detachment of an added water molecule, show a double maximum in their angular distribution (<xref ref-type="fig" rid="F9">Figure 9B</xref>). One maximum in the range <inline-formula id="inf27">
<mml:math>
<mml:mrow>
<mml:mn>0.8</mml:mn>
<mml:mo>&#x3c;</mml:mo>
<mml:mtext>cos</mml:mtext>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> shows OH groups pointing more or less perpendicularly away from the surface, whereas another set of OH groups points away from the surface at an oblique angle. The OH groups formed from the transferred water proton and a surface oxygen atom only take on the almost perpendicular configurations (<xref ref-type="fig" rid="F9">Figure 9C</xref>), which forms when an oxygen bridge in the surface layers is broken up by protolysis. With increasing water coverage, the average orientation deviates increasingly from the normal direction, probably to facilitate the incorporation of these OH groups into the growing hydrogen bonding network. From the differences in <xref ref-type="fig" rid="F9">Figure 9B</xref> (relative to c) one may arrive at the conclusion that, in addition to OH groups perpendicular to the surface, there is another set of OH that solvates cobalt ions in a more sideways fashion. From an inspection of <xref ref-type="fig" rid="F6">Figure 6</xref> one may conclude that these ions are predominantly Co<sup>2&#x2b;</sup> ions. This is further supported by the fact that, with increasing water content, the relative abundance of the tilted configurations increases, when simultaneously the Co<sup>2&#x2b;</sup> ions become more strongly solvated (see above). The majority of OH groups that point perpendicularly away from the surface bridge, as a detailed analysis shows, two Co<sup>3&#x2b;</sup> sites and are preferentially located close to the sites that do not contain Co<sup>2&#x2b;</sup> ions in layer L0. These bridging OH groups form many of the orange traces in <xref ref-type="fig" rid="F6">Figure 6</xref> that bridge two purple traces at a right angle, where no green trace of Co<sup>2&#x2b;</sup> can be found.</p>
</sec>
<sec id="s3-6">
<title>B-Termination</title>
<p>On the B-termination only Co<sup>3&#x2b;</sup> sites are available for adsorption and 16 potential adsorption sites are available in the 2 &#xd7; 2 supercell. The adsorption of 16, 24, and 32 water molecules thus corresponds to 1, 1.5, and 2&#xa0;ML coverage, respectively. As seen from <xref ref-type="fig" rid="F4">Figures 4D-F</xref>, water films are partially dissociated in the contact layers at all coverages, and the dissociation degree is 25% on average. Proton transfer occurs here also on unsaturated surface oxygen atoms that are not connected to an adjacent Co<sup>2&#x2b;</sup> in the lower (L2) layer, but less frequently than on those connected to the Co<sup>2&#x2b;</sup> in the top (L0) layer in the case of the A-termination. As is also visible in <xref ref-type="fig" rid="F6">Figure 6</xref>, at 1&#xa0;ML the lateral distributions in the (x, y) plane of topmost Co<sup>3&#x2b;</sup> ions (L1) show mostly water molecules in an epitactic layer binding on top of the Co<sup>3&#x2b;</sup> ions, while a few OH groups bridge them in a manner similarly as on the A-terminated surface.</p>
<p>
<xref ref-type="fig" rid="F5">Figure 5B</xref> indicates the formation of an almost flat water layer at 1&#xa0;ML as illustrated by the oxygen density peaks at <inline-formula id="inf28">
<mml:math>
<mml:mrow>
<mml:mn>2.1</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> and the hydrogen peaks at <inline-formula id="inf29">
<mml:math>
<mml:mrow>
<mml:mn>2.3</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>. This feature is still present at 1.5&#xa0;ML; at 2&#xa0;ML the onset of the formation of a diffuse region can be seen as the appearance of tails at larger distances. Transferred protons which are bound to surface oxygen atoms are also present, like on the A- termination, but in lower proportion, as illustrated by the small peaks centered at <inline-formula id="inf30">
<mml:math>
<mml:mrow>
<mml:mn>1.1</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>The radial distribution function of Co<sup>3&#x2b;</sup>-O<sub>w</sub> (<xref ref-type="fig" rid="F8">Figure 8C</xref>) shows a first peak centered around <inline-formula id="inf31">
<mml:math>
<mml:mrow>
<mml:mn>1.9</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> that extends up to <inline-formula id="inf32">
<mml:math>
<mml:mrow>
<mml:mn>2.2</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#xc5;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>, which illustrates the presence of intact water molecules and OH groups that bind via a strong bond. At higher coverage, the width of the peak does not change, which indicates that the rigidity of the contact layer is almost independent of the water content in the range studied here.</p>
<p>Water molecules at small surface coverage are oriented almost parallel to the surface (<xref ref-type="fig" rid="F9">Figure 9D</xref>). With increasing water content additional configurations occur where one OH points towards the surface while the other points at a small angle away from the surface. This behaviour is indicative of a growing hydrogen bond network. Hydroxide ions created from water molecules by proton transfer to the surface are strongly tilted (<xref ref-type="fig" rid="F9">Figure 9E</xref>), with the OH vector almost parallel to the surface at low water coverage. With increasing coverage, the OH direction is more parallel to the surface normal. This behavior may also be associated with the fact that there is substantial proton exchange between OH species and water molecules (see red and grey spots in <xref ref-type="fig" rid="F6">Figure 6</xref>). Most OH groups originating from a surface oxygen atom (<xref ref-type="fig" rid="F9">Figure 9F</xref>), however, point almost upward, similarly to the A-terminated surface, although a few OH groups are also substantially tilted (at <inline-formula id="inf33">
<mml:math>
<mml:mrow>
<mml:mtext>cos</mml:mtext>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2248;</mml:mo>
<mml:mn>0.7</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>). Please note that the distribution functions in <xref ref-type="fig" rid="F9">Figure 9</xref> are not normalized to one, but the area is proportional to the number of OH vectors contributing to the distribution functions.</p>
</sec>
<sec id="s3-7">
<title>Proton Transfer</title>
<p>Water dissociation into OH species and subsequent proton transfer to the surface is a frequent scenario on both surface terminations. Both surfaces are strongly hydroxylated. On the A-terminated surface, on average 5.5 of the 16 water molecules at the 0.5&#xa0;ML coverage remain undissociated, while about 21 OH groups exist on average. With increasing water content the average number of OH groups increases to 23.5 and 28.0 for the simulations at 0.75 and 1.0&#xa0;ML, respectively. The number of water molecules also increases from about 5.5 to 12 to 18 for 0.5, 0.75, and 1&#xa0;ML, respectively. For the B-terminated surface, the number of OH groups increases only moderately from 6.0 to 6.2 and 10.1 and the number of water molecules from 13 to 20.9 to 26.8 (at 1, 1.5 and 2&#xa0;ML, respectively).</p>
</sec>
</sec>
<sec sec-type="conclusion" id="s4">
<title>Conclusion</title>
<p>
<italic>Ab initio</italic> molecular dynamics simulations have shown that water adsorbs on both terminations of the Co<sub>3</sub>O<sub>4</sub> (001) surface via a partially dissociative mode. The dissociation degree on the A-termination in contact with a water layer is significantly higher than on the B-termination.</p>
<p>On the A-terminated surface, water fully dissociates on or near all sites and thereby breaks up bridging oxygen atoms, mostly between Co<sup>3&#x2b;</sup> sites, but also near Co<sup>2&#x2b;</sup>. One can even observe coadsorption of intact and dissociated water molecules on the same Co<sup>2&#x2b;</sup> ion. OH<sup>&#x2212;</sup> adsorption promotes a bulk-like geometry of the reconstructed clean surface. The degree of dissociation on the adjacent Co<sup>3&#x2b;</sup> sites is 50% and remains constant with increasing water coverage. The interface is dynamic with frequent, localised proton transfers between water molecules and hydroxide species via a Grotthuss mechanism. Nevertheless, the oxygen arrangement in the water layer is quite rigid.</p>
<p>On the B-termination, a similar behaviour is observed. However, the interface is less reactive. Here, water forms epitactic layers on top of Co<sup>3&#x2b;</sup> sites and the dissociation degree is less than 25% in the contact layer at all coverages. Proton transfer within the water film and towards the surface is also observed, but with lower probability than on the A-terminated surface.</p>
<p>Both surfaces are strongly hydroxylated with rigid oxygen layers but mobile and disordered proton subsystems. Overall, water solvates the A-terminated surface more strongly. It is also noteworthy that on the A-terminated surface, if water is present, the surface reconstruction observed on the clean surface is released again, due to the ability of water to increase the coordination number of the Co<sup>2&#x2b;</sup> ions substantially <italic>without</italic> moving them along the rows.</p>
<p>In summary, the larger degree of water dissociation on the A-terminated surface and the higher local abundance of OH groups (relative to water molecules) around Co<sup>2&#x2b;</sup> ions as compared to Co<sup>3&#x2b;</sup> ions is consistent with thermal desorption experiments which show high-temperature desorption peaks most probably due to OH groups. These difficult to displace tightly bound OH groups may be at the origin of the observed low oxidation activity on the A-terminated (001) surface.</p>
</sec>
<sec id="s5">
<title>Data Availability Statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s6">
<title>Author Contributions</title>
<p>We have contributed equally to the work: planning, discussions and writting.</p>
</sec>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
</body>
<back>
<ack>
<p>This study was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) &#x2013; 388390466 &#x2013; TRR 247 within the work of Project A6. The authors gratefully acknowledge computing time granted by the Center for Computational Sciences and Simulation (CCSS) of the Universit&#xe4;t of Duisburg-Essen and provided on the supercomputer magnitUDE (DFG grants INST 20876/209-1 FUGG, INST 20876/243-1 FUGG) at the Zentrum f&#xfc;r Informations- und Mediendienste (ZIM). We would also like to thank Rossitza Pentcheva and her group (TRR 247 Project B4) for helpful discussions and optimised oxide structures.</p>
</ack>
<sec id="s8">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="https://www.frontiersin.org/articles/10.3389/fenrg.2020.604799/full&#x23;supplementary-material">https://www.frontiersin.org/articles/10.3389/fenrg.2020.604799/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="datasheet1.pdf" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Anke</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Bendt</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Sinev</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Hajiyani</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Antoni</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zegkinoglou</surname>
<given-names>I.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Selective 2-propanol oxidation over unsupported Co3O4 spinel nanoparticles: mechanistic insights into aerobic oxidation of alcohols</article-title>. <source>ACS Catal.</source> <volume>9</volume>(<issue>7</issue>):<fpage>5974</fpage>&#x2013;<lpage>5985</lpage> <pub-id pub-id-type="doi">10.1063/1.5053729</pub-id>.</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cai</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>J.</given-names>
</name>
</person-group>. (<year>2019</year>). <article-title>Ultrathin, polycrystalline, two-dimensional Co3O4 for low-temperature co oxidation</article-title>. <source>ACS Catal.</source>, <volume>9</volume>(<issue>3</issue>):<fpage>2558</fpage>&#x2013;<lpage>2567</lpage> <pub-id pub-id-type="doi">10.1021/acscatal.8b04064</pub-id>.</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen.</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Selloni</surname>
<given-names>A.</given-names>
</name>
</person-group>. <article-title>Water adsorption and oxidation at the Co3O4 (110) surface</article-title>. (<year>2012</year>). <source>J. Phys. Chem. Lett.</source> <volume>3</volume>(<issue>19</issue>):<fpage>2808</fpage>&#x2013;<lpage>2814</lpage>. <pub-id pub-id-type="doi">10.1016/j.cplett.2010.12.015</pub-id>.</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Selloni</surname>
<given-names>A.</given-names>
</name>
</person-group>. (<year>2011</year>). <article-title>Electronic structure and bonding properties of cobalt oxide in the spinel structure</article-title>. <source>Phys. Rev. B</source>. <volume>83</volume>:<fpage>245204</fpage>. <pub-id pub-id-type="doi">10.1103/physrevlett.77.3865</pub-id>.</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Creazzo</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Ruth Galimberti</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Pezzotti</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Gaigeot</surname>
<given-names>M.-P.</given-names>
</name>
</person-group>. (<year>2019</year>). <article-title>DFT-MD of the (110)-Co3O4 cobalt oxide semiconductor in contact with liquid water, preliminary chemical and physical insights into the electrochemical environment</article-title>. <source>J. Chem. Phys.</source> <volume>150</volume>(<issue>4</issue>):<fpage>041721</fpage>. <pub-id pub-id-type="doi">10.1063/1.5053729</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Deng</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>T&#xfc;ys&#xfc;z</surname>
<given-names>H.</given-names>
</name>
</person-group>. (<year>2014</year>). <article-title>Cobalt-oxide-based materials as water oxidation catalyst: recent progress and challenges</article-title>. <source>ACS Catal.</source> <volume>4</volume>(<issue>10</issue>):<fpage>3701</fpage>&#x2013;<lpage>3714</lpage>. <pub-id pub-id-type="doi">10.1021/acscatal.9b01048</pub-id>.</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Grimme</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Antony</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ehrlich</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Krieg</surname>
<given-names>H.</given-names>
</name>
</person-group>. (<year>2010</year>). <article-title>A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu</article-title>. <source>J. Chem. Phys.</source>, <volume>132</volume>(<issue>15</issue>):<fpage>154104</fpage>.</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hubbard</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Flowers</surname>
<given-names>B. H.</given-names>
</name>
</person-group> <article-title>Electron correlations in narrow energy bands</article-title>. (<year>1963</year>). <source>Proc. Roy. Soc. Lond. Math. Phys. Sci.</source> <volume>276</volume>(<issue>1365</issue>):<fpage>238</fpage>&#x2013;<lpage>257</lpage>. <pub-id pub-id-type="doi">10.1021/jp109264b</pub-id>.</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>K. J.</given-names>
</name>
</person-group> and <person-group person-group-type="author">
<name>
<surname>Park</surname>
<given-names>Y. R.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Optical investigation of charge-transfer transitions in spinel Co3O4</article-title>. <source>Solid State Commun.</source>, <volume>127</volume>(<issue>1</issue>):<fpage>25</fpage>&#x2013;<lpage>28</lpage>. <pub-id pub-id-type="doi">10.1021/jp200581s</pub-id>.</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Llorca</surname>
<given-names>J. P.</given-names>
</name>
<name>
<surname>Dalmon</surname>
<given-names>J.-A.</given-names>
</name>
<name>
<surname>Homs</surname>
<given-names>N.</given-names>
</name>
</person-group>. (<year>2004</year>). <article-title>Transformation of Co3O4 during ethanol steam-re-forming. activation process for hydrogen production</article-title>. <source>Chem. Mater.</source> <volume>16</volume>(<issue>18</issue>):<fpage>3573</fpage>&#x2013;<lpage>3578</lpage>. <pub-id pub-id-type="doi">10.1021/cs500713d</pub-id>.</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ma</surname>
<given-names>C. Y.</given-names>
</name>
<name>
<surname>Mu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J. J.</given-names>
</name>
<name>
<surname>Jin</surname>
<given-names>Y. G.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>G. Q.</given-names>
</name>
<etal/>
</person-group>. (<year>2010</year>). <article-title>Mesoporous Co3O4 and au/Co3O4 catalysts for low-temperature oxidation of trace ethylene</article-title>. <source>J. Am. Chem. Soc.</source>, <volume>132</volume>(<issue>8</issue>):<fpage>2608</fpage>&#x2013;<lpage>2613</lpage>. <pub-id pub-id-type="doi">10.1021/cm049311p</pub-id>.</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Montoya</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Haynes</surname>
<given-names>B. S.</given-names>
</name>
</person-group>. (<year>2011</year>). <article-title>Periodic density functional study of Co3O4 surfaces</article-title>. <source>Chem. Phys. Lett.</source> <volume>502</volume>(<issue>1</issue>):<fpage>63</fpage>&#x2013;<lpage>68</lpage>. <pub-id pub-id-type="doi">10.1063/1.3382344</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Perdew</surname>
<given-names>J. P.</given-names>
</name>
<name>
<surname>Burke</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Ernzerhof</surname>
<given-names>M.</given-names>
</name>
</person-group>. (<year>1996</year>). <article-title>Generalized gradient approximation made simple</article-title>. <source>Phys. Rev. Lett.</source> <volume>77</volume>:<fpage>3865</fpage>&#x2013;<lpage>3868</lpage>. <pub-id pub-id-type="doi">10.1039/c3ta12960h</pub-id>.</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shinde</surname>
<given-names>V. R.</given-names>
</name>
<name>
<surname>Mahadik</surname>
<given-names>S. B.</given-names>
</name>
<name>
<surname>Gujar</surname>
<given-names>T. P.</given-names>
</name>
<name>
<surname>Lokhande</surname>
<given-names>C. D.</given-names>
</name>
</person-group>. (<year>2006</year>). <article-title>Supercapacitive cobalt oxide (Co3O4) thin films by spray pyrolysis</article-title>. <source>Appl. Surf. Sci.</source> <volume>252</volume>(<issue>20</issue>):<fpage>7487</fpage>&#x2013;<lpage>7492</lpage>. <pub-id pub-id-type="doi">10.1016/j.apsusc.2005.09.004</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Ang</surname>
<given-names>H. M.</given-names>
</name>
<name>
<surname>Tad&#xe9;</surname>
<given-names>M. O.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
</person-group>. (<year>2013</year>). <article-title>Co3O4 nanocrystals with predominantly exposed facets: synthesis, environmental and energy applications</article-title>. <source>J. Mater. Chem.</source> <volume>1</volume>:<fpage>14427</fpage>&#x2013;<lpage>14442</lpage>. <pub-id pub-id-type="doi">10.1039/C3TA12960H</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="web">
<collab>The CP2K developers group</collab>. <comment>CP2K is freely available from:</comment> <ext-link ext-link-type="uri" xlink:href="http://www.cp2k.org/">http://www.cp2k.org/</ext-link>, (<year>2016</year>). </citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Ding</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Meng</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Qiu</surname>
<given-names>J.</given-names>
</name>
</person-group>. (<year>2011</year>). <article-title>Novel hydrodesulfurization nano-catalysts derived from Co3O4 nanocrystals with different shapes</article-title>. <source>Catal. Today</source> <volume>175</volume>(<issue>1</issue>):<fpage>509</fpage>&#x2013;<lpage>514</lpage>. <pub-id pub-id-type="doi">10.1016/j.cattod.2011.02.052</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname>
<given-names>X. L.</given-names>
</name>
</person-group>, and <person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>J. Q.</given-names>
</name>
</person-group>. (<year>2011</year>). <article-title>DFT studies on H2O adsorption and its effect on CO oxidation over spinel Co3O4 (110) surface</article-title>. <source>Surf. Sci.</source> <volume>605</volume>(<issue>23</issue>):<fpage>1962</fpage>&#x2013;<lpage>1967</lpage>. <pub-id pub-id-type="doi">10.1021/ja906274t</pub-id>.</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zasada</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Piskorz</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Cristol</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Paul</surname>
<given-names>J.-F.</given-names>
</name>
<name>
<surname>Kotarba</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Sojka</surname>
<given-names>Z.</given-names>
</name>
</person-group>. (<year>2010</year>). <article-title>Periodic density functional theory and atomistic thermodynamic studies of cobalt spinel nanocrystals in wet environment: molecular interpretation of water adsorption equilibria</article-title>. <source>J. Phys. Chem. C</source>. <volume>114</volume>(<issue>50</issue>):<fpage>22245</fpage>&#x2013;<lpage>22253</lpage>. <pub-id pub-id-type="doi">10.1021/jp109264b</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zasada</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Piskorz</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Stelmachowski</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Kotarba</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Paul</surname>
<given-names>J.-F.</given-names>
</name>
<name>
<surname>P&#x142;oci&#x144;ski</surname>
<given-names>T.</given-names>
</name>
<etal/>
</person-group>. (<year>2011</year>). <article-title>Periodic DFT and HR-STEM studies of surface structure and morphology of cobalt spinel nanocrystals. retrieving 3D shapes from 2D images</article-title>. <source>J. Phys. Chem. C</source>. <volume>115</volume>(<issue>14</issue>):<fpage>6423</fpage>&#x2013;<lpage>6432</lpage>. <pub-id pub-id-type="doi">10.1021/jp109264b</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>
