On Catalytic Behavior of Bulk Mo 2 C in the Hydrodenitrogenation of Indole over a Wide Range of Conversion Thereof

: The catalytic activity of bulk molybdenum carbide (Mo 2 C) in the hydrodenitrogenation (HDN) of indole was studied. The catalyst was synthesized using a temperature-programmed reaction of the respective oxide precursor (MoO 3 ) with the carburizing gas mixture of 10 vol.% CH 4 / H 2 . The resultant material was characterized using X-ray di ﬀ raction, CO chemisorption, and nitrogen adsorption. The catalytic activity was studied in the HDN of indole over a wide range of conversion thereof and in the presence of a low amount of sulfur (50 ppm), which was used to simulate the processing of real petroleum intermediates. The molybdenum carbide has shown high activity under the tested operating conditions. Apparently, the bulk molybdenum carbide turned out to be selective towards the formation of aromatic products such as ethylbenzene, toluene, and benzene. The main products of HDN were ethylbenzene and ethylcyclohexane. After 99% conversion of indole HDN was reached (i.e., lack of N-containing compounds in the products was observed), the hydrogenation of ethylbenzene to ethylcyclohexane took place. Thus, the catalytic behavior of bulk molybdenum carbide for the HDN of indole is completely di ﬀ erent compared to previously studied sulﬁde-based systems.


Introduction
The petroleum industry is under constantly increasing pressure exerted by legislators in terms of improving the quality of diesel fuel towards reducing the toxicity of its exhaust emissions. In recent decades, the Environmental Protection Agency (EPA) has established several regulations targeted at the lowering of sulfur emission in diesel fuel from 50 down to 10 ppm [1][2][3]. Indeed, the decrease in sulfur content in diesel fuels has an additional advantage-it allows the use of the catalytic converters in the engines so that nitrogen oxide emissions may also be reduced. However, it is well known that organic compounds containing nitrogen are inhibitors of the hydrodesulfurization (HDS) process. That is why the permissible content of such N-compounds should also be reduced to a level below

Characterization of the Catalyst
The specific surface area (S g ) of the fresh passivated Mo 2 C was determined to be 33 m 2 g −1 .
The XRD patterns (Figure 1) collected for a passivated catalyst before and after HDN reaction may be assigned to the lattice planes of the hexagonal compact phase (hcp) β-Mo 2 C structures.
The sharp peak at 2θ = 36 • of a weak intensity may be assigned to the characteristic diffraction line of the residues of silicon carbide SiC (the reflection ascribed to (1 1 1) Miller index), which was used to dilute the catalyst in the catalytic bed (see Section 3.4). As seen, no alteration of the crystallographic structure of the Mo 2 C catalyst was observed after the test of HDN of indole. The experimental CO uptake of the fresh sample was equal to 142 µmol g −1 and decreased down to 7.69 µmol g −1 after the HDN reaction. This change may be due to the coking of the catalyst's surface. On the other hand, it should be noted that the catalyst was exposed to contact with air after it was removed from the reactor for the purpose of studying the phase composition (XRD). This, in turn, could lead to passivation and a significant reduction in the number of active sites adsorbing CO. For this reason, care must be taken when analyzing the result of CO sorption for the spent catalyst. More details are provided in our former work [7]. Catalysts 2020, 10 The sharp peak at 2θ = 36° of a weak intensity may be assigned to the characteristic diffraction line of the residues of silicon carbide SiC (the reflection ascribed to (1 1 1) Miller index), which was used to dilute the catalyst in the catalytic bed (see Section 3.4). As seen, no alteration of the crystallographic structure of the Mo2C catalyst was observed after the test of HDN of indole. The experimental CO uptake of the fresh sample was equal to 142 μmol g −1 and decreased down to 7.69 μmol g −1 after the HDN reaction. This change may be due to the coking of the catalyst's surface. On the other hand, it should be noted that the catalyst was exposed to contact with air after it was removed from the reactor for the purpose of studying the phase composition (XRD). This, in turn, could lead to passivation and a significant reduction in the number of active sites adsorbing CO. For this reason, care must be taken when analyzing the result of CO sorption for the spent catalyst. More details are provided in our former work [7].

HDN of Indole over Mo2C
The possible reaction scenarios of indole HDN over the bulk Mo2C with all detected compounds are presented in Scheme 1.

HDN of Indole over Mo 2 C
The possible reaction scenarios of indole HDN over the bulk Mo 2 C with all detected compounds are presented in Scheme 1. The sharp peak at 2θ = 36° of a weak intensity may be assigned to the characteristic diffraction line of the residues of silicon carbide SiC (the reflection ascribed to (1 1 1) Miller index), which was used to dilute the catalyst in the catalytic bed (see Section 3.4). As seen, no alteration of the crystallographic structure of the Mo2C catalyst was observed after the test of HDN of indole. The experimental CO uptake of the fresh sample was equal to 142 μmol g −1 and decreased down to 7.69 μmol g −1 after the HDN reaction. This change may be due to the coking of the catalyst's surface. On the other hand, it should be noted that the catalyst was exposed to contact with air after it was removed from the reactor for the purpose of studying the phase composition (XRD). This, in turn, could lead to passivation and a significant reduction in the number of active sites adsorbing CO. For this reason, care must be taken when analyzing the result of CO sorption for the spent catalyst. More details are provided in our former work [7].

HDN of Indole over Mo2C
The possible reaction scenarios of indole HDN over the bulk Mo2C with all detected compounds are presented in Scheme 1.   As seen, notwithstanding the further reaction pathways, the first step thereof is the hydrogenation of indole towards indoline, followed by the opening of the N-containing ring to o-ethylaniline (OEA). The latter one can undergo either direct denitrogenation (called DDN route or route #1 in this work) to the ethylbenzene or full hydrogenation of the ring before denitrogenation (called HYD or route #2) to o-ethylaniline. The third possibility is OEA dealkylation to aniline (route #3).
It is pertinent to mention that in our research, no intermediate compound such as octahydroindole (OHI) was detected in the reaction product mixture. This is in agreement with the results presented by Bunch et al. [53] for the HDN of indole over sulfides as hydrotreating catalysts. They stated that in the presence of H 2 S, indole undergoes transformation to OEA because full hydrogenation of indoline into OHI was not observed. It may be concluded that in case of using molybdenum carbide, the preferred route is the formation of OEA. Benzene present in the products originates from the subsequent dealkylation of ethylbenzene (EB) and toluene (T), as well as from the hydrogenolysis of C (sp2) -N of aniline. Similar results were reported by Li et al. [8], who carried out the indole HDN process over Mo carbides and nitrides.
The presence of traces of ethylcyclohexene (ECHen) for the tests with short contact time undoubtedly indicates that it was formed by the extraction of a hydrogen atom from the carbon atom in the β position of the o-ethylcyclohexylamine molecule (OECHA). Indeed, the hydrogenolysis of the C (sp3) -N bond of OECHA is significantly easier compared with the C (sp2) -N bond scission of OEA.
The conversion of OEA occurs following two main competitive reaction routes. The first one (route #1) results in the formation of ethylbenzene (EB), toluene (T), and benzene (B) (hydrogenolysis of C (sp2) -N bond resulting in the formation of ammonia and EB). Route #2 leads to the formation of OECHA, ECHen, and ethylcyclohexane (ECH) (hydrogenation of aromatic ring to OECHA) followed by hydrogenolysis of C (sp3) -N bond resulting in the elimination of NH 3 and formation of ethylcyclohexene (ECHen) and, ultimately, ethylcyclohexane (ECH) [47,56]. The aforementioned hydrogenolysis of C-N bond in OEA molecules occurs according to the β-elimination mechanism (Hoffman type), although it may also follow the alternative mechanism of nucleophilic substitution that would result in the same products [63]. The formation of small amounts of benzene and toluene is the result of hydrogenolysis of the C-C bond, which takes place on the metallic centers of molybdenum carbide [47].
The hydrogenation of EB towards ECH is inhibited because of the presence of nitrogen-containing compounds being stronger adsorbed on the active sites. Since measurable amounts of aniline (A) were also observed in the reaction products, an additional pathway of OEA dealkylation (route #3 in Scheme 1) could not be excluded. Aniline is formed after the hydrogenolysis of the C-C bond between the ethyl and phenyl group in the OEA molecule. However, for the sake of simplicity and taking into account the absence of OHI and the negligible amounts of aniline (<0.5%), only two pathways for the HDN of indole were considered (routes #1 and #2) as shown below and discussed hereinafter. Table 1 displays the shares of indole, all the intermediates and final products, as well as total conversion, and HDN conversion over bulk Mo 2 C for a selected contact time of 0.37 s.
As seen, the main reaction products are ethylbenzene, ethylcyclohexane, and compounds resulting from ethylbenzene dealkylation, i.e., toluene and benzene. Moreover, nitrogen-containing compounds such as indoline, o-ethylaniline, and aniline were also present in the products mixture. Additionally, traces of ethylcyclohexene have been detected. The data collected in Table 1 as well as in Figure 2 indicate that molybdenum carbide shows relatively high conversion under the tested operating conditions. Indeed, for a contact time of 0.37 s, the HDN conversion (i.e., degree of HDN) over Mo 2 C reached 79%, while the total conversion of indole was equal to 99%.
The selectivity towards route #1 (DDN route) and #2 (HYD route) was calculated for different conversions and is presented in Table 2. Interestingly, the favored reaction pathway strongly depends on the conversion of indole (and, consequently, on contact time). Namely, in the case of low conversions (<33%), the selectivity data indicate that the HYD route (route #2) is favored, whereas for the highest conversion, the DDN route (route #1) takes over with a selectivity of 1.4 at 99% conversion. With this, it appears evident that for the highest conversion, the DDN scenario is the preferred one. As seen, the main reaction products are ethylbenzene, ethylcyclohexane, and compou resulting from ethylbenzene dealkylation, i.e., toluene and benzene. Moreover, nitrogen-contain compounds such as indoline, o-ethylaniline, and aniline were also present in the products mixt Additionally, traces of ethylcyclohexene have been detected. The data collected in Table 1 as we in Figure 2 indicate that molybdenum carbide shows relatively high conversion under the te operating conditions. Indeed, for a contact time of 0.37 s, the HDN conversion (i.e., degree of HD over Mo2C reached 79%, while the total conversion of indole was equal to 99%. The selectivity towards route #1 (DDN route) and #2 (HYD route) was calculated for diffe conversions and is presented in Table 2. Interestingly, the favored reaction pathway strongly depe on the conversion of indole (and, consequently, on contact time). Namely, in the case of conversions (<33%), the selectivity data indicate that the HYD route (route #2) is favored, whereas the highest conversion, the DDN route (route #1) takes over with a selectivity of 1.4 at 99% convers With this, it appears evident that for the highest conversion, the DDN scenario is the preferred o      Apparently, three regions can be distinguished. In the first one (region I: tc < 0.27 s), the HYD route prevails with the presence of indoline and OEA to some extent. The second region (region II: 0.27 < tc < 0.37 s) is characterized by a sudden decrease in indole mole fraction along with the increase in OEA mole fraction as well as a significant increase in EB and ECH concentrations. The linear drop in indole content in regions I and II indicates that the indole extinction process is of zero order. In the second region, the real HDN reaction of OEA occurred (in this region, OEA is clearly the most abundant intermediate).
When the contact time is longer than 0.37 s (region III), indole is consumed almost quantitatively (i.e., the total conversion reaches 99%). The mole fraction of OEA decreases and takes over the role of reactants leading to HDN products such as EB and ECH (Figure 3). The simultaneous increase in both EB and ECH shares indicates that the reaction occurs according to both pathways (route #1 and #2); however, it is worth noting that route #1 is the favored one. These results are contradict the results reported by Ozkan et al. [54]. Therein, the authors have shown that the presence of H2S could invert the selectivity towards route #2. In our work, the presence of sulfur (dimethyl disulfide (DMDS) as a source of H2S), did not switch the reaction pathway. The alteration in selectivity between ECH and EB observed for tc > 0.8 s (Figure 3) is only due to the progressive hydrogenation of EB to ECH, which can occur without appreciable limitations because of the absence of N-containing compounds as shown elsewhere [57,[64][65][66][67][68][69][70]. To the best of our knowledge, such specific behavior (i.e., the privileged pathway of aromatics formation) is a unique feature of the molybdenum carbide catalyst.
Over Mo2C, aromatic products are favored because the nucleophilic substitution SN2 is the dominant mechanism governing the HDN reactions [65]. On the other hand, it is well known that the presence of indole and indoline strongly inhibits the conversion of ethylaniline [66] and, as a consequence, it inhibits the whole HDN process [67]. The formation of ethylbenzene as the dominant HDN product of indole was also reported by Piskorz et al. [68]. Based on the density functional theory (DFT) simulations, they found out that the selectivity of the reaction depends on the manner of the H: attack on the aromatic ring or the amino group of the OEA intermediate. The difference in the energy barriers between the DDN and HYD molecular paths calculated by the authors was equal to 33.5 kJ mol −1 (8 kcal mol −1 ) [68]. This justifies the higher amounts of the aromatic compounds (i.e., EB) compared to aliphatic ones (such as ECH). Apparently, three regions can be distinguished. In the first one (region I: t c < 0.27 s), the HYD route prevails with the presence of indoline and OEA to some extent. The second region (region II: 0.27 < t c < 0.37 s) is characterized by a sudden decrease in indole mole fraction along with the increase in OEA mole fraction as well as a significant increase in EB and ECH concentrations. The linear drop in indole content in regions I and II indicates that the indole extinction process is of zero order. In the second region, the real HDN reaction of OEA occurred (in this region, OEA is clearly the most abundant intermediate).
When the contact time is longer than 0.37 s (region III), indole is consumed almost quantitatively (i.e., the total conversion reaches 99%). The mole fraction of OEA decreases and takes over the role of reactants leading to HDN products such as EB and ECH (Figure 3). The simultaneous increase in both EB and ECH shares indicates that the reaction occurs according to both pathways (route #1 and #2); however, it is worth noting that route #1 is the favored one. These results are contradict the results reported by Ozkan et al. [54]. Therein, the authors have shown that the presence of H 2 S could invert the selectivity towards route #2. In our work, the presence of sulfur (dimethyl disulfide (DMDS) as a source of H 2 S), did not switch the reaction pathway. The alteration in selectivity between ECH and EB observed for t c > 0.8 s (Figure 3) is only due to the progressive hydrogenation of EB to ECH, which can occur without appreciable limitations because of the absence of N-containing compounds as shown elsewhere [57,[64][65][66][67][68][69][70]. To the best of our knowledge, such specific behavior (i.e., the privileged pathway of aromatics formation) is a unique feature of the molybdenum carbide catalyst.
Over Mo 2 C, aromatic products are favored because the nucleophilic substitution S N2 is the dominant mechanism governing the HDN reactions [65]. On the other hand, it is well known that the presence of indole and indoline strongly inhibits the conversion of ethylaniline [66] and, as a consequence, it inhibits the whole HDN process [67]. The formation of ethylbenzene as the dominant HDN product of indole was also reported by Piskorz et al. [68]. Based on the density functional theory (DFT) simulations, they found out that the selectivity of the reaction depends on the manner of the H: attack on the aromatic ring or the amino group of the OEA intermediate. The difference in the energy barriers between the DDN and HYD molecular paths calculated by the authors was equal to 33.5 kJ mol −1 (8 kcal mol −1 ) [68]. This justifies the higher amounts of the aromatic compounds (i.e., EB) compared to aliphatic ones (such as ECH).
At 99% indole conversion, the dominating products are EB and ECH. When the contact time is longer than 0.54 s, there are no more nitrogen-containing products. Interestingly, a change can be seen in the preferred route of the reaction-ECH is produced with increased efficiency, while the amount of formed EB decreases. While a significant number of nitrogen-compounds are still present in the feed, hydrogenation of the unsaturated products is inhibited.
In the case of the contact time t c > 0.53 s (no N-containing compounds), the effect of the total change in the selectivity was observed. This may be associated with the hydrogenation of EB to ECH. This is in accordance with the results reported formerly by Mordenti [69], who indicated that OEA (as the main N-containing intermediate product comprising nitrogen) is responsible for blocking the catalytic centers. Consequently, the hydrogenation of ethylbenzene to ethylcyclohexane is impeded. Similar research has been published previously by Sayag [70], who presented in his work evidence supporting the proposed inhibitory effect of o-propylaniline (OPA) in isopropyl benzene hydrogenation to propylcyclohexane.
The same researcher with co-workers [57] studied the HDN of indole over Mo 2 C supported on activated carbon black composites (Mo 2 C/CBC). They reported a maximum HDN conversion (100%) reached over all the studied catalysts notwithstanding the contact time. The authors found out that the increase in saturated products comes only from the hydrogenation of the unsaturated ones. Indeed, while N-compounds are still present in the feed in abundant amounts, hydrogenation of the unsaturated products is inhibited because of a competitive strong adsorption of these molecules onto the Mo 2 C surface [57,67].
A similar phenomenon was observed by Da Costa et al. [64], who studied the influence of 4,6-dimetylodibenzotiophene (4,6-DMDBT) on hydrogenation of biphenyl (BPh) over Mo 2 C supported on Al 2 O 3 . When 4,6-DMDBT was present in the reaction system, the route of hydrogenation of BPh to cyclohexylbenzene (CHB) was impeded because of strong adsorption of S-compounds on the surface of the catalyst. Furthermore, a similar effect was observed in the case of HDS of dibenzothiophene (DBT) [16]; the CHB was produced from tetrahydrodibenzothiophene because of the strong adsorption of S-compounds. Furthermore, Mamede et al. [71] have observed a similar behavior when carrying out the tests of hydrogenation of toluene over Mo 2 C in the presence of thiophene.
For the sake of comparison, Table 3 reports the chosen literature data on HDN of indole over carbides and sulfides. Sulfide catalysts feature a high selectivity towards the saturated product (ECH), in opposition to observation over molybdenum carbide. This effect is due to the enhancement of the C (sp3) -N bond cleavage in the presence of sulfur [72] instead of the C (sp2) -N bond cleavage. The selectivity of HDN of indole over Mo 2 C found in majority of works is shifted towards the aromatic product (i.e., EB). This is consistent with our results (see last line in Table 3). It can be concluded that, despite the presence of 50 ppm of S, the behavior of Mo 2 C is completely different than that of the sulfide catalysts. The undoubted advantage of Mo 2 C is the possibility of significant reduction in hydrogen consumption for the process. As such, bulk Mo 2 C may be a promising alternative for commercial HDN catalysts used on a technological scale.

Synthesis of Bulk Mo 2 C
Molybdenum carbide (Mo 2 C) was prepared using a modified procedure reported by Volpe and Boudart [46]. The temperature profile used for the synthesis is presented in Figure 4.

Synthesis of Bulk Mo2C
Molybdenum carbide (Mo2C) was prepared using a modified procedure reported by Volpe and Boudart [46]. The temperature profile used for the synthesis is presented in Figure 4.  Molybdenum trioxide was put on a sintered holder placed in a quartz reactor and was subjected to carburization with a mixture of 10 vol. % CH4/H2 at a flow rate of 13.2 L g −1 h −1 . The temperature was raised linearly from 300 to 973 K at a heating rate of 53 K h −1 and then the material was held at 973 K for the next 1 h. Eventually, the sample was cooled to room temperature (RT) under the constant flow of the carburizing mixture. In the case when exposure to air was needed, the sample was passivated using 1 vol. % O2/He (pulse method) in order to form a protective oxide layer on its surface. Molybdenum trioxide was put on a sintered holder placed in a quartz reactor and was subjected to carburization with a mixture of 10 vol.% CH 4 /H 2 at a flow rate of 13.2 L g −1 h −1 . The temperature was raised linearly from 300 to 973 K at a heating rate of 53 K h −1 and then the material was held at 973 K for the next 1 h. Eventually, the sample was cooled to room temperature (RT) under the constant flow of the carburizing mixture. In the case when exposure to air was needed, the sample was passivated using 1 vol.% O 2 /He (pulse method) in order to form a protective oxide layer on its surface.

Characterization Techniques
Structural characterization of the Mo 2 C catalyst was carried out by means of X-ray diffraction (XRD) using a SIEMENS D-500 diffractometer (Siemens, Berlin, Germany) working with the CuKα monochromatized radiation (λ = 1.54184 Å). The XRD patterns were analyzed using HighScore (Malvern Panalytical, United Kingdom) software. The degrees of crystallinity of the catalyst were estimated based on the intensity of the reflections recorded in the angular range of 2θ = 10-90 • .
A Quantachrome-Quantasorb Jr apparatus (Quantachrome Instruments, Boynton Beach, FL, USA) was employed for the S g measurement in dynamic conditions. The S g of the catalyst was calculated from the N 2 adsorption/desorption at different partial pressures of nitrogen. The standard S g measurement using the Brunauer-Emmett-Teller (BET) model is based on a three-point analysis based on the adsorption of mixtures containing 10, 20, and 30 vol.% of N 2 in He. Prior to the measurement, the material (100 mg) was outgassed under flowing nitrogen at 623 K for 2 h.
CO adsorption is a classic technique commonly used for estimating the number of surface metallic sites. It is also employed to titrate the active centers on the surfaces of the transition metal carbides or nitrides [47,[60][61][62]. The CO titration was performed in situ, i.e., in the synthesis reactor without exposing the fresh material to air or to the passivating mixture after synthesis. The pulses of given quantities of CO (17 µmol) were injected to the reactor containing Mo 2 C with regular intervals at RT under He flow (40 mL min −1 ) purified before by an oxygen trap (Oxisorb, Messer Griesheim GmbH, Bad Soden am Taunus, Germany). After each pulse, the number of CO molecules not chemisorbed was analyzed using a conventional detection system equipped with a TCD. The CO injections were continued until the probe molecule quantitatively saturated the catalyst's surface. The acquired data were processed and the number of micromoles of the chemisorbed CO per gram of the sample was calculated.

Kinetic Study of HDN of Indole
The HDN tests of indole were carried out using a downflow fixed-bed microreactor working in a high-pressure flow system. The catalytic tests were performed over 0.8 g of catalyst (diluted with SiC at a mass ratio SiC/Mo 2 C = 5:1) in the presence of sulfur (50 ppm S), evolved by the decomposition of 70 ppm dimethyl disulfide (DMDS) in the feed. The addition of sulfur was employed in order to simulate the HDN process of real petroleum intermediates that always contain certain amounts of S-compounds. The liquid feed containing 2.5 wt. % of indole dissolved in decaline was supplied to the flow manifold by a high-pressure piston pump (Gilson model 307, Artisan Technology Group, Mercury DriveIL, Champaign, IL, USA). Hydrogen flow was controlled by a mass flow controller (Brooks 5850TR, Brooks Instruments, PA, USA), while the total pressure was adjusted by a back-pressure regulator (Brooks 5866, Brooks 5850TR, Brooks Instruments, PA, USA). The temperature in the reactor oven was adjusted using a temperature controller (Hanyoung MX9, Hanyyoung, Vietnam).
The kinetic study of the HDN process was performed at a total pressure of 60 bar, and a temperature of 613 K using the H 2 /feed volume ratio of 600 and contact times (t c ) between 0.13 and 1.07 s. The contact time (t c ) was calculated as follows: The liquid products were collected with constant intervals of 1 h and analyzed via manual injections (0.8 µL) to a gas chromatograph (HP 4890, Agilent Technologies, Santa Clara, CA, USA) equipped with a capillary column (DB1: 30 m × 0.25 mm × 0.25 µm, Agilent Technologies, Santa Clara, CA, USA) and a FID detector. The temperature of the injector was kept at 523 K, and the split ratio was 1:50. High-purity argon was used as the carrier gas. The oven temperature program was the following: the temperature was increased from 323 K to 453 K at a rate of 5 K min −1 , and then held at 453 K for 1 min. Next, the temperature was increased from 453 K to 493 K at 10 K min −1 and held at this temperature for 2 min. The total time of the analysis was 33 min.
The carbide catalyst was stabilized in situ for 10 h using a mixture of 70 ppm of DMDS solution in decaline in hydrogen (60 mL min −1 ) as a carrier gas, at the processing temperature under the total pressure of 60 bar.
The total conversion of indole, degree of HDN, and selectivity were calculated according to Equations (2)-(4): Conv. of indole [%] = N N N N + N indole ·100%; ( Degree of HDN [%] = N H N N + N indole ·100%; (3) where: N N is the mole % of all products of decomposition of indole (including N-containing compounds except indole); N indole is the mole % of indole in the liquid product after reaction; N H is the mole % of non-nitrogen products from the decomposition of indole; EB + T + B are mole % of ethylbenzene, toluene, and benzene, respectively; ECH means mole % of ethylcyclohexane.

Conclusions
Molybdenum carbide was found to be active in HDN of indole in the presence of a low amount of sulfur (50 ppm). It was found that the preferred HDN reaction path depends on the indole conversion. Namely, in the case of low conversion, direct denitrogenation is the favored reaction route with EB being the main product. Contrarily, in the case of the highest conversion, the privileged pathway is the full hydrogenation towards ethylcyclohexane. The preferential direct denitrogenation route is opposite to the mechanism, which was hitherto observed over sulfide catalyst, which in general leads mainly to the hydrogenated products. Thus, the indole HDN mechanism over the bulk molybdenum carbide is different compared to that of sulfide catalysts. Over the studied Mo 2 C material, the hydrogenated products (such as ECH) become the main products only in the case where N-containing products are eradicated from the reactant mixture. Simultaneously, a minor dealkylation route towards aniline was observed. This brings certain significant advantages, i.e., the consumption of hydrogen in HDN processes may be significantly reduced. To the best of our knowledge, herein we report on this peculiar behavior of the Mo 2 C catalyst for the first time.