Recently, sustainable and CO2-netual biomass has attracted much attention because it is the only renewable energy source that can be converted to liquid transportation fuels and chemicals [1]. Biomass can be readily converted to a liquid product, namely bio-oil, via a well-developed and commercially viable process, i.e., fast pyrolysis [2]. However, bio-oil cannot be directly used as fuel because its high oxygen content (usually 20-40 wt%) [2] results in undesired properties such as low thermal and chemical stabilities, a low heating value, high acidity, and immiscibility with fossil fuels [3]. Bio-oil therefore needs to be upgraded by removing oxygen to yield hydrocarbons compatible with petroleum-derived fuels [4, 5]. In particular, lignin-derived bio-oil, which consists mainly of phenolic compounds (e.g., phenol, anisole, guaiacol, cresol, and syringol), can be converted to benzene, toluene, and xylene (BTX) via hydrodeoxygenation (HDO) [4, 5]. BTX consists of basic chemicals that are currently mainly derived from petroleum and coal and also have high octane numbers as gasoline components [6, 7].
Selective cleavage of CAr-O bonds (CAr denotes a benzene ring carbon atom), i.e., direct deoxygenation (DDO), is necessary for HDO of phenolic compounds to yield aromatic hydrocarbons (especially BTX) [8], and the benzene ring should remain intact. A low H2 pressure and high temperature are thermodynamically favorable for producing BTX [9, 10], but the use of catalysts is crucial for obtaining high BTX yields.
Conventional sulfide catalysts have been used for HDO of phenolic compounds [11]. CoMo/Al2O3 is more favorable than NiMo/Al2O3 for DDO [12]. However, to prevent sulfide catalyst deactivation because of sulfur loss, sulfur-containing reagents (e.g., H2S) must be added to the reaction system, and products can be contaminated with sulfur during upgrading [4, 13]. The development of non-sulfide catalysts is therefore important. Several researchers have reported that molybdenum carbide [14], FeMoP [15], and Fe/SiO2 [16] show high selectivities for the DDO pathway in HDO of phenolic compounds. Many studies of HDO of phenolic compounds have focused on noble-metal and Ni catalysts. Noble-metal catalysts are highly active in HDO, but they facilitate saturation hydrogenation of the benzene ring [17, 18]. Moreover, their applications are limited because of their high cost and scarce resources. Low-cost monometallic Ni also shows excellent activity in HDO of phenolic compounds [9], but serious saturation hydrogenation of benzene rings occurs. Moreover, metallic Ni is highly active in C-C bond hydrogenolysis [19], resulting in a decreased carbon yield and increased H2 consumption. Much research has been performed on suppressing deep hydrogenation of benzene rings on noble metal and metallic Ni catalysts. Strategies based on the effect of metal-support interactions [20] and modification with a second metal promoter have been investigated.
Metal promoters can be divided into two types, according to their oxophilicity, i.e., oxophilic metals (e.g., Fe, Zn, Re, and Mo) and less oxophilic metals (e.g., Cu). Chia et al. [21] developed a new class of bifunctional heterogeneous catalysts, i.e., highly reducible Rh metal combined with the oxophilic metal Re, which facilitate selective C-O hydrogenolysis of biomass-derived feedstocks (such as cyclic ethers and polyols). Density functional theory (DFT) calculations have shown that in HDO of m-cresol on a bimetallic Pt-Mo catalyst [17], oxophilic Mo atoms favor geometries in which aromatic surface species have an upright orientation through favorable Mo-O interactions at the surface; this reduces the barrier to CAr-O bond scission. The addition of Fe to Pd [22] and Ni [23, 24] also helps to retain the benzene ring structure by selective cleavage of CAr-O bonds because the interactions between Fe and oxygen facilitate CAr-O bond activation. A DFT calculation [23] indicated that the presence of oxophilic Fe causes repulsion by the π-electron system of the benzene ring. Notably, a monometallic Ni catalyst modified with a less oxophilic metal, i.e., Cu, still had high activity in saturation hydrogenation of benzene rings [25, 26]; however, the presence of Cu suppressed C-C bond hydrogenolysis because of the geometric effect of Cu [27, 28].
Metallic In, which is much less active than metallic Ni in hydrogenation, has been reported to be an efficient catalyst for selective hydrogenation of C=O bonds to form C-OH [29-31]. Our previous work [32] showed that Ni-In alloys and intermetallic compounds are much less active than metallic Ni in the hydrogenation of ethylene and C-C bond hydrogenolysis. This is attributed to the geometric and electronic effects of In. On the basis of these results, we speculate that a metallic-In-modified Ni catalyst could preferentially yield BTX in HDO of phenolic compounds. To the best of our knowledge, there have been no reports of HDO of phenolic compounds on Ni-In bimetallic catalysts.
In the present work, we prepared supported bimetallic Ni-In catalysts with different Ni/In ratios and Ni contents. SiO2 was used as the support because of its inertness. Anisole, which has a typical CCH3-O-CAr phenolic structure [33, 34], was selected as a model compound for testing the HDO performances of bimetallic Ni-In/SiO2 catalysts. In thermodynamic terms, decreasing the reaction pressure is beneficial for BTX production [9, 35], therefore the catalytic performances were evaluated at atmospheric pressure in the present work. We found that the geometric and electronic effects of In in the bimetallic catalysts increased the selectivity for BTX and suppressed saturation hydrogenation of benzene rings, C-C bond hydrogenolysis, and methanation. This enhanced the carbon yield and reduced H2 consumption.
SiO2-supported monometallic Ni and Ni-In bimetallic catalysts were prepared using an incipient impregnation method. SiO2 was purchased from the Qingdao Haiyang Chemicals Co., Ltd. (Qingdao, China). First, SiO2 particles (diameter 150-250 μm) were incipiently impregnated with an aqueous solution of Ni(NO3)2 or a mixed aqueous solution of Ni(NO3)2 and In(NO3)3. The resulting samples were kept at room temperature for about 48 h, followed by drying at 120 ℃ for 12 h in a drying oven and calcination at 500 ℃ for 4 h in a muffle furnace to obtain the catalyst precursors in an oxidized state. Two impregnations were used to prepare catalyst precursors with Ni mass contents greater than 30%. Secondly, the precursors were reduced in a H2 flow on a quartz fixed-bed reactor at 450 ℃ for 3 h to prepare monometallic Ni and bimetallic Ni-In catalysts. Before the reduced catalysts were exposed to air for characterization using X-ray diffraction (XRD), transmission electron microscopy (TEM), X-ray photoelectron spectroscopy (XPS), and N2 sorption, they were passivated using a 0.5% O2/99.5% N2 flow (320 mL/min) at room temperature for 4 h.
Two series of catalysts were prepared using this procedure: (1) a SiO2-supported Ni catalyst with a Ni mass content of 40% (denoted by Ni/SiO2), and SiO2-supported bimetallic Ni-In catalysts with different Ni/In atomic ratios (y) and a Ni mass content of 40% (denoted by NiyIn/SiO2); (2) SiO2-supported bimetallic Ni-In catalysts with different Ni mass contents (x) and a Ni/In atomic ratio of 30:1, denoted by xNi30In/SiO2. In all the catalysts, the Ni mass content was set according to the equation mNi/(mNi + mSiO2).
H2 temperature-programmed reduction (H2-TPR) was performed using laboratory-made apparatus. The catalyst precursor (50 mg) was loaded into a U-tube quartz reactor (diameter 4.0 mm). The sample was kept at 50 ℃ in a 10% H2/90% N2 flow (60 mL/min) until the thermal conductivity detector (TCD) stabilized. The precursor was then reduced at a heating rate of 10 ℃ /min. The H2 consumption was determined using a TCD.
NH3 temperature-programmed desorption (NH3-TPD) was performed using the same apparatus as for H2-TPR. The catalyst precursor (150 mg) was loaded into the reactor and reduced by a H2 flow (60 mL/min) at 450 ℃ for 1 h. After the system had cooled to 100 ℃, NH3 adsorption was performed for 30 min. Physically adsorbed NH3 on the catalyst surface was purged using a He flow (40 mL/min) at 100 ℃; NH3-TPD was conducted at a heating rate of 15 ℃ /min in the He flow. The NH3 desorption signal was recorded using a TCD.
H2 chemisorption was also performed using the same apparatus as for H2-TPR. The catalyst precursor (100 mg) was loaded into the reactor and reduced in situ by a H2 flow (60 mL/min) at 450 ℃ for 1 h, followed by purging with a N2 flow (60 mL/min) at the same temperature (450 ℃) for 1 h to remove H2 adsorbed on the catalyst surface. The U-tube quartz reactor was then cooled to 30 ℃ in the N2 flow. When the TCD was stable, H2 (50 μL) was pulsed into the N2 flow until the effluent areas of consecutive pulses were constant. The total dynamic H2 uptake was calculated.
XRD patterns were recorded using a D8-Focus X-ray diffractometer (Bruker, Germany) with Cu Kα radiation (λ = 0.154 nm), operated at 40 kV and 40 mA. Powder samples were scanned from 10° to 90° at a rate of 8°/min.
N2 adsorption/desorption isotherms were recorded using a Quantachrome Quadrasorb SI instrument at -196 ℃ in liquid N2. Passivated catalysts were used. The specific surface area (ABET) was calculated using the multi-point Brunauer-Emmett-Teller (BET) method, and the pore volume (Vp) and average pore diameter (dp) were determined using the Barret-Joyner-Halenda method.
XPS was performed using a PHI5000 VersaProbe instrument with Al Kα radiation (1486.6 eV). The binding energy (BE) was corrected based on adventitious carbon (C 1s at 284.8 eV). The oxide layer on the passivated catalyst was removed by sputtering with an Ar+ ion beam (4 kV, 3.0 MPa, 1.5 min).
TEM and high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) mapping and line scanning were performed using a JEM-2100F instrument (operated at 2000 kV) equipped with an energy-dispersive X-ray spectroscopy (EDS) facility. The passivated catalyst was ground to a powder, ultrasonically dispersed in anhydrous ethanol, and then deposited on a Cu grid for characterization.
The catalytic performances in HDO of anisole were evaluated using a fixed-bed quartz-tube reactor (diameter 12 mm). An isothermal catalyst bed was maintained by diluting the catalyst precursor (400 mg) with quartz sand (3 g) of a uniform diameter. Quartz sand (diameter 380-830 μm, 2 g) was placed on the catalyst layer for feedstock preheating and dispersion. The catalyst precursor was reduced in situ by a H2 flow (100 mL/min) at 450 ℃ for 3 h. The catalyst was then cooled to 300 ℃. The H2 flow rate was adjusted to 13.81 mL/min and an octane solution containing 8 wt% anisole was fed into the reactor at a rate of 0.047 mL/min using a micro pump. The weight hourly space velocity (WHSV) of anisole and the H2/anisole molar ratio were 0.4 h-1 and 25, respectively. The liquid products were analyzed using an SP-3420 gas chromatography (GC) system equipped with a flame ionization detector (FID) and a commercial SE-30 capillary column (50 m × 0.32 mm × 3.0 μm); ethylbenzene was used as an internal standard. The C1 gaseous products (i.e., CO, CO2, and CH4) were detected using an on-line 102 GC equipped with a TCD and a TDX-101 packed column; N2 was used as an internal standard. The gaseous light hydrocarbon products (C1-C5) were detected using an on-line SP-3420 GC system equipped with a FID and an HP-AL/S capillary column (50 m × 0.535 mm × 15 μm).
The performances of the catalysts were described by the anisole conversion (X), selectivity for product i (Si), and BTX yield (YBTX), which were calculated using the following equations:
where nAnisole, 0 and nAnisole are the number of moles of anisole in the feedstock and products, respectively, ni is the number of moles of product i, and SBTX is the selectivity for BTX.
The turnover frequency (TOF, s-1) of anisole was calculated as
where X is the anisole conversion, W is the mass of catalyst used (g), F is the anisole feed rate (μmol/s), and M is the H2 uptake (μmol/g). This equation was used because the conversions in the present work were mostly higher than 50%, and even close to 100%; these are far from differential conditions, therefore an integral analysis was required [36]. This equation, in which -ln(1-X) substitutes for X, assumes a pseudo-first-order reaction, which is reasonable because of the large excess of H2 and the isothermal catalyst bed.
DFT calculations were performed using Gaussian 09. The Lee-Yang-Parr correlation functional (B3LYP) and the 6-31g(d, p) basis set were used; the molecular structure of anisole was optimized and the interactions between anisole and single Ni and In atoms via adsorption of the oxygen atom in anisole were calculated.
The adsorption energies of anisole on metallic Ni and In (M) atoms were calculated using the equation
where EM and EAnisole are the total energies of the M atom and free anisole, respectively, and EAnisole/M is the total energy of anisole bound to the M atom through an oxygen atom.
Fig. 1 shows the H2-TPR profiles for the Ni/SiO2 and Ni-In bimetallic catalyst precursors. As Fig. 1(a) shows, the Ni/SiO2 and Ni50In/SiO2 catalyst precursors gave similar curves, i.e., a main peak (at ~355 ℃) with a shoulder at higher temperature, attributed to reduction of NiO and nickel silicate, respectively [37, 38]. The main peak shifted from 355 to 460 ℃, and the shoulder also shifted to higher temperature, with decreasing Ni/In atomic ratio from 50:1 to 10:1, i.e., the introduction of In suppressed the reduction of Ni species. This could be because reduction of indium oxide is more difficult than NiO reduction [32]; it also indicates interactions between Ni and In species. The presence of In species could retard the diffusion of hydrogen species during H2-TPR because metallic In interacts less with hydrogen than metallic Ni does [29], making reduction difficult. The results also indicate that there were interactions between Ni and In species. It is worth noting that the initially reduced Ni species can promote the reduction of indium oxides via spilt-over hydrogen [32]. The co-reduction of Ni and In species could occur during TPR, which would lead to strong interactions between Ni and In in the reduced catalysts.
Fig. 1(b) shows that for the xNi30In/SiO2 catalyst precursors the main peak area increased but the shoulder peak area did not change obviously with increasing Ni content, i.e., an increase in the Ni content did not promote the formation of nickel silicate in the xNi30In/SiO2 precursors. This can be explained by nickel silicate formation on the SiO2 surface via a reaction between Ni2+ species and Si-OH groups [37], with the number of surface Si-OH groups on SiO2 remaining constant.
Fig. 2 shows the XRD patterns for the Ni/SiO2 and bimetallic Ni-In catalysts. In all the patterns, three diffraction peaks are visible, at ~44.5°, 51.8°, and 76.8°, assigned to the (111), (200), and (220) lattice planes, respectively, of metallic Ni (PDF 04-0850). For NiyIn/SiO2, the peaks from metallic Ni shift to slightly lower angles with decreasing Ni/In ratio (Fig. 2(a)). This shift to lower angles is not obvious because the Ni/In atomic ratios (50-10) are very high (i.e., the In content is very low). Our previous work indicated that this shift was more obvious when the Ni/In atomic ratio decreased to 6:1 and 4:1 [32]. This shift indicates incorporation of In atoms into the metallic Ni lattice, which is reasonable because the In atom (radius 0.166 nm) is larger than the Ni one (radius 0.124 nm) and In incorporation expands the metallic Ni lattice. The Ni-In bimetallic crystallite sizes for NiyIn/SiO2 were about 8.7 nm, regardless of the Ni/In ratio (Table 1), similar to the metallic Ni crystallite size (8.8 nm) for Ni/SiO2. This shows that introducing a small amount of In into Ni/SiO2 has a negligible effect on the crystallite size. For xNi30In/SiO2 (Fig. 2(b)), the diffraction peaks became stronger and sharper with increasing Ni content, and the crystallite size increased from 7.5 to 10.0 nm with increasing Ni content from 16 wt% to 50 wt% (Table 2).
Tables 1 and 2 show the textural properties of the catalysts. The NiyIn/SiO2 catalysts with different Ni/In atomic ratios had similar pore diameters, and the presence of In in NiyIn/SiO2 (apart from Ni40In/SiO2) led to a slight increase in the ABET and pore volume compared with those for Ni/SiO2 (Table 1). For xNi30In/SiO2, the ABET and pore volume decreased with increasing Ni content, but all the catalysts had similar average pore diameters (Table 2).
TEM was used to determine the distributions of Ni and In in Ni10In/SiO2. The TEM image in Fig. 3(a) shows that the Ni-In bimetallic particle size was mainly in the range 6-10 nm, with an average size of 8.5 nm; this is consistent with the calculated value (8.5 nm) obtained using the Scherrer equation. The ED spectrum of a random region (Fig. 3(a)) shows that the Ni/In atom ratio was 11.1, close to the nominal value of 10:1, indicating uniform distributions of Ni and In. This was confirmed by a HAADF-STEM mapping image of a random region (Fig. 3(b)). HAADF-STEM-EDS line scanning further confirmed that Ni and In were uniformly distributed in the bimetallic particles (Fig. 3(c)), indicating intimate contact between Ni and In. Because of this intimate contact, In atoms would isolate Ni ones, leading to a decrease in the size of Ni atom ensembles. This is also shown by the H2 uptakes.
H2 chemisorption was used to determine the surface density of Ni. As shown in Table 1, the density for Ni50In/SiO2, which had the lowest In content (i.e., the highest Ni/In ratio, 50:1), was much lower (38 μmol/g) than that for Ni/SiO2 (200 μmol/g). The H2 uptake decreased almost linearly from 38 to 6.7 μmol/g (Fig. 4) with decreasing Ni/In ratio from 50:1 to 10:1. For xNi30In/SiO2, the H2 uptake increased slightly from 21 to 25 μmol/g with increasing Ni content from 16 wt% to 50 wt%. Clearly, the introduction of In significantly decreased the surface Ni site density. The surface energy of metallic In (~560 mJ/m2 [39]) is lower than that of metallic Ni (~2700 mJ/m2 [40]), therefore In enrichment of the surfaces of Ni-In bimetallic crystallites can occur. The resulting dilution of Ni atoms by In ones limits the capacity of Ni to adsorb and activate H2. This indicates that In metal atoms are inactive in H2 adsorption, consistent with the inertness of metallic In in hydrogenation [29, 32]. The Ni content had almost no effect on the H2 uptakes by the xNi30In/SiO2 catalysts because they had the same Ni/In atomic ratios. We suggest that consecutive Ni atoms were separated by In atoms in the Ni-In bimetallic catalysts and the ensembles of surface Ni atoms became smaller with decreasing Ni/In ratio.
Fig. 5 shows the XPS spectra of Ni/SiO2, Ni30In/SiO2, and Ni10In/SiO2 in the Ni 2p3/2 and In 3d regions. As shown in Fig. 5(a), the Ni 2p3/2 BEs for Ni/SiO2, Ni30In/SiO2, and Ni10In/SiO2 are similar (~852.6 eV), and are ascribed to metallic Ni0 [41, 42]. This differs from a report that charge transfer from In to Ni occurred in Ni-In intermetallic compounds, shown by extended X-ray absorption fine structure results [29]. Charge transfer from In to Ni should be reflected in a decrease in the Ni 2p3/2 BE. This discrepancy may arise because of the small amount of metallic In in the Ni-In bimetallic catalysts in our study. Our previous work shows that the BEs of Ni 2p3/2 for metallic Ni and Ni2In are 852.8 and 852.3 eV, respectively [32]. The In 3d BEs provide information on the electronic interactions between metallic Ni and In. It has been reported that the In 3d5/2 and In 3d3/2 BEs for metallic In0 are 444.0 and 451.4 eV [43, 44], respectively. Here, the In 3d5/2 (444.5-444.6 eV) and In 3d3/2 (452.0-452.2 eV) BEs for Ni30In/SiO2 and Ni10In/SiO2 are higher than those for metallic In0 (Fig. 5(b)), indicating charge transfer from In to Ni. This is in agreement with the higher electronegativity of Ni (1.91) compared with that of In (1.78) [29]. There is therefore a higher electron density of Ni atoms in the Ni-In bimetallic catalysts than in the monometallic Ni one, and the In species in the Ni-In bimetallic catalysts are positively charged. Charge transfer from In to Ni also indicates strong interactions between Ni and In on the catalyst surface. Additionally, the surface Ni/In ratios for Ni30In/SiO2 and Ni10In/SiO2 were 8.2 and 3.3, respectively, i.e., lower than the nominal ones (30:1 and 10:1, respectively) in the bulk catalysts. This also indicates In enrichment on the catalyst surface, consistent with the H2 chemisorption results.
It has been reported that only Lewis acidic sites are present on Ni/SiO2 and Ni-In/SiO2 [32]. For Ni/SiO2, unreduced Ni species (especially nickel silicate) contribute to the Lewis acidity [45, 46]; In species with a small positive charge can also act as Lewis acidic sites.
Fig. 6 shows the NH3-TPD profiles for the catalysts. Ni/SiO2 and NiyIn/SiO2 (Fig. 6(a)) gave two main NH3 desorption peaks, at 180 and above 330 ℃, and a shoulder at higher temperature. The second peak shifted to higher temperature with increasing In content (i.e., decreasing Ni/In ratio). There was also another weak peak at ~700 ℃ for NiyIn/SiO2 when the Ni/In atomic ratio was less than 20:1. xNi30In/SiO2 also gave two main NH3 desorption peaks at around 180 and 330 ℃, and the shoulder became more obvious with increasing Ni content (Fig. 6(b)). The total number of acidic sites on NiyIn/SiO2 (especially strong acidic sites) was lower than that on Ni/SiO2, and this became more obvious with decreasing Ni/In ratio. The amount of acidic sites on xNi30In/SiO2 increased with increasing Ni content. This is reasonable because Lewis acidic sites are mainly derived from Ni species, therefore a higher Ni content leads to an increased number of acidic sites.
In HDO of anisole on different catalysts, apart from anisole and n-octane (i.e., the solvent), the compounds detected in the liquid effluent included benzene, toluene, xylene, cyclohexane, phenol, cresol, and n-hexane. No cyclohexanol, cyclohexanone, and methoxycyclohexane were found. In the gaseous effluent, C1-C5 hydrocarbons were detected. According to the literature [4, 25, 47-49], two parallel routes can occur in HDO of anisole (Scheme 1). One is DDO, which involves direct cleavage of CAr-O bonds (the area marked by a dashed line in Scheme 1). Demethoxylation of anisole can occur to form benzene; demethylation and methyl transfer can also initially occur to yield phenol, cresol, and xylenol, which can be converted to benzene, toluene, and xylene, respectively, via CAr-O bond cleavage. Benzene can be further hydrogenated to cyclohexane. The other route is hydrogenation-deoxygenation (HYD), in which benzene ring hydrogenation yields methoxycyclohexane and cyclohexanol, followed by demethoxylation-hydrogenation and dehydration-hydrogenation to form cyclohexane. Ring opening of cyclohexane via C-C bond hydrogenolysis gives n-hexane, and then further C-C bond hydrogenolysis yields C1-C5 hydrocarbons. As discussed below, in the present work, benzene was the main product. Given that methoxycyclohexane and cyclohexanol were not detected, the DDO route may be more favorable than HYD at atmospheric pressure [35]. The CArO-CH3, CAr-OH, and CAr-OCH3 bond dissociation energies are 339, 468, and 422 kJ/mol, respectively [4, 35], therefore the CArO-CH3 bond is the most easily broken to produce phenol, followed by the CAr-OCH3 and CAr-OH bonds. In thermodynamic terms, demethoxylation of anisole to yield benzene is more favorable than conversion of phenol to benzene.
Additionally, as indicated below, n-heptane was predominantly generated by hydrogenolysis of n-octane (i.e., the solvent).
As shown in Fig. 7(a), the introduction of In led to a decrease in anisole conversion. The anisole conversion on Ni/SiO2 was 99.3%, and it gradually decreased to 55.7% as the Ni/In atomic ratio decreased to 10:1. This is primarily ascribed to the large decrease in the surface density of Ni, as indicated by the H2 uptake results (Table 1). Metallic In is inert in hydrogenation [29, 32]. The decrease in the surface density of Ni with decreasing Ni/In atomic ratio (Table 1) is not favorable for anisole conversion. The presence of In enhances the anisole TOF. The data in Table 1 show that the TOF for Ni/SiO2 was 1.3 × 10-2 s-1, whereas the values were (5.3-7.2) × 10-2 s-1 for NiyIn/SiO2 with Ni/In ratios between 50:1 and 10:1. This suggests that although metallic In is inert in hydrogenation, it may be active in anisole conversion.
Here, the target products are BTX, and the introduction of a suitable amount of In enhanced the selectivity for BTX. The selectivity for BTX was 52.0% on Ni/SiO2, and increased to 64.0% with decreasing Ni/In atomic ratio to 30:1, followed by a drop to 49.8% at a Ni/In ratio of 10:1. The BTX yield reached the maximum, i.e., 60.4%, at a Ni/In ratio of 40:1 (Fig. 7(a)). The total selectivity for phenol and cresol was low, and increased from 0.8% to 3.1% when the Ni/In atomic ratio decreased from 1:0 to 10:1.
Benzene was the main BTX component produced, and the trend in its selectivity changes was the same as that for BTX (45.9%, 58.6%, 60.0%, 62.2%, 60.0%, and 47.6% at Ni/In atomic ratios of 1:0, 50:1, 40:1, 30:1, 20:1, and 10:1, respectively). Toluene and xylene were formed by methyl transfer during the reaction. The total selectivity for toluene and xylene decreased from 6.0% to 2.1% with decreasing Ni/In ratio from 1:0 to 10:1, showing that the presence of In restrained the methyl-transfer reaction. This is also clearly shown by the molar ratio of benzene to toluene and xylene combined (denoted by nB/nTX). The nB/nTX ratio was 27.7-21.2 at Ni/In ratios between 50:1 and 10:1, much larger than that (7.7) for Ni/SiO2. The reduced activity in methyl transfer can be explained by a decrease in the total amount of acidic sites because these sites are involved in methyl transfer [10, 36, 50].
Fig. 7(b) shows the selectivities for cyclohexane and n-hexane. The selectivity for cyclohexane slowly increased from 1.6% to 8.4% with decreasing Ni/In atom ratio from 1:0 to 10:1; conversely, the selectivity for n-hexane decreased from 22.5% to 0.11%. As mentioned above, n-hexane and n-heptane were both detected in the liquid effluent. The selectivity for n-heptane (calculated on the basis of consumed anisole) and total selectivity for C6-C8 compounds also decreased with decreasing Ni/In atomic ratio. The total selectivities for C6-C8 compounds exceeded 100% on Ni/SiO2 and Ni50In/SiO2, which is impossible. Reasonably, n-hexane can be derived from the ring opening of cyclohexane. n-Heptane could be generated by hydrogenolysis of specific C-C bonds in methylcyclohexane; however, no methylcyclohexane was detected in the liquid effluent. Given the high activity of metallic Ni in hydrogenolysis of C-C bonds [9, 19], n-heptane might be produced by hydrogenolysis of n-octane (i.e., the solvent). To verify this, n-octane alone was pumped into the reactor to evaluate the activities of Ni/SiO2, Ni30In/SiO2, and Ni10In/SiO2; the n-octane conversions were 42.4%, 24.5%, and 6.89%, respectively, and C1-C7 hydrocarbons were detected. This shows that hydrogenolysis of n-octane occurred on Ni/SiO2 and NiyIn/SiO2, and the activity of metallic Ni in C-C bond hydrogenolysis was significantly suppressed by the introduction of In. Because n-hexane could also be produced by hydrogenolysis of the solvent (n-octane), the selectivity for n-hexane calculated on the basis of the anisole conversion is not an accurate value. The selectivity for n-heptane is also not the real value because it was calculated on the basis of anisole consumption. Nevertheless, the selectivity for n-hexane can be used to reflect the catalytic activity in C-C bond hydrogenolysis, i.e., a higher selectivity for n-hexane indicates a higher activity in C-C bond hydrogenolysis. In short, the activity of Ni/SiO2 in the hydrogenolysis of C-C bonds was suppressed by the introduction of metallic In, and the higher the In content was, the greater the suppression. This low activity in C-C bond hydrogenolysis increases the selectivity for cyclohexane.
The above results indicate that the introduction of In into Ni/SiO2 increased the selectivities for BTX and cyclohexane. Given the large difference between the anisole conversions on Ni/SiO2 and NiyIn/SiO2 under the same reaction conditions, this comparison of product selectivities may be unsatisfactory. We therefore also compared the selectivities for BTX and cyclohexane at similar anisole conversions. Similar anisole conversions of ~60.0% on Ni/SiO2 and NiyIn/SiO2 were obtained by changing the anisole WHSV. The selectivity for BTX was 71.3% on Ni/SiO2, about 2%-5% lower than those on NiyIn/SiO2 (apart from Ni10In/SiO2). We suggest that the presence of In suppressed benzene ring hydrogenation, resulting in higher selectivities for BTX on NiyIn/SiO2. This is confirmed by the lower activity of NiyIn/SiO2 compared with that of Ni/SiO2 in the hydrogenation of benzene. Under the same conditions (300 ℃, 0.1 MPa, H2/benzene molar ratio 17.4, and benzene WHSV 0.55 h-1), Ni/SiO2, Ni30In/SiO2, Ni20In/SiO2, and Ni10In/SiO2 gave benzene conversions of 65.4%, 39.3%, 22.9%, and 17.0%, respectively. Clearly, the activity of NiyIn/SiO2 decreased with increasing In content. Additionally, at similar anisole conversions, i.e., ~60%, the selectivity for cyclohexane on Ni/SiO2 was 0.7, and it increased from 3.2 to 9.9 with decreasing Ni/In ratio from 50:1 to 10:1. This is consistent with a lower activity of NiyIn/SiO2 in C-C bond hydrogenolysis, as indicated above.
The molar ratio of CH4 to converted anisole (denoted by nCH4/nΔAnisole) also reflects the catalytic activity in C-C bond hydrogenolysis. CH4 was the main gaseous product on NiyIn/SiO2 and can be derived from the methoxy group in anisole via the following pathways: (1) ArO-CH3 bond hydrogenolysis; (2) demethoxylation (i.e., Ar-OCH3 bond cleavage) to form methanol followed by hydrogenation; and (3) methanation of CO derived from methanol decomposition. CH4 can also be produced by C-C hydrogenolysis of n-hexane (from cyclohexane) and n-octane (i.e., the solvent). If C-C hydrogenolysis does not occur, the nCH4/nΔAnisole ratio should be ≤ 1.0. As shown in Fig. 7(c), Ni/SiO2 gave an nCH4/nΔAnisole ratio of 13.3, and In introduction significantly decreased the nCH4/nΔAnisole ratio. The nCH4/nΔAnisole ratio decreased from 3.6 to 0.5 with decreasing Ni/In ratio from 50:1 to 10:1. CO was not detected on Ni/SiO2, but was generated on NiyIn/SiO2. The molar ratio of CO to consumed anisole (nCO/nΔAnisole) increased from 0.1 to 0.3 with decreasing Ni/In atomic ratio from 50:1 to 10:1 (Fig. 7(c)), i.e., the Ni-In bimetallic catalysts were less active than the monometallic Ni catalyst in methanation of CO.
The effects of the Ni content on the performances of Ni-In bimetallic catalysts in HDO of anisole were investigated using xNi30In/SiO2 catalysts with various Ni contents.
Fig. 8(a) shows that the anisole conversion increased with increasing Ni content, and exceeded 90% on xNi30In/SiO2 with a Ni content ≥ 30 wt%. This is because of the increased surface density of Ni (Table 2). The TOFs for anisole on xNi30In/SiO2 were also calculated. The data in Table 2 show that the TOF was 2.8 × 102 s-1 on 16%Ni30In/SiO2, and ~5.0 × 10-2 s-1 on xNi30In/SiO2 with Ni contents between (24 and 50) wt%. The Ni-In bimetallic crystallite size did not affect anisole conversion when the Ni-In bimetallic crystallites were larger than 8.0 nm.
The most abundant products were still BTX; their selectivities and yields on xNi30In/SiO2 are shown in Fig. 8(a). For xNi30In/SiO2, the BTX selectivity remained at ~ 61%, and the BTX yield increased from 42.6% to 59.0% with increasing Ni content from 16 to 40 wt%. The Ni content therefore did not affect the selectivity for BTX, indicating that the BTX selectivity was determined by the Ni/In atomic ratio rather than the Ni content. The BTX yield increased with increasing Ni content because of the increase in the anisole conversion. Small amounts of phenol and cresol were detected in the liquid effluent, and the total selectivity for phenol and cresol was low (ranging from 0% to 2%).
Benzene was the main BTX component produced on xNi30In/SiO2. The selectivity for benzene on xNi30In/SiO2 did not change significantly, and ranged from 59% to 62%. The total selectivity for toluene and xylene on xNi30In/SiO2 was low (ranging from 1.5% to 1.9%), indicating that xNi30In/SiO2 had low activity in the methyl-transfer reaction. This is also clearly shown by the nB/nTX ratio on xNi30In/SiO2; the ratio remained at ~35 with increasing Ni content from 16 wt% to 40 wt%, which is much higher than that (7.7) on Ni/SiO2. Consistent with the above results, for NiyIn/SiO2, the introduction of In has a negative effect on the methyl-transfer reaction, mainly because of the decreased amount of acidic sites.
Fig. 8(a) shows that the selectivity for cyclohexane on xNi30In/SiO2 was maintained at around 5.0%. n-Hexane and n-heptane were always detected in the liquid products. The selectivity for n-hexane was ~2% on xNi30In/SiO2. CH4 was the main gaseous product on xNi30In/SiO2. Fig. 8(b) shows that the nCH4/nΔAnisole ratios were between 1.2 and 1.5 on xNi30In/SiO2, much lower than that (13.3) on Ni/SiO2, further indicating much lower activities in C-C bond hydrogenolysis on the Ni-In bimetallic catalysts.
As mentioned above, CO can be formed by methanol decomposition. CO was always detected on xNi30In/SiO2; nCO/nΔAnisole was 0.14-0.18 on xNi30In/SiO2 (Fig. 8(b)). Clearly, the Ni content had little effect on the Ni-In bimetallic catalytic activity in CO methanation.
We also investigated the performances of Ni/SiO2 with various Ni contents (denoted by xNi/SiO2). The conversion was 90.6% at a Ni mass content of 16% and ~99.5% for Ni mass contents between 24% and 50%. At the same Ni content, the anisole conversion on xNi30In/SiO2 was lower than that on xNi/SiO2. This is because of the smaller H2 uptake (21-25 μmol/g) by xNi30In/SiO2). For xNi/SiO2, the H2 uptakes were 44, 94, 138, 200, and 244 μmol/g for Ni mass contents of 16%, 24%, 30%, 40%, and 50%, respectively. However, xNi30In/SiO2 gave higher selectivities for BTX and cyclohexane in all cases. The selectivity for BTX on xNi/SiO2 ranged from 51.5% to 59.1%, and that for cyclohexane ranged from 1.1% to 1.6%.
In summary, the introduction of a small amount of metallic In suppressed the activity of Ni/SiO2 in benzene ring hydrogenation, C-C bond hydrogenolysis, and methanation. This is important in terms of increasing the carbon yield and reducing H2 consumption. Additionally, the effect of the Ni/In atomic ratio on the product selectivities was greater than that of the Ni content.
As mentioned above, the introduction of In into Ni/SiO2 has the following effects. (1) It decreases the catalytic activity in anisole conversion and increases the selectivity for BTX (mostly benzene), i.e., it decreases the activity in benzene ring hydrogenation. (2) It increases the selectivity for cyclohexane and decreases the selectivity for n-hexane, i.e., it decreases the catalytic activity in C-C bond hydrogenolysis. (3) It decreases the catalytic activity in methanation. All these effects are related to geometric and electronic modification of metallic Ni by metallic In. Specifically, ensembles of Ni atoms become smaller, and the electron density on Ni is increased by charge transfer from In.
The H2 uptakes indicate that the surface Ni site density decreased significantly when In was introduced into Ni/SiO2 (Tables 1 and 2). This is unfavorable for anisole conversion because metallic Ni is active in hydrogenation, whereas metallic In is relatively inert [29]. However, the introduction of metallic In enhanced the selectivity for BTX. We speculate that there are two reasons for these results. One is that the selectivities for CAr-O bond cleavage of the Ni-In bimetallic catalysts could be higher than that of the monometallic Ni catalyst. The other is that Ni-In bimetallic catalysts might have lower activities in benzene ring hydrogenation. In the first case, given that reduction of indium oxide is more difficult than NiO reduction [29], metallic In is more oxophilic than metallic Ni. It has been reported that when more oxophilic Fe and Zn were introduced into Ni [51], Pd [22], or Pt [52] catalysts, they facilitated selective cleavage of CAr-O bonds to produce benzene series hydrocarbons in HDO of m-cresol, guaiacol, and anisole; this is because of the strong interactions between Fe or Zn and oxygen. Additionally, Ni-In intermetallic compounds have high selectivities for hydrogenation of C=O bonds in α, β-unsaturated aldehydes [29] and fatty acids [30, 31], producing C-OH groups. In this study, DFT calculations using Gaussian 09 were used to determine the interactions between metallic Ni or In atoms and anisole via oxygen adsorption. The results show that the adsorption energies of anisole on metallic Ni and In atoms were -86.9 and 5.80 kJ/mol, respectively, i.e., there are strong interactions between metallic Ni and anisole via adsorption of oxygen, but spontaneous adsorption of anisole on In atoms via the anisole oxygen is not possible. The data in Table 3 show that the CAr-O and CCH3-O (CCH3 denotes the carbon in methyl) bonds in anisole both lengthen when anisole is adsorbed on a Ni atom, whereas their lengths are almost unchanged when anisole is adsorbed on an In atom. These results indicate that the presence of In does not favor anisole adsorption via oxygen, and activation and cleavage of the CAr-O bond, and that metallic Ni atoms are the predominant active sites for CAr-O bond cleavage. The higher selectivities for BTX on Ni-In bimetallic catalysts compared with that on Ni/SiO2 can therefore be mainly attributed to suppression of benzene ring hydrogenation, as discussed below.
Metallic Ni is highly active in benzene ring hydrogenation because of its high affinity for both H2 and benzene rings [23]. Ni-In bimetallic catalysts have lower H2 adsorption abilities as a result of In incorporation (Table 1). Charge transfer from In to Ni leads to increased electron density on the Ni atoms, which is unfavorable for the adsorption of electron-rich benzene rings [53-55]. The dilution of surface Ni atoms by In atoms may not favor adsorption and hydrogenation of benzene rings. As discussed in section 3.2.1, the activities of the Ni-In bimetallic catalysts in benzene hydrogenation were lower than that of Ni/SiO2, and the lower the Ni/In atomic ratio was, the lower the activity was. Additionally, it has been reported that, compared with monometallic Ni/SiO2, Ni-In/SiO2 bimetallic catalysts suppress the adsorption and hydrogenation of electron-rich ethylene [32, 56].
The introduction of In suppressed benzene hydrogenation, therefore the higher selectivity for cyclohexane on the In-modified catalysts compared with that on Ni/SiO2 is puzzling. This can be explained by greater C-C bond hydrogenolysis on monometallic Ni than on the Ni-In bimetallic catalysts, i.e., ring opening of cyclohexane is easier on metallic Ni, therefore cyclohexane is consumed. The lower activities in C-C bond hydrogenolysis of the Ni-In bimetallic catalysts compared with that of monometallic Ni is mainly ascribed to a decrease in the size of Ni atom ensembles. It has been widely reported that C-C bond cleavage requires a large ensemble of contiguous Ni atoms [25, 27, 57]. The mechanism of ethane hydrogenolysis suggests that co-adsorption of two adjacent atoms is needed to activate and facilitate C-C bond cleavage [57]. Here, although the monometallic Ni and Ni-In bimetallic crystallites in Ni/SiO2 and NiyIn/SiO2 were similar in size (Tables 1 and 2), the incorporation of In into the metallic Ni lattice decreased the size of the ensembles of Ni atoms. The H2 uptakes (Tables 1 and 2, and Fig. 4), suggest that the Ni atom ensembles became smaller with decreasing Ni/In ratio. Consequently, the selectivity for cyclohexane increased with decreasing Ni/In ratio, and the selectivity for n-hexane and the nCH4/nΔAnisole ratio decreased (Fig. 7). Low activities of Ni-In alloys and intermetallic compounds in C-C bond hydrogenolysis have also been reported for HDO of acetic acid and selective hydrogenation of acetylene [31, 32].
In anisole HDO, hydrogenolysis of CAr-OCH3 bonds (i.e., demethoxylation) produces methanol, which can be either hydrogenated to CH4 or decomposed to CO. Given that the Ni-In bimetallic catalysts were less active than monometallic Ni, methanol hydrogenation would be suppressed by the presence of In. It is reasonable that CO was not formed on Ni/SiO2 because metallic Ni is highly active in methanation [58]. CO was detected on the Ni-In bimetallic catalysts, indicating that bimetallic Ni-In was less active than monometallic Ni in methanation. Charge transfer from In to Ni may favor methanation of CO because the increased electron density on Ni increases the strength of CO adsorption and therefore activation of CO [59, 60]. However, the decrease in the size of Ni atom ensembles is unfavorable for methanation of CO [61]. Here, the geometric effect of In may be dominant in the methanation performance.
For NiyIn/SiO2 catalysts, the H2 uptake decreased linearly from 38 to 6.7 μmol/g with decreasing Ni/In atomic ratio from 50:1 to 10:1; these values are much lower than that for Ni/SiO2, namely 200 μmol/g. The H2 uptake increased slightly from 22 to 25 μmol/g for xNi30In/SiO2 with increasing Ni content from 16 wt% to 50 wt%. The monometallic Ni crystallites in Ni/SiO2 and the Ni-In bimetallic crystallites in NiyIn/SiO2 were similar in size (~8.7 nm), therefore dilution of Ni atoms by In atoms greatly suppressed H2 chemisorption, and a higher In content led to a lower H2 uptake. XPS showed that charge transfer from In to Ni occurred in the Ni-In bimetallic catalysts. There were strong interactions between metallic Ni and In, and metallic In geometrically and electronically modified metallic Ni in the Ni-In bimetallic catalysts. Additionally, the number of acidic sites on the Ni-In bimetallic catalysts was lower than that on Ni/SiO2.
In anisole HDO, although the anisole conversions on the Ni-In bimetallic catalysts were lower than that on monometallic Ni, because of the lower density of surface Ni sites on the bimetallic catalysts, the selectivities for BTX (mainly benzene) and cyclohexane were higher as a result of suppression of benzene ring hydrogenation and C-C bond hydrogenolysis. Additionally, bimetallic Ni-In was less active than monometallic Ni in methyl transfer because of the smaller number of acidic sites. The selectivity for BTX increased with decreasing Ni/In ratio and was maximum at a Ni/In ratio of 30:1. However, the Ni content had little effect on the selectivity for BTX. At 300 ℃, 0.1 MPa, an anisole WHSV of 0.4 h-1, and a H2/anisole ratio of 25, a BTX yield of 60.4% was obtained on Ni40In/SiO2 with a Ni content of 40%, i.e., higher than that (51.6%) obtained on Ni/SiO2.
We greatly appreciate Dr. Yifei Chen and Prof. Haoxi Jiang for the DFT calculation and discussion.