Abstract

More than half of the total world oil reserves are heavy oil, extra heavy oil, and bitumen; however their catalytic conversion to more valuable products is challenging. The use of submicronic particles or nanoparticles of catalysts suspended in the feedstock may be a viable alternative to the conversion of heavy oils at refinery level or downhole (in situ upgrading). In the present work, molybdenum sulfide (MoS2) particles with varying diameters (10000–10 nm) were prepared using polyvinylpyrrolidone as capping agent. The prepared particles were characterized by DLS, TEM, XRD, and XPS and tested in the hydrodesulfurization (HDS) of a vacuum gas oil (VGO). A correlation between particle size and activity is presented. It was found that particles with diameters around 13 nm show double the HDS activity compared with the material with micrometric particle sizes (diameter ≈ 10,000 nm).

1. Introduction

The depletion of light oils supplies along with the expected increase in energy demand has shifted the attention to the exploitation and conversion of heavy oils and heavy petroleum fractions. More than half of the total world oil reserves are heavy oil, extra heavy oil, and bitumen, with vast deposits in Canada (Alberta province) and Venezuela (Orinoco region). Even though catalytic hydroprocessing is a well-established process in the refining industry the conversion of such heavy oil fractions is challenging, from the catalytic point of view, due to metal sulfides and coke deposition in the pore mouth of the support in the commonly used Ni(Co)Mo(W)/alumina catalysts, which leads to the reduction of the catalysts effective life. One alternative to avoid the problem posed by the use of supported porous catalysts is the application of ultradispersed catalysts, which are very small particles suspended in the feed to be processed [18]. Also, in situ (or downhole) processing has been proposed as a way of producing synthetic oils from heavy oils and bitumen in the reservoir, thus reducing energy demands for field hydroprocessing and decreasing emissions [911]. For this application the particle size of the catalyst should be as small as possible (preferably in the nanorange, i.e., <100 nm) in order to increase their efficiency and to be able to navigate in the porous media without causing any plugging or impairment of the reservoir [12].

MoS2 is the base of hydroprocessing catalysts; however there are no reports on the effect of the particle size on the HDS activity for real feeds using unsupported (bulk) MoS2. Kuriki et al. [13] obtained particles of different diameters (2600–40 nm) by milling commercial MoS2 and tested them in the hydrogenation of 1-methylnaphthalene. From the data presented it can be inferred that smaller particles possess higher specific activity but a correlation between conversion and particle size is not presented.

Even though smaller particles are expected to have higher activity, this may not always be the case [14]. Also, the kind of feed and the morphology of the particles could have an important effect on their activity [15].

Preparation of MoS2 particles has been performed by different methodologies including the sulfidation of oxides [16, 17], the decomposition of precursors [18], ultrasonic irradiation [1921], the use of emulsions [2227], or using Layered Double Hydroxide (LDH) as hard template [28]. Most of these methodologies do not produce nanoparticles or allow controlling the particle diameter. Additionally, in some cases, the yields obtained from such preparations are very low and further testing the activity of the particles prepared is not possible.

Synthesis of nanoparticles assisted by capping agents has been reported in the literature for the preparation of defined shapes and sizes of nanostructures, such as nanorods and nanowires [2931]. The preparation of particles coated with a thin protective polymer layer is a good approach since they are more stable against aggregation (due to the large decrease in surface energy). This is important since it is known that when approaching the nanometric range the high surface energy of inorganic nanoparticles makes them very unstable, prone to undergo chemical reactions with the environment and also to self-aggregate [29].

In the present work, MoS2 nanoparticles with different diameters were prepared in presence of polyvinylpyrrolidone (PVP) as capping agent using a simple methodology suitable for preparing homogeneous nanoparticles with different diameters and reasonable quantities.

The particles were characterized by Dynamic Light Scattering (DLS), Transmission Electron Microscopy (TEM), X-ray Diffraction (XRD), and X-ray Photoelectron Spectroscopy (XPS). The hydrodesulfurization (HDS) of a vacuum gas oil (VGO) in a batch reactor using 3.44 MPa of hydrogen and 653 K was used for catalytic activity determinations.

2. Materials and Methods

2.1. Materials

Crystalline ammonium tetrathiomolybdate (99.97%) and polyvinylpyrrolidone (99%) were obtained from Aldrich and used without any further treatment. Sulfuric acid (98%) and methanol of analytical grade were also used in the preparations of the particles. A vacuum gas oil (VGO), whose characteristics are shown in Table 1, was used as feedstock for the hydrodesulfurization (HDS) reactions.

2.2. Particles Preparation

The preparation of the MoS2 was made based on a methodology adapted from the one reported by Amma and collaborators [31] for the preparation of cadmium sulfides nanoparticles using different capping organic agents. A simple preparation methodology was developed by precipitating MoS3 (which is converted to MoS2 by thermal treatment), acidifying a solution of ammonium tetrathiomolybdate (ATTM) and the capping agent according to the following chemical reaction:Thus, to 30 cm3 of a 0.01 M solution of ATTM in methanol : water (1 : 1; v/v), 15 cm3 of PVP in methanol solution (0.5, 0.1, 0.02, or 0.01 mM) was added and mixed well. Then 0.5 cm3 of a 0.5 M sulfuric acid solution was added drop by drop with intense agitation, which immediately turned the color of the solution from light orange to dark brown. The solutions were then centrifuged at 10000 rpm for periods of 20 minutes and washed three times with warm water and twice with methanol to eliminate (as much as possible) any remaining capping agent from the surface. The separation of the solids prepared using PVP concentration of 0.5, 0.1, and 0.02 mM was very difficult and the amount of solid obtained was low representing a yield of around 50%. For the other two preparations (0.01 and 0 mM of PVP) the separation was easy and the yield was about 75%. Finally, the obtained solids were air and vacuum dried (at 353 K) overnight to produce the as-prepared particles.

2.3. Particle Characterization

The size of the particles was determined by DLS using a Zetasizer Nano S (Malvern Instruments) equipped with a 4.0 mW He-Ne laser (633 nm) operated at an angle of 173° and a temperature of 298 K to obtain the average particle diameter (Z-Ave) of the suspended particles. The particles were redispersed in acetone and sonicated in a laboratory bath for 15 minutes prior to the size determination.

TEM images were obtained on a HITACHI 7650 microscope (point-to-point resolution using an accelerating voltage of 80 kV). Solid samples (a tip of a small spatula) were added to a vial containing 3 cm3 of a 1: 1 ethanol/water mixture and sonicated in a bath until the solution appeared homogeneous. Then one drop of the suspension was placed on the carbon side of a copper grid and let to dry in a fume hood. Once dried the samples were taken to the microscope for analysis.

XRD patterns were obtained with a Rigaku X-ray diffractometer (Model Ultima III) over the 2 theta range of 5 to 90 degrees, scanning speed of 0.2 degrees/second, and a step size of 0.02 degrees (2 theta) with Cu Kα radiation and 40 KV. The XRD patterns were referenced to the powder diffraction files (International Centre for Powder Diffraction PDF 2005) using the pattern processing analysis software JADE from Materials Data Inc.

A PHI VersaProbe 5000-XPS was used to record XPS spectra. The spectra were taken on samples pressed on double sided tape, using monochromatic Al source, 1486.6 eV, with 49.3 W of power and beam diameter of 200.0 μm. Double neutralization, that is, a low energy electron beam and low energy Ar+ beam, was used during spectrum taking to avoid nonuniform charging. Binding energies relative to at 284.8 eV are reported. Data analysis was performed using MultiPak software developed by Physical Electronics and provided with the XPS instrument.

2.4. Catalytic Activity

The particles produced were tested for the HDS of a VGO in a batch Parr reactor under 3.44 MPa of hydrogen at 653 K during 6 hours. The reactor was first loaded with the mixture of particles (≈0.0580 g.) suspended in VGO (≈20.000 g.). The reactor was tested for leaks and purge with hydrogen, before heating was started at 20 K min−1 and the stirring set at 500 rpm. Hydrogen inlet was kept open during the test in order to avoid pressure dropping due to hydrogen consumption. Once the reaction time was completed, the reactor was cooled down and opened and the liquid products were separated from the catalysts by centrifugation. The liquid products were dissolved in toluene and the sulfur content, of the original VGO and reaction products, was determined on an Antek 9000 NS elemental analyzer.

3. Results and Discussion

Particle sizes obtained by DLS, for the as-prepared solids, are presented in Table 2 along with the corresponding standard deviation. Solid are named using the following nomenclature: -PVP(X), where value in parenthesis indicates the concentration of PVP used in the preparation. Each value is the average of 3 preparations except for the solid -PVP(0.01), where only particles from two different preparations were measured, and the solid made without PVP from which only one preparation was done. The sample with no PVP was obtained following the same methodology as the samples with PVP, but with no addition of the capping agent solution.

As can be seen, particles with diameters as small as 10 nm were produced. The use of PVP concentrations higher than 0.1 mM does not produce significant changes in the particle size (see Table 2). A similar effect (i.e., increasing the concentration of the capping agent produces only minor changes in particle diameter) has been reported for the preparation of CdS nanoparticles using glucose as capping agent [32] and ZnS with mercaptoethanol [33].

The solids -PVP(0.5), -PVP(0.02), -PVP(0.01), and -PVP(0) were chosen for the preparation of the particles to study the catalytic activity dependence on particle size.

In order to confirm the particle size obtained by DLS, TEM analysis was also performed (however a statistical analysis was not done) and micrographs of -PVP(0.5), -PVP(0.02), and -PVP(0.01) are shown in Figure 1, where small spherical particles with diameters around 17 nm and 118 nm are observed for -PVP(0.5) and -PVP(0.02), respectively. These values estimated by TEM are in good agreement with those obtained by DLS (see Table 2). For the solid -PVP(0.01) particles with the shape of small rods were obtained. Why the produce particles presented that shape is not clear. The rods are around 366 nm long and 102 nm thick. The solid prepared without the use of the capping agent was observed by light microscopy (Figure 2), and very large particles with diameters of around 10 μm, forming agglomerates, were observed.

The particles before reaction were also analyzed by XRD (Figure 3). The distinctive amorphous pattern of MoS3 [34], which is the expected product according to reaction , is clearly observed for the samples prepared with PVP.

XRD patterns for the solids after reaction (Figure 3) show very broad lines characteristic of amorphous materials, consistent with the low particle diameter observed by DLS and TEM. The (002) basal plane peaks (line at around 15 2θ°) show very low intensity indicating low or no stacking at all and highly disordered packing of MoS2 slabs. Other broad peaks at higher angle (2θ > 30°), also characteristic of MoS2, are observed. TEM image obtained for the -PVP(0.02) after reaction (not shown) shows the solid to be mainly composed of single slabs, characteristic of the MoS2 structure, in a much disorganized pattern.

According to the geometrical model proposed by Kasztelan et al. [35], the active site on a supported active phase is a distribution of very small two-dimensional patches of MoS2. Eijsbouts et al. [36] reported a linear correlation between the hydrodenitrogenation activity (HDN) and the number of accessible corners and edges in the MoS2 for supported NiMo/Al2O3 commercial catalysts. Also, it has been suggested [37] that the activity of the particles is related not only to the amount of molybdenum atoms found at the edges, where sulfur vacancies can be formed easier, but also to the stacking number. Those results seem to indicate that very small (15–10 nm) nonstacked nanoparticles, like the ones obtained in this work, with a higher density of edge sites should produce higher HDS activity.

X-ray photoelectron spectra are shown in Figures 4 and 5. MoS3 synthesized by acidification of ATTM [34] solution is also included. This sample was prepared in order to compare the corresponding binding energies with those of the precipitated nanoparticles before reaction since the binding energy values reported in the literature vary widely [3747].

Figure 4(a) shows the Mo3d signal for the precipitated MoS3. A main doublet located at 228.77 eV (Mo3d5/2) and 232.00 eV (Mo3d3/2) that should correspond to MoS3 is observed, and the binding energies are close to the ones reported by Stevens and Edmonds [43] and Vrubel et al. [46] for this compound. Another less intense doublet is observed at 230.13 eV for Mo3d5/2 and 233.54 eV for Mo3d3/2 which can be assigned to Mo oxysulfide phase () [39] or to Mo2S5 [8, 44]. Mo2S5 could have been formed while the sample was dried at 353 K under vacuum overnight. Finally, a third minor doublet (231.81 eV for Mo3d5/2 and 235.09 eV for Mo3d3/2) can be assigned to MoVI oxide phase. Bands at lower binding energy that correspond to S2s photoemission with contribution for S2− and at 226.13 and 227.40, respectively, are also obtained.

The high resolution X-ray photoelectron spectrum for the sample -PVP(0.5) before reaction (Figure 4(b)) also shows three 3d5/2-3d3/2 spin-orbit coupling signals. Similar doublets at about the same binding energies as the ones obtained for MoS3 are seen (see Table 3); however, the signal due to or Mo2S5 (230.11 eV for Mo3d5/2 and 233.39 eV for Mo3d3/2) is more intense in this case, suggesting that the prepared MoS3 nanoparticles are less pure than the ones prepared by precipitation. These results indicate that the nanoparticles, before reaction, are mainly constituted by MoS3, which is also supported by the XRD results (Figure 3). The third spin-orbit doublet for Mo3d in Figure 4(b) (231.93 eV for Mo3d5/2 and 235.13 eV for Mo3d3/2) may correspond to MoO3 due to partial oxidation of the sample during manipulation.

The spectrum for the sample after reaction has a simpler structure than the Mo3d signal of the one before reaction, showing a couple of doublets. The first one has binding energy values of 229.52 and 232.66 eV, characteristic of MoS2 [48], and a minor second doublet assigned to MoO3, possibly formed as a result of sample handling. Also signal at 226.8 eV, which could be assigned to S2s for MoS2, is observed.

Figure 5(a) shows the XPS S2p signal for MoS3. Two doublets corresponding to the S2p3/2 and S2p1/2, are observed. The first one with the 2p3/2 peak centered at 161.45 eV can be assigned to S2−, and the second one with a binding energy of 162.68 eV can be attributed to polysulfide (). The first signal shows higher intensity and the ratio (S2−)/() is 1.4. Similar signals have been reported in the literature for MoS3 with (S2−)/() ratios from 0.5 to 2.4 [3943, 46]. There is also a small signal at higher binding energy which can be due to SO4− formed by oxidation of the sample during manipulation.

For the solids before reaction (Figure 5(b)) two S2p3/2-S2p1/2 are also observed with similar binding energies; however, the (S2−)/() is 0.88 and the difference could be due to the presence of some other sulfides or oxysulfides as evidenced in the Mo3d spectrum. After reaction, the solid shows S2p spectrum with one main doublet characteristic of S2− in MoS2 with another small one, possibly due to impurities.

Table 3 presents a summary of the XPS binding energies for both samples and the assignment for the Mo3d and S2p signals. It is important to say that the S/Mo ratio of 1.88 for sample after reaction is similar to previously reported S/Mo ratios for bulk MoS2 [37], while the ratio for the sample before reaction is 2.97, which additionally supports the presence of MoS3.

From all the characterization shown, it is clear that during the reaction almost all (if not all) Mo is converted to MoS2. The formation of MoS2 from an intermediate MoS3 phase has been reported to easily happens when MoS3 is heated above 623 K [29]. The XPS and XRD results support the formation of MoS2 on sample -PVP(0.5) after the reaction, which was performed at 653 K.

The catalytic activity of the particles produced is shown in Table 4.

Important differences in the activity of these catalytic particles are observed in the HDS. A trend can be observed for the HDS activity of the solids with respect to the particle diameter (see Figure 6). The -PVP(0.5) solid with the smallest particle size is the one showing highest HDS conversion. The higher activity of the smaller particles could be the result of an increase in the number of available surface atoms located at the edges of the molybdenum sulfide clusters which have been proposed to be the active sites for HDS [35]. Figure 6 presents a plot of HDS conversion versus the particle diameter (in logarithmic scale) where a clear trend can be observed; however -PVP(0.02) and -PVP(0.01) show similar conversions. This effect might be due to the different morphology evidenced by TEM for these two materials. In fact, it has been reported that the morphology of MoS2 nanoparticles can lead to variations in the hydroconversion of vacuum residue [15]. Molybdenum sulfide particles obtained with diameters in the order of 13 nm show double the HDS conversion compared with the material prepared without addition of capping agent.

4. Conclusions

The present study shows the feasibility of preparing molybdenum sulfide nanoparticles with controlled diameters from 10 to 400 nm and very low stacking, as evidence by DLS, SEM, TEM, and XRD.

A trend was observed for the relationship between particle size and HDS conversion. For particles with diameters around 13 nm, the % of HDS was double that of bulk molybdenum sulfide and also significantly higher when compared with particles with diameters around 100–400 nm.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The authors want to acknowledge the Natural Sciences and Engineering Research Council of Canada (NSERC), Nexen-CNOOC Ltd., and Alberta Innovates-Energy and Environment Solutions (AIEES) for the financial support provided through the NSERC/NEXEN/AIEES Industrial Research Chair in Catalysis for Bitumen Upgrading. Also, the contribution of facilities from the Canada Foundation for Innovation, the Institute for Sustainable Energy, Environment and Economy, and the Schulich School of Engineering at the University of Calgary are greatly appreciated. Special thanks are also extended to Dr. Tobias Fürstenhaupt for his help with the Electron Microscopy determinations and Ms. Jodi David for revision of the manuscript.