Origin of Nb2 O5 Lewis Acid Catalysis for Activation of Carboxylic Acids in the Presence of a Hard Base.

The Nb2 O5 surface catalyzes the amidation of carboxylic acids with amines through Nb5+ Lewis acid activation of the C=O group. In this work, DFT calculations were applied to theoretically investigate the C=O bond activation of a model carboxylic acid (acetic acid) on θ-Al2 O3 (110), anatase TiO2 (101), and T-Nb2 O5 (100) surfaces. The adsorption sites, adsorption energies, reaction energy barriers, electronic properties, and vibrational frequency of acetic acid were examined in detail. It was found that the bond activation of the carbonyl group is most efficient on Nb2 O5 , although the adsorption energy is larger on Al2 O3 and TiO2 . The most efficient C=O bond activation on Nb2 O5 results in the lowest energy barrier of C-N bond formation during amidation. The Nb2 O5 surface also shows larger tolerance to methylamine and water molecules than Al2 O3 and TiO2 surfaces. These crucial factors contribute to the highest amidation catalytic reactivity on Nb2 O5 . Furthermore, the position of the mean density of states of the d-conduction band of the active metal site relative to the Fermi energy level correlates well with the efficiency in the C=O bond activation and, consequently, the catalytic activity for amidation. These results suggest that, unlike a classical understanding of strong acid sites of metal oxide surfaces, interaction of a carbonyl HOMO with an unoccupied metal d-orbital, or, in other words, covalent-like interaction between a carbonyl group and metal adsorption site, is relevant to the present system.

Introduction

Lewis acids can promote nucleophilic substitution reactions of carbonyl compounds, but not when the nucleophiles have strong basicity. In principle, Lewis acids are deactivated by strong coordination of basic substrates (such as amines). In condensation reactions, water (coproduct) has a greater affinity for Lewis acids than carboxylic acids and deactivates the Lewis acids. Amide synthesis from carboxylic acids (and their derivatives) and amines is a typical example of nucleophilic substitution reactions of carbonyl compounds with a strong base. Hence, amide condensation by Lewis acid catalysts remains a challenging reaction. Conventionally, amides and imides are synthesized from (di)carboxylic acids and amines via activated carboxylic acid derivatives, such as carboxylic acid anhydrides or acyl chlorides, or using stoichiometric amounts of condensation reagents,1 which generates large amounts of unwanted coproducts.

Recently reported methods for the amidation of n class="Chemical">carboxylic acids or its derivatives with amines, using homogeneous Lewis acid catalysts, are of importance due to their high atom efficiency.2 However, these methods suffer from low catalytic efficiency, and difficulties associated with catalyst reuse and catalyst/product separation. To develop a practical Lewis acid catalyst for amidation reactions, one should design a heterogeneous Lewis acid catalyst that preferentially interacts with the carbonyl group of a carboxylic acid in the presence of a strong base (amine and water), as proposed in a recent study.3

Recently, one of the authors of this paper (Shimizu), together with other coworkers, reported a new heterogeneous Lewis acid catn class="Chemical">alyst, Nb2O5, which exhibits high catalytic performance and good catalyst reusability for four types of amidation reactions: synthesis of amides from amines with carboxylic acids (reaction 1) or esters (reaction 2) and synthesis of cyclic imides from amines with dicarboxylic acids (reaction 3) or dicarboxylic anhydrides (reaction 4) (see Figure 1).4, 5, 6, 7 The selected results of catalyst screening tests for these reactions are also shown in Figure 1. Interestingly, the catalytic activity of the metal oxides in all the reactions changes in the following order: Nb2O5>TiO2>θ‐Al2O3. It should be noted that Nb2O5 shows the best performance among these catalysts even if the activities are normalized by their BET specific surface areas (54 m2 g−1 for Nb2O5, 45 m2 g−1 for TiO2, and 74 m2 g−1 for θ‐Al2O3). The activity of Nb2O5 is even higher than that of Sc(OTf)3, a well‐known homogeneous water‐tolerant Lewis acid catalyst.8,9

Figure 1

Catalytic results for amidation‐type reactions (1–4): (1) amide yield from carboxylic acid and amine (navy bar),4 (2) amide yield from ester and amine (blue bar),5 (3) cyclic imide yield from dicarboxylic acid and amine (yellow bar),6 and (4) cyclic imide yield from dicarboxylic anhydride and amine (gray bar).7

Catalyn class="Chemical">tic results for amidation‐type reactions (1–4): (1) amide yield from carboxylic acid and amine (navy bar),4 (2) amide yield from ester and amine (blue bar),5 (3) cyclic imide yield from dicarboxylic acid and amine (yellow bar),6 and (4) cyclic imide yield from dicarboxylic anhydride and amine (gray bar).7

To explain the reasons for the activity order, they carried out kinetic and infrared (IR) experiments. IR results of pyridine adsorption on these metal oxides showed a predominantly Lewis acidic nature of these catalysts. IR spectra of a model carboxylic acid (acetic acid) at −75 °C showed that the C=O stretching band of the adsorbed acetic acid (ν C=O) on Nb2O5 was located at lower wavenumber (1686 cm−1) than that those of Al2O3 (1697 cm−1) and TiO2 (1695 cm−1).4 Results indicate that the interaction between the Lewis acid site (Nb5+ cation) on the surface of Nb2O5 and the carbonyl group results in weaker carbonyl C=O bond than the sites of other Lewis acidic oxides (Al2O3 and TiO2). The acid−base interaction between the Nb5+ Lewis acid site and carbonyl group can increase the electrophilicity of the carbonyl group, which results in high reactivity of the carboxyl acid (or ester) with a nucleophile, amine. Kinetic studies suggest that Lewis acid sites of Nb2O5 are less negatively affected by basic molecules (amines, water) in solution than other Lewis acidic oxides (Al2O3 and TiO2) and homogeneous Lewis acid catalysts.4 Experimental evidence suggests that the high activity of Nb2O5 for amidation‐type reactions is due to its base‐tolerant Lewis acidity, as well as suitable acid−base interaction between the Nb5+ Lewis acid site and carbonyl oxygen, which increases the reactivity of the carbonyl group with amines. To verify this hypothesis, a theoretical study is required to achieve an atomic level understanding of the selective activation mechanism.

In this study, we apply density functionn class="Chemical">al theory (DFT) calculations to examine the adsorption of CH3COOH, CH3NH2, and H2O and the energy barrier for amidation on θ‐Al2O3, TiO2, and Nb2O5 surfaces. The catalytic activity trend among these metal oxides can be described by the energy barrier, adsorption energy, the density of states (DOS) of catalytic active sites, the adsorbed geometric bond lengths, molecular vibrational frequency, and the electronic charges. These properties enable us to describe the origin of the carbonyl C=O bond activation, the nature of the interaction between catalyst and the carbonyl oxygen atom, and base‐ and water‐tolerance properties, leading to an understanding of the relatively high catalytic activity of the Nb2O5 catalyst.

Models and Computational Details

Bulk Structure

The phase of θ‐Al2O3, anatase n class="Chemical">TiO2 and T‐phase of Nb2O5 were chosen because they were applied in the experiment.4,5 The bulk structure of θ‐Al2O3 has a monoclinic structure (Figure S1 in the Supporting Information). A unit cell consists of a total of twenty atoms, with eight Al atoms and twelve O atoms, having four formula units per unit cell. A 6×6×4 Monkhorst−Pack k‐point mesh was applied. The calculated cell parameters and those obtained from the experiment12 are shown in Table S1. The calculated parameters are in excellent agreement with what has been reported previously13 and the experimental values; deviation is in the order of ∼1 %. These parameters are also sufficient for the total energy to converge within 0.001 eV.

The bulk structure of the T‐phase of Nb2O5 is based on the reported structure.14 The cn class="Chemical">alculations employed 8×8×3 k‐points Monkhorst−Pack mesh sampling in the surface Brillouin zone for the TiO2 structure and 6×3×6 k‐points mesh sampling for the T‐Nb2O5 structure. The calculated lattice parameters of anatase TiO2 and T‐Nb2O5 are in excellent agreement with the experimental values.14,17 They are reported in Table S1.

Surface Structure

The θ‐Al2O3(110) surface is represented by a slab model with similar geometric details as applied in our previous work.18 The (110) surface of the face‐centered cubic (fcc) n class="Chemical">oxygen sublattice corresponds to (010) of the surface index of a monoclinic unit cell structure of θ‐Al2O3.19 Previous studies showed that the preferential exposed surface of θ‐Al2O3 is the (110) surface of the fcc oxygen sublattice.20,21 We employed the 2×1 supercell with a three θ‐Al2O3 cell thickness slab containing 180 atoms and a vacuum region of ∼15 Å in the z direction. It was previously found that the θ‐Al2O3 charge‐neutral slab was achieved with a three to four θ‐Al2O3 cell thickness.19 The θ‐Al2O3(110) surface is shown in Figure 2(a). The θ‐Al2O3(110) surface calculations employed a 2×2×1 k‐points Monkhorst−Pack mesh sampling. All of the Al and O atoms, and adsorbate molecules, were fully optimized, except for the Al and O atoms on the bottom atomic layer of Al2O3, which were fixed to the calculated bulk distances.

Figure 2

Top view of the clean surface structures and the possible active sites of (a) θ‐Al2O3(110), (b) TiO2(101), and (c) T‐Nb2O5(100). AlO3 and AlO4 refer to an Al surface atom having atomic coordination numbers 3 and 4, respectively. TiO5 and TiO6 refer to a Ti surface atom having atomic coordination numbers 5 and 6, respectively. NbO5 and NbO6 refer to a Nb surface atom having atomic coordination numbers 5 and 6, respectively.

Top view of the clean surface structures and the possible active sites of (a) θ‐n class="Chemical">Al2O3(110), (b) TiO2(101), and (c) T‐Nb2O5(100). AlO3 and AlO4 refer to an Al surface atom having atomic coordination numbers 3 and 4, respectively. TiO5 and TiO6 refer to a Ti surface atom having atomic coordination numbers 5 and 6, respectively. NbO5 and NbO6 refer to a Nb surface atom having atomic coordination numbers 5 and 6, respectively.

The 3×1 supercell of the TiO2(101) surface consisn class="Chemical">ting of six Ti atomic layers and twelve oxygen atomic layers containing a total of 108 atoms and a vacuum region of ∼15 Å in the z direction was applied. The TiO2(101) surface structure is shown in Figure 2(b). The calculations employed 2×2×1 k‐points Monkhorst−Pack mesh sampling. The bottom two atomic layers of Ti atoms and four atomic layers of O atoms were fixed at their bulk optimized positions, while the rest were fully relaxed.

The T‐Nb2O5(100) surface model was constructed from the cn class="Chemical">alculated T‐Nb2O5 bulk structure. The Nb terminated (100) surface was applied in this work because it was found to have the lowest surface energy.22 The 1×2 supercell of T‐Nb2O5(100) surface is applied in this work and contains a total of 168 atoms (48 Nb and 120 O atoms). The top view of T‐Nb2O5(100) surface structure is shown in Figure 2c. The Nb surface atoms consist of 5 and 6 atomic coordination numbers (Figure 2c). The side view of T‐Nb2O5(100) surface structure is shown in the Figure S2. The slab thickness of approximately 8.5 Å was employed and a vacuum region of ∼15 Å was added in the z direction. Adsorbates and all atomic layers from the surface down to ∼3 Å in the z direction were fully relaxed and the rest in the slab were fixed at their bulk positions. The reciprocal space was only sampled at the Γ k‐point of the Brillouin zone due to the relatively large size of the slab unit cell. The similar T‐Nb2O5(100) surface model was previously applied for glycerol dehydration investigation.22 It should be noted that deactivation of Nb2O5 catalyst was not experimentally found at least for several recycling tests.4, 5, 6, 7 The X‐ray diffraction (XRD) measurements of Nb2O5 catalyst structure both before and after the amidation‐type reactions is shown in Figure S3. The XRD patterns remained essentially the same after the reactions, indicating that the Nb2O5 catalyst is stable during the reactions.

General Computational Settings

The specific calculations parameters for each bulk and surface structures are described above and the general settings are described in this section. The fully periodic plane‐wave DFT calculations as implemented in the Vienna ab initio simulation program (VASP)23,24 were employed. Spin‐polarized DFT calculations were performed using the Perdew−Burke−Ernzerhof (PBE) exchange−correlation functional10,11 described within the generalized gradient approximation (GGA) implemented with the projector augmented‐wave function (PAW)15,16 method for representing the nonvalence core electrons. The pseudopotentials were used to model the semi‐core p states of Ti and Nb. The plane‐wave cutoff energy was optimized at 400 eV. For the inclusion of long‐range dispersion, the semiempirical dispersion potential correction method described by Grimme and coworkers25,26 was applied. Gaussian broadening27 was employed with a smearing width of 0.1 eV. The dipole corrections in all three spatial directions were included for isolated molecules calculations and only in z direction was included for surface calculations. The results were checked for convergence with respect to numbers of k‐point and energy cutoff. The convergence criteria for electronic self‐consistent iteration was set to 1.5×10−7 eV, and the ionic relaxation loop was limited for all forces smaller than 0.03 eV/Å for free atoms. Bader charge analysis was performed using VASP‐VTST.28, 29, 30 The transition state structures were located using nudged elastic band (NEB)31 and dimer32 methods. The transition state structures were characterized by a normal mode analysis to ensure that it has one imaginary frequency.

Results and Discussion

Carboxylic Acid, Methylamine, and Water adsorption

CH3COOH Adsorption

The adsorption energy can be used as a preliminary probe of n class="Chemical">Lewis acidity of metal oxide surfaces. The adsorption energy ( ) was calculated by [Eq. (1)]:

where is the energy of the system upon adsorption, is the energy of the clean surface, and is the energy of isolated adsorbate in gas phase. The cn class="Chemical">alculated adsorption energies of the reactants on three metal oxide surfaces at the most stable structure in each site are summarized in Table 1. The adsorption energies of CH3COOH on TiO2(101) surface are in good agreement with those previously reported as shown in Table S2.

Table 1

The strongest adsorption energies of CH3COOH, CH3NH2, and H2O on θ‐Al2O3(110), TiO2(101), and Nb2O5(100) surfaces (in eV).[a]

Adsorbate	θ‐Al2O3(110)		TiO2(101)[c]		Nb2O5(100)
	AlO3	AlO4	TiO5‐1	TiO5‐2	NbO5	NbO6
CH3COOH	−3.93[b]	−1.21	−1.23	−0.98	−0.83	−0.66
CH3NH2	−2.40	−1.85	−1.49	–	−1.14	−1.09
H2O	−1.92	−1.56	−0.87	–	−0.64	−0.55

[a] The corresponding structures are shown in Figures 3, 4, and S5, for CH3COOH, CH3NH2, and H2O, respectively. [b] Dissociative adsorption. [c] The stable adsorption of CH3COOH at TiO6 site was not found.

The strongest adsorption energies of n class="Chemical">CH3COOH, CH3NH2, and H2O on θ‐Al2O3(110), TiO2(101), and Nb2O5(100) surfaces (in eV).[a]

[a] The corresponding structures are shown in Figures 3, 4, and S5, for CH3COOH, n class="Chemical">CH3NH2, and H2O, respectively. [b] Dissociative adsorption. [c] The stable adsorption of CH3COOH at TiO6 site was not found.

The amidation reaction on metal oxide surfaces commences with the Lewis acid activation of CH3COOH, followed by nucleophilic amine attack on the carbonyl carbon of the adsorbed CH3COOH. Thus, the C=O bond activation of CH3COOH is a key factor in the reaction. Therefore, we now discuss the CH3COOH adsorption. The adsorption of CH3NH2 and H2O is discussed later in sections 3.1.2 and 3.1.3. The surface reactivity varies at different surface sites. This variation can be ascribed to the different coordination environments of the surface atoms. The adsorption energies for CH3COOH were calculated at all metal cation sites with different coordination numbers. Different configurations are found on the three surfaces, where the OH group in CH3COOH is dissociated or nondissociated. The dissociative adsorption of CH3COOH shows that a proton adsorbs on the oxygen site of the surface and the CH3COO− species on a metal site. The CH3COO− adsorbs as bidentate chelating, in which both oxygen atoms interact to the same metal atom site, or as bridging, in which the oxygen atoms bind to different metal atom sites. For all surfaces shown in Figure 3, the most stable adsorption structures show CH3COOH interacting with the metal sites through the oxygen atom of the C=O group.

Figure 3

Configurations of (a) CH3COOH adsorption on θ‐Al2O3(110) surface at AlO3 (left) and AlO4 (right) sites, (b) on TiO2(101) surface at TiO5 site with hydrogen bond formed with O 2 atomic coordination number (left) and O 3 atomic coordination number (right) sites, and (c) on Nb2O5(100) surface at NbO5 (left) and NbO6 (right) sites.

Configurations of (a) n class="Chemical">CH3COOH adsorption on θ‐Al2O3(110) surface at AlO3 (left) and AlO4 (right) sites, (b) on TiO2(101) surface at TiO5 site with hydrogen bond formed with O 2 atomic coordination number (left) and O 3 atomic coordination number (right) sites, and (c) on Nb2O5(100) surface at NbO5 (left) and NbO6 (right) sites.

The adsorption of n class="Chemical">CH3COOH on the AlO3 site was found to be dissociative (Figure 3(a), left). Both oxygen atoms of CH3COO− bind to the AlO3 sites (bridging on different rows). However, the O−H dissociative structure is assumed to be too stable for further amidation reaction. We examined other possible adsorption configurations, yet we could not locate the nondissociative adsorption of CH3COOH on the AlO3 site. The dissociative adsorptions of CH3COOH on TiO2(101) and Nb2O5(100) surfaces were also found but were significantly less stable than the nondissociative adsorptions reported in Table 1 by 0.81 and 0.47 eV, respectively.

The adsorption structures of n class="Chemical">CH3COOH on TiO2(101) and Nb2O5(100) surfaces include hydrogen bond formation between a hydrogen atom of the carboxylic OH group and an oxygen of the metal oxide surface, in which the bond distances are in the range 1.51–1.78 Å (Figure 3b and 3c). The structure of θ‐Al2O3(110) at the AlO4 site shows that the CH3 group, instead of the carboxylic OH group, forms very weak hydrogen bonds (bond lengths are 2.07 and 2.31 Å) with the oxygen of the Al2O3 surface (Figure 3(a), right).

The adsorption energy trend indicates the weakest adsorption of CH3COOH on the Nb2O5 surface and the strongest adsorption on the Al2O3 surface (Table 1). The metal−O(carbonyl) (M−OC=O) bond length trend correlates well with the adsorption energy trend. The M−OC=O bond length is in the order Al2O3 (1.94 Å)<TiO2 (2.11 Å)<Nb2O5 (2.24 Å). However, a different trend was found in the C=O bond length. Upon adsorption, the carboxylic C=O bond length increases in the order Al2O3 (1.24 Å)<TiO2 (1.25 Å)<Nb2O5 (1.26 Å) (free molecule: 1.22 Å). This suggests that upon carboxylic acid adsorption, the C=O bond is most activated on the Nb2O5 surface. The trend of this bond activation was also found in the lowest vibrational frequency of the C=O stretching mode of the adsorbed CH3COOH on the Nb2O5 surface and the highest one on the Al2O3 surface. The calculated vibrational frequencies of C=O stretching mode for the most stable adsorbed structure are 1685, 1623, and 1616 cm−1 on Al2O3(110), TiO2(101), and Nb2O5(100), respectively. The qualitative trend agrees well with the respective experimental values of 1697, 1695, and 1686 cm−1.4

The partin class="Chemical">al charges of the adsorbed CH3COOH and metal oxide surface are shown in Table 2. Charges of the carbonyl oxygen are more negative upon adsorption on all metal oxide surfaces, with the highest negative charge (−1.34|e|) on θ‐Al2O3(110), the least negative charge (−1.15|e|) on Nb2O5(100), and a moderate charge of −1.19|e| on TiO2(101). The change in metal atom charge due to CH3COOH adsorption is very small (0.03–0.07|e|) on all metal oxide surfaces. Charges of the carbonyl oxygen and metal atom on Nb2O5(100) show the smallest changes upon adsorption compared to those on θ‐Al2O3(110) and TiO2(101). This clearly indicates that the nonelectrostatic or covalent‐like interaction is relevant, especially for the Nb−OC=O bond.

Table 2

Bader charge (values given in |e|) analysis of metal sites and CH3COOH upon adsorption on Al2O3(110), TiO2(101), and Nb2O5(100).

θ‐Al2O3(110)			TiO2(101)				Nb2O5(100)
	clean surface	AlO4		clean surface	TiO5‐1	TiO5‐2		clean surface	NbO5	NbO6
Al3 [a]	+2.44	+2.47	Ti3	+2.13	+2.20	+2.19	Nb11	+2.62	+2.66	+2.61
O2	−1.57	−1.57	O6	−0.97	−1.07	−0.97	Nb1	+2.67	+2.67	+2.70
O4	−1.57	−1.56	O9	−1.15	−1.13	−1.19	O9	−1.04	−1.09	−1.04
							O3	−0.98	−0.99	−1.04

[a] The atom numbers correspond to the adsorption geometries shown in Figure 3.

Bader charge (values given in |e|) ann class="Chemical">alysis of metal sites and CH3COOH upon adsorption on Al2O3(110), TiO2(101), and Nb2O5(100).

Here, we focus on the adsorption energy of n class="Chemical">carbxylic acid and its activation on three surfaces. The higher activities of Nb2O5 have experimentally been validated in our previous works.4, 5, 6, 7 The results of in situ FT‐IR studies4 support the activation of C=O bond in the carboxylic acids. The adsorption of reactants is a fundamental step for surface‐catalyzed reactions. Various characteristic properties of the adsorbed species such as bond distances, bond orders, vibrational frequencies are also frequently used as descriptors to rationalize and predict the catalytic activities.33,34 Therefore, the adsorption properties are very informative and lead to understanding activation of reactants and eventually the catalytic activities.

Base Tolerance

Basic molecules such as amines (reagents in the amidation reactions discussed here) generally suppress the activity of Lewis acid catalysts by strong adsorption on the acid sites of the catalysts. Based on earlier experiments carried out on the amidation reaction that we investigate in this paper,4,5 the adsorption of the other reactant, CH3NH2, on the three metal oxide surfaces was now examined here. The calculated adsorption energies of CH3NH2 on θ‐Al2O3(110), TiO2(101), and Nb2O5(100) surfaces are summarized in Table 1 and their adsorption configurations are shown in Figure 4.

Figure 4

Structures of CH3NH2 adsorption on (a) θ‐Al2O3(110) surface at AlO3 (left) and AlO4 (right) sites, (b) TiO2(101) surface at TiO5 site, and (c) Nb2O5(100) surface at NbO5 (left) and NbO6 (right) sites.

Structures of CH3NH2 adsorpn class="Chemical">tion on (a) θ‐Al2O3(110) surface at AlO3 (left) and AlO4 (right) sites, (b) TiO2(101) surface at TiO5 site, and (c) Nb2O5(100) surface at NbO5 (left) and NbO6 (right) sites.

Table 1 shows that the adsorption of n class="Chemical">CH3NH2 is weakest on the Nb2O5 surface (−1.16 eV) and strongest on the Al2O3 surface (−2.42 eV). The CH3NH2 adsorption is non‐dissociative on the three surfaces. The atomic coordination environment is responsible for differences in adsorption energy, especially on the Al2O3 surface; the adsorption of CH3NH2 on the AlO3 site is considerably stronger than that on the AlO4 site. The trend in bond distance between the metal and nitrogen atom (M−N) naturally correlates to the trend in adsorption energy, in which the stronger adsorption results in the shorter M−N bond length (Figure 4).

The adsorption of n class="Chemical">CH3NH2 is stronger than that of CH3COOH on the three surfaces. However, the adsorption of CH3NH2 is significantly weak on the Nb2O5 surface compared to that on Al2O3 and TiO2 surfaces. This suggests that base poisoning on the Nb2O5 surface is less likely to occur compared to Al2O3 and TiO2 surfaces, which allows the metal active sites to be available for the reactant, CH3COOH, resulting in higher catalytic activity of the Nb2O5 catalyst.

The results of Bader charge analysis of n class="Chemical">CH3NH2 adsorption are given in the Supporting Information (Table S3). Upon adsorption on the Nb2O5 surface, the partial charge of the N atom varies by only −0.02|e|, while it changes by −0.17 and −0.23|e| upon adsorption on the Al2O3 and TiO2 surfaces, respectively. The change in metal charge due to amine adsorption is very small on all surfaces; the change is in the range +0.00–0.04|e|. This reveals that the interaction between the N atom of the amine and the Nb atom induces very small charge transfer between the N and the NbO5 site, compared to between N and the AlO3, AlO4, and TiO5 sites, which results in weak adsorption of CH3NH2 on the Nb2O5 surface.

Water Tolerance

In condensation reactions in which the water molecule is a coproduct, water produced during such reactions can suppress catalytic activity by strong adsorption on acid sites of the catalysts. Substantially weak interaction between water molecules and acidic active sites is required to maintain the high catalytic activity of the catalyst.

Table 1 shows that the calculated n class="Chemical">H2O adsorption energy trend is similar to that of CH3NH2 adsorption; Nb2O5 shows the weakest H2O adsorption. The Nb2O5 catalyst shows the better water tolerance compared to the Al2O3 and TiO2 surfaces. Thus, water, which is a reaction coproduct, is not likely to poison the catalytic active sites, hence leading to the higher catalytic activity on the Nb2O5 catalyst. This is in agreement with many previous findings, i. e., that Al2O3 is quite reactive toward H2O. It was found that water molecules dissociate and strongly adsorb on Al2O3 and TiO2 surfaces, with the formation of hydroxylated sites.35, 36, 37, 38 Also, our recent experimental work demonstrated water‐tolerant Lewis acid catalysis of the related material, Nb2O5 ⋅ nH2O.39 The adsorption geometries of water molecule are shown in Figure S5.

Metal Oxides: Density of States

Generally, at the n class="Chemical">Lewis acid site, the charge‐transfer or electron‐donating interaction occurs between the highest occupied molecular orbital (HOMO) of the adsorbate and the states in the conduction band of the surface. To understand surface−adsorbate interaction in the present system on the three metal oxide surfaces, we analyzed the mean of the partial density of states (PDOS) for the active‐site surface atoms referenced to the Fermi level. For example, the mean PDOS of s‐conduction ( ) and s‐valence ( ) bands is defined by Equations (2) and (2):40

where is the Fermi level, is the s PDOS, and is the energy. By definition, the HOMO of the adsorbate is located at the Fermi level. , , , and were cn class="Chemical">alculated in the same manner, applying the integration over and . The standard method to align the energy scales of PDOS is applied by shifting the energy scales of PDOS with respect to their Fermi level (E F=0). The calculated mean PDOS values are tabulated in Tables S5 and S6.

Figure 5 shows the PDOSs of metal surface atoms (possessing different n class="Disease">atomic coordination numbers) before adsorption, and the mean PDOSs in the conduction and valence bands (vertical dotted lines). The mean PDOSs of s‐ and p‐conduction bands of the AlO3 site are located closer to the Fermi level (closer by ∼0.6 eV), than those of the AlO4 site. Thus, the Al atom at the AlO3 site tends to accept electrons from the HOMO of the adsorbate more favorably than the Al atom at the AlO4 site. This suggests that the AlO3 site could be more active as a Lewis acid than the AlO4 site, leading to strong interaction at the AlO3 site, where the spontaneous dissociative adsorption is thermodynamically favorable. Correlation of the mean PDOSs of the p‐conduction band with the activity of an Al surface atom as Lewis acid has also been reported earlier.40 Furthermore, the d electrons are dominant for the reactivity of metal sites on TiO2 and Nb2O5 surfaces; therefore, we mainly consider d electrons for the Ti and Nb sites. However, the descriptor of the mean PDOSs of d‐conduction bands does not correlate to the adsorption energy difference between TiO5 and TiO6 sites and that between NbO5 and NbO6 sites (see Figure 5). The NbO5 and NbO6 sites show the mean PDOSs of d‐conduction bands differences of 0.05 eV, while the CH3COOH adsorption was significantly stronger at the NbO5 site, by 0.18 eV (Table 1).

Figure 5

PDOSs of s‐, p‐, and d‐electrons for the AlO3 and AlO4 sites of Al2O3, Ti surface atom at TiO5 and TiO6 sites, and the Nb surface atom at NbO5 and NbO6 sites. The vertical dotted lines show the mean PDOSs of s (black), p (blue), and d (red) electrons in the conduction and valence bands. The energy is referenced to the Fermi level and the Fermi level is located at .

PDOSs of s‐, p‐, and d‐electrons for the AlO3 and n class="Chemical">AlO4 sites of Al2O3, Ti surface atom at TiO5 and TiO6 sites, and the Nb surface atom at NbO5 and NbO6 sites. The vertical dotted lines show the mean PDOSs of s (black), p (blue), and d (red) electrons in the conduction and valence bands. The energy is referenced to the Fermi level and the Fermi level is located at .

Charge transfer upon carboxylic acid adsorption differs between the Al site and the transition metal sites of Ti and Nb atoms. As seen in Table 2, there is a significant charge change upon adsorption of OC=O on the AlO4 site (0.23|e|), while adsorptions on the TiO5 and NbO5 sites are very small (<0.08|e|). This difference is also evident in the change of PDOS at the adsorption site before and after adsorption. Figure 6(a) shows that, upon adsorption, the mean PDOSs of the s‐ and p‐conduction bands of the AlO4 site show large energy shifts to higher energies, by 1.1 eV and 0.66 eV, respectively. Therefore, the charge transfer from the HOMO of adsorbates to s‐ and p‐ conduction states of the AlO4 site (Lewis acid) is seen on the Al2O3 surface. However, at the TiO5, NbO5, and NbO6 sites, the OC=O charge change due to adsorption is very small (<0.08|e|, Table 2) and the shifts of the mean d‐conduction states upon adsorption are very small compared with the mean PDOS of the p‐conduction band of the AlO4 site; TiO5=0.09 eV and NbO5=0.22 eV (Figure 6(b and c)). The shifts of the mean PDOSs of the s‐ and p‐conduction bands of TiO5 and NbO5 upon adsorption are also very small (<0.2 eV) (Figure S6, and Tables S5 and S6). This indicates that the same descriptor (the mean PDOS of the p‐conduction band) used to probe both the electron acceptance ability as a Lewis acid and adsorbate binding strength of the Al2O3 surface could be different from that used to probe that of the TiO2 and Nb2O5 surfaces. It should be noted that the presence of the shifts of the mean d‐conduction states upon adsorption on the TiO5 and NbO5 sites suggests charge transfer from HOMO of adsorbates to the d‐conduction states on the TiO2(101) and Nb2O5(100) surfaces as well; however, the shifts are much less pronounced than on the Al2O3(110) surface.

Figure 6

PDOSs of metal active sites before and after CH3COOH adsorption and PDOSs of OC=O and CC=O of CH3COOH after adsorption on (a) Al2O3(110), (b) TiO2(101), and (c) Nb2O5(100). The top and second panels show PDOSs of s, p, and d electrons of the metal active sites on clean surfaces and after CH3COOH adsorption, respectively. The bottom panel shows PDOS of p electrons of OC=O and CC=O of adsorbed CH3COOH. The dotted vertical lines show the mean PDOSs of s, p, and d electrons in the conduction and valence bands. The energy is referenced on the Fermi level and the Fermi level is located at .

The interaction n class="Gene">between the metal d‐valence states and the orbitals of adsorbates could be an important factor in understanding the adsorption strength on metal surfaces due to back‐donation interaction.41 The mean PDOS of the d‐valence band (the d‐band center theory) of the Ti atom on a clean TiO2 surface is at higher energy and closer to the Fermi level than that of the Nb atom on a clean Nb2O5 surface. Thus, the mean PDOS of the d‐valence band suggests the stronger adsorption of CH3COOH on the TiO2 surface than on the Nb2O5 surface.

The plot in the bottom panel of Figure 6 shows the OC=O(p) upon adsorption of CH3COOH on surfaces. Upon adsorption, the LUMO of OC=O(p) is shifted to the energy level so that it can hybridize with the states in the d‐conduction bands of the Ti and Nb atoms. Also, the HOMO of OC=O(p) is shifted to the energy level so that it can interact with the states in the d‐valence bands of the Ti and Nb atoms. Then the onset of the occupied orbitals of OC=O(p) is at approximately the same energy level with the onset of the d‐valence energy level of the TiO5 and NbO5 sites. In contrast, upon adsorption on the Al2O3 surface, the LUMO of the O(p) carbonyl does not shift to hybridize with most states in the s‐ and p‐conduction bands of the AlO4 site. Thus, the PDOSs of OC=O(p) (Figure 6) also provide interesting insight into the difference in the interaction between OC=O(p) with Al(p) and between OC=O(p) with Ti(d) and Nb(d) sites. This indicates that the covalent‐like interaction is relevant for the TiO2 and Nb2O5.

For the Al2O3 surface, the mean PDOS of the p‐conduction band of the AlO4 site and the adsorption energy (E ad) are correlated in the way that the conventional Lewis acid characteristic should possess. The conventional Lewis acid catalyst is anticipated to show the stronger adsorption, with the mean PDOS of the conduction band of the atomic metal located at lower energy, which is closer to its Fermi level. The mean PDOS of the p‐conduction band of the AlO3 site located at lower energy and closer to the Fermi energy level clearly identifies the stronger Lewis acidity of the AlO3 site than that of the AlO4 site. However, this relationship does not simply hold in the case of the Ti and Nb sites of TiO2 and Nb2O5. The well‐known d‐band center descriptor41 can apply to the adsorption strength trend at the TiO5 and NbO5 sites. The d‐band center of the TiO5 site locates closer to the Fermi energy level than that of the NbO5 site, due to the stronger adsorption on the TiO5 site.

The Nb2O5(100) surface shows the weakest n class="Chemical">CH3COOH adsorption, resulting in the longest M‐OC=O bond length. However, the adsorbed carboxylic C=O on Nb2O5(100) is activated the most, the adsorbed carboxylic C=O bond length increases the most, and the vibrational frequency of the C=O stretching mode is the weakest. The relatively small atomic charge change of OC=O and the hybridization between the PDOS of OC=O(p) and PDOSs of d‐conduction bands of Ti and Nb upon adsorption of CH3COOH indicate that the covalent‐like interaction could be relevant to the higher degree of C=O bond activation on the TiO2 and Nb2O5 surfaces than that on the Al2O3 surface. Furthermore, analysis suggests that the location of the mean PDOSs of d‐conduction bands of TiO2 and Nb2O5 and the p‐conduction band of Al2O3 can be correlated to the degree of C=O bond activation: Nb2O5>TiO2>Al2O3.

Carboxylic Acid Amidation Pathway

The possible reaction mechanism of amidation reaction between carboxylic acid and amine over Nb2O5 is shown in Scheme 1.

Scheme 1

Possible reaction mechanism of amidation reaction between carboxylic acid and amine over Nb2O5.

The produced H+ and OH− species on the Nb2O5 surface are consumed to form n class="Chemical">water so that the active sites would not be occupied by these species. The crucial step of amidation reaction is the C−N bond formation between acid and amine (Step 1). Thus, the energy barrier of step (1) were calculated on Al2O3(110), TiO2(101), and Nb2O5(100) surfaces and shown in Table 3.

Table 3

Energy barrier (eV) of step (1) C−N bond formation and calculated C=O vibrational frequency (ν C=O, cm−1) of adsorbed carboxylic acid on Al2O3(110), TiO2(101), and Nb2O5(100) surfaces.

Surface	Energy barrier [eV][a]	ν C=O [cm−1]
Al2O3(110)	1.11	1685
TiO2(101)	0.92	1623
Nb2O5(100)	0.84	1616

[a] The corresponding transition state structures are shown in Figure S7.

Energy barrier (eV) of step (1) C−N bond formation and cn class="Chemical">alculated C=O vibrational frequency (ν C=O, cm−1) of adsorbed carboxylic acid on Al2O3(110), TiO2(101), and Nb2O5(100) surfaces.

Nb2O5 surface yields the lowest energy barrier of 0.84 eV for C−N bond formation compared to Al2O3 and TiO2 surfaces. The lowest energy barrier of the C−N bond formation on Nb2O5 surface could result from that the C=O bond is most activated on Nb2O5 showing the C=O longest bond length and the lowest vibration frequency of C=O stretching mode upon adsorption as discussed in section 3.1.1. The trend of C−N bond formation energy barrier, (Nb2O5<TiO2<Al2O3) agrees with the trend of C=O bond activation which is evident by the vibration frequency C=O stretching mode.

Conclusions

To understand the performance of the recently developed direct amide bond formation on metal oxide surfaces, we investigated the C=O bond activation of carboxylic acid on Al2O3, TiO2, and Nb2O5 surfaces. DFT calculations were performed for the θ‐Al2O3(110), anatase TiO2(101), and T‐Nb2O5(100) surfaces. The reactive adsorption sites of carboxylic acid were found to be at the AlO4, TiO5, and NbO5 sites. The adsorption energies of carboxylic acid were in the order Al2O3>TiO2>Nb2O5, while the bond activation of C=O and its vibrational frequency showed the opposite trend. This means that the C=O bond activation of the carbonyl group is most efficient on Nb2O5(100) surface, although the adsorption energy is larger on the Al2O3 and TiO2 surfaces. The base and water tolerance of the metal oxide surfaces is also a crucial factor for the catalysts of direct amidation reactions. The adsorption of CH3NH2 is significantly weak on Nb2O5 compared to that on Al2O3 and TiO2 surfaces. Water molecule is also adsorbed substantially weak on Nb2O5 surface. These suggest that Nb2O5 surface possess better base and water tolerance property than Al2O3 and TiO2 surfaces. Furthermore, the most activation of the C=O bond on Nb2O5 leads to the lowest energy barrier of C−N bond formation which is the key step of the amidation. These crucial factors are the origin for the high activity of the Nb2O5 surface in the direct amidation reaction.

The relevant covalent‐like interaction found in Nb2O5, which is seen in the hybridization between the PDOS of the Nb d‐conduction band and the PDOS of OC=O(p), accompanied by adsorption of CH3COOH, indicates the important role it plays in carbonyl C=O bond activation. Therefore, unlike the standard interpretation of the Lewis acid site of metal oxide surfaces, the covalent‐like interaction between carbonyl group of carboxylic acid and metal site is shown to be relevant in carbonyl C=O bond activation and, consequently, in the reactivity of the present direct amidation reaction of carboxylic acids and amines. The position of the mean PDOS of the d‐conduction band demonstrates a good correlation with the degree of carbonyl C=O bond activation; namely, the donative interaction between carbonyl HOMO and metal unoccupied d‐orbital. The closer the mean PDOS of the d‐conduction band to the Fermi level, the higher the degree of carbonyl C=O bond activation on the Lewis acid site of the metal oxide. This new understanding of the Lewis acid site elucidates the reason for Nb2O5 being the most effective heterogeneous catalyst for the direct amidation of carboxylic acids among the metal oxide surfaces examined here. Furthermore, it is useful for the further development of heterogeneous catalysts.

Associated Content

Supporting Information

Bulk structures, XRD patterns of the Nb2O5 catn class="Chemical">alyst, Charge density differences, H2O adsorption structures, Bader charges with the adsorption of CH3NH2 and H2O, and Transition state structures.

Conflict of interest

The authors declare no conflict of interest.

As a service to our authors and readers, this journal provides supporn class="Chemical">ting information supplied by the authors. Such materials are peer reviewed and may be re‐organized for online delivery, but are not copy‐edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors.

Supplementary

Click here for additionn class="Chemical">al data file.