DFT Study of Azole Corrosion Inhibitors on Cu 2 O Model of Oxidized Copper Surfaces : II . Lateral Interactions and Thermodynamic Stability

Abstract: The adsorption of imidazole, triazole, and tetrazole—used as simple models of azole corrosion inhibitors—on various Cu2O(111)and Cu2O(110)-type surfaces was characterized using density functional theory (DFT) calculations with the focus on lateral intermolecular interactions and the thermodynamic stability of various adsorption structures. To this end, an ab initio thermodynamics approach was used to construct two-dimensional phase diagrams for all three molecules. The impact of van der Waals dispersion interactions on molecular adsorption bonding was also addressed. Lateral intermolecular interactions were found to be the most repulsive for imidazole and the least for tetrazole, for which they are usually even slightly attractive. Both non-dissociative and dissociative adsorption modes were considered and although dissociated molecules bind to surfaces more strongly, none of the considered structures that involve dissociated molecules appear on the phase diagrams. Our results show that the three azole molecules display a strong tendency to preferentially adsorb at reactive coordinatively unsaturated (CUS) Cu surface sites and stabilize them. According to the calculated phase diagrams for Cu2O(111)-type surfaces, the three azole molecules adsorb to specific CUS sites, designated as CuCUS, under all conditions at which molecular adsorption is stable. This tentatively suggests that their corrosion inhibition capability may stem, at least in part, from their ability to passivate reactive surface sites. We further comment on a specific drawback due to neglect of configurational entropy that is usually utilized within the ab initio thermodynamics approach. We analyze the issue for Langmuir and Frumkin adsorption models and show that when configurational entropy is neglected, the ab initio thermodynamics approach is too hasty to predict phase-transition like behavior.


Introduction
This is the second part of the two-part series of articles about the adsorption of three simple azole molecules-imidazole, triazole, tetrazole-on various sites of several Cu 2 O(111)-and Cu 2 O(110)-type surfaces, including Cu and O vacancies, characterized by means of density functional theory (DFT) calculations.In Part I [1] we focused on the azole-surface bonding described at lower coverage without reference to lateral intermolecular interactions-these are addressed in the current Part II-and compared it to the Cl-surface bonding.While the three azole molecules (their molecular structures are shown in Scheme 1) can be seen as models of azole corrosion inhibitors, chloride is a prototypical corrosion activator.As for molecular adsorption, both non-dissociative and dissociative adsorption modes were considered-the latter involves N1-H bond cleavage-and we showed that intact molecules bind to the surface considerably more weakly in comparison to dissociated molecules, which adsorb strong enough to rival the Cl-surface bonds, although even intact molecules can bind rather strongly to specific coordinatively unsaturated (CUS) Cu sites.Indeed, specific coordinatively unsaturated Cu sites bind adsorbates much more strongly than saturated sites, which is precisely what one would expect.Among them, the O vacancy sites, which consist of a triplet of unsaturated Cu ions, were found particularly reactive for the dissociation of the N1-H bond.While non-dissociative adsorption energies were found similar for all three molecules, significant difference between them appeared for the dissociative adsorption mode, i.e., dissociated triazole and tetrazole bind considerably stronger than dissociated imidazole, because the former two can form two strong N-Cu bonds, but imidazole cannot due to its incompatible molecular geometry.Dissociative adsorption is consequently favorable only for triazole and tetrazole, but only at the aforementioned oxygen vacancy sites, where it proceeds without barrier (or nearly so).
imidazole triazole tetrazole In the current paper, we first address the lateral intermolecular interactions as a function of molecular surface coverage and then consider thermodynamic stability of various adsorption structures using the method known as ab initio thermodynamics [2].Due to obvious modeling reasons the adsorption is considered at a solid/vacuum interface in the current two-part series of articles, although in the context of corrosion inhibition it would be more appropriate to consider adsorption at a solid/water interface.We further comment on a specific drawback concerning the way the ab initio thermodynamics approach is usually employed, which is due to neglect of configurational entropy and present an analysis of how this neglect affects the adsorption isotherms in the case of Langmuir and Frumkin adsorption models.
This paper is organized as follows: the next section presents technical details, which include computational method, description of the surface models and surface sites thereon, definition of various adsorption modes as well as the utilized thermodynamic framework.Results are presented in Section 3, whereas in Section 4 we discuss our results and analyze the above-mentioned drawback related to neglect of configurational entropy in approximate ab initio thermodynamics treatment.Finally, Section 5 contains our conclusions.

Computational Details
DFT calculations were performed with the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) [3].The reason for this choice is motivated by the observation that the GGA+U method, which is often used for the description of transition-metal oxides, does not substantially improve the band-gap of Cu 2 O [4,5].Nevertheless, the effect of the +U correction on the adsorption properties of current molecules on Cu 2 O surfaces was evaluated to some extent and the comparison between the PBE and PBE+U methods is provided in Appendix A of Part I [1].
To better describe lateral intermolecular interactions, we also utilized the PBE-D functional, which includes a reparametrized empirical dispersion correction of Grimme [6,7] that consists of a damped C 6 R −6 like energy term on top of PBE.The double prime in the PBE-D label indicates the reparametrization of the original method.The reason for the reparametrization is that the original PBE-D overestimates molecular bonding tocopper surfaces [8][9][10][11], which can be attributed to a too large C 6 value of the Cu atom [10].The PBE-D is reparametrized so as to match the experimental adsorption energy of a flat lying benzene on Cu(111) (for details, see our previous publication [12]).In particular, the C 6 parameter of Cu was set to the value of 140 Ry/Bohr 6 (the original value is 375 Ry/Bohr 6 ).
The pseudopotential method with ultrasoft pseudopotentials was used [13,14].The Kohn-Sham orbitals were expanded in a plane-wave basis set up to a kinetic energy cutoff of 30 Ry (240 Ry for the charge density cutoff).Brillouin zone (BZ) integrations were performed employing the special-point technique [15] using Marzari-Vanderbilt cold smearing [16] of 0.01 Ry.All calculations were done using the PWscf code from the QUANTUM ESPRESSO distribution [17,18].

Cu 2 O as a Model of Oxidized Copper Surface and Description of Considered Surface Models
Cu 2 O slabs without a metal support underneath were used as models of oxidized copper surfaces; the justification of this choice is provided in Part I [1].We use the following surface models (cf. Figure 1; for the more detailed description and graphical representation of considered surfaces, see Figures 1 and 2 in Part I): (i) stoichiometric Cu 2 O(111) -this surface contains two distinct copper ions: a coordinatively saturated (CSA) and coordinatively unsaturated (CUS), labeled as Cu CSA and Cu CUS , respectively.It also contains two distinct oxygen ions: O up and O dn located above (up) and below (dn) the surface Cu layer, respectively.Note that the high-symmetry Cu 2 O(111) is not stable and when the symmetry is broken the Cu CUS ions relax laterally toward adjacent Cu CSA ions [19]; the corresponding structure is designated as Cu 2 O(111)r, where "r" stands for "relaxed".(ii) Cu 2 O(111)-w/o-Cu CUS  Note that among the five models only the first one is stoichiometric, while all the others are non-stoichiometric.In contrast, the label with the postfix specification of the surface supercell, such as Cu 2 O(111) w/o +1Cu CUS -(2 × 2) designates that there is only one Cu CUS ion per specified surface supercell.With this designation the number of Cu CUS ions is not related to the number of adsorbed molecules.

Adsorption Calculations
Adsorption calculations on Cu 2 O(111)-type surfaces were modeled with slabs consisting of four O-Cu-O trilayers, whereas Cu 2 O(110):CuO was modeled with slabs consisting of five layers, that is, two CuO-Cu bilayers plus an extra CuO layer at the bottom of the slab.For Cu 2 O(111)-type surfaces the bottom trilayer and for Cu 2 O(110):CuO the bottom bilayer were constrained to bulk positions, while all other degrees of freedom were relaxed.Molecules were adsorbed on the top side of the slab and the thickness of the vacuum region-the distance between the top of the ad-molecule and the adjacent slab-was set to about 20 Å.The dipole correction of Bengtsson [20] was applied to cancel an artificial electric field that develops along the direction normal to the slab due to periodic boundary conditions imposed on the electrostatic potential.In contrast to adsorption calculations, surface energies of clean surfaces were estimated with thicker and fully relaxed symmetric slabs as described below in Section 2.4.
Adsorption properties were calculated with a variety of different supercells so as to model adsorption at different coverages and we used the following uniformly shifted k-meshes for the BZ integrations:

Definition of Surface Coverage
Given the size of the current molecules (and assuming adsorption onto discrete sites), one molecule per (1 × 1)-Cu 2 O(111) or (1 × 1)-Cu 2 O(110):CuO can be reasonably considered as the largest viable coverage before the lateral intermolecular Pauli repulsion sets in.A relative coverage of 1 monolayer (ML) is therefore defined as one molecule per (1 × 1)-Cu 2 O(111) or (1 × 1)-Cu 2 O(110):CuO unit cell; relative coverage will be also designated in percents as, e.g., 1 ML ≡ r100%, where r is used as a mnemonic for "relative".

Non-Dissociative, Dissociative, Mixed-Site, and Mixed-Mode Adsorption
We consider both non-dissociative and dissociative adsorption of azole molecules.Non-dissociative adsorption means that no molecular bonds are broken upon adsorption, i.e., molecule adsorbs in intact form, whereas dissociative adsorption involves the N1-H bond cleavage upon adsorption (the N1 atom is indicated in Scheme 1).
The term "mixed-site" adsorption implies that in a given adsorption structure the molecules are adsorbed at different types of surface sites, whereas the term "mixed-mode" adsorption designates that some molecules are adsorbed non-dissociatively and others dissociatively.In the current study, mixed-mode adsorption is typically also mixed-site adsorption, because dissociated and intact molecules adsorb to different sites.Hence the term mixed-mode adsorption will also implicitly imply mixed-site adsorption, whereas the term mixed-site adsorption will implicitly pertain to non-dissociative adsorption modes, unless stated otherwise.

Adsorption Equations
The reaction for non-dissociative adsorption of azole molecules can be written as: whereas dissociative adsorption, which involves N1-H bond cleavage (cf.Scheme 1), can be described as: where the label MolH stands for an isolated intact molecule, standalone * designates a free adsorption site, while MolH * , Mol * , and H * denote adsorbed species.Analogously to the MolH labeling, the specific labels for intact imidazole, triazole, and tetrazole are ImiH, TriH, and TetH, whereas dissociated molecules are designated by Imi, Tri, and Tet labels, respectively.The respective non-dissociative and dissociative adsorption energies (per adsorbed molecule) are calculated as: and where m and n are the numbers of non-dissociatively and dissociatively adsorbed molecules per supercell, respectively.E MolH , E slab , E mMolH/slab , and E n(Mol+H)/slab are the total energies of isolated intact MolH molecule, clean slab, slab with m adsorbed molecules per supercell, and slab with n coadsorbed Mol+H species per supercell, respectively.Mixed-mode adsorption is described as: and the respective average adsorption energy per molecule is calculated as: To estimate how strongly a given species binds to the surface, we will utilize the adsorption bond strength (D), which by convention is positive and hence opposite to the binding energy (E b ); the two quantities are calculated as: where A stands for adsorbate (MolH or Mol) and the meaning of the energy terms is analogous to those defined above.By this definition the E b of Mol is calculated with respect to the isolated radical (Mol • ) and not the anion (Mol − ); isolated radicals were calculated with spin-polarized calculations.Note that for non-dissociative adsorption at low coverage the binding energy is equivalent to the adsorption energy, i.e., The relative stability of molecular vs. dissociative adsorption is evaluated by considering the dissociation reaction on the surface, MolH * + * −→ Mol * + H * , and the respective dissociation energy is calculated as (for simplicity of notation we consider one molecule per supercell): and dissociative adsorption is favored over the non-dissociative adsorption when ∆E < 0.

Thermodynamic Stability
The thermodynamic stability of molecule-surface systems is evaluated by following the treatment from our previous publication [12].In particular, adsorption surface free energy ( γads ) is considered to be a two-dimensional function of chemical potentials of oxygen (µ O ) and azole molecule (µ MolH ): where γ surf (µ O ) is the surface free energy of the pristine surface and it is a function of the oxygen chemical potential, while γ ads (µ MolH ) is the term due to molecular adsorption and it is a function of the molecular chemical potential.Note that the dependence of γads on the temperature and partial pressures of oxygen and azole molecules is implicitly taken into account by µ O and µ MolH parameters; the mapping of µ O and µ MolH into temperature and partial pressures is described in Appendix A.
To a first approximation, the term due to molecular adsorption (i.e., the adsorption free energy per unit area) can be written as: where A is the area of the supercell and n is the number of adsorbed molecules per supercell (for mixed-mode adsorption n ≡ m + n, cf.Equations ( 5) and ( 6 We should remark that the current use of the symbol ε ads is different from the one used in our previous publication [12], i.e., currently ε ads designates adsorption energy per molecule, while in ref [12] it designated total adsorption energy per unit area. In contrast to molecular adsorption, surface free energies were calculated with symmetric slabs where both surfaces are equivalent.For the (111)-type structures slabs consisting of 6, 7, and 9 O-Cu-O trilayers were used, whereas for the Cu 2 O(110):CuO surface, slabs consisting of 4, 5, and 6 CuO-Cu bilayers plus a terminating CuO layer were used.The surface free energy as a function of oxygen chemical potential is evaluated as: is the number of excess O atoms in the slab and accounts for its non-stochiometry; it is given by: where N O and N Cu are the numbers of O and Cu ions in the slab (for stoichiometric slabs ∆N stoich O = 0).We consider oxygen and copper chemical potentials to be interdependent via the relation: where E bulk Cu 2 O is the total energy of the Cu 2 O bulk per formula unit.Note that µ Cu and µ O cannot physically vary without bounds (e.g., if the µ O becomes too low the Cu 2 O would decompose into Cu-bulk and oxygen gas), hence the range of viable µ O is defined to be between oxygen-lean (O lean ) and oxygen-rich (O rich ) limits, which are chosen according to Equations (16a) and (16b), respectively: where E bulk Cu is the total energy of Cu atom in the Cu-bulk and E O 2 is the total energy of isolated O 2 molecule, respectively.Half of the total energy of the O 2 molecule is chosen as the zero reference for µ O , i.e., Hence at the O lean limit ∆µ O = −1.27eV (calculated value) and at the O rich limit ∆µ O = 0 eV.

The Thermodynamic Stability of Bare Cu 2 O Surfaces
The stability of the currently considered bare Cu 2 O surfaces is evaluated in Figure 2 by means of surface free energies as a function of oxygen chemical potential, γ surf (µ O ) of Equation ( 13 • reconstructed surface (the difference between the two models is that model B lacks Cu CUS ions) is currently not considered for the following two reasons: (i) according to DFT calculations the surface free energy of model A is considerably higher than that of model B (see Figure A1 in the Appendix B and also ref [24]) and (ii) no new types of Cu and O ions are introduced in model A. Finally, let us comment on the stability of Cu 2 O(111) w/o +1Cu CUS -(N × N) surface models, which lack all but a single Cu CUS ion per supercell, that is, they lack (N 2 − 1) Cu CUS ions.These models can be seen as a linear combination of Cu 2 O(111)r and Cu 2 O(111)-w/o-Cu CUS , hence their surface free energies should be in between the values displayed by the latter two.Given that the considered Cu 2 O(111) w/o +1Cu CUS -(2 × 2) and Cu 2 O(111) w/o +1Cu CUS -(3 × 3) models contain only 25% and 11% of the Cu CUS ions, respectively, their surface free energies should be closer to Cu 2 O(111)-w/o-Cu CUS than to Cu 2 O(111)r, which is indeed the case as seen in Figure 2; for the comparison between the estimated (via linear combination) and actually calculated surface free energies of the two Cu 2 O(111) w/o +1Cu CUS models, see Figure A2 in the Appendix B.
Considering the surface free energies, shown in Figure 2, we should make the following comment.While our results are in good agreement with previous GGA studies, Nilius et al. [25] recently argued that the GGA predicted cost of non-stoichiometry is too small and correspondingly GGA overestimates the stability of non-stoichiometric surfaces.Using the HSE hybrid functional they showed, in agreement with their STM experiments, that stoichiometric Cu 2 O(111) becomes more stable than Cu 2 O(111)-w/o-Cu CUS under oxygen-lean conditions.

Molecular Adsorption at Lower Coverage
The top row of Figure 3 summarizes the results presented in Part I [1], that is, adsorption bond-strengths for MolH * and Mol * adsorption modes of imidazole, triazole, and tetrazole at various considered Cu sites on Cu 2 O surfaces and the respective reaction energies for MolH * + * → Mol * + H * dissociation.These results were obtained for a coverage, which is low enough that it allows for focussing predominantly on the molecule-surface interactions (see Figure 3 in Part I); currently this coverage can be seen as intermediate, because herein we consider also lower and higher coverages.While in Part I, we focused on the molecule-surface interactions and presented only the PBE results, currently we mainly focus on the lateral intermolecular interactions, hence we utilize a dispersion correction to better describe them (PBE-D results, shown in the bottom row of Figure 3).The comparison between PBE and PBE-D results presented in Figure 3 shows that the two sets of results give the same trends; the main difference is that due to attractive dispersion interactions the PBE-D adsorption bond-strengths are by about 0.3 to 0.5 eV stronger than the PBE ones.@Cu CSA @Cu CUS @Cu surf @Cu Ovac @Cu CSA @Cu CUS @Cu surf @Cu Ovac @Cu CSA @Cu CUS @Cu surf @Cu Ovac @Cu CSA @Cu CUS @Cu surf @Cu Ovac  bind adsorbates more strongly than saturated Cu sites (Cu CSA and Cu surf ), (iii) on a given site non-dissociative adsorption energies are similar for all three molecules, while dissociated triazole and tetrazole bind much stronger than dissociated imidazole, hence (iv) molecular dissociation is favored just for triazole and tetrazole, but only at Cu Ovac oxygen vacancy sites.Note that after dissociation H binds to a nearby lattice O ion thus forming an OH group.Figure 3c therefore reveals that non-dissociative adsorption is usually favored over dissociative adsorption; the only two exceptions are triazole and tetrazole adsorbed at Cu Ovac on Cu 2 O(111)-recon-( • , for which dissociative adsorption is favored.For benzotriazole, the dissociative adsorption was also found favourable only at oxygen vacancy sites [26] and furthermore oxygen vacancy sites were found reactive for the cleavage of the O-H bond of methanol [27].

Lateral and Dispersion Interactions: PBE vs. PBE-D
Azoles are polar molecules with rather large permanent dipole moment-PBE calculated values are 3.79, 4.42, and 5.45 D for imidazole, triazole, and tetrazole, respectively [28]-which can result in rather long ranged dipole-dipole intermolecular interactions between adsorbed molecules.Depending on the orientation of the dipole moment of adsorbed molecules, the interactions can be either attractive (for dipoles pointing parallel to the surface) or repulsive (for perpendicular dipoles) [12,28,29].In our previous publication [12], we found that lateral interactions are repulsive for imidazole adsorbed at Cu CUS or Cu CSA sites on Cu 2 O(111), whereas for triazole and tetrazole they are either less pronounced or even attractive; these dependencies are consistent with the orientation of molecular dipoles in the adsorbed states [12]., for which dissociative adsorption is the stablest (cf. Figure 3c).However, when all the Cu Ovac sites are occupied on Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • , which happens at a relative coverage of 1/3 ML, the molecules start to adsorb non-dissociatively; this implies a mixed-mode adsorption for coverages larger than 1/3 ML.The transition from dissociative to mixed-mode adsorption is seen by the large upward "jump" in adsorption energy as going 1/3 to 2/3 ML coverage; coverages at which the mixed-mode adsorption takes place are indicated by a yellowish band.Bottom row: the difference between PBE-D and PBE adsorption energies.
On Cu 2 O(111) based models the adsorbed molecules form centered hexagonal patterns (see Figure 3 in Part I), where each molecule has six neighbors at R nn .However, such a high symmetry pattern is not possible for mixed-site and mixed-mode adsorption on Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • and for Cu 2 O(110):CuO, hence the lateral dependence is plotted against the effective nearest-neighbor intermolecular distance, R eff nn , which is defined as: where A mol is the surface area per molecule, i.e., A mol = A/n, where A is the area of the supercell and n is the number of adsorbed molecules per supercell.This definition stems from the relation between the R nn and A mol on Cu 2 O(111) based surface models, where A mol = √ 3R 2 nn /2. Figure 4 reveals that-similarly to the above described lateral dependencies on Cu CSA and Cu CUS sites of Cu 2 O(111)-based models-imidazole displays repulsive interactions also on Cu 2 O(110):CuO, while the lateral interactions for triazole and tetrazole are less pronounced thereon.On Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • , however, the situation is more complicated due to the presence of 1/3 ML of oxygen vacancies that result in Cu Ovac sites and for triazole and tetrazole also due to the mixed-mode adsorption.Notice the pronounced upward shift of the average adsorption energy when the coverage becomes higher than 1/3 ML (or R eff nn < 20 Bohrs; region marked with a yellowish band).This upshift is, however, to a large extent due to mixed-mode adsorption rather than the repulsive intermolecular interactions.Namely, triazole and tetrazole prefer to adsorb dissociatively onto Cu Ovac sites, but at a relative coverage of 1/3 ML all Cu Ovac sites are occupied.Beyond this coverage the molecules continue to adsorb non-dissociatively onto already singly occupied Cu Ovac sites, though in hindered geometry, up to the coverage of 2/3 ML after which imidazole adsorbs to uncovered parts of the surface via the N1-H• • • O up hydrogen bond, whereas triazole and tetrazole at 1 ML coverage prefer to singly occupy the Cu Ovac sites and the remaining molecules adsorb via N1-H• • • O up hydrogen bonds.
As for the role of dispersion interactions, it can be observed from Figure 4 that dispersion enhances the molecule-surface interactions by about 0.3 to 0.5 eV, depending on the specific case (for the PBE-D vs. PBE differences see the bottom row plots of Figure 4).Dispersion interactions also stabilize high coverage configurations (in relative sense) due to attractive intermolecular dispersion interactions.This effect is the most clearly seen on Cu 2 O(110):CuO, where the largest number of structures is considered: notice that the PBE-D vs. PBE difference changes from −0.3 to −0.4 eV as going from R eff nn of 16 Bohrs down to about 10 Bohrs.We should comment on the anomalous PBE-D vs. PBE difference for tetrazole on Cu 2 O(111)-w/o-Cu CUS , where the PBE-D apparently destabilizes the high coverage (1 × 1) phase for tetrazole: the reason is that in this specific case PBE and PBE-D predict different tetrazole configurations to be the most stable and the two configurations display different lateral dependence.If instead the same configurations are considered for both functionals, then the expected behavior is recovered (see Figure A3 in the Appendix B).The PBE-D vs. PBE difference also appears anomalous on Cu 2 O(111)-recon-( The reason is that on this surface various sites are occupied at high coverage and the PBE-D vs. PBE difference reflects dispersion contributions not only to the lateral intermolecular interactions but also to the molecule-surface interactions.
The corollary of the results presented in Figure 4 is that dispersion interactions stabilize adsorption structures and that imidazole displays the most repulsive lateral interactions among the three considered molecules; note that for triazole or tetrazole on Cu 2 O(111)-recon-( • the large reduction of average adsorption energy magnitude at higher coverages is due to the mixed-mode adsorption rather than intermolecular repulsion.

Thermodynamic Stability
The fact that we evaluated on the order of a few hundred different adsorption structures implies that their thermodynamic (TD) stability has to be considered in a systematic way.Hence, we will first consider adsorption free energies as a one-dimensional function of molecular chemical potential, γ ads (µ MolH ) of Equation (11).In this way we can characterize the TD stability of different coverages on a given surface; the condition is that the surface itself is kept the same for all the coverages.This poses a problem for Cu 2 O(111) w/o +1Cu CUS model, for which the fraction of Cu CUS ions changes with molecular coverage, i.e., the fraction of Cu CUS ions is the same as the molecular coverage.Hence, we characterize instead the TD stability of adsorption structures on two "fixed" surface models: Cu 2 O(111) and Cu 2 O(111) w/o +1Cu CUS -(2 × 2).The first contains 100% and the latter 25% of Cu CUS ions, irrespective of the molecular coverage.On the former model the molecules adsorb strictly at Cu CUS sites, whereas on the latter model the molecules adsorb also at Cu CSA sites after all the Cu CUS sites are occupied, which happens at molecular coverage of r25%.
PBE and PBE-D calculated adsorption surface free energies as functions of the molecular chemical potential are shown in Appendix B: Figure A4 The most stable structure of a given molecule on a given surface at a given molecular chemical potential is the one that displays the lowest γ ads .Figures A4 and A5 thus reveal that high coverage r100% configurations usually dominate, i.e., they are the stablest over the broadest range of the molecular chemical potential.There are only two exceptions for which lower coverage phases are the most stable over a broad range of low molecular chemical potential, i.e., (i) the r33% phases of triazole and tetrazole on Cu 2 O(111)-recon-( • and (ii) the r25% phases of all three molecules on Cu 2 O(111) w/o +1Cu CUS -(2 × 2).According to Figures A4 and A5 the largest viable stabilization of surface free energies due to molecular adsorption corresponds to high coverage r100% phases.This reasoning is corroborated by Figure 5, which plots the total adsorption energy per unit area, nε ads /A, for the three molecules on the five considered surfaces.Notice that for all the cases the nε ads /A monotonically decreases as coverage increases, thus being the most exothermic for the r100% phases.Current results reveal that the molecules bind considerably more strongly to coordinatively unsaturated Cu sites (cf.Figures 3 and 4), which raises the question whether stronger molecular bonding can stabilize them.This issue can be addressed by considering the stabilization of surface free energies due to molecular adsorption, which according to Figure 5 is the strongest for r100% phases, hence Figure 6 compares the surface free energies of pristine surfaces to those that are fully covered with imidazole, triazole, or tetrazole; for this analysis the molecular chemical potential was fixed at ∆µ MolH = 0.It can be seen that molecular adsorption alters the surface stability trend, because it stabilizes the stoichiometric Cu 2 O(111) surface to such an extent that it becomes the stablest under oxygen-lean conditions (prior to molecular adsorption the Cu 2 O(111)-w/o-Cu CUS is the most stable).In contrast, under oxygen-rich conditions the Cu 2 O(110):CuO is the most stable prior and after molecular adsorption.Figure 6 further reveals that molecular adsorption does not stabilize sufficiently the Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • surface (this surface is never the stablest), despite the fact that the Cu Ovac sites are very reactive toward molecular adsorption and are able to dissociate triazole and tetrazole by cleaving the N1-H bond.This issue will be further commented in the discussion Section 4.1.

Two-Dimensional Phase Diagrams
The above one-dimensional treatments (cf.Figures 6, A4 and A5) can be extended by considering the adsorption surface free energy as a two-dimensional (2D) function of µ MolH and µ O via Equation (10).Such a 2D treatment has already been considered in our previous publication [12], but only for Cu   8.In these figures the dependence of each chemical potential was also recast into a temperature scale at partial pressure p = 1 atm and into a pressure scale at T = 300 K.Note that each azole molecule should have its own p and T scales, but the scales for the three molecules are so similar that they are visually indistinguishable.We first comment on the 2D phase diagrams for (111)-type surfaces only (Figure 7).Apart from one exception (i.e., high coverage r100% phase of triazole on Cu 2 O(111)-w/o-Cu CUS that appears only on PBE-D calculated phase-diagram at the oxygen-rich limit), only two high coverage molecular phases appear on the top part of the phase diagrams: r100% @ Cu 2 O(111) on the left (oxygen-lean conditions) and r100% @ Cu 2 O(111) w/o +1Cu CUS -(2 × 2) on the right (oxygen-rich conditions).Below these high coverage molecular phases, two diagonal stripe like regions appear: the r25% @ Cu 2 O(111) w/o +1Cu CUS -(2 × 2) and beneath it the r11% @ Cu 2 O(111) w/o +1Cu CUS -(3 × 3).The domains of these two phases are the largest for imidazole and the smallest for tetrazole.In some cases very small regions of r50% and r75% @ Cu 2 O(111) w/o +1Cu CUS -(2 × 2) phases also appear in between the respective r25% and r100% phases.Finally, for sufficiently low molecular chemical potentials the bare Cu 2 O(111)-w/o-Cu CUS surface is the stablest (bottom part of the phase diagrams).
Comparing the current 2D phase diagrams of only Cu 2 O(111)-type surfaces (Figure 7) to those published in our previous publication [12], we note that apart from one exception none of the molecular phases on Cu 2 O(111)-w/o-Cu CUS appear on the current phase diagrams.Instead the high coverage r100% phase on Cu 2 O(111) w/o +1Cu CUS -(2 × 2) appears on the top-right part of the phase diagrams.The reason for the discrepancy between the current 2D phase-diagrams and those previously published [12] is that we did not consider the Cu 2 O(111) w/o +1Cu CUS -(N × N) surfaces in our previous publication.
When both (111)-and (110)-type structures are considered together (Figure 8) then Cu 2 O(111)-type structures remain on the left part of phase diagrams (oxygen-lean conditions), whereas Cu 2 O(110):CuO structures appear on the right (oxygen-rich conditions).As for the latter, high coverage r100% molecular-phase appears at the top-right and bare Cu 2 O(110):CuO surface at the bottom-right.The transition between the two-i.e., from high coverage to bare surface-is the broadest and the most continuous-like for imidazole.Namely, one passes from high coverage to lower and lower coverages as ∆µ MolH decreases: note the various blue band like regions on the PBE phase diagram indicating the r50%, r33%, r17%, r11%, and r8% phases.For triazole this transition is sharper and for tetrazole the transition from r100% to bare Cu 2 O(110):CuO is sharp on the PBE-D diagram, whereas on the PBE diagram there is only a thin dark blue band region corresponding to r11% phase.This trend is due to the intermolecular lateral-interactions, which are the most repulsive for imidazole and slightly attractive for tetrazole (cf. Figure 4).
The comparison between the PBE and PBE-D calculated phase diagrams shows that the effect of dispersion interactions is mainly twofold: (i) the boundary between molecularly covered and bare surfaces shifts to lower values of ∆µ MolH due to dispersion enhanced adsorption bonding; (ii) the transition from bare Cu 2 O(110):CuO surface to high coverage r100% molecular phase thereon is sharper due to attractive dispersion interactions between molecules at higher coverage (by "sharper" we mean that intermediate coverage blue bands are either narrower (imidazole and triazole) or even disappear (tetrazole) on phase diagrams).

On the Role of Cu CUS and Cu Ovac Sites
According to the phase diagrams for Cu 2 O(111)-type surfaces (Figure 7), the three azole molecules are able to stabilize the Cu CUS sites under all conditions at which molecular adsorption is stable (the only exception appears for triazole on the PBE-D phase diagram at the oxygen-rich limit, where the r100% @ Cu 2 O(111)-w/o-Cu CUS is the stablest).Cu CUS sites stabilized by molecular adsorption are particularly interesting in view of the observations made recently by Nilius et al. [25], who drew attention to the well-known drawback of GGA related to the underestimation of band-gaps and incorrect predictions of absolute positions of band-edges.Consequently, they showed that GGA too strongly favors non-stoichiometric Cu 2 O(111)-w/o-Cu CUS against Cu 2 O(111).However, despite the fact that GGA overestimates the instability of Cu CUS sites, it predicts that azole molecules are able to stabilize them (note that currently utilized PBE and PBE-D are GGA).In this respect we can anticipate two corollaries: (i) GGA may overestimate molecular bonding to Cu CUS sites, which leads to error cancellation, i.e., overestimated Cu CUS instability cancels with overestimated molecular bonding.(ii) GGA may reasonably predict the molecular bonding to Cu CUS sites; in this case the molecular bonding to Cu CUS sites would dominate to an even larger extent as currently predicted.It is worth noting that the comparison between GGA and GGA+U calculations (see Appendix A in Part I [1]) seems to indicate that GGA does not overestimate molecular bonding to Cu CUS sites, because GGA+U predicts slightly stronger bonding.In both aforementioned cases, it is clear that molecular adsorption stabilizes the thermodynamic deficiency of coordinatively unsaturated Cu CUS ions.Furthermore, azole molecules also bind more strongly to undercoordinated defects on metallic Cu surfaces [30][31][32].Apparently, azoles display a strong tendency to preferentially adsorb at reactive undercoordinated or unsaturated surface sites and to stabilize them, which in turn tentatively suggests that their corrosion inhibition capability may, at least in part, stem from their ability to passivate reactive surface sites.
The observation that the Cu CUS sites are stabilized by molecular adsorption is relevant in light of the fact that adsorption data, obtained from combined experimental-computational investigations [27,33], were recently used to evaluate various surface models of Cu 2 O surfaces.While such a procedure is elegant and appropriate for inferring the structure of surfaces covered with molecules, current results clearly show that it may not always be possible to correctly deduce the stability of bare surfaces on the basis of adsorption data.
Another notable observation from Figure 7 is that the Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • surface is not sufficiently stabilized by the strong bonding of dissociated molecules to Cu Ovac sites thereon, because no Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • structure appears on the phase diagrams.Nevertheless, oxygen vacancy Cu Ovac sites were observed experimentally [23].Furthermore, it is in fact very difficult to prepare surfaces with no oxygen vacancies in vacuum [34].Here several issues should be mentioned.First, it is obvious that oxygen vacancies and other surface defects are (to some extent) stabilized by configurational entropy, which is currently neglected.Secondly, oxygen vacancies are currently considered only at relatively large concentration of 33% using the Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • model; according to calculations, the surface with such a large concentration of vacancies is too unstable (cf. Figure 2), hence it seems that a lower concentration of vacancies would be more realistic.Thirdly, the comparison between GGA and GGA+U calculations (see Appendix A in Part I [1]) indicates that GGA probably underestimates the dissociation tendency and correspondingly also the stabilization of the Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • surface.It can be therefore reasonably anticipated that oxygen vacancies may be stabilized at lower concentration due to both configurational entropy effects and strong molecular adsorption.Further support to the presence of oxygen vacancy sites on Cu 2 O(111) surfaces was recently provided by means of combined experimental-computational studies of methanol dehydrogenation and adsorption of SO 2 [27,33].Finally, one also needs to consider non-equilibrium and kinetic aspects, which may be very important for corrosion and in particular for corrosion inhibition.Under such circumstances the reactive oxygen vacancy Cu Ovac sites may be passivated by the strongly adsorbed dissociated azole molecules.

The Relevance of Current Results in the Context of Corrosion Inhibition
The corollary of the above discussion about the inaccuracies introduced by the utilized approximations is that the phase diagrams presented in Figures 7 and 8 should be considered qualitatively rather than quantitatively (note that differences appear already between PBE vs PBE-D phase diagrams).The two most important points are that the Cu CUS sites are stabilized by molecular adsorption and that azole molecules are able to cover Cu 2 O surfaces in high coverage, unless the molecular chemical potential is too low (at a given temperature low molecular chemical potential corresponds to low molecular partial pressure in the gas-phase or low concentration in solution, i.e., the smaller the concentration the smaller the chemical potential).Note that high coverage of the inhibitor on the surface is typically deemed very important to achieve efficient corrosion inhibition, because it is often assumed that inhibitor coverage (θ) is directly related to corrosion inhibition efficiency (η), i.e., η ≈ θ [35].However, it should be noted that the current phase diagrams pertain to noncompetitive adsorption, that is, adsorption of a single molecular species at the solid/vacuum interface; this interface was chosen due to obvious modeling reasons.However, in the context of corrosion the solid/water interface is relevant.Adsorption at the latter interface is competitive, because the surface is always covered with solvent molecules and other species, such as hydroxyls and ions from the solution.Correspondingly, a given molecule adsorbs only if its adsorption is competitive enough to substitute other species from the surface.In contrast, for noncompetitive adsorption at the solid/vacuum interface the molecule adsorbs readily, provided that its interaction with the surface is attractive.This point has to be kept in mind when considering the phase diagrams presented above.
The presented phase diagrams therefore reveal that azole molecules are able to cover Cu 2 O surfaces in high coverage only in principle.Whether they will really do under given conditions for a given corrosion medium is an open-question that cannot be answered on the basis of the current results.It is known, though, that the currently considered plain azole molecules are not very good inhibitors [36][37][38], one of the reasons being their too high solubility [36].However, their functionalized forms are known to be excellent inhibitors, e.g., molecules where the benzene ring is fused with the five-membered azole ring (benzimidazole and benzotriazole) or their mercapto derivatives (e.g., mercapto-benzimidazole).

Drawbacks of the Approximate Ab Initio Thermodynamics Approach
We would like to comment on two drawbacks of the approximate ab initio thermodynamics ("abTD") approach used to construct the phase diagrams presented in Figures 6-8, A4 and A5.The first is general and not specific to "abTD", i.e., it is due to the enormous landscape of potential structures.Namely, a given phase diagram can only contain structures that are explicitly considered.If a given structure is not considered-for example, in calculations one can easily fail to identify some very stable structure-then it cannot appear on the phase diagram, even though it may be the most stable.This is why the current phase diagrams for (111)-type surfaces, presented in Figure 7, are different from those reported previously by us [12].Namely, in our previous publication we did not consider the Cu 2 O(111) w/o +1Cu CUS surfaces.
The second drawback is specific to the approximate "abTD" approach, in particular, it is the drawback associated with Equation (11).In this equation the adsorption energy, ε ads , is used to represent the chemical potential of adsorbed species, whereas entropic effects, in particular the configurational entropy, are neglected.The consequence is that in some cases the transition between the clean and fully covered r100% phases is discrete, which indicates a phase-transition (see for example the PBE-D phase diagram of tetrazole in Figure 8).Below we show that such a phase-transition can be an artifact due to the neglect of configurational entropy.
Let us consider a very simple case where the adsorption energy ε ads is independent from the coverage, i.e., ε ads = constant (Langmuir adsorption model).Figure 9a plots the adsorption surface free energies as a function of molecular chemical potential as given by the approximate "abTD" approach, Equation (11).Three regimes are evident, in particular: (i) for ∆µ MolH < ε ads the clean surface is the most stable (black horizontal line), (ii) at ∆µ MolH = ε ads all lines intersect indicating the degeneracy of θ ∈ [0, 1] coverages, and (iii) at ∆µ MolH > ε ads the 100% covered surface is the most stable.This implies that there is a discrete jump (phase-transition) from clean to full covered surface at ∆µ MolH = ε ads (Figure 9b).The resulting isotherm can be stated as: This phase-transition is, however, an artifact due to the neglect of configurational entropy.By taking configuration entropy into account one obtains the Langmuir adsorption isotherm, which can be expressed as a function of molecular chemical potential: where T is the temperature and k is the Boltzmann constant.Consequently, the transition from the clean to the fully-covered surface is smooth.Figure 9c plots the Langmuir adsorption isotherm as a function of molecular chemical potential at two different temperatures (100 and 300 K); it can be seen that when the temperature is lower, the curve is steeper; the "abTD" curve therefore corresponds to Langmuir adsorption isotherm at T = 0 K.This treatment therefore clearly demonstrates that when the ε ads is independent from the coverage, the phase-transition, predicted by the "abTD" approach (Figure 9b), is an artifact due to the neglect of configurational entropy.However, for current molecules the ε ads depends on the coverage (cf. Figure 4), hence we also consider how the "abTD" approach performs in such cases.We again consider a very simple model, that is, the Frumkin model, where the adsorption energy linearly depends on the coverage (θ) [39]: where ε o ads is the adsorption energy at zero coverage and ε lat is the lateral intermolecular energy per molecule at full coverage; ε lat < 0 for attractive and ε lat > 0 for repulsive lateral interactions.We will take a relatively small magnitude of |ε lat | = 0.05 eV, which corresponds to about 2kT at 300 K.By approximating the configurational entropy with the ideal one, which assumes that all configurations at a given θ are equally probable (this is currently not the case due to ε lat ), the corresponding relation between the coverage and ∆µ MolH (Frumkin adsorption isotherm) is: We are only interested in comparing the "abTD" isotherm to the Frumkin isotherm, hence instead of solving this equation for θ, which is non-trivial, we can plot the Frumkin isotherm parametrically.The origin of the factor 2 in the 2θε lat term in the above equation is the following.The total adsorption energy in the Frumkin model is: where N is the number of adsorbed molecules and A is the surface area.The chemical potential of the adsorbed molecule is therefore: where the pV term was neglected, because this term is negligible for solids at ambient pressures.
In the case of attractive interactions, the "abTD" approach again predicts the phase-transition from clean to fully-covered surface at ∆µ MolH = ε o ads + ε lat .However, this time the Frumkin isotherm is much closer to the "abTD" isotherm (Figure 9d) than in the Langmuir case (Figure 9b,c).For T = 300 K and ε lat = −0.05eV, the Frumkin isotherm almost coincides with the "abTD" isotherm for θ between 0.25 and 0.75, but below and above these coverages the Frumkin isotherm goes smoothly toward clean and fully-covered surfaces, respectively.
Only in the case of repulsive interactions does "abTD" approach not predict the phase-transition.Instead, the variation of coverage with ∆µ MolH can be described as: that is, coverage increases linearly with ∆µ MolH when ∆µ MolH ∈ (ε ads , ε ads + 2ε lat ).However, even in this case the "abTD" predicted transition from clean to fully-covered surface is sharper than that predicted by Frumkin adsorption isotherm (Figure 9e).It should be noted, though, that the Frumkin isotherms plotted in Figure 9d,e are too gradual, because it can be reasonably anticipated that the actual configurational entropy is smaller than the ideal configurational entropy that was used in Equation (23).This implies that the actual isotherms would be between those of "abTD" and Frumkin isotherms, when the adsorption energy is linearly dependent on the coverage.interactions were found the most repulsive for imidazole and the least repulsive for tetrazole, for which they are usually even slightly attractive.The impact of van der Waals dispersion interactions on molecular adsorption bonding was also addressed; they enhance the molecule-surface bonding by about 0.4 eV/molecule.The thermodynamic analysis reveals that although dissociated azole molecules bind more strongly to surfaces than intact molecules and although at oxygen vacancy sites the dissociative adsorption is preferred over the non-dissociative, none of the considered structures that involve dissociated molecules appear on the phase diagrams, i.e., whatever the oxygen and azole molecular chemical potentials there is always some other structure that is more stable and does not involve oxygen vacancies and dissociated molecules.Nevertheless, it is known from experiments that oxygen vacancies exist on Cu 2 O surface and that it is in fact very difficult to prepare surfaces with no oxygen vacancies in vacuum [34].The reason for this discord between the calculated phase diagrams and experimental observations is threefold: (i) we neglect the configurational entropy effects, which definitely favor the formation of vacancies, (ii) we consider oxygen vacancies only at relatively large concentration of 33% using the Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • model; according to calculations the surface with such a large concentration of vacancies is thermodynamically unstable (it seems that a lower concentration of vacancies is more realistic), and (iii) GGA may underestimate the molecular dissociation tendency (as inferred from comparison with GGA+U results) and correspondingly the stabilization of the Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • surface.Therefore, it can be reasonably anticipated that oxygen vacancies may be stabilized thermodynamically at lower concentration due to both configurational entropy effects and strong molecular adsorption.

Lateral intermolecular interactions
The thermodynamic analysis further reveals that molecular adsorption can alter the relative stability of surfaces, i.e., molecular adsorption stabilizes some surfaces more than others.In particular, our results show that the three azole molecules adsorb to otherwise thermodynamically deficient Cu CUS sites under all conditions for which molecular adsorption is stable on Cu 2 O(111)-type surfaces.This clearly indicates that molecular adsorption stabilizes the thermodynamic deficiency of coordinatively unsaturated Cu CUS ions.Indeed, azole molecules display a strong tendency to preferentially adsorb at reactive undercoordinated surface sites, such as currently considered Cu CUS and oxygen vacancy sites, and to stabilize them.This in turn tentatively suggests that the ability of these azole molecules to inhibit corrosion may stem, at least in part, from their ability to passivate reactive surface sites.There is, however, an important point to keep in mind when interpreting the current results in the context of corrosion inhibition.Namely, phase diagrams presented in the current paper pertain to noncompetitive adsorption (i.e., adsorption of a single molecular species) at the solid/vacuum interface, but in the context of corrosion the solid/water interface is relevant and adsorption thereon is competitive.Current phase diagrams therefore reveal that azole molecules are able to cover Cu 2 O surfaces in high coverage only in principle (note that inhibitor surface-coverage is often assumed to be directly related to its corrosion inhibition efficiency); whether they actually do so under given conditions for a given corrosion medium is an open-question that cannot be answered on the basis of current results.
We further analyzed a specific drawback of the approximate "abTD" approach that was used for the thermodynamic analysis and the construction of phase diagrams.In the "abTD" approach configurational entropy effects are often neglected and this results in a too hasty prediction of phase-transition like behavior between different coverages.While such a phase-transition can be reasonably expected when intermolecular lateral interactions are significant, the "abTD" approach with neglected configurational entropy also predicts a phase-transition like behavior from the clean to the fully covered surface when the adsorption energy is independent from the coverage.

Figure 1 .
Figure 1.Top-views of considered Cu 2 O surface structures with various Cu (bigger brown balls) and O (smaller red balls) ions labeled graphically.(a)-(d) four different Cu 2 O(111)-type surfaces and (e) Cu 2 O(110)-type surface.The unit cell of each surface structure is indicated with a white parallelogram.For a more detailed graphical representation of these models, see Figures 1 and 2 in Part I [1].

× √ 3 )
); these surfaces are shown graphically in Figure 1.Our calculations reasonably reproduce the results of Soon et al. [21] for Cu 2 O(111), Cu 2 O(111)-w/o-Cu CUS , and Cu 2 O(110):CuO as well as those of Li et al. [22] for the Cu 2 O(111), Cu 2 O(111)-w/o-Cu CUS , and Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • surfaces.The stoichiometric high-symmetry Cu 2 O(111) and the low-symmetry Cu 2 O(111)r display the highest surface free energies and are thus thermodynamically the least stable among the considered surfaces in the whole viable range of ∆µ O .Under oxygen-lean conditions the non-stoichiometric Cu 2 O(111)-w/o-Cu CUS is the stablest, whereas under oxygen-rich conditions Cu 2 O(110):CuO surface is preferred.As for the reconstructed Cu 2 O(111)-recon-( √ 3 R30 • , observed experimentally by Önsten et al. [23] (their model B), its surface free energy is lower than that of the stoichiometric Cu 2 O(111)r, but also significantly higher in comparison to the two non-stoichiometric surfaces, Cu 2 O(111)-w/o-Cu CUS and Cu 2 O(110):CuO.Önsten's model A of the (

Figure 2 .
Figure 2. PBE calculated surface free energies for the considered Cu 2 O surfaces as a function of the oxygen chemical potential, Equation (13).The O lean and O rich labels mark the oxygen-lean and oxygen-rich limits, respectively.The corresponding surface structures are depicted in Figure1; see also Figures1 and 2in Part I[1].

Figure 3
Figure 3  reveals that (i) intact molecules (MolH * ) bind more weakly to the surface in comparison to dissociated molecules (Mol * ), (ii) coordinatively unsaturated Cu sites (Cu CUS and Cu Ovac ) bind adsorbates more strongly than saturated Cu sites (Cu CSA and Cu surf ), (iii) on a given site non-dissociative adsorption energies are similar for all three molecules, while dissociated triazole and tetrazole bind much stronger than dissociated imidazole, hence (iv) molecular dissociation is favored just for triazole and tetrazole, but only at Cu Ovac oxygen vacancy sites.Note that after dissociation H binds to a nearby lattice O ion thus forming an OH group.Figure3ctherefore reveals that non-dissociative adsorption is usually favored over dissociative adsorption; the only two exceptions are triazole and tetrazole adsorbed at Cu Ovac on Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • , for which dissociative adsorption is favored.For benzotriazole, the dissociative adsorption was also found favourable only at oxygen vacancy sites[26] and furthermore oxygen vacancy sites were found reactive for the cleavage of the O-H bond of methanol[27].

Figure 5 .× √ 3 )
Figure 5. PBE calculated total adsorption energy per unit area, nε ads /A, as a function of coverage (the PBE-D results display the same trend and are thus not shown).Except for triazole and tetrazole on Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • , the molecules prefer to adsorb non-dissociatively.The yellowish regions on graphs for Cu 2 O(111) w/o +1Cu CUS -(2 × 2) and Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • surfaces indicate the regions of mixed-site and mixed-mode adsorption, i.e., on the former the molecules adsorb onto Cu CUS sites up to the coverage of r25%, at which all the Cu CUS sites are occupied, and above it the molecules adsorb onto remaining surface sites, whereas on Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • triazole and tetrazole adsorb dissociatively onto Cu Ovac sites up to the coverage of r33% and above it non-dissociatively.

Figure 6 .
Figure 6.Surface free energies (top row) as a function of oxygen chemical potential, γ surf (µ O ), for five considered surfaces and their stabilization (two middle rows) due to high coverage r100% adsorption, γ surf (µ O ) + nε ads /A, of imidazole (left), triazole (middle), and tetrazole (right) as calculated with PBE and PBE-D functionals; note that ordinate axes are different for each row but within the row they are the same.The bottom row shows top-view snapshots of high coverage triazole structures on the five considered surfaces; surface Cu atoms are shown as bigger brown balls and O atoms as smaller red balls, whereas molecular atoms are colored as in Scheme 1.Notice the mixed-site and mixed-mode adsorption on Cu 2 O(111) w/o +1Cu CUS -(2 × 2) and Cu 2 O(111)-recon-( √ 3 × √ 3)R30 • , respectively.In particular, on the former one molecule adsorbs at Cu CUS and three molecules at Cu CSA sites per supercell, whereas on the latter one molecule adsorbs dissociatively and two molecules non-dissociatively (labeled as H+Mol & 2MolH); dissociated molecules are encircled.
Figure 7 shows 2D phase diagrams for only Cu 2 O(111)-type surfaces, whereas in Figure 8 the Cu 2 O(110):CuO structures are superposed with Cu 2 O(111)-type structures; snapshots of the structures that appear on these phase diagrams are shown in the middle of Figure Oxygen chemical potential, Δµ O (eV)

Figure 7 .
Figure 7. Two-dimensional phase diagrams as a function of azole and oxygen chemical potentials, γads (µ O , µ MolH ), for imidazole (left), triazole (middle), and tetrazole (right) on Cu 2 O(111)-type surfaces, calculated with PBE (top row) and PBE-D (bottom row) functionals.Top-view snapshots of structures appearing on the phase diagrams are shown in Figure 8.The dependence on each chemical potential was also recast into a temperature scale at partial pressure p = 1 atm and into a partial pressure scale at T = 300 K (shown by two additional axes for each chemical potential).

Figure 8 .
Figure 8.As in Figure 7, but in addition to Cu 2 O(111)-type surfaces also the Cu 2 O(110):CuO surface is considered.Notice that Cu 2 O(111)-type structures appear on the left part of phase diagrams, whereas Cu 2 O(110):CuO structures appear on the right.Top-view snapshots and naming of structures that appear on the phase diagrams are shown for imidazole.

Figure 9 .
Figure 9.Comparison of approximate ab initio thermodynamics (abTD) approach, Equation (11), compared to Langmuir and Frumkin adsorption isotherms.(a) Adsorption free energy per unit area as a function of molecular chemical potential, calculated by Equation (11), for the case when the adsorption energy is independent from the coverage (i.e., Langmuir adsorption model) and the corresponding (b) "abTD" isotherm.(c) Langmuir adsorption isotherms at T = 100 and 300 K. Comparison between "abTD" and Frumkin isotherms at T = 300 K for (d) attractive lateral interactions, ε lat = −0.05eV and (e) repulsive lateral interactions, ε lat = +0.05eV.
-this model is the Cu 2 O(111) that lacks all the Cu CUS ions; notation "-w/o-Cu CUS " stands for "Cu 2 O(111) without Cu CUS ions".(iii) Cu 2 O(111) w/o +1Cu CUS -this model is the Cu 2 O(111) that lacks all but a single Cu R30 • pattern.Each O-vacancy is surrounded by three coordinatively unsaturated Cu ions, labeled as Cu Ovac .The Cu Ovac ions therefore always appear in triplets and we use the term "Cu Ovac site" to designate the site composed of these three ions.(v) Cu 2 O(110):CuO -this is the Cu 2 O(110) model terminated by a CuO layer on both sides of the slab.Note that the stoichiometric Cu 2 O(110) slab consists of CuO-Cu bilayers and is thus polar (terminated with a CuO layer on one side and a Cu layer on the other side).Hence, a symmetric slab terminated with the CuO layer on both sides is used instead; thus the denotation Cu CUSion either per adsorbed molecule or per supercell; the subscript "w/o" is a shorthand for "w/o-Cu CUS " and the suffix "+1Cu CUS " conveys that one Cu CUS ion is retained.2 O(110):CuO.The surface CuO layer of Cu 2 O(110):CuO consists of coordinatively saturated Cu ions, labeled as Cu surf , and O ions labeled as O surf .
)), ∆µ MolH is the molecular chemical potential measured with respect to the total energy of the isolated molecule, ∆µ MolH = µ MolH − E MolH , while ε ads stands for by fitting this equation for several values of N L , where N L is the number of either O-Cu-O trilayers in a (111) slab or CuO-Cu bilayers in a (110) slab, E L is the total energy of a single trilayer or bilayer in the bulk (both trilayer or bilayer are stoichiometric and hence multiples of the Cu 2 O formula unit), A is the area spanned by the supercell (factor 2 is due to the two equivalent surfaces of symmetric slab), and ∆N stoich O imidazole, triazole, and tetrazole on considered Cu 2 O surfaces; note that the ordinate axes are different for each surface.Lines are drawn to guide the eye.Molecules preferentially bind to the following sites: Cu CUS on Cu 2 O(111) w/o +1Cu CUS , Cu CSA on Cu 2 O(111)-w/o-Cu CUS , Cu Ovac on Only the most stable configurations at each site are considered, which-apart from two exceptions-correspond to non-dissociative adsorption mode.The two exceptions are triazole and tetrazole at Cu Ovac sites on Cu 2 O(111)-recon-( Figure 4. Top row: dependence of PBE (solid symbols and lines) and PBE-D (open symbols and dashed lines) calculated average adsorption energies on either the nearest-neighbor intermolecular distance (R nn ) or the effective nearest-neighbor intermolecular distance (R eff nn of Equation (18)) for • , and Cu surf on Cu 2 O(110):CuO.