Simultaneous Adsorption of 4,6-Dimethyldibenzothiophene and Quinoline over Nickel and Boron Modiﬁed Gamma-Al 2 O 3 Adsorbent

: The simultaneous adsorption of quinoline and 4,6-dimethyldibenzothiophene over adsorbents, based on alumina modiﬁed with boron and nickel under ambient temperature and pressure, was studied. The adsorbents were characterized by BET speciﬁc surface area, a potentiometric method for the determination of acid strength, electrophoretic migration, and X-ray di ﬀ raction. The results showed that the adsorbent containing nickel had better adsorption capacity than the adsorbent modiﬁed with nickel and boron, which was attributed to its greater acidity and ability to generate π -complexation between the adsorbent and the molecules. In terms of selectivity, quinoline was more adsorbed than 4,6-dimethyldibenzothiophene in all systems, due to the basic nature of quinoline. The experimental data in all cases were adjusted by three kinetic models (Yoon–Nelson, Yan and Thomas), and the regression coe ﬃ cients in all models were close to one. Finally, the values of the kinetic constant obtained by the Thomas model were used to relate the adsorption capacity results.


Introduction
Diesel is a complex mixture of hydrocarbons formed by organosulfur and organonitrogen molecules, among other compounds. The combustion of these molecules is responsible for environmental pollution, acid rain and health problems. Also, the sulfur compounds can poison the catalysts used in catalytic converters [1], while the nitrogen compounds inhibit the removal of sulfur compounds [2]. To eliminate these pollutants, different hydrotreating (HT) reactions, specifically hydrodesulfurization (HDS) and hydrodenitrogenation (HDN), are carried out during the refining of petroleum. These processes use Co-Mo/γ-Al 2 O 3 , Ni-Mo/γ-Al 2 O 3 or Ni-W/γ-Al 2 O 3 catalysts at high temperatures (300-400 • C) and under elevated hydrogen pressures (20-100 atm).
However, due to depletion of continental crude oil deposits, the deeper reservoirs of petroleum contain a large number of sulfur and nitrogen polynuclear aromatic molecules that are refractory to the conventional methodologies. To eliminate these compounds, new costly processes of HT that produce considerable amounts of carbon must be applied. In order to economize the desulfurization or denitrogenation processes, other milder methods must be developed. In this regard, several research groups have studied the elimination of these types of molecules by adsorption [3][4][5][6][7], which, unlike HT reactions, is carried out at ambient temperature and pressure. Despite there being a wide variety

Preparation of the Adsorbents
Gamma alumina (BASF D10-10) was modified with 0.4 wt% of boron (B(0.4) γ-Al 2 O 3 ) and with 4.0 wt% of nickel (Ni(4.0) γ-Al 2 O 3 ), using the wet impregnation method. The incorporation of boron was realized with H 3 BO 3 (Sigma Aldrich p.a.) using the required amount to obtain a 0.4 wt% of boron. Then, the adsorbent was dried at 378 K overnight and then calcined in air at 823 K for 4.5 h. The Ni-modified alumina adsorbent was obtained with the required amount of Ni(NO 3 ) 2 ·6H 2 O to obtain a 4.0 wt% of nickel. This adsorbent was dried at 378 K overnight and calcined in air at 823 K for 4.5 h. The Ni and boron modified alumina adsorbent (Ni(4.0)B(0.4) γ-Al 2 O 3 ) was obtained by wet impregnation of Ni(NO 3 ) 2 ·6H 2 O over B(0.4) γ-Al 2 O 3 , then were dried at 378 K overnight and calcined in air at 823 K for 4.5 h.

Characterizations of the Adsorbents
The determination of the specific surface area (S BET ) [34] for each adsorbent was done in a Micromeritics ASAP-2010 apparatus for volumetric nitrogen adsorption-desorption, using 0.2-0.4 g of each sample, which were degassed at 120 • C for 18 h, reaching a final pressure of 1 × 10 −3 mmHg. Subsequently, the analysis was performed within a glass slide at 77 K. The determination of the Isoelectric Point (IEP) was realized by electrophoretic migration, as in previous studies [15]. The surface acidity of the adsorbents was measured by potentiometric titration [35]. The X-ray diffraction patterns of the powder samples were made in a Panalytical X'Pert diffractometer equipped with an Xcelerator detector with automatic divergence slits and Cu K α1/α2 radiation (40 kV, 40 mA). The K-beta radiation of Cu was excluded using nickel filters and the K-alpha2 radiation was removed arithmetically using the Panalytical HighScore Plus Software. The positions of the peaks and the profiles were adjusted by the Pseudo-Voigt function and the WinXPow software. The identification of the phases was carried out using the PDF-2 database of the International Center of Diffraction Data (ICDD).

Adsorption Experiments
Adsorption experiments were performed at ambient temperature and pressure using a vertical Pyrex reactor equipped with a supporting glass porous disk. A bed of 500 mg of adsorbent was placed in the reactor. Prior to the determination of the adsorption capacity (q e ), a surface cleaning process was performed for each adsorbent, with a 20 mL min −1 flow of pure Ar at 378 K for 1 h and subsequent cooling to room temperature. For the adsorption experiments, the liquid flow was driven into the reactor by means of a variable flow peristaltic mini-pump. The experiment was performed using a mixture of 250 ppmw of QN and 250 ppmw 4,6-DMDBT in isooctane at a feed rate (Q) of 0.5 mL min −1 . Samples were collected every 10 min, until saturation of the adsorbent was reached; the total time and number of samples depended on the adsorption capacity of the adsorbent.
The samples were analyzed by gas chromatography with a flame ionization detector (GC-FID), using a Shimadzu GC-2010 equipped with an SPB-5 capillary column (L 30 m, I.D. 0.25 mm, Film 0.25 mm), and the following conditions: detector at 553 K, injector at 523 K, and carrier flow, 30 mL min −1 (N 2 ). The column temperature was set to increase from 433 K to 513 K at a rate of 20 K min −1 ; 1.0 µL of the sample volume was injected, using an autosampler, for each GC-FID run.
The adsorption results for the different systems studied were represented as adsorption curves (adsorbed fraction vs eluted volume), where the adsorbed fraction is expressed as 1-(C t C 0 −1 ), where C t corresponds to the concentration of QN at time t (after adsorption) and C 0 is the feed concentration. The QN and 4,6-DMDBT adsorbed per unit mass of the adsorbent (q e ), and was calculated using Equation (1). The selectivity was studied with the relative selectivity factor used by Kim et al. [4] (Equation (2)), where q ei and q en are the adsorption capacity for QN and 4,6-DMDBT, respectively.

Adsorption Kinetics
To estimate the kinetic parameters that govern the adsorption process in a fixed bed system, several mathematical models have been studied by different research groups, and many of them are summarized in a review by Xu et al. [36]. Each of these models describes the dynamic behavior of the adsorption curves and attempts to relate experimental data with those calculated by the prediction model. For the estimation of the parameters of each model, the non-linear analysis of the curves (C t C 0 −1 vs. t) was performed, using the nonlinear curve fit available in Origin 8.0 Software.

Yoon-Nelson Model
The Yoon-Nelson model is extremely concise in its form for the prediction of mass transfer phenomena. This model assumes the adsorption rate is proportional to the adsorption capacity of the solid and the initial concentration of feed. The nonlinear form of this model is presented in Equation (3).
where k YN correspond to the Yoon-Nelson velocity constant (min −1 ), and τ (min) denotes the adsorption time when 50% has been reached in adsorbate extraction.

Yan Model
This model is a variation of the Dose-Response model, initially developed for pharmacological studies. This model has recently been used to describe the kinetics of metal adsorption [32], however it is interesting to study its behavior in the study of adsorption of organic molecules. This model is represented by Equation (4).
where q e is the maximum concentration of the adsorbate on the surface of the adsorbent (mg g −1 ), and a corresponds to the dimensionless constant of the model.

Thomas Model
This model is used to estimate the adsorption capacity (q e ) of an adsorbent and predict adsorption curves, assuming second order reversible kinetics. This model is based on the Langmuir isotherms, so the adsorption process is not limited only by the reaction on the surface, but also by the phenomena of mass transfer in the liquid-solid interface. This model is described by Equation (5).
where C t corresponds to the concentration of the adsorbate at time t, k th corresponds the Thomas velocity constant (mL min −1 mg −1 ), Q is the feed flow (mL min −1 ), q e is the maximum concentration of the adsorbate on the surface of the adsorbent (mg g −1 ), m is the mass of adsorbent (g), and C 0 corresponds to the initial concentration of the adsorbate (mg mL −1 ).

Results and Discussion
The results of characterization by BET method, acid strength and isoelectric point of adsorbents are shown in Table 1. Initially, it can be observed that the specific surface area of γ-alumina decreases from 213 to 170 and 209 m 2 g −1 for B(0.4) γ-Al 2 O 3 and Ni(4.0) γ-Al 2 O 3 , respectively, when boron and nickel are incorporated into this support. For the nickel modified adsorbent, the slight decreases in the specific surface area could be related to the pore obstruction of support by the presence of the metal. For the boron modified adsorbent, despite the decrease in the specific surface area values, an increase in the pore diameter values was observed. This behavior can be due to the formation of new pores in the γ-Al 2 O 3 by the incorporation of boron or for the obstruction of pores with minor diameter, which was also observed by Engelso et al. [37] when studying the boron incorporation on alumina. These authors suggest that with boron incorporation, the diameters of pores < 10 nm will be diminished, and new pores > 50 nm will appear.
On the other hand, the acid strength values listed in the Table 1, shows an increase in the values of B(0.4) γ-Al 2 O 3 (−33.1) and Ni(4.0) γ-Al 2 O 3 (−6.8) adsorbents with respect to pure alumina (−61.7). This behavior can be related with the formation of new Lewis and Brönsted sites by the incorporation of nickel to the alumina [25,26], and the formation of new Brönsted sites caused by the incorporation of boron.
For the case of electrophoretic migration measures, a decrease in the values of IEP of alumina (8.02) to 7.93 and 6.80 is observed, when boron and nickel were incorporated to this support, respectively (see Table 1). In the case of nickel, this decrease can be due to the formation of nickel aluminate type (NiAl 2 O 4 ) species with lower values of IEP near to 2.7 [38]. For boron incorporation, the decrease in the IEP value can be explained by the formation of negative surface complexes with boric acid [39]. On the other hand, the acid strength values listed in the Table 1, shows an increase in the values of B(0.4) γ-Al2O3 (−33.1) and Ni(4.0) γ-Al2O3 (−6.8) adsorbents with respect to pure alumina (−61.7). This behavior can be related with the formation of new Lewis and Brönsted sites by the incorporation of nickel to the alumina [25,26], and the formation of new Brönsted sites caused by the incorporation of boron.
For the case of electrophoretic migration measures, a decrease in the values of IEP of alumina (8.02) to 7.93 and 6.80 is observed, when boron and nickel were incorporated to this support, respectively (see Table 1). In the case of nickel, this decrease can be due to the formation of nickel aluminate type (NiAl2O4) species with lower values of IEP near to 2.7 [38]. For boron incorporation, the decrease in the IEP value can be explained by the formation of negative surface complexes with boric acid [39]. For the Ni  In the case of alumina, the characteristic diffraction peaks of γ-Al2O3 (included in Figure 1 as reference) can be clearly seen. When boron is incorporated into alumina, the diffractogram of the B(0.4) γ-Al2O3 adsorbent does not show significant variation with respect to pure γ-Al2O3, presenting the same diffraction peaks, and no crystalline boron oxide species (B2O3) were detected. Therefore, in the case of the B(0.4) γ-Al2O3 adsorbent, it could be noted that there are no crystalline or agglomerated B species, and it is proposed that B species are highly dispersed on the surface of γ-Al2O3, possibly forming a thin coating or monolayer, which cannot be detected by XRD. The above can be based on the low loading of B (0.4 wt %) with respect to the high specific surface area of the alumina (see Table  1). In addition, in the diffractogram of the B(0.4) γ-Al2O3 adsorbent, it is possible to observe a slight increase in the intensity of the diffraction peaks, that is, a slight increase in the crystallinity of γ-Al2O3. One reason for this increase in the crystallinity of γ-Al2O3 could be related to the calcination temperature of the B(0.4) γ-Al2O3 adsorbent (823 K for 4.5 h).
In the case of the Ni(4.0) γ-Al2O3 adsorbent, it can be seen in Figure 1 that its diffractogram is similar to that obtained in pure γ-Al2O3, but unlike that with B(0.4) γ-Al2O3 adsorbent, in this adsorbent the incorporation of nickel on the alumina produces a notable decrease in intensity of the diffraction peaks corresponding to γ-Al2O3, as seen in Figure 1, which represents a loss of crystallinity in Ni(4.0) γ-Al2O3 adsorbent. The above behavior could be explained considering the higher Ni load than B load (4.0 > 0.4 wt%) on the alumina support. However, this 4.0 wt% of Ni loading is not so high as to exceed the dispersion capacity of alumina support, remembering that its specific surface area is 213 m 2 g −1 (see Table 1). In fact, the Ni(4.0) γ-Al2O3 adsorbent has a specific surface area of 209 m 2 g −1 , almost the same as the specific surface area of γ-Al2O3. Therefore, Ni loading (4.0 wt%) does not affects alumina specific surface area but does affect its crystallinity as mentioned above. In addition, in the diffractogram of the Ni(4.0) γ-Al2O3 adsorbent, it is not possible to observe nickel In the case of alumina, the characteristic diffraction peaks of γ-Al 2 O 3 (included in Figure 1 as reference) can be clearly seen. When boron is incorporated into alumina, the diffractogram of the B(0.4) γ-Al 2 O 3 adsorbent does not show significant variation with respect to pure γ-Al 2 O 3 , presenting the same diffraction peaks, and no crystalline boron oxide species (B 2 O 3 ) were detected. Therefore, in the case of the B(0.4) γ-Al 2 O 3 adsorbent, it could be noted that there are no crystalline or agglomerated B species, and it is proposed that B species are highly dispersed on the surface of γ-Al 2 O 3 , possibly forming a thin coating or monolayer, which cannot be detected by XRD. The above can be based on the low loading of B (0.4 wt %) with respect to the high specific surface area of the alumina (see Table 1). In addition, in the diffractogram of the B(0.4) γ-Al 2 O 3 adsorbent, it is possible to observe a slight increase in the intensity of the diffraction peaks, that is, a slight increase in the crystallinity of γ-Al 2 O 3 . One reason for this increase in the crystallinity of γ-Al 2 O 3 could be related to the calcination temperature of the B(0.4) γ-Al 2 O 3 adsorbent (823 K for 4.5 h).
In the case of the Ni(4.0) γ-Al 2 O 3 adsorbent, it can be seen in Figure 1 that its diffractogram is similar to that obtained in pure γ-Al 2 O 3 , but unlike that with B(0.4) γ-Al 2 O 3 adsorbent, in this adsorbent the incorporation of nickel on the alumina produces a notable decrease in intensity of the diffraction peaks corresponding to γ-Al 2 O 3 , as seen in Figure 1, which represents a loss of crystallinity in Ni(4.0) γ-Al 2 O 3 adsorbent. The above behavior could be explained considering the higher Ni load than B load (4.0 > 0.4 wt%) on the alumina support. However, this 4.0 wt% of Ni loading is not so high as to exceed the dispersion capacity of alumina support, remembering that its specific surface area is 213 m 2 g −1 (see Table 1). In fact, the Ni(4.0) γ-Al 2 O 3 adsorbent has a specific surface area of 209 m 2 g −1 , almost the same as the specific surface area of γ-Al 2 O 3 . Therefore, Ni loading (4.0 wt%) does not affects alumina specific surface area but does affect its crystallinity as mentioned above. In addition, in the diffractogram of the Ni(4.0) γ-Al 2 O 3 adsorbent, it is not possible to observe nickel oxides species (NiO, Ni 2 O 3 ), so it is likely that Ni species formed on the alumina are highly dispersed, cannot be detected by XRD, or Ni species are in an agglomerate with low crystallinity or amorphous state, but this last explanation is less probable. Finally, in the case of the diffractogram of the Ni(4.0) B(0.4) γ-Al 2 O 3 adsorbent, it presents a similar result to the Ni(4.0) γ-Al 2 O 3 adsorbent (see Figure 1), i.e., a decrease in the intensity of the diffraction peaks corresponding to γ-Al 2 O 3 . The XRD analysis of the Ni(4.0) B(0.4) γ-Al 2 O 3 adsorbent could be summarized as a weighted combination of the B(0.4) γ-Al 2 O 3 and Ni(4.0) γ-Al 2 O 3 adsorbents, considering that the Ni(4.0) B(0.4) γ-Al 2 O 3 adsorbent presents different metals loads on alumina (Ni loading is 10 times more than B loading), and as it was noted above, the incorporation of B on γ-Al 2 O 3 produces an increase in the crystallinity of alumina, the subsequent incorporation of Ni on B(0.4) γ-Al 2 O 3 produces a loss in its crystallinity.
In addition, in the diffractogram of the Ni(4.0) B(0.4) γ-Al 2 O 3 adsorbent, it is not possible to identify the presence of boron and nickel oxides, which would indicate that the presence of B at the concentration used, does not affect the dispersion of Ni. It is worth recalling the BET area results of these adsorbents (see Table 1), where it is shown that the incorporation of B on alumina produce a decrease in its specific surface area, from 213 to 170 m 2 g −1 , unlike the incorporation of Ni that does not produce changes in the specific surface area of γ-Al 2 O 3 , from 213 to 209 m 2 g −1 . It is interesting to note that when Ni is incorporated on boron-modified alumina, which corresponds to Ni(4.0) B(0.4) γ-Al 2 O 3 adsorbent, there is an increase in its specific surface area, from 170 to 197 m 2 g −1 (see Table 1).
This last result with the XRD results (see Figure 1), allows us to propose that Ni species would be interacting with surface B species, forming highly dispersed Ni-B species, since they are not detected by XRD, but at the same time this Ni-B species could be responsible of increasing the specific surface area of the Ni(4.0) B(0.4) γ-Al 2 O 3 adsorbent with respect to the B(0.4) γ-Al 2 O 3 adsorbent. The possible existence of Ni-B species interacting on the surface of the alumina will be discussed later.

Adsorption Tests
The breakthrough curves, obtained for the adsorption of mixed 4,6-DMDBT and QN on the different adsorbents, are shown in Figure 2. The experimental adsorption capacity (q e ) and selectivity (α i-n ), determined by Equations (1) and (2), respectively for 4,6-DMDBT and QN are summarized in Table 2. The results of this table show that all systems studied adsorbed higher amounts of QN than 4,6-DMDBT. This behavior can be attributed to the strong interaction between the Al(III) sites on γ-Al 2 O 3 with the electronic pair of the nitrogen atom of QN that acts as a Lewis base [40]. In the case of 4,6-DMDBT, the steric hindrance effect of the methyl groups of this molecule generate a decrease in the interaction between OH functional groups of γ-Al 2 O 3 with the aromatic rings of 4,6-DMDBT, producing lower adsorption.  The comparison of the adsorption of 4,6-DMDBT and QN on the other supports used in hydrotreating (HT) reactions indicate that the adsorption capacity is dependent on factors such as specific surface area, adsorption pH, density and acid strength, among others [41,42]. For example, we studied the adsorption of 4,6-dimethyldibenzothiophene (4,6-DMDBT) on different supports (SiO2, Al2O3, zeolite 4A, zeolite 13X, natural zeolite, and activated carbon). Despite that the higher adsorption was observed on active carbon due to the high specific surface area, the alumina also exhibits considerable adsorption capacity, which would be related to the presence of OHfunctional groups, which may also have a certain adsorption capacity [43].
On the other hand, the incorporation of boron or nickel to γ-Al2O3 increases the qe values of QN and 4,6-DMDBT. The increases of adsorption values of QN and 4,6-DMDBT by the B(0.4) γ-Al2O3 system can be attributed to the increase of the acidity of the adsorbent, produced by the incorporation of boron to alumina [44]. In the case of the Ni (4.0) γ-Al2O3 system, the better adsorption of both molecules should be due to the proven ability of Ni to interact with nitrogen or sulfur organic rings through the π-complexation mechanism. In this mechanism, the metallic cation can form the usual σ bonds with their s-orbitals, while their d-orbitals can back-donate electron density to the antibonding π-orbitals of the sulfur or nitrogen rings [7]. The comparison of the adsorption of 4,6-DMDBT and QN on the other supports used in hydrotreating (HT) reactions indicate that the adsorption capacity is dependent on factors such as specific surface area, adsorption pH, density and acid strength, among others [41,42]. For example, we studied the adsorption of 4,6-dimethyldibenzothiophene (4,6-DMDBT) on different supports (SiO 2 , Al 2 O 3 , zeolite 4A, zeolite 13X, natural zeolite, and activated carbon). Despite that the higher adsorption was observed on active carbon due to the high specific surface area, the alumina also exhibits considerable adsorption capacity, which would be related to the presence of OHfunctional groups, which may also have a certain adsorption capacity [43].
On the other hand, the incorporation of boron or nickel to γ-Al 2 O 3 increases the q e values of QN and 4,6-DMDBT. The increases of adsorption values of QN and 4,6-DMDBT by the B(0.4) γ-Al 2 O 3 system can be attributed to the increase of the acidity of the adsorbent, produced by the incorporation of boron to alumina [44]. In the case of the Ni (4.0) γ-Al 2 O 3 system, the better adsorption of both molecules should be due to the proven ability of Ni to interact with nitrogen or sulfur organic rings through the π-complexation mechanism. In this mechanism, the metallic cation can form the usual σ bonds with their s-orbitals, while their d-orbitals can back-donate electron density to the antibonding π-orbitals of the sulfur or nitrogen rings [7].
When the nickel was incorporated into B(0.4) γ-Al 2 O 3 , a decrease in the values of q e of QN and 4,6-DMDBT was observed. These results demonstrated that the boron did not improve the dispersion of nickel on the alumina to the concentrations of metals used, as was evidenced in the other similar works [45]. The change in the acidity and the IEP values of Ni(4.0) B(0.4) γ-Al 2 O 3 adsorbent in comparison with the Ni (4.0) γ-Al 2 O 3 could indicates a possible formation of a mixed species of nickel-boron, with lower adsorption capacity of QN and 4,6-DMDBT. In this context, Graff et al. [46], studied the formation of triborate nickel complex in an experiment based on electrochemical reactions and pH monitoring, in which nickel ions were gradually formed by the oxidation of a nickel metal electrode in a solution of boric acid. The authors demonstrated that nickel cations could react with the triborate, forming a complex between both elements. However, the possible formation of a boron nickel complex on the surface of the alumina must be probed with additional characterization techniques.
In terms of selectivity (α i-n ), it is possible to observe in all cases that QN is favored in the adsorption with respect to 4,6-DMDBT. However, this selectivity varies depending on the adsorbent used. This behavior have been observed by diverse authors, establishing that the π-complexation adsorbents exhibit high capacity, but low selectivity for sulfur compounds as the result of competitive adsorption of the aromatic compounds [6,[47][48][49][50][51][52][53]. The selectivity to QN for B(0.4) γ-Al 2 O 3 adsorbent decreases with respect to unmodified γ-Al 2 O 3 , possibly due to a decrease in the Lewis sites of Al(III) when boron is incorporated. Subsequently, the incorporation of nickel generates an increase in the selectivity of QN, probably due to the increase in the acidity of the adsorbents, which can be attributed to the generation of Lewis sites of Ni(II) on the surface of the alumina. Finally, in a mixed system of boron and nickel, the selectivity decreases, due to a combined effect of both elements. Similar results of selectivity were obtained by authors that studied the simultaneous adsorption of sulfur and nitrogen molecules. In this respect, Thaligari et al. studied the simultaneous removal of the dibenzothiophene and quinoline by using zinc impregnated granular activated carbon (Zn-GAC). The study, which used a Taguchi's experimental design methodology, demonstrated a higher adsorption of quinoline in comparison to dibenzothiophene [54]. Likewise, Mohammadian et al. found that the adsorption selectivity of all adsorbents for nitrogen was slightly higher than for sulfur, when studying the desulfurization and denitrogenation of a model fuel containing benzothiophene (BT), dibenzothiophene (DBT), quinoline and carbazole, using MSU-S, CeO 2 -MSU-S, and Cu 2 O-MSU-S adsorbents [22].

Adsorption Kinetics
The kinetics parameters obtained by the three models studied for 4,6-DMDBT and QN adsorption, are shown in Tables 3 and 4, respectively. Firstly, the values of the regression coefficient (r 2 ) obtained for all systems studied are very close to one, which demonstrates that the three models agree with the experimental data.
By using the Yoon-Nelson model, it can be noted that the adsorption capacity is directly related to the values of tau, which corresponds to the time in which 50% of the adsorption has been reached. Therefore, this model can be used to estimate how high the adsorption capacity of the material will be, and how long the saturation of the adsorbent takes. Table 3. Kinetics parameters obtained by the three models studied for 4,6-DMDBT.

Exp
Yoon-Nelson Yan Thomas  Table 4. Kinetics parameters obtained by the three models studied for QN. In the case of the Yan model, although this model as mentioned above is used in pharmacological studies, when applied in the adsorption of the different molecules under study, the r 2 values fits well with the experimental data, so this model could be used to also predict the adsorption curves of these compounds.

Yoon-Nelson Yan Thomas
As can see in Moreover, the Thomas model delivers more information about the adsorption process. Taking into account the considerations of the model, it is possible to establish that the adsorption process of the molecules follows a reversible kinetic of second order, where all adsorption sites are equivalent and each of these sites can contain a molecule adsorbed which does not interact with adjacent sites. Likewise, the information obtained allows the determination of the kinetic constant of the adsorption process of 4,6-DMDBT and QN on the studied adsorbents. On this matter, it is possible to establish that the kinetic constant obtained by this model is inversely proportional to the adsorption capacity of the different materials studied, due probably to the longer contact time used in the interaction between the molecule and the adsorbent. This is seen more clearly in Figure 3.

Conclusions
Among the adsorbents studied in this work, Ni(4.0) γ-Al2O3 has the better performance for the simultaneous adsorption of quinoline and 4,6-DMDBT, this due to its higher acidity and capacity to generate π-complexation between the molecule and the adsorbent. The Ni(4.0) B(0.4) γ-Al2O3 adsorbent has a lower adsorption capacity of both molecules, presumably due to formation of a mixed species of nickel-boron. In terms of selectivity, QN is favored in the adsorption with respect to 4,6-DMDBT, due to its basic nature. In all adsorbent systems studied, the adsorption curves fitted well with the three models, which show high r 2 values. However, the Thomas model predicts better the adsorption capacity values, and also gives more information about the adsorption process.

Conclusions
Among the adsorbents studied in this work, Ni(4.0) γ-Al 2 O 3 has the better performance for the simultaneous adsorption of quinoline and 4,6-DMDBT, this due to its higher acidity and capacity to generate π-complexation between the molecule and the adsorbent. The Ni(4.0) B(0.4) γ-Al 2 O 3 adsorbent has a lower adsorption capacity of both molecules, presumably due to formation of a mixed species of nickel-boron. In terms of selectivity, QN is favored in the adsorption with respect to