Effects of pore size, mesostructure and aluminum modification on FDU-12 supported NiMo catalysts for hydrodesulfurization

A series of NiMo/FDU-12 catalysts with tunable pore diameters and mesostructures have been controllably synthesized by adjusting the synthetic hydrothermal temperature and applied for the hydrodesulfurization of dibenzothiophene and its derivative. The state-of-the-art electron tomography revealed that the pore sizes of FDU-12 supports were enlarged with the increase in the hydrothermal temperature and the mesostructures were transformed from ordered cage-type pores to locally disordered channels. Meanwhile, the MoS2 morphology altered from small straight bar to semibending arc to spherical shape and finally to larger straight bar with the change of support structures. Among them, FDU-12 hydrothermally treated at 150 °C possessed appropriate pore diameter and connected pore structure and was favorable for the formation of highly active MoS2 with curved morphology; thus, its corresponding catalyst exhibited the best HDS activity. Furthermore, it was indicated that the isomerization pathway could be significantly improved for HDS of 4,6-dimethyldibenzothiophene after the addition of aluminum, which was expected to be applied to the removal of the macromolecular sulfur compounds. Our study sheds lights on the relationship between support effect, active sites morphology and HDS performance, and also provides a guidance for the development of highly active HDS catalysts.


Introduction
Nowadays, the demand of low-sulfur, high-quality transportation fuels is growing rapidly by reason of the necessity to solve pollution problems induced by SO x emissions and to meet the increasingly stringent environmental protection laws and regulations (Rajendran et al. 2020). The main industrial processes to realize sulfur removal include hydrodesulfurization (HDS), oxidative desulfurization Wu et al. 2020;Wei et al. 2020;Jiang et al. 2020) and adsorptive desulfurization (Dong et al. 2020;Liu et al. 2019;Wu et al. 2018a;Li et al. 2018a). During the HDS process, catalysts play an important role in determining the efficiency. Therefore, the development of novel HDS catalysts with improved desulfurization efficiency has attracted ever-increasing attention and become a hot issue needed to be solved urgently. Support is one of the most important parts for HDS catalyst, which plays a vital role in bearing active components, improving their dispersion as well as providing a space for the catalytic reactions. In addition, the support can affect the adsorption, diffusion and the accessibility of reactants to the active sites. Therefore, the properties of the support significantly influence the HDS catalytic activity. There are many tunable factors for support, such as composition, pore size and pore structure. Hence, the investigation of the influence Handling Editor: Jian Liu

Edited by Xiu-Qiu Peng
Electronic supplementary material The online version of this article (https ://doi.org/10.1007/s1218 2-020-00502 -5) contains supplementary material, which is available to authorized users. of support structural characteristics on HDS performance is of great importance and is the foundation for guiding the development of highly active HDS catalysts.
In recent years, mesoporous materials with highly ordered pore structure, large pore volume and specific surface area have drawn increasing interests in the field of catalysis and exhibited advantageous performance in the removal of aromatic sulfide molecules. Different kinds of mesoporous materials have been applied as the supports of HDS catalysts, such as MCM (Kaluza et al. 2019;Mendez et al. 2017;Wang et al. 2005), HMS (Vutolkina et al. 2019;Palcheva et al. 2016), SBA Thanh Tung and Qian 2018;Yuan et al. 2014) and KIT (Meng et al. 2020;Soni et al. 2010) series. The results showed that such supports favored the formation of type II active phases with good dispersion which exhibited higher HDS activity than γ-Al 2 O 3 supported catalysts. Meanwhile, both the pore diameter and structure of catalyst supports showed important effects on the HDS activity (Zhou et al. 2020;Meng et al. 2020;Zhang et al. 2019). Song et al. (2016) synthesized a highly ordered cubic mesoporous silica material (SBA-16@hexane) with a pore diameter of 15.1 nm, the largest one among all the reported pore sizes of SBA-16 materials. The prepared catalyst NiMo/Al-SBA16@hexane showed the highest HDS activity, which is attributed to the superior diffusivity of the novel SBA-16@hexane with ultra-large pore size. Soni et al. (2010) reported that, compared to SBA-15 supported catalyst with a 2D hexagonal structure, NiMo/KIT-6 with a 3D bicontinuous structure exhibited a much higher HDS performance since the 3D mesostructured KIT-6 could afford a better dispersion of active species and faster diffusion of reactants and products.
FDU-12 is a type of 3D mesoporous material with a cage-like face-centered cubic mesostructure which has been attracted increasing attention as supports for HDS in recent years (Hu et al. 2020;Meng et al. 2018;Wu et al. 2018b;Li et al. 2018bLi et al. , 2019Liu et al. 2016Liu et al. , 2017Wang et al. 2017;Zhang et al. 2017;Du et al. 2016). For example, FDU-12based micro-and mesoporous composites (Wu et al. 2018b;Li et al. 2018b;Du et al. 2016), Ti-or Zr-modified FDU-12 (Meng et al. 2018;Cao et al. 2014), etc. were successively developed, showing a high HDS efficiency due to the appropriate synergistic effect of pore properties and acid property. However, FDU-12 has a cage-like structure with large cage pore size and small entrance size, which restricts the entry of the active metals and limits the utilization of the large pore. Meanwhile, the effect of FDU-12 structure on the morphology of the active phases and HDS activity have not been reported so far to the best of our knowledge.
In this work, a series of FDU-12 with different pore diameters and structures have been controllably synthesized by adjusting the hydrothermal temperature from 100 to 210 °C. The properties of corresponding NiMo supported catalysts are comprehensively characterized, and the effects of the pore size, mesostructure and aluminum modification on the MoS 2 morphology and their catalytic performance for HDS are systematically explored. It is found that the morphology of the MoS 2 active site changes with the change of the support structure, further affecting the HDS performance.

Preparation of FDU-12 supports with different structures
The FDU-12 materials were synthesized according to the method reported previously with a minor modification (Yang et al. 2013). Firstly, F127 (2.0 g) and KCl (10.0 g) were completely dissolved in HCl (2 M, 120 mL) at 15 °C; then, TMB (2.4 g) was added into the mixture and kept stirring for 6 h. Afterward, TEOS (8.32 g) was added into the above mixture. After stirring for 24 h at 15 °C, the mixtures were directly transferred into an autoclave for hydrothermal treatment at a given temperature (100, 130, 150, 170, 190 and 210 °C) for 24 h. The solid product was collected by filtration, washed with deionized water, dried and then calcined at 550 °C for 5 h to remove the templates. For convenience, the final samples were marked as F-(T), in which T represents the hydrothermal temperature.

Preparation of the modified supports Al-FDU-12
The support F-150 was selected to modify with aluminum by the post-synthesis method (Klimova et al. 2008). The catalysts were prepared as follows: 1.0 g of F-150 support was dissolved in 100 mL of dry ethanol containing the required amounts of aluminum isopropoxide the desired Si/Al molar ratios of 100, 60, 40, 30 and 20. This mixture was kept for 24 h at 60 °C with magnetic stirring. Then, the solid material was collected by a vacuum filtered device, washed several times with ethanol and deionized water to eliminate chloride anions. Finally, the sample dried at room temperature was calcined at 550 °C for 5 h. The Al-modified materials are named as Al-F(x), in which x represents the theoretical Si/ Al molar ratio.

Catalyst preparation
NiMo catalysts supported on F-(T) and Al-F(x) were prepared by a standard incipient wetness co-impregnation technique reported elsewhere by using appropriate concentration of ammonium molybdate tetrahydrate and nickel nitrate hexahydrate as Mo and Ni sources, respectively. The nominal compositions of MoO 3 and NiO in all prepared catalysts were 12 wt% and 3 wt%, respectively, corresponding to the molar ratio of Ni:Mo = 3:7. After impregnation, the materials were dried at 80 °C for 12 h in an air atmosphere, followed by calcination at 500 °C for 4 h, denoted NiMo/F-(T) and NiMo/Al-F(x), respectively.

Catalyst characterization
Small-angle X-ray scattering (SAXS) profiles of the support and corresponding catalysts were obtained from Bruker NanoSTAR with a 2-D detector and X-ray beam pinhole collimated (40 kV, 30 mA). X-ray diffraction (XRD) experiments were collected on a Philips X'Pert diffractometer equipped with Cu Kα radiation at a wavelength of 1.5406 Å in the range of 2θ = 8°-70°. N 2 adsorption-desorption isotherms were determined with a Micromeritics ASAP 2002 instrument (USA). The specific surface areas were calculated using the Brunauer-Emmett-Teller (BET) method. The samples to be measured were firstly degassed in the preparation station at 180 °C with a vacuum of 10 −5 Torr for 15 h and then switched to the analysis station for adsorption-desorption experiment at − 196 °C. UV-Vis diffuse reflection spectroscopy (UV-Vis DRS) experiments were carried out on a Hitachi U-4100 UV-Vis spectrophotometer with the integration sphere diffuse reflectance attachment in the region of 200-800 nm. The standard support reflectance was used as the baseline for the corresponding catalyst measurement. Temperatureprogrammed reduction analyses with hydrogen (H 2 -TPR) for the oxidic catalysts were performed on a home-built apparatus. The sample (100 mg) needed to be pretreated in an Ar stream at 450 °C for 2 h and then cooled down to room temperature. Then, the Ar flow was switched to a 10% H 2 /Ar flow, and the sample was heated to 1000 °C at a rate of 10 °C/min. Finally, the H 2 consumption to reduce the metal oxides was detected by a thermal conductivity detector (TCD). Temperature-programmed desorption with ammonia (NH 3 -TPD) measurements using the same apparatus described for the H 2 -TPR experiments. Firstly, the sample (100 mg) was heated in a pure He flow from room temperature to 600 °C at a rate of 10 °C/min and then cooled to 100 °C. Secondly, NH 3 was adsorbed for 20 min, and the sample was subsequently purged by a flowing He stream to remove excessive and physically adsorbed NH 3 at 100 °C for 1 h. Thirdly, the sample was heated to 600 °C at a rate of 10 °C/min in a pure He flow and desorption pattern was recorded. Pyridine FT-IR adsorption experiments were measured on self-supported wafers in an in situ cell. The samples were dehydrated at 500 °C for 5 h under a vacuum of 1.33 × 10 −3 Pa followed by the adsorption of pure pyridine vapor at room temperature for 20 min. The system was evacuated at 200 and 350 °C, respectively, to record FTIR spectra. X-ray photoelectron spectroscopy (XPS) was carried out on a VG ESCA Lab 250 spectrometer using Al Kα radiation for the sulfided catalysts. C1s of adventitious carbon (285.0 eV) was used as reference to calibrate all binding energies (BE), and the obtained XPS spectra were fitted by using XPSPEAK software in order to quantify the contents of Mo and Ni species. Highresolution transmission electron microscopy (HRTEM) was obtained on a Philips Tecnai G2 F20 transmission electron microscope operated at an accelerating voltage of 200 kV, and the images were recorded on a TVIPS 1 k × 1 k CCD camera.

Catalytic activity
Prior to the catalytic activity testing, the catalysts were firstly presulfided for 3 h at 360 °C by using CS 2 (3 wt%) in cyclohexane. The HDS activity was evaluated in a continuously flowing tubular fixed-bed microreactor with the internal diameter of 10.0 mm and length of 500 mm. The assessment condition is 360 °C for NiMo/F-(T) and 340 °C for NiMo/Al-F(x) with the liquid hourly space velocity (LHSV) of 20 1/h, total pressure of 4.0 MPa and H 2 /hydrocarbon volumetric ratio of 300. 1 wt% DBT/4,6-DMDBT in heptane was used as a model compound. The catalyst (0.5 g) was diluted to a constant volume of 2.0 mL by using quartz particles (0.25 mm) before being put into the reactor. After the reaction was stabilized for 4 h, the products were collected and the sulfur content of original solution and products was analyzed by a WK-2C-type microcoulomb meter. The products were analyzed by an offline Finnigan Trace GC/MS instrument installed with an HP-5MS capillary column (30 m × 0.25 mm × 0.25 μm) and a pulsed flame photometric detector (O. I. Co., PFPD 5380). Assuming that the HDS reaction of DBT is a quasifirst-order reaction, the HDS activity of the catalyst can be obtained from the following formula (Klimova et al. 2008): where F is the number of moles of DBT entering the reactor per second, unit mol/s; m is the loading mass of the catalyst, unit g; τ is the total conversion of DBT hydrodesulfurization; and k HDS is the HDS reaction rate constant, unit mol/g·s.

Characterization of F-(T) supports and corresponding NiMo catalysts
FDU-(T) supports prepared with different hydrothermal temperatures were characterized by the SAXS and N 2 sorption, and the results are displayed in Fig. S1. Five well-resolved diffraction peaks assigned to the 111, 220, 311, 420, 511 and 531 reflections can be clearly observed in F-100 support ( Fig. S1a), indicating highly ordered large cage-type Fm3m mesoporous silica FDU-12 has been successfully synthesized (Yu et al. 2006). Although some of the diffraction peaks disappear with further increase in the hydrothermal temperature from 100 to 210 °C, the primary diffraction peaks indexed as 111, 220 and 311 are always present, indicating the preservation of the ordered mesostructure. It is found that the 111 diffraction peak gradually shifts to lower q value with the increase in the hydrothermal temperature, and accordingly, the unit cell parameters (a 0 ) increase from 35.6 to 39.6 nm (Table 1). Figure S1b shows the N 2 adsorption-desorption isotherms of FDU-(T) supports, and they all exhibit type IV isotherm. However, the shape of hysteresis loop changes from H2 to H1 (F-210) as increasing the hydrothermal temperature, indicating that the structure of FDU-(T) supports may transforms from 3D to 2D or to analogous 2D. Figures S1c and 1d give the pore size and entrance size distribution curves of FDU-(T) supports calculated from the adsorption and desorption branch isotherms, respectively. The entrance sizes of F-100 and F-130 supports are less than 3 nm, implying that the pores are rarely connected or even partially isolated. With the increase in the hydrothermal temperature from 150 °C to 190 °C, the pore sizes increase from 17.4 to 23.1 nm and the entrance sizes increase from 8.8 to 13.2 nm, indicating the pore connection is more full. The pore size (26.5 nm) and the entrance size (18.1 nm) of F-210 support are gradually approaching with further increase in the hydrothermal temperature and the entrance of the cage-type hole is mutually linked completely, suggesting the support structure has been transformed to homologous 2D straight channel. Figure 1 shows the electron tomography slice images of the supports under different hydrothermal temperatures. The results indicate that the supports (F-100 and F-130) obtained   -190 (d) at lower hydrothermal temperature have relatively independent cage-type holes with small connection in between, while the shape of cavity is relatively regular. As the hydrothermal temperature increases (F-150), the cavity size and window aperture increase significantly, and the cage-type holes are connected by a larger window. With the hydrothermal temperature further increases to 190 °C, the cavity and entrance sizes of the support are gradually approaching, while the cages are interpenetrated, which accords well with the N 2 sorption results, indicating that suitable hydrothermal temperature can effectively improve the cavity and entrance size. The SAXS patterns, N 2 adsorption-desorption isotherms and pore size distribution of NiMo/F-(T) catalysts are presented in Fig. 2. The SAXS patterns (Fig. 2a) of NiMo/F-(T) catalysts all exhibit two or three well-resolved peaks which are the characteristic of FDU-12 materials, suggesting that ordered large cage-type structure is maintained in all catalysts after deposition of Ni and Mo species. The N 2 adsorption-desorption isotherms of NiMo/F-(T) catalysts ( Fig. 2b) all correspond to a type IV, and the shape of hysteresis loop is almost unchanged after impregnation of NiMo species. In addition, the position of the capillary condensation is shifted to higher P/P 0 as the rise of hydrothermal temperature, implying the enlargement of the pore size, which is confirmed by the pore size and entrance size distribution curves calculated from adsorption (Fig. 2c) and desorption branches (Fig. 2d) of NiMo/F-(T) catalysts, respectively.
The textural and structural characteristics (specific surface area S BET , pore volume V P , pore size, entrance size and unit cell parameter a 0 ) of FDU-(T) supports and corresponding catalysts are summarized in Table 1. It can be observed that the specific surface areas of FDU-(T) supports decrease, and the pore volumes, pore sizes and entrance sizes increase gradually with an increase in the hydrothermal treatment temperature. The specific surface area of F-210 support (224.2 m 2 /g) is less than half of F-100 support (541.1 m 2 /g); nevertheless, the pore size and entrance size of the former (26.5 nm and 18.1 nm) are over twice and six times as large as the latter (10.8 nm and < 3 nm), respectively. F127 is one of the most studied amphiphilic surfactants that are commonly used as a template for synthesizing mesoporous materials, and it is believed that F127 can drive the formation of micelles in which the hydrophobic PPO blocks can act as the core and the hydrophilic PEO blocks as the corona (Mayanovic et al. 2014). It has been known that the hydrophilicity of the PEO blocks decreases with increased hydrothermal temperature and PEO blocks can be drawn back into the micelles which increases the hydrophobic segment volume of block copolymer micelles, consequently resulting in the enlargement of the pore size and entrance size and the concomitant decline of the surface area. Compared to FDU-(T) supports, the structural parameters of NiMo/F-(T) catalysts are all slightly decreased, which should be related to the increase in the samples' density after incorporation of Ni and Mo precursors or due to some blockage of pores (Wang et al. 2005). Figure 3a shows the powder XRD patterns of NiMo/F-(T) catalysts, and they all present a very broad peak at around 23° which is typical for amorphous silica (Cao et al. 2014). It can be seen that no signal of any crystalline phase is detected in NiMo/F-100 and NiMo/F-130 catalysts, suggesting that the metallic particles are well dispersed on the surface of the supports with high specific surface areas (Wang et al. 2002). However, weak diffraction peaks are observed for NiMo/F-150 catalyst at around 26.7°, which are corresponding to the reflections of MoO 3 in orthorhombic phase (JCPDS No. 09-0209). Furthermore, the intensity of MoO 3 diffraction gradually strengthens with the decline of the specific surface areas, indicating the dispersion of active particles decreases for NiMo supported catalysts. The active components with the same loading amount are easier to aggregate in the supports with lower specific surface areas and form large oxide particles after the high-temperature calcination.
To discern the chemical state and obtain more information about the dispersion of Mo species in NiMo/F-(T) catalysts, UV-Vis DRS spectra were recorded after subtracting the spectra of the corresponding supports as displayed in Fig. 3b. The adsorption bands represent the ligand-to-metal charge transfer O 2− → Mo 6+ , and the exact positions of these bands are affected by the aggregation and coordination state of the Mo 6+ species. The isolated molybdate species in tetrahedral coordination, Mo(T d ), show a characteristic absorption band at ~ 230 nm, whereas polymolybdate octahedral Mo(O h ) species absorb at 280-320 nm depending on the aggregation degree of these species (Weber 1995). As shown in the spectra (Fig. 3b), a mixture of Mo species in tetrahedral and octahedral coordination is present in all NiMo/F-(T) catalysts, but the Mo absorption edges produce a slight red shift with increased hydrothermal temperature, indicating the occurrence of larger MoO 3 clusters, which is in agreement with the XRD result (Fig. 3a).
To further compare the degree of sulfidation and the dispersion of metal species on the surface of FDU-(T) supports, the freshly sulfided catalysts were characterized by XPS. The Mo 3d XPS spectra of NiMo/F-(T) catalysts containing the peaks of S 2s and three Mo 3d doublets were deconvoluted as shown in Intensity, a.u. Moreover, S 2s level which the binding energy is at about 225.9 eV in the spectra must be subtracted from the total spectra of Mo 3d in the concentration calculation of the different Mo species (Escobar et al. 2009). The sulfidation degree of Mo species, Mo sulfidation , is expressed as the ratio of Mo 4+ /M TOT (M TOT = Mo 4+ + Mo 5+ + Mo 6+ ) (Gao et al. 2018;Han et al. 2012), and the corresponding results are given in Table 2. It is shown that the highest sulfidation degree is obtained from the NiMo/F-100 catalyst (51.3%). When the T is 150, the Mo sulfidation of the corresponding catalyst generates only a slight drop to 49.0%. However, further increase in the hydrothermal temperature results in a rapid decline of the sulfidation degree. It is decreased to 38.4% for the NiMo/F-210 catalyst probably owing to the formation of bulk MoO 3 crystallites. The surface Ni/(Ni + Mo) and S/Mo atomic ratios were also determined by XPS and the results are given in Table 2. The NiMo catalysts correspond to theoretical Ni/(Ni + Mo) atomic ratio of 0.3 and S/Mo atomic ratio of 2. In our experiment, the experimental Ni/(Ni + Mo) and S/Mo atomic ratios at lower hydrothermal temperature (100 and 130 °C) are almost around the theoretical value, indicating that Ni is evenly dispersed on the surface of catalysts. However, experimental Ni/(Ni + Mo) and S/Mo atomic ratios of the catalysts NiMo/F-190 and NiMo/F-210 gravely deviate from the theoretical value due to the poor dispersion of Mo species.
The Ni 2p XPS spectra of the sulfided NiMo/FDU-12 catalysts are presented in Fig. S3, and the Ni 2p envelope is decomposed by considering three major contributions corresponding to NiS x sulfided phase (which could arise from Ni 2 S 3 , Ni 9 S 8 or NiS, with binding energy between 853.5 and 853.8 eV), NiO x (binding energy of 856.9 eV) and the so-called NiMoS phase (with binding energy between 854.0 and 855.0 eV) (Ninh et al. 2011;Escobar et al. 2009). Two parameters of promotion rate (PR) and promoter ratio (see the corresponding formulas in the footnote of Table 2) are closely related to the formation of NiMoS active phases over the sulfided NiMo catalysts. The former represents the proportion of total deposited nickel that is actually engaged in the NiMoS phase formation; meanwhile, the latter is related to the amount of that highly active species (Cao et al. 2014;Merida-Robles et al. 1999). It is clear that the PR and (Ni/ Mo) slabs values of NiMo/FDU-12 catalysts also decrease with the increase in the hydrothermal temperature, which consists with the changing trend of the Mo sulfidation values.
HRTEM characterization of the sulfided NiMo/FDU-(T) catalysts was performed to obtain more information about the morphology and dispersion of MoS 2 crystallites induced by the support structures. Representative micrographs selected for each one of the catalysts are shown in Fig. 4. The typical lattice fringes of MoS 2 crystallites with a spacing of 0.61 nm (corresponds to the 002 interplanar distances) and some features of the FDU-12 mesopore structure are observed in micrographs of the sulfided catalysts (Cao et al. 2014). The straight MoS 2 slabs are dispersed well over the sulfided NiMo/F-100 and NiMo/F-130 catalysts (Fig. 4a, b). With an increase in hydrothermal temperature, the entrance size of the support increases significantly and massive active phases migrate into the cage-like holes and form curved MoS 2 crystallites (Fig. 4c) due to the confinement effect of the spherical pores. The entrance size and pore size of the supports further increase with the continue raise of the hydrothermal temperature (T = 170 and 190) and more active phases enter the cell. The curved MoS 2 crystallites keep growing along with the eyeball until forming a closed-shell structure with a much more stacked arrangement and poor dispersion of MoS 2 on NiMo/F-190 catalyst (Fig. 4e). The closed-shell MoS 2 structure consists of slabs running all along the periphery of the unit cells, while a quite small proportion of flat slabs is detected inside the periphery region. With further increase in the hydrothermal temperature, the entrance size of the support enlarges obviously and approaches the cell size. The support structure transforms to homologous 2D straight channel and straight MoS 2 slabs are formed on the NiMo/F-210 catalyst again (Fig. 4f).
From the above characterizations, the mesoporous material F-150 keeps ordered 3D cage-type structure with high specific surface area and large entrance size, and the corresponding NiMo catalyst shows relatively high Mo sulfidation , PR (Ni/Mo) slabs and better MoS 2 morphology. Therefore, we chose F-150 as support to carry out Al modification with different molar ratios of Si:Al by chemical grafting method and investigate the effect of alumina content on the HDS activity and selectivity of the corresponding NiMo catalysts.

Characterization of Al-F(x) supports and corresponding NiMo catalysts
F-150 support was modified with different amounts of Al by using aluminum isopropoxide as alumina source through the post-synthesis method to obtain Al-F(x) supports. Figure S4 shows the N 2 adsorption-desorption isotherms of Al-F(x) supports, and they all display a typical type IV isotherm with sharp inflection in the relative pressure range of 0.7-0.9. H2-type hysteresis loop confirms the ordered nature of the mesoporous F-150. No change is detected in the form of the isotherms or in the shape of the hysteresis loop, indicating the preservation of the support pore structure after Mo and Ni deposition (Fig. 5a). Analyzing from the adsorption and desorption branches of the isotherms, the maximum pore size and entrance size are located at 17.0 nm and 8.3 nm, respectively (Fig. 5b, c). The parameters about surface areas, pore volumes, pore sizes and entrance sizes of NiMo/F-150 and NiMo/Al-F(x) catalysts are given in Table 3. It can be observed that the alumination of F-150 by reacting with Al 3+ ions only leads to a slight drop in surface area and pore volume which might be attributed to the increase in the sample density after incorporation of aluminum (Klimova et al. 2008). This drop becomes more significant with the increase in the Al loading (x = 20) due to pore blockage caused by a low dispersion of the metal phases (Mo and Ni). However, the pore sizes and entrance sizes do not suffer from a significant change after Al grafting. In order to study the effect of incorporation of Al atoms on the dispersion of Mo and Ni species, powder XRD was performed and the results are displayed in Fig. 6a. As shown in the spectra, the intensity of the diffraction corresponding to reflections of MoO 3 with orthorhombic phase (JCPDS No. 09-0209) decreases with the increase in the Al loading in the F-150 support, indicating the dispersion of the active components is improved. However, there is a peak located at 2θ = 23.4° whose intensity increases evidently with excess Al loading, implying the formation of aluminum molybdate Al 2 (MoO 4 ) 3 (JCPDS No. 23-0764) in the NiMo/Al-F(20) catalyst. Klimova (Gutierrez-Alejandre et al. 2015) reported that part of excess extra-framework Al 3+ species which were in the form of octahedral coordination was easy to form Anderson-type heteropolymolybdate and then transformed to Al 2 (MoO 4 ) 3 after calcinations (Klimova et al. 2008). Therefore, Al incorporation in the F-150 support, on the one hand, enhances the dispersion of Mo oxide species. On the other hand, it induces the formation of Al 2 (MoO 4 ) 3 . The DRS spectra of NiMo catalysts (Fig. 6b) show that Mo species in tetrahedral and octahedral coordination all present and also confirms the improved dispersion of Mo species with Al addition, and the NiMo/Al-F(x) catalysts produce a blue shift of the Mo absorption edge with decreasing x.
H 2 -TPR measurements of NiMo/F-150 and NiMo/Al-F(x) catalysts were taken to explore the reducibility of the Ni and Mo species present in the oxidic precursors of the different catalysts (Fig. 6c). The reduction profiles of all the NiMo catalysts show hydrogen consumption in a broad temperature interval (between 300 and 900 °C) with two areas along with many reduction peaks. The low-temperature peaks (300-600 °C) can be attributed to the first-step reduction of polymeric octahedral Mo species from Mo 6+ to Mo 4+ . The peaks at high temperature (600-900 °C) can be ascribed to the second-step reduction from Mo 4+ to Mo 0 of the polymeric octahedral Mo species with different aggregation degrees and the reduction of tetrahedral Mo species from Mo 6+ to Mo 0 (Qu et al. 2003). It is observed that both the low-and high-temperature peaks shift to higher temperature with the decrease in the Si/Al molar ratio, indicating the enhancement of the metal-support interaction after the incorporation of Al. The effect of Al modification on the acidity of catalysts was investigated by NH 3 -TPD (Fig. 6d) and pyridine FT-IR (Fig. S5). The strength of the acid sites can be determined by the desorption temperature of the adsorbed NH 3 . Based on the desorption temperature, the acid sites are classified as weak (100-250 °C), medium (250-400 °C) and strong (T > 400 °C) (Li et al. 2012). It can be seen that all the catalysts display a broad peak in the temperature range from 75 to 400 °C, revealing the presence of both weak and medium acid sites. In addition, the intensity of the desorption peaks shows that the number of acid sites in the NiMo/Al-F(100) catalyst is much higher than that in the NiMo/F-150 catalyst. Moreover, the increase in the Al content leads to a higher intensity of the desorption peaks, which means that Al modification can improve the acidity.   The pyridine adsorption IR spectroscopy measurements were taken after evacuation at 200 and 350 °C in the wave number range 1600-400 cm −1 , and the results are shown in Fig. S5. The peak at 1450 cm −1 is assigned to pyridine adsorbed onto L acid sites and the peak at 1540 cm −1 is ascribed to pyridine adsorbed onto B acid sites, while the peak at 1490 cm −1 is attributed to pyridine adsorbed onto both B and L acid sites (Kataoka and Dumesic 1988). An attempt has also been made to quantitatively estimate the number of B and L acid sites for the five catalysts using pyridine adsorption followed by degassing at 200 and 350 °C, corresponding to the total and medium/strong acid sites, respectively. The acid strength distributions and acid quantities of NiMo/F-150 and NiMo/Al-F(x) catalysts are listed in Table 4. After degassing at 200 °C, the amount of B acid sites of the NiMo/Al-F(100) catalyst is about 3.6 times as high as that of the NiMo/F-150 catalyst, while this difference is slightly smaller (about two times) for the amounts of medium and strong acid sites after degassing at 350 °C, and the amounts of B acid sites of the NiMo/Al-F(x) catalysts all increase with increasing Al content. The L acid sites of the NiMo catalysts present a similar changing trend with B acid sites. In addition, the B/L ratio of the NiMo catalysts gradually increases and then decreases with the increase in the Al loading amount and peaks at the Si/Al molar ratio of 30. The Mo 3d and Ni 2p XPS spectra of the sulfided NiMo catalysts accompanied by their deconvolutions are displayed in Figs. S6 and S7. The corresponding parameters of the different contributions Mo 3d and Ni 2p are concluded in Table 5. It is shown that the values of Mo sulfidation , S/Mo and Ni/(Ni + Mo) of Al-incorporated catalysts are higher than that of the NiMo/F-150 catalyst and they all increase with the decrease in the Si/Al molar ratio from 100 to 30, indicating the incorporation of an appropriate amount of Al can improve the dispersion of Mo and Ni species. However, the above three values of the NiMo/Al-F(20) catalyst decrease in comparison with the NiMo/Al-F(30) catalyst, which could be attributed to the appearance of aluminum molybdate and the decrease in the surface area.
The HRTEM images of the sulfided catalysts NiMo/ Al-F(100) and NiMo/Al-F(40) are shown in Fig. S8, and it can be observed that the morphology of MoS 2 crystallites changes depending on the support. As described previously, in the pure silica F-150 supported catalyst, MoS 2 crystallites with length between 10 and 15 nm and stacking layers from five to twelve are found (Fig. 4c). Al incorporation in F-150 as support results in a better dispersion of MoS 2 phase.

HDS catalytic activity
Catalytic behaviors of the NiMo/F-(T) catalysts were evaluated in HDS of DBT at 360 °C, and the results are shown in Fig. 7. As the hydrothermal temperature increases, the corresponding catalyst desulfurization rate increases first and then decreases, and NiMo/F-150 shows the maximum HDS rate of 94.5%. Although the supports (F-100 and F-130) prepared at lower hydrothermal temperature have larger specific surface area and the corresponding catalysts have the highest proportion of sulfidation degree and NiMoS active phase, the obtained support entrance sizes are too small and restrict   1 3 the diffusion of materials. In addition, the small entrance sizes of the catalysts obtained at a lower hydrothermal temperature not only limit the diffusion of reactants and product molecules during the reaction (Table S1), but also prevent the sufficient contact between the active phases and the reactants, thereby reducing the HDS performance. As the hydrothermal temperature increases, the cavity and entrance sizes of the support are increased significantly, and the desulfurization rate of the catalyst NiMo/F-150 reached the maximum. This is because the larger pore sizes can facilitate the rapid transport of the materials and promote the reactants to contact sufficiently with the active phases. Further increasing the hydrothermal temperature, the specific surface area of prepared materials decreases significantly and the diffusion resistance is small, resulting in a deterioration of the dispersion of the active component so that some of the active phases do not function during the reaction. In addition, the proportion of the sulfidation degree of the catalyst and the active phases of NiMoS are also greatly reduced, which seriously affect the HDS reaction. The values of the HDS ratio (performed at 340 °C) and pseudo-first-order rate constant k HDS for the catalysts with Al incorporation are summarized in Table 6. It is shown that all the catalysts exhibit higher catalytic activity than that without Al. The HDS ratio increases when Si/Al molar ratio decreases from 100 to 30. NiMo/Al-F(30) catalyst possesses the highest HDS ratio (97.8%) and k HDS (1.27) values, which are about 1.25 times and 3.43 times as high as NiMo/F-150 catalyst (HDS ratio of 78.1% and k HDS of 0.37), respectively. Further increase in Al loading (Si/Al molar ratio of 20) leads to a decrease in the catalytic activity, and the obtained HDS ratio and k HDS are 88.7% and 0.81 for NiMo/Al-F(20), respectively.
HDS of DBT can be realized through two parallel reaction routes (as shown in Fig. S9). One is the direct desulfurization route (DDS) to yield biphenyl (BP) through the hydrogenolysis of C-S bonds and then slowly generate cyclohexylbenzene (CHB) by hydrogenation, which can be negligible. The other one is the hydrogenation pathway (HYD) to produce tetrahydrodibenzothiophene (THDBT) and hexahydrodibenzothiophene (HHDBT) intermediates first and then CHB through the hydrogenation of the aromatic rings of DBT molecule. To compare the contributions of the DDS and HYD pathways to the overall HDS activity of the five prepared catalysts, the HYD/DDS selectivity is calculated by dividing the productivity of CHB, THDBT and HHDBT by the yield of BP at the similar DBT conversion (about 50%), and the results are presented in Table 6 ( Pena et al. 2014;Klimova et al. 2009). It is clearly seen that the HDS of DBT proceeds predominantly through the DDS route and the HYD/DDS ratio is positively correlated with the specific activities of the catalysts, indicating that the total HDS activities of DBT of the sulfided NiMo/F-150 and NiMo/Al-F(x) catalysts depend mostly on the conversion rate to BP (Nava et al. 2007). The HYD/DDS ratios follow the order: NiMo/F-150 < NiMo/Al-F(100) < NiMo/ Al-F(20) < NiMo/Al-F(60) < NiMo/Al-F(40) < NiMo/ Al-F(30). The HYD/DDS ratios of the Al-incorporated catalysts are higher than that of the NiMo/F-150 catalyst, indicating that the desulfurization of DBT over the Al-modified supported catalysts promotes the HYD route more than the  pure silica supported catalyst. The effect of Al incorporation can be attributed to the structure and specific electronic properties of the active phases formed on the Al-modified support. The acceleration of the HYD pathway must be associated with a good dispersion of the active phase from the above analysis results. Moreover, the amount of B acid sites has also some relationship to HYD route. The NiMo/ Al-F(30) catalyst presents the highest HYD/DDS ratio, while the ratio decreases in the NiMo/Al-F(20) catalyst which may be explained by the fact that the excessive Al loading causes a decrease in surface area and the appearance of Al 2 (MoO 4 ) 3 , leading to a decrease in the amount of active phases. NiMo/Al-F(30), the best desulfurization catalyst, was evaluated with the model compound 4,6-DMDBT at the reaction temperature of 340 °C, the weight hourly space velocity (WHSV) of 15 l/h, the H 2 /oil volume ratio of 300 and the pressure of 4.0 MPa. We speculated the reaction path (Fig. 8) of 4,6-DMDBT according to the product distribution analyzed using GC-MS. It presents three reaction paths of DDS, HYD and isomerization (ISO) in which the first two paths were reported in most of the studies (Zhang et al. 2017). In the DDS pathway, 3,3′-dimethylbiphenyl (3,3′-DMBP) is obtained by hydrogenolysis of 4,6-DMDBT molecules directly. In the HYD pathway, a benzene ring is first hydrogenated to form 4,6-tetrahydrodimethyldibenzothiophene (4,6-THDMDBT) or 4,6-hexahydrodimethyldibenzothiophene (4,6-HHDMDBT), and after desulfurization, 3,3′-dimethylcyclohexenebenzene (3,3′-DMCHEB) and 3,3′-dimethylbicyclohexane (3,3′-DMBCH) are formed. Partially isomerized products are produced ascribing to the increase in B acid content after Al modification. It promotes methyl transfer and produces 3,7-dimethyldibenzothiophene (3,7-DMDBT) and 2,7-dimethyldibenzo thiophene (2,8-DMDBT). Afterward, 4,4′-dimethylbiphenyl (4,4′-DMBP) and 4,5′-dimethylbiphenyl (4,5′-DMBP) are formed after desulfurization. Further hydrogenation generates a small amount of 4,4′-dimethylbicyclohexane (4,4′-DMBCH) and 4,5′-dimethylbicyclohexane (4,5′-DMBCH). The specific contents are shown in Fig. 8. It is known that the 4,6-DMDBT molecule is difficult to be desulfurized because of the strong steric hindrance which makes the active sites difficult to directly contact with sulfur atoms. But herein, the methyl transfer occurring in the isomerization can greatly reduce this steric effect, and additionally, the sizes of the entrance and cavity of the NiMo/Al-F(30) catalyst are large, which reduces the diffusion resistance of the 4,6-DMDBT molecule in the reaction and is expected to be beneficial for the hydrodesulfurization.

Conclusion
A series of NiMo/FDU-12 catalysts have been comprehensively characterized to explore the effects of the pore size, mesostructure and aluminum modification on the MoS 2 morphology and their catalytic performance for HDS. The support structures transform from 3D to homologous 2D straight channel with the increase in the hydrothermal temperature, while the MoS 2 morphology alters from straight bar to semibending arc then to closed shell and finally to straight bar on the corresponding catalysts. It is found that MoS 2 with bended morphology has more active sites with high activity. Thus, F-150 support possessing appropriate pore diameter, 3D structure and curved MoS 2 exhibits the best HDS performance. The incorporation of aluminum facilitates the dispersion of Mo and Ni species. With an appropriate Al incorporation, Mo and Ni species are easier to be sulfided to form MoS 2 and NiMoS phases. HDS tests show that the catalytic activities gradually increase and then decrease with the increase in the Al content, and the NiMo/ Al-F(30) catalyst exhibits the highest HDS efficiency due to its relatively high surface area, good dispersion of the active phases and high acid content. Our study demonstrates that the support structure plays an important role on active sites morphology and HDS performance, which might help to develop highly active HDS catalysts for macromolecular sulfur compounds.