Hydrolysis of Acetamide on Low-Index CeO2 Surfaces: Ceria as a Deamidation and General De-esterification Catalyst

Using DFT calculations and acetamide as the main example, we show that ceria is a potential catalyst for the hydrolysis of amide and similar bonds. The overall reaction is endergonic in the gas phase, yielding acetic acid and ammonia, but is slightly exergonic in the aqueous phase, which facilitates ionization of the products (CH3COO– and NH4+). Neighboring Ce and O sites on the CeO2(111), (110), and (100) facets are conducive to the formation of an activated metastable tetrahedral intermediate (TI) complex, followed by C–N bond scission. With van der Waals and solvation effects taken into account, the overall reaction energetics is found to be most favorable on the (111) facet as desorption of acetic acid is much more uphill energetically on (110) and (100). We further suggest that the Ce–O–Ce sites on ceria surfaces can activate X(=Y)–Z type bonds in amides, amidines, and carboxylate and phosphate esters, among many others that we term “generalized esters”. A Brønsted-Evans–Polanyi relationship is identified correlating the stability of the transition and final states of the X–Z generalized ester bond scission. A simple descriptor (ΣΔχ) based on the electronegativity of the atoms that constitute the bond (X, Y, Z) versus those of the catalytic site (O, Ce, Ce) captures the trend in the stability of the transition state of generalized ester bond scission and suggests a direction for modifying ceria for targeting specific organic substrates.


INTRODUCTION
The stability of peptides, proteins, antibodies, and some polymers is largely attributed to their amide (C(�O)−N) bonds, which retain a partial double bond character via resonance stabilization. 1 Large rotational barriers associated with the amide C−N bonds help preserve the 3D structure of amino acid chains in proteins. 2 Scission of peptide bonds in peptides and proteins (deamidation), and of C−N bonds in related C(�X)−N compounds, such as amidines, plays a fundamental role in proteomics, therapeutics, chemical genetics, analytical biochemistry, and biotechnology. In environmental and biochemical settings, the reaction proceeds through the formation of a tetrahedral intermediate (TI) following nucleophilic attack at the carbonyl C by water. 3,4 A base-catalyzed mechanism, in which a hydroxide group functions as the nucleophile, appears to be operative even in neutral aqueous solutions. 5−7 Oxygen exchange between amides and water is observed, which is attributed to reversible formation of the TI. The dissociation of the C−N bond in the TI is the rate-limiting step, in contrast to the dissociation of the C−O bond in the hydrolysis of carboxylate esters where the formation of the TI is rate-limiting because alkoxides (−OR) are better leaving groups than amides (−NR 2 ). 4,8 At neutral pH and ambient temperature, peptide bonds have a long half-life of 250−600 years, with respect to nonenzymatic hydrolysis. 9,10 Enzymes such as peptidase and deaminase catalyze the hydrolysis of amide C−N bonds in a narrow range of pH and temperature. Many peptidases do not contain metals, while a particular subclass, metallopeptidase, does use metal cations to catalyze the reaction. Representative members of this class include carboxypeptidase A and thermolysin, both of which contain a Zn 2+ at their active sites. How Zn 2+ accelerates amide C−N bond scission has long been debated. 11−13 The role of Zn 2+ is currently understood as positioning the carbonyl group by coordinating to its O atom, stabilizing a water molecule or hydroxide group that attacks the carbonyl C, and stabilizing one or both negatively charged C−O groups. 13,14 There is also evidence that intramolecular nucleophilic attack by an hydroxyl or thiol group that is a part of certain amino acids (e.g., serine, threonine, and cysteine) on the central C helps break the amide C−N bond under unfavorable external pH conditions, which is referred to as N−O/S acyl rearrangement. 15,16 Interest in artificial metallopeptidase has been growing. 14,17,18 Based on the understanding of the action of metallopeptidase, complexes of Cu 2+ , 19,20 Pd 2+ , 21 Pt 2+ , 22 and Co 3+ (ref 23) have been investigated for catalytic hydrolysis of amides and peptides under mild conditions. Kita et al. reported Zn acetate and Zn triflate to have higher deamidation activity than Co, Mn, and Cu acetates, while Pd, Ni, and Ag acetates have no activity. 24 Even the inorganic salt ZnCl 2 has noticeable deamidation activity, but Zn dust, which is likely oxidized Zn, has little activity. 24,25 Alcohols can be used in lieu of water for the alcoholysis of amides. 24 Solid heterogeneous catalysts have certain advantages over homogeneous catalysts with respect to recovery and separation from reaction media. Since the first reports in 1950s, 26,27 complexes of Ce 4+ have been investigated for catalyzing the hydrolysis of small peptides 28,29 and proteins 30−33 in homogeneous solutions. Higher conversions were reported for Ce 4+ complexes than Ce 3+ and other metal ion species, which is attributed to the large coordination number of Ce 4+ and its ability to form a more stable reduced Ce 3+ state by withdrawing electrons from the amide linkage. This has led to investigation into solid ceria (CeO 2 ) as a possible deamidation catalyst. Over the past decade, Shimizu and co-workers demonstrated efficient alcoholysis and subsequent esterification of amides over high loadings of CeO 2 in batch mode in hot, boiling alcohols. 34−37 Recently, the same group investigated the hydrolysis of acetamide under similar conditions of boiling solvents and high catalyst loadings using various solid oxides including Nb 2 O 5 , ZrO 2 , CeO 2 , TiO 2 , SiO 2 , and Al 2 O 3 , among others as catalysts. 38 The first three oxides were reported to achieve conversion of acetamide to acetic acid over 45%, including 100% conversion over Nb 2 O 5 . The authors proposed mechanisms for the hydrolysis and alcoholysis of amides based on the concerted actions of Lewis-acidic metal cation and basic O anion that attack the carbonyl O and carbonyl C, respectively, and supported them with density functional theory (DFT) calculations on CeO 2 (111). 35,38 They also identified a correlation between the increasing basicity of lattice oxygen in metal oxides, as represented by the O 1s binding energy, and higher initial reaction rates for several oxides. It remains unclear why CeO 2 is more effective than Nb 2 O 5 in alcoholysis 37 but not in hydrolysis. 38 Ceria-catalyzed hydrolysis and alcoholysis reported thus far in the literature have been based on commercially available polycrystalline ceria powder. 34−38 It is not clear which lowindex facet or what surface structure, if any, plays a dominant role in the observed deamidation activity. Surface structure sensitivity has been observed for many reactions on ceria. 39−42 Mullins and co-workers showed that small organic molecules, including water, alcohols, acetic acid, and acetaldehyde, exhibit qualitatively different temperature-programmed desorption (TPD) results on CeO 2 (111) vs CeO 2 (100) thin films. 43 Ceria nanocubes and nanorods, which predominantly expose (100) and/or (110) facets, and nano-octahedra, which predominantly expose (111) facets, exhibit different activities or product distributions when catalyzing a wide range of reactions including CO oxidation, 44 ethanol oxidation, 45 carbamoylation, 46 CO 2 −methanol coupling, 47 and aqueousphase dephosphorylation. 48 The differences are variously attributed to factors including different binding strength, basicity of lattice oxygen, and reducibility of the facets. 49,50 In the present work, we perform periodic DFT calculations to examine in detail the mechanism of hydrolysis of the amide C−N bond, using acetamide as the main model compound, on the three low-index facets of CeO 2 : (111), (110), and (100). The basic reaction pathway, as will be presented and discussed below, is outlined in Figure 1. Because of the significance of deamidation reactions in biochemical contexts, the effects of solvation and dispersion interactions are taken into consideration using an implicit solvation model and a van der Waals functional, respectively. We further consider C−N bond scission in benzamide, N-methylacetamide, acetamidine, and adenine on CeO 2 (111) to highlight the ability of the Ce−O− Ce site ensemble to adsorb and activate C(�X)−N bonds, which is representative of the reactivity of ceria toward an even broader class of heteroatomic organic compounds containing the motif of X(�Y)−Z, including organic amides, phosphates, and carboxylates, all of which we term generalized esters (GEs). Our results shed light on the fundamental properties that endow ceria with versatile reactivity toward biologically and environmentally relevant organic molecules, including peptides and nucleobases. They make ceria a promising enzyme-mimetic material for the development of novel diagnostic, pharmaceutical, and therapeutic technologies, among others. 51−53 The main challenge for ambient temperature hydrolysis applications is identified to be the desorption of the carboxylate product.

METHODS
Periodic DFT calculations were performed within the generalized gradient approximation (GGA-PW91 54 ) and with van der Waals corrections (optB86b-vdW), 55−58 using the Vienna Ab initio Simulation Package (VASP) 59−62 and Quantum Espresso (QE). 63,64 The potentials due to core electrons were described using the projector-augmented wave (PAW) method, 65,66 and the Kohn−Sham valence states were expanded in a plane wave basis set with a cutoff energy of 400 eV in VASP. PAW potentials for Ce, O, and C, and ultrasoft pseudopotentials (USPP) for H and N, all taken from the standard solid-state pseudopotentials (SSSP) library v1.1.2, 67 were used in QE. The (111), (110), and (100) facets of CeO 2 were modeled as p(3 × 3) slabs ( Figure 2) of 9, 4, and 9 atomic layers (corresponding to 3, 4, and 4 O−Ce−O trilayers), with periodic images of the slabs separated by 12−16 Å of vacuum and dipole decoupling applied in the z-direction. A Γ-centered Monkhorst−Pack 2 × 2 × 1 k-point grid was used to sample the surface Brillouin zone. 68 Adsorption of a single molecule per surface unit cell resulted in 1/9 monolayer (ML) coverage. The top three, two, and five atomic layers of the (111), (110), and (100) slabs were relaxed, respectively, while the remainder of each slab was held fixed at bulk coordinates. Oxygen vacancies were not considered, because those in the surface and near-surface regions are readily blocked or annihilated by the oxygenate species involved in this reaction system such as water and acetate, 69−71 and more generally, by O 2 and water occupying and dissociating in them in ambient settings. Geometry optimization was performed until the maximum residual force was 0.03 eV/Å or less in each relaxed degree of freedom.
The adsorption energy (ΔE ads ) was evaluated as ΔE ads = E total − E slab − E gas , where E total , E slab , and E gas refer to the total energies of a slab with an adsorbate, the slab without any adsorbate, and the isolated adsorbate in gas phase, respectively. A more negative value of ΔE ads indicates stronger bonding between an adsorbate and a slab. Gas-phase species were optimized in a (15 Å) 3 simulation cell, with dipole decoupling applied in all directions in VASP 72 and with the Martyna− Tuckerman correction applied in QE. 73 The elementary steps with associated transition states (TS) in the catalytic C−N bond scission and hydrolysis of acetamide and other GEs were determined using the climbing-image nudged elastic band method 74,75 and dimer method 76,77 until the maximum residual force converged to 0.03 eV/Å or less in all relaxed degrees of freedom. The reaction energy (ΔE rxn ) was evaluated as ΔE rxn = E FS − E IS , and the activation barrier (E a ) was evaluated as E a = E TS − E IS , where E IS , E TS , and E FS are the total energies of the initial, transition, and final states, respectively. Each TS was verified to possess only one vibrational mode with a negative curvature in the direction of the reaction under consideration. Vibrational modes and frequencies were calculated using finite difference approximation of the dynamical matrix with a displacement of 0.01 Å. Infrared spectra of adsorbed states were simulated using Atomic Simulation Environment (ASE). 78 The DFT+U approach of Dudarev et al. 79 was employed to partially rectify the delocalization of Ce 4f states resulting from self-interaction error. 50 A U eff value of 2 eV was used in this study based on better agreements between GGA+U predictions obtained at U eff = 2−3 eV and available experimental measurements, including the bulk reduction energy of CeO 2 , 80,81 chemisorption energy of CO on CeO 2 (110), 82 and peak desorption temperatures of AcH. 83 The calculated equilibrium lattice constants of 5.476 Å (GGA-PW91) and 5.452 Å (optB86b-vdW) 84 are in good agreement with the experimental value of 5.41 Å. 85,86 The effects of solvation by water were modeled implicitly by treating the solvent as a dielectric continuum. While VASPsol 87,88 has been popularly used to study solvation effects for a variety of molecules and extended structures, selfconsistent field (SCF) cycles failed to converge for the CeO 2 facets. Indeed, there has been no study reporting the application of VASPsol to modeling solvation of CeO 2 surfaces, to the best of our knowledge. Therefore, VASPsol was applied only to isolated molecular species at default settings. For periodic CeO 2 slabs, we used the self-consistent continuum solvation (SCCS) model as implemented in Environ for QE, 89,90 with permittivity set to 78.3 and surface tension and pressure both set to 0. A higher cutoff energy of 35 Ry, or 476 eV, was used to improve SCF convergence. While the maximum residual forces converged to below 0.04 eV/Å for molecules, ions, and the clean CeO 2 slabs, they could not be reduced to below 0.3, 0.1, and 0.1 eV/Å for some of the proposed reaction intermediates adsorbed on the (111), (110), and (100) slabs, respectively.

RESULTS AND DISCUSSION
3.1. Adsorption of Acetamide. Two stable molecular adsorption configurations can be found for acetamide on the stoichiometric CeO 2 facets: the η 1 state (Figures 3a−c), in which acetamide is located above a 3fc or 4fc site through C� O with one of the amine H atoms simultaneously interacting with a neighboring O atom in the oxide surface (denoted O latt ; see Figure 2); and a metastable state, in which acetamide binds to an O latt through its central carbonyl C, with the carbonyl O and the amine N each coordinated to an adjacent Ce (Figures 3d−f). Strong Lewis acid−base interaction between the Ce cation and the carbonyl oxygen of acetamide is previously concluded by Kamachi et al. 35 Key bond distances in the minimum-energy structures shown in Figure 3 are reported in Table 1.
Adsorption in the η 1 state causes minimal changes in the geometry of acetamide, which appears to be a nonactivated process. The η 1 state, because of its polarized C�O and C−N bonds, is open to nucleophilic attack on the central C by a nucleophile, e.g., O latt in the surface. It must overcome an activation barrier as the sp 2 hybridization of the >C�O group transitions to sp 3 hybridization that is seen in the TI in the hydrolysis of amide and peptide bonds by biochemists. 4,7,11,23,91,92 In metal ion-catalyzed hydrolysis, binding of the carbonyl O to a metal ion, followed by nucleophilic attack by water or OH, results in the formation of the TI complex. In contrast, here the TI is formed directly on CeO 2 surfaces without an additional nucleophile due to the presence of adjacent acid−base sites. The calculated E a values for this conversion (step b → c, Table 2) are 0.65, 0.57, and 0.80 eV on the (111), (110), and (100) facets, respectively. For comparison, Kamachi et al. calculated this E a value to be 0.73 eV for acetamide on (111) using GGA-PBE and an atomic basis set. 35 Similar complexes have been previously predicted theoretically for dimethylcarbonate, 46,93 benzyl acetate, 35 and phosphate monoesters on CeO 2 (111). 94,95 The C�O and C− N bonds lengthen significantly from the η 1 state to the TI state (Table 1), suggesting the latter to be a precursor state to the dissociation of the C−N bond.
The adsorption energies (ΔE ads ) of the molecular species adsorbed on the three CeO 2 facets in this reaction calculated using GGA-PW91 and optB86b-vdW are listed in Table S1 in the Supporting Information, together with available literature values. ΔE ads consistently follows the trend (111) > (110) > (100) (i.e., adsorption most stable on (100)), according to both GGA-PW91 and optB86b-vdW, including acetamide in the η 1 and TI states. The trend reflects the fact that the combined local coordination of the topmost Ce and O latt atoms is the lowest on (100) and highest on (111). 96 Stronger binding of the products (acetic acid and NH 3 ) than the reactants is evident.
Kamachi et al. reported Fourier transform infrared (FTIR) spectra of gas-phase acetamide exposed to polycrystalline CeO 2 at 120°C. 35 Three prominent peaks were found at 1430, 1554, and 1656 cm −1 . The 1656 cm −1 peak was attributed to the ν(C�O) mode of η 1 acetamide, the intensity of which decreased appreciably on the order of minutes. This assignment agrees with our simulated IR spectra for acetamide (see Figure S1 in the Supporting Information), which put this mode at 1620−1640 cm −1 on the three lowest-index facets of ceria. The mode is absent in the TI, nor is IR detection of it expected, because of the metastable nature of the TI and facile C−N bond scission (see the next section). The other two peaks in the FTIR increased in intensity with time and were attributed to acetate species. The peaks at 1430 and 1554 cm −1 agree with previous IR studies of acetic acid adsorbed on CeO 2 (111) 71 and polycrystalline ceria. 97,98 They are consistent with our analysis below, which suggests acetate to be desorption-limited.

Deamidation and Hydrolysis of Acetamide.
Based on the literature, we propose and investigate the following mechanism for the deamidation and hydrolysis of acetamide, which represents a direct reaction pathway to the products, acetic acid and ammonia. The steps include (1) adsorption of acetamide in the η 1 state and then the TI; (2) C−N scission forming acetyl and NH 2 ; (3) adsorption of H 2 O; (4) hydrogenation of NH 2 forming NH 3 ; (5) attack of acetyl by OH forming acetate; (6) desorption of NH 3 and acetic acid. The reaction parameters for each of these steps on the three low-index CeO 2 facets calculated using GGA-PW91 are tabulated in Table 2. The snapshots of the intermediate states on the (111) facet are shown in Figure 4, and those on the (110) and (100) Table 2 are mapped out in Figure 5. Since free-energy contributions at the gas/solid  interface are expected to be minor for reactions of small organic compounds on ceria under ambient conditions, 84 they are omitted from Figure 5.
A comparison of the calculated reaction energy profiles for the same overall reaction in gas phase, AcD (g) + H 2 O (g) → NH 3(g) + CH 3 COOH (g) (black lines, Figures 5a−c), reveals a progressively deepening energy landscape from (111) to (110) and to (100). The most stable state on each facet is either h (acetyl + NH 3 + OH) or i, l (NH 3 , acetate + H), meaning that the reaction is limited by product desorption, at least in the gas phase.
The E a value for C−N bond scission in the TI (step c → d) producing acetyl and NH 2 is 0.73 eV on (111), and less on (100) and (110). The TS of C−N scission is higher in energy than gas-phase acetamide (i.e., zero on ΔE axis) on (111), meaning that it is more likely for the molecule to desorb than to undergo decomposition on this facet, according to GGA-PW91. The step is endothermic by +0.62 eV on (111) (110) and (111) facets (cf. Figure 2) each provide a pair of Ce 4+ cation sites at a suitable distance and angle with respect to O latt (the nucleophile attacking the carbonyl C) that can coordinate to both the O and N atoms in acetamide simultaneously in the TS. A parallel situation is found in ester decomposition on metal surfaces, where a triangular ensemble of three metal atoms activates the C(� O)−O bond. 99 The situation is different on the (100) facet where Ce 4+ and O latt atoms are linearly aligned, so that in the TI and TS of C− N bond scission, the entire dioxo complex must shift to a nearest br site to better accommodate the rigid O−C−N angle (see Figures S3c and S3c  ‡ d). A much larger barrier exists for the TI to enolize at its methyl end (1.60 eV on (111)), so that reaction pathway is not considered further.
Akin to what has been reported for coadsorbed pairs of open-shell molecules on oxides that have large band gaps, 100−102 cooperative interaction is notable between the molecular fragments produced from C−N bond scission, acetyl, and NH 2 , on all three facets, being −0.88, −1.26, and −1.27 eV on the (111), (110), and (100) facets (see Table S2 in the Supporting Information). Therefore, the two fragments are expected to remain in close vicinity of each other at mild temperatures, and NH 2 hydrogenation to NH 3 is investigated in the presence of acetyl.
In agreement with previous reports, 103−105 we find water dissociation by itself to have small barriers on all three facets ( Table 2, step e → f). On (111), coadsorption with water lowers the E a of C−N bond scission by ca. 0.1 eV, and a coadsorbed water molecule nearly spontaneously transfers a hydrogen atom to NH 2 forming NH 3 and OH so hydrogen transfer is included as part of step c + f → h in Table 2. A water molecule in close vicinity of acetamide dissociates spontaneously on (110) and (100), although C−N bond scission that occurs in the presence of a pair of OH and H (step c + g → d + g) on these two facets does not differ much from C−N bond scission in the absence of water (cf. step c → d). C−N bond scission is followed by NH 2 accepting the dissociated H atom to forming NH 3 (step d + g → h), which has an smaller E a still than the preceding step. The energetic difference between the NH 2 and NH 3 states is the smallest on the (100) facet ( Figure 5), making it more likely for the coverage of NH 2 to build up and become detectable on (100).
The acetyl and OH groups stabilize each other strongly, like acetyl and NH 2 ( Table S2). The OH group attacks acetyl to form an η 1 acetate coordinated to an H atom (step k → l), i.e., CH 3 COO−HO latt , which amounts to a partially dissociated acetic acid. This state is iso-energetic to the μ acetate coadsorbed with an H atom on an adjacent O latt reported in ref 71. This hydroxylation step has a small E a (0.17 and 0.23 eV) and exothermic ΔE rxn on (111) and (110), but is endothermic with a sizable E a (0.78 eV) on (100). The minimum-energy path for this step on all three facets involves the OH group displacing the O latt atom to which acetyl is bonded, that is, in a sense, the OH group pushes an acetate group (i.e., acetyl−O latt ) out of the lattice site that it occupies (cf. Figure 4k) and heals the oxygen vacancy that is left behind. Both NH 3 and acetic acid adsorb on the ceria surfaces strongly (cf. Table S1) so that the desorption of the product molecules to the gas phase, particularly acetic acid, (l → m, Table 2) is more rate-limiting than all of the surface elementary steps. This is more so on the (110) and (100) facets than on (111).
We further observe that the reverse of steps l → m, k → l, f → g, and e → f constitutes a reaction channel for acetic acid to reduce the ceria surfaces, which is energetically competitive   with acetic acid desorption (Table S3 in the Supporting Information). By removing OH and H as water, a mixture of acetyl and acetate (the other species in step k → l), and not solely acetate, would populate the surface. This surface dehydration channel is not expected to be favorable in aqueous phase based on Le Chatelier's principle, but may become operative when the availability of water is low or nill. This may explain the evolution of water from CeO 2 (111) at near-ambient temperature following acetic acid adsorption in the TPD experiments of Calaza et al. 71 Although the main reaction considered here is the hydrolysis of acetamide, it is worthwhile to remark that an alternate, dehydration pathway exists for acetamide also; see discussion in the Supporting Information. The dehydration of acetamide is energetically not competitive with the deamidation and hydrolysis of acetamide but likewise could become favorable under dry conditions at elevated temperatures.

Effects of vdW Interaction and Solvation.
Our prior work 52,84 suggests that vdW interactions contribute to the stability of organic molecules as small as acetaldehyde when they are adsorbed on ceria. Thus, we have recalculated the minimum-energy reaction energy profiles on the three lowindex facets using the optB86b-vdW functional (red dashed traces, Figures 5a−c). Each profile is deepened by vdW interactions on the ΔE axis, particularly for the states in which two molecules are coadsorbed.
Biochemically, deamidation reactions occur in aqueous environments. For ceria functioning as mimics of natural enzymes, an additional possibility of ionized products (NH 4 + and CH 3 COO − ) exists since ammonia and acetic acid are mildly basic and acidic, respectively. We have calculated the free energy of hydration (ΔG hyd ) for the reactant and product species, using VASP/VASPsol and QE/Environ. The results are tabulated in Table S4 in the Supporting Information. It can be seen that the two implicit solvation models agree closely with each other, irrespective of the exchange-correlation functional used. For neutral, closed-shell species, our results also agree closely with the values calculated based on the Langevin dipoles solvation model by Warshel and co-workers and with experimental measurements. 106,107 Noticeable deviations are seen for the ionic species, with the theoretical models, particularly VASPsol and Environ, underpredicting solvation for the acetate anion and overpredicting it for the ammonium cation. Since ionic species carry localized charges that are strongly solvated by water, a lack of explicitly description of localized hydrogen bonding by the implicit solvation models is expected to lead to some errors. While the overall deamidation reaction is slightly endothermic and endergonic in the gas phase, the reaction AcD(aq) + H 2 O(aq) → NH 4 + (aq) + CH 3 COO − (aq) is computed to be mildly exergonic at −0.18 eV, which is 0.31 eV lower than the ΔG rxn for the gas phase reaction (see Table S5 in the Supporting Information). It suggests that aqueous phase is a thermodynamically necessary condition for the deamidation reaction to occur under ambient conditions. The solvated ionic final product state is indicated by a hollow triangle in Figure 5.
We have then attempted to estimate the effect of hydration on the stability of the surface intermediates in combination with vdW interactions, using the SCCS implicit solvation model implemented in QE/Environ. Solvation by water is predicted to reduce the surface energy of the clean ceria facets by 11, 21, and 25 meV/Å 2 for the (111), (110), and (100) facets, respectively, which is large, compared to, e.g. Pt (111), for which the surface energy is reduced by implicit solvation by water by 2 meV/Å 2 . 87 The sizable stabilization is consistent with the fact that ceria surfaces exhibit arrays of localized charge centers, which elicit a strong dielectric response by the solvent, independent of specific chemical or hydrogen bonding interactions between the aqueous phase and the surface.
The adsorption of the reactant and product molecules generally interferes with the interfacial dielectric response of aqueous phase to the ceria surfaces. The solvation of the surfaces is reduced (i.e., positive ΔΔG hyd , Table 3) particularly for (100). The exceptions include H 2 O on (111), which has no effect on surface solvation, and acetic acid on (111) and (110), which enhances surface solvation. Viewed from a different angle, the adsorption strengths of all the molecules are reduced (i.e., positive ΔΔE ads , Table 3), meaning more facile desorption in the aqueous phase than in the gas phase.
Overall, the proposed pathway has a shallower reaction energy profile on (111) than on (110) and (100)   (c + g ‡ d + g) TS of C−N scission + OH + H; (d + g) acetyl + NH 2 + OH + H; (d + g ‡ h) TS of NH 2 hydrogenation; (h) acetyl + NH 3 + OH; (i, k ‡ l) NH 3 , TS of OH attack; (i, l) NH 3 , acetate + H, (j, l) NH 3 desorbed; and (j, m) acetic acid desorbed. Labels correspond to those in Table 2. Hollow blue triangle in each panel corresponds to aqueous NH 4 + and CH 3 COO − . The optB86b-vdW profile, when calculated using VASP, differs from QE (red squares) by 0.11 eV or less. Acetic acid is present as an acetate + H pair on the three facets. Differences in optB86b-vdW ΔE ads calculated using VASP (reported in Table S1) vs QE are 0.11 eV or less in all cases. ΔG hyd (aq) = free energy of hydration for molecules in aqueous phase; ΔΔG hyd = change in free energy of hydration for surface with adsorbate vs bare surface; ΔΔE ads = decrease in adsorption strength due to hydration, i.e., amount by which desorption becomes more facile; ΔΔE ads = ΔΔG hyd − ΔG hyd (aq) = ΔE ads SCCS − ΔE ads repeated with vdW interaction and solvation effects included due to convergence difficulties. Instead, we make the assumption that the vdW and solvation effects on a TS are the average of these effects on the IS and FS of an elementary step. This approach would yield errors of 0−0.15 eV for the stability of TS based on the results of our previous study on aldol condensation of acetaldehyde on CeO 2 (111). 84 Under this assumption, conversion from the η 1 state to the TI is seen to be slightly outcompeted by acetamide desorption on all three facets, of which the (110) facet is the most efficient at directing a given molecular flux toward acetamide deamidation versus acetamide desorption, followed by (100), and then (111). On the other hand, the adsorption energy of acetic acid (ΔE ads SCCS , Table 3) remains at −1 eV on (111) and lower on (110) and (100) (ΔE ads SCCS , Table 3), suggesting that product desorption continues to be limiting, particularly on the two more open facets, even when assisted by an aqueous phase. A strategy that modifies the C−N bond scission activity and acetate adsorption strength differently will be needed to further optimize an oxide-based catalyst including ceria for the reaction. Once desorbed, acetic acid and NH 3 can undergo acid−base neutralization by ionizing, which provides the thermodynamic driving force for the overall reaction.

Decomposition of Amides, Amidines, and Generalized Esters on Ceria.
The activation and dissociation of the C−N bond over Ce−O latt −Ce sites are not limited to acetamide. We demonstrate this for several other amine compounds, including benzamide, N-methylacetamide, acetamidine, and adenine. The reaction energy profiles up to C−N bond scission for these amine compounds on CeO 2 (111) are plotted in Figure 6 for comparison with acetamide. The immediate fragments that result are benzoyl + NH 2 (for benzamide), acetyl + CH 3 NH (for N-methylacetamide), CH 3 CNH + NH 2 (for acetamidine), and purine + NH 2 (for adenine). These fragments also mutually stabilize on the surface (Table S2). Secondary or tertiary amides (e.g., Nmethylacetamide) produce amines instead of ammonia.
Benzamide and N-methylacetamide can be viewed as derivatives of acetamide. Benzamide has the methyl group in acetamide replaced with a phenyl group, whereas Nmethylacetamide is obtained by replacing one of the amine hydrogen atoms with a methyl group, making it perhaps the simplest model for peptide bonds in proteins. 108 Acetamidine is a simple amidine, an imine analog of an amide with the C� O group replaced with a C�NH group. These molecules are likewise open to nucleophilic attack by O latt of ceria at the carbonyl C to form TI and then undergo C−N bond scission, with E a similar to that for acetamide ( Figure 6) except for one.
Previously, we studied the adsorption of a primary nucleobase, adenine, on CeO 2 (111) as a model system for organic/inorganic interfaces. 52 Adenine can be viewed as a cyclic amidine in which the C�NH group is a part of an aromatic ring. In nature, adenine deaminase catalyzes the conversion of adenine to NH 3 and hypoxanthine. The enzyme contains a binuclear metal center of two Fe 2+ , which stabilizes a hydroxide group that attacks the C6 position of adenine, the carbon atom to which −NH 2 is attached, and replaces the amine group. 92 On CeO 2 (111), because of the aromaticity of the purine group, nucleophilic attack by O latt at the C6 position forming TI has a notably higher E a value (1.16 eV) than in the other compounds, but is followed by a more modest E a value for C−N bond scission than in acetamide (0.48 eV).
The C−N bond scission steps for acetamide, benzamide, Nmethylacetamide, acetamidine, and adenine are analyzed in terms of a linear energy relationship between the TS and the FS. 109,110 To this dataset, we also add the results for several other examples of generalized esters (GEs), which all contain the X(�Y)−Z moiety, where Y and Z are more electronegative atoms than X, resulting in a polarized X−Z generalized ester bond (Figure 7). The additional GEs include methyl formate and methyl acetate (X�C, Y�Z�O), and phenyl phosphate, para-nitrophenyl phosphate, para-chlorophenyl phosphate, and chloromethyl phosphate (X�P, Y� Z�O) that are taken from our previous work on phosphate monoesters on CeO 2 (111). 94,95 The results for acetamide, acetamidine, and methyl acetate on the (110) and (100) facets are also included, for which the dissociated fragments of RX� Y and ZR′ are notably more stable than on (111). When the stability of the TS of the X−Z bond scission is plotted against that of the FS for this diverse group of GEs, a linear Brønsted− Evans−Polanyi relationship is obtained (Figure 8a). The TS for the carboxylates and phosphates on (111), and for all the compounds considered on (110) and (100), lie below 0 eV versus the gas phase, which suggests ceria to be a potentially effective catalyst for de-esterification of such compounds. This is consistent with previous theoretical studies that report low E a for the scission of ester bonds in other examples of GEs including dimethylcarbonate 93 and dimethyl methylphosphonate. 111 Using this linear relationship, one may quickly

ACS Catalysis pubs.acs.org/acscatalysis
Research Article estimate the stability of the TS of the scission of the generalized ester bond in GEs. The data in Figure 8a cluster in separate groups for phosphates, carboxylates, and amides/amidines on a given facet. This clustering hints at underlying, systematic factors that influence the stability of the surface complexes that these organic compounds form on ceria. For acetamide and the related compounds, the carbonyl C−O latt bond obviously plays an important role in the stability of the TS of C−N bond scission on the ceria surfaces. At the same time, the carbonyl O and the amine N atoms are coordinated to adjacent Ce sites, but how strong the O(N)−Ce interaction is versus the C−O latt bond is unclear. We surmise that the nature of the X, Y, and Z atoms in a GE that correspond to the ensemble of Ce−O latt − Ce surface sites determines the stability of the TS for a given GE on CeO 2 (111). Inspired by Capdevila-Cortada et al., 96 we propose a descriptor to measure this interaction, which is based on the electronegativity (χ i ) of the atoms involved: We use the Pauling electronegativity for the elements involved, i.e., χ Ce = 1.12, χ C = 2.55, χ N = 3.04, χ O = 3.44, and χ P = 2.19. a is an adjustable parameter to represent different weighting of the X−O latt bonding relative to the interaction between the Y and Z atoms and Ce. The stability of the TS for the scission of the X−Z generalized ester bond on CeO 2 (111) included in Figure 8a is replotted against ∑Δχ in Figure 8b. While this simple descriptor does not account for the chemical environment beyond the X(�Y)−Z moiety or the catalytic site structure, Figure 8b suggests that it captures the trend in the stability of the TS, which supports our hypothesis. The optimized value of a is 1.83, indicating that the X−O latt bond is notably stronger than the interaction between the Y or Z atom and Ce 4+ . The different coordination environment and geometry of the Ce−O latt −Ce ensemble would account for the different reactivity of the facets toward the decomposition of GEs.
Of the two fragments that are produced upon scission of the generalized ester bond, the ZR′ fragment is hydrogenated into an alcohol or amine/NH 3 , while the RX(�Y) fragment is hydroxylated into RX(�Y)OH that is tautomeric with RX(� O)YH. Conceivably, if RX(�Y)OH is not tied up in an aromatic moiety (e.g., adenine), it could undergo hydrolysis again until all X−Y bonds (when Y ≠ O) are replaced by X−O bonds to produce the corresponding acid (RXOOH); e.g., acetamidine fully hydrolyzed to acetic acid. Based on the examples used in Figure 8, we further compare the ΔE ads values of several carboxylic acids, including formic, acetic, and benzoic acids, on CeO 2 (111). The difference in ΔE ads values turns out to be noticeable in the gas phase but less pronounced at the water/solid interface (see Figure S5 in the Supporting Information).

CONCLUSIONS
Ceria has been regarded as enzyme-mimicking because it exhibits catalytic activity under ambient conditions for the decomposition via hydrolysis of a range of organic compounds including carboxylates, phosphates, and amides, which is reminiscent of the action of natural metallohydrolases, including phosphatase and peptidase, that use metal ions to catalyze the hydrolysis of many biological compounds of the same types. We have performed periodic DFT calculations to investigate the mechanism and energetics of deamidation via hydrolysis on the three low-index CeO 2 facets, (111), (110), and (100), using acetamide as the main model compound.
The reaction shares similar features on the three facets, of which (111) is the least reactive and binds reaction intermediates least strongly while (100) is the most reactive. Acetamide adsorbs molecularly with its carbonyl O coordinated to a Ce cation and one of its amine H coordinated to an adjacent lattice O anion. The amide C−N bond scission is preceded by the formation of a tetrahedral intermediate (TI) state, formed via nucleophilic attack by lattice O on the carbonyl C, which has an activation barrier of E a = 0.65, 0.57, and 0.80 eV on the (111), (110), and (100) facets, respectively. In the presence of a water molecule, C−N  (6), methyl formate (7), chloromethyl phosphate (8), para-chlorophenyl phosphate (9), and phenyl phosphate (10), and para-nitrophenyl phosphate (11). Unmarked numbers indicate results on (111), single prime symbol (′) indicates results on (110), and double prime symbol (′′) indicates results on (100). The TS and FS of the phosphates are those that do not involve spontaneous dissociation of an acidic proton. Data values are listed in Table S6 in the Supporting Information, and snapshots of the TS are shown in Figure S4  scission in the TI occurs with a more modest barrier (E a = 0.65, 0.46, and 0.47 eV on (111), (110), and (100)), producing acetyl and NH 2 fragments. NH 2 readily extracts a hydrogen atom from water to form NH 3 . The remaining OH group attacks acetyl to form an acetate + H pair, which is facile on (111) (E a = 0.17 eV) but has a more substantial barrier on (100) (E a = 0.78 eV). The overall reaction is desorptionlimited in the gas phase, particularly for acetic acid on the two more open facets. The calculated desorption energies for NH 3 and acetic acid are 0.74 and 0.93 eV on (111), 0.76 and 1.51 eV on (110), and 1.09 and 1.70 eV on (100). These values are larger than the calculated E a values for all the surface reaction steps in the minimum-energy reaction mechanism on the respective facets.
Since biologically and environmentally relevant deamidation reactions occur in the aqueous phase, we have further used the SCCS implicit solvent model in combination with van der Waals corrections to estimate the reaction energetics at the aqueous interface. While vdW interactions stabilize the surface intermediates, the implicit solvent has the opposite effect and makes the reaction energy profiles less corrugated than in the gas phase. Desorption of all reactant and product species is predicted to be enhanced by the aqueous phase, although the desorption of acetic acid is expected to remain kinetically hindered at ambient temperature.
A survey of the literature and our own studies suggest that the Ce−O latt −Ce site ensembles present on ceria surfaces are well-suited to activating organic moieties of the X(�Y)−Z type, in which Y and Z are more electronegative atoms than X, resulting in a polarized X−Z bond. We term such compounds generalized esters (GEs), which includes carboxylates, phosphates, and amides among many others. A simple descriptor (∑Δχ) is proposed that takes into account the electronegativity of both the atoms constituting the bond (X, Z) and those of the catalytic site (Ce, O). It captures a linear trend in the stability of the TS of X−Z bond scission on CeO 2 (111), revealing a connection between the composition of the oxide surface and the catalytic activity for the decomposition of GEs.
Our work identifies certain factors that will need to be addressed in future catalyst design based on ceria and other metal compounds for catalytic decomposition of a wide range of organic and biologically active compounds, such as those controlling the catalytic activity for cleaving generalized ester bonds and the rate-limiting desorption of carboxylate products. The findings are also relevant to ceria-and oxide-catalyzed generalized esterification and trans-esterification reactions, which involve the formation instead of scission of generalized ester bonds. 46,93,112 ■ ASSOCIATED CONTENT
Comparison of adsorption energies of reaction reactants and products on CeO 2 (111), (110), and (100) with literature values; interaction energies of reaction intermediates on CeO 2 (111), (110), and (100); comparison of reaction energies for the desorption of water compared to desorption of acetic acid on CeO 2 (111), (110), and (100); comparison of free energies of hydration of reactant and product species with literature values; overall total and free energies for acetamide hydrolysis; tabulated data used in Figure 8; simulated vibrational spectra of acetamide adsorbed on CeO 2 (111), (110), and (100); snapshots of stable intermediates and transition states of the elementary steps in deamidation and hydrolysis of acetamide on stoichiometric CeO 2 (110) and (100) facets; snapshots of transition states of X−Z bond scission in GEs on CeO 2 (111), (110), and (100); adsorption energies of several carboxylic acids on CeO 2 (111); results and discussion on acetamide dehydration on CeO 2 (111), (110), and (100) (PDF)