Studying Rotational Mobility of V O Complexes in Atmospheric

Apr 5, 2017 - Shell buys stake in solar developer Silicon Ranch. Shell will buy a 43.8% stake in the U.S. solar developer Silicon Ranch from private i...
0 downloads 11 Views 2MB Size
Subscriber access provided by University of Newcastle, Australia

Article

Studying Rotational Mobility of V=O Complexes in Atmospheric Residues, Their Resins and Asphaltenes by ESR Qingyan Cui, Koji Nakabayashi, Xiaoliang Ma, Jin Miyawaki, Keiko Ideta, Yoshika Tennichi, Morio Ueda, Adel Al-Mutairi, Abdulazim MJ Marafi, Joo-Il Park, Seong-Ho Yoon, and Isao Mochida Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.6b03279 • Publication Date (Web): 05 Apr 2017 Downloaded from http://pubs.acs.org on April 7, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Studying Rotational Mobility of V=O Complexes in Atmospheric Residues, Their Resins and Asphaltenes by ESR Qingyan Cui,b Koji Nakabayashi,b Xiaoliang Ma,c Jin Miyawaki,b Keiko Ideta,b Yoshika Tennichi,a Morio Ueda,a Adel Al-Mutairi,c Abdulazim MJ Marafi,c Joo-Il Park,c Seong-Ho Yoon b and Isao Mochida*,a a b

Kyushu Environmental Evaluation Association, Fukuoka, Japan Institute for Materials Chemistry and Engineering, Kyushu University, Kasuga,

Fukuoka, Japan c

Petroleum Research Center, Kuwait Institute for Scientific Research, Safat, Kuwait

ABSTRACT Behaviors of V=O complexes in atmospheric residues (AR) from two Kuwait crudes, their resins and asphaltenes were studied using the electron spin resonance (ESR) at 20 and 100 oC to examine the effects of the surrounding matrix, concentration in solvent and temperature on V=O rotational mobility. The results show that the surrounding molecules in the petroleum fractions constrain the V=O rotational mobility significantly. The constraint on the V=O complexes by the surrounding matrix in different environments increases in the order of resin < AR < asphaltene. Less constraint by the AR than by the asphaltene can be ascribed to the solvent role of the lighter components in the AR. The higher measurement temperature of 100 oC significantly decreases the constraint on V=O complexes by the matrix, while a higher sample concentration in the solvent shows stronger constraints. However, the constraint on the V=O complexes in the asphaltene is hardly moderated, even when dissolved in toluene at 100 oC. Additionally, Kuwait Export Crude atmospheric residues (KEC-AR) and its resin dissolved in toluene show weaker constraints on the V=O complexes than Lower Fars 1

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

atmospheric residues (LF-AR) and its resin, possibly due to higher dispersibility or lower aggregation of KEC-AR and its resin. KEY WORDS: Electron spin resonance, atmospheric residue, asphaltene, resin, V=O complex, rotational mobility. 1. INTRODUCTION Metal species in petroleum have been targeted extensively for removal in the refining process, because they may otherwise be deposited on catalysts, transfer tubes and filters in the downstream steps of petroleum refining.1 Metal species are usually removed through hydrodemetallization (HDM) over Mo or NiMo alumina catalyst, where the formed V2S6 and NiS deposit on the surface of catalysts, which ultimately results in deactivation of the catalyst.2-5 A deeper extent and larger capacity for metal removal have been explored for longer-term continuous operations. Fully understanding the structure and state of the metal complexes, as well as their relationships with surrounding molecules, is important in designing a better HDM catalyst and a more efficient HDM process. Metal species are present mainly in the heavy fractions of petroleum, such as resin and especially asphaltene, which consists of condensed polynuclear aromatics and contains a high concentration of heteroatoms. The analysis of asphaltene has been reported through fractionation with a binary solvent.6-8 Additionally, asphaltene has been studied widely using nuclear magnetic resonance (NMR), X-ray diffraction (XRD) and transmission electron microscopy (TEM) to demonstrate its structure and properties.9-13 It is important to understand the state of the metal complexes in asphaltene and petroleum, which influences the deposition of the metal species together with carbonaceous substances on the catalyst during the HDM process. 2

ACS Paragon Plus Environment

Page 2 of 35

Page 3 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

V and Ni porphyrins are common metal species in petroleum, although V and Ni non-porphyrin, iron sulfide and naphthenic alkali earth metal salts also exist in petroleum.14 The structures of the metal porphyrins have been determined.14,15 Recently, high-resolution Fourier transform ion cyclotron resonance (FT-ICR-MS) with positiveand negative- ion electrospray16-18 and atmospheric pressure photoionization19,20 procedures detected V=O and Ni porphyrin species in crude oil and residues, which may provide helpful information regarding the molecular structure of the metal species. High- temperature gas chromatography with atomic emission detection (GC-AED)

21

has also provided the molecular information on the metal species. Both analyses provided information on isolated metal species. This is a strong basis for understanding metal species at a molecular level. However, V and Ni species co-exist with organic components in petroleum, where there are interactions between the metal species and surrounding organic components through non-covalent, hydrogen bond, and aromatic π-π interactions. Thus, it is also important to study the environment of the metal species and their relationships with surrounding organic components. The structure of metal species in petroleum can be described as a four-tiered structure, as follows. 1) A tetradentate core ligand, where some core ligands can contain atoms other than nitrogen. 2) Pendants on the core ligand, where some of the pendants are bound covalently to large species around the tetradentate core ligand. 3) Axial ligands to the complex, where the ligands can be present in surrounding molecules, which carry the coordinative functions to central metal ions in an axial direction. 4) Physical bonding of the complex with surrounding molecules or molecular groups, to form micelles or stack into the large aromatic sheets.22, 23 High-performance liquid chromatography (HPLC)24-26 and gel permeation 3

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

chromatography (size exclusion chromatography) inductively coupled plasma GPC(SEC)-ICP27-30 have been used to detect metal complexes in the fractions, such as AR, resin and asphaltene, where the metal complexes are bound to matrix molecules or dissolved in the matrix. It is not yet fully understood how the metal complexes exist in petroleum, or how the surrounding matrix interacts with them to change their mobility, although such an understanding is important for improving the metal removal performance of the HDM process. The electron spin resonance (ESR) technique has been used to analyze vanadium (Ⅱ) porphyrin (V=O) complexes in petroleum. Dickson31 et al. proposed using ESR to identify the environment of V=O complexes in petroleum, and to examine g values and hyperfine coupling constants (A) for several square pyramidal environments around vanadyl complexes. Additionally, vanadyl porphyrin concentrations have been estimated using ESR.21,32 However, there are few reports examining the effects of ligands and the surrounding matrix on V=O rotational mobility in petroleum. Campbell and Freed33 measured a series of ESR spectra of V=O complexes dissolved in solvents at different temperatures and simulated them on the basis of the tumbling rates and τ values of the complexes at a given temperature. Wong and Yen34 studied the mobility of V=O complexes in asphaltene by ESR through distinguishing anisotropic and isotropic spectra, and illustrated the mobility through liberating the V=O complexes to be mobile in the asphaltene dissolved in some solvents in ESR measurement. These studies indicate that ESR can be useful for studying the mobility of V=O complexes in different environments. In the present study, the effects of solvents (toluene and tetrahydrofuran (THF)), measurement temperature, concentration of the samples in solution and surrounding 4

ACS Paragon Plus Environment

Page 4 of 35

Page 5 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

organic matrixes on the mobility of V=O complexes were studied by distinguishing the ESR anisotropic and isotropic spectra, extending the approaches of Campbell and Freed33 and Wong and Yen.34 ESR spectra of Etio and tetraphenyl porphyrin V=O complexes dissolved in toluene, were determined at a series of temperatures as a frame of reference to classify the profiles of ESR spectra of the petroleum fractions according to the tumbling rate, as reported by Campbell and Freed.33 Such mobility information may provide insight into the interaction between V=O complexes and their surrounding matrixes, which may be significant in improving the metal removal performance of the HDM process. 2. EXPERIMENTAL SECTION 2.1. Sample Preparation Two kinds of Kuwait atmospheric residues were from Lower Fars (LF) crude and Kuwait Export Crude (KEC). The carbon (C), hydrogen (H), and nitrogen (N) contents were obtained using an elemental analyzer (model EA-1110, CE Instruments, Milan, Italy). The sulfur (S) contents were measured by X-ray fluorescence XRF (XGT-1700WR, Horiba, Kyoto, Japan). The vanadium (V) and nickel (Ni) contents were determined by inductively coupled plasma mass spectrometry (ICP-MS) (Agilent 7700 series, Agilent Technologies, USA). The separation of saturate, aromatic, resin and asphaltene fractions was performed as follows. AR was dissolved in n-heptane at a ratio of 1/50 (g/g) with stirring at 60 oC for 5 h, and the mixture was then filtered. The insoluble fraction was extracted with toluene in a Soxhlet apparatus and dried to obtain asphaltene. The n-heptane solution fraction (maltene) was eluted consecutively with n-heptane, toluene and toluene/methanol (9:1, v/v) at a ratio of solvent to maltene of 250 ml to 1 g with a glass column packed with activated neutral alumina to yield the 5

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

saturate, aromatic and resin fractions, respectively. Their properties and compositions are summarized in Table 1. 2.2. Electron Spin Resonance Analysis The electron spin resonance (ESR) analysis was carried out using a JESFA200 ESR spectrometer with an X-band bridge (JEOL Ltd., Tokyo, Japan) with standard 100-kHz field modulation. The sample measurement was performed at a microwave frequency of 9.0 GHz and power of 1.0 mW using a cylindrical TE cavity, where the quartz tube (5 mm o.d., 4 mm i.d.) for holding the sample was inserted along the cylindrical axis of the cavity. The magnetic field was calibrated with an ESR marker (Mn2+ powder). A temperature accessory (DVT controller, JEOL) was used to control the measurement temperature of the samples in the cavity, and liquid nitrogen gas was used to cool the sample when it was measured at a temperature below 0 oC. The V etioporphyrin (Etio) (Sigma-Aldrich, St. Louis, MO, USA) and V 5, 10, 15, 20-tetraphenylporphyrin (TPP) (Wako Chemicals, Tokyo Japan) samples were dissolved in toluene at 200 ppm. The ARs, resin and asphaltene were diluted with toluene (Wako Chemicals, Tokyo Japan) and THF (Wako Chemicals, Tokyo Japan) at concentrations from 0.1 to 80 wt % to study the V=O tumbling stage. A typical V isotropic spectrum obtained with a TPP V=O complex in toluene measured at 20 oC shows eight perpendicular lines as shown in Figure 1a, indicating that the V complex can tumble rapidly, without constraint. A typical V anisotropic spectrum obtained by TPP V=O complex in toluene measured at -120 oC consists of 16 partially overlapped components of the hyperfine structure in Figure 1b, 8 lines in parallel and 8 lines perpendicular because of I=7/2, suggesting that the V complex is strongly constrained, and in a rigid state.35-37 6

ACS Paragon Plus Environment

Page 6 of 35

Page 7 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

The B parameter, which has been reported as a sensitive indicator of the tetragonal distortion that occurs with a change in the V=O bond length and distance of four nitrogen ligands in the basal plane,38 can be derived as follows: B = ∆g‫ ׀׀‬/∆g┴ = 4  /    ∆g ‫ = ׀׀‬g ‫ ׀׀‬− g  = −8  /  ∆g ┴ = g ┴ − g  = −2  / 

where g ‫ ׀׀‬and g ┴ are the magnetic field values, and ‫ ׀׀‬and ┴ are the hyperfine coupling tensors. The g and A values are related to the electronic structure of the vanadyl species. That g ┴ > g ‫׀׀‬

in the electronic ground state suggests tetragonal

distortion. g ‫ ׀׀‬and g ┴ are used as parameters for the vanadyl complexes to clarify the vanadium electronic structure. ge is the free electron g value of 2.0023, and λ is the spin-orbit coupling constant of the free ion. The method to obtain these parameter values from an ESR analysis has been described.38,39 The parameter values of ‫ ׀׀‬, ┴, g ‫ ׀׀‬and g ┴ are obtained using anisotropic simulation software, and the deviations of ‫ ׀׀‬, ┴, g ‫ ׀׀‬, g┴ and B are 0.03, 0.01, 0.0003, 0.0001 and 0.05, respectively. 3. RESULTS 3.1. ESR Spectra of Etio and TPP V=O Complexes Dissolved in Toluene Figure 2 shows a series of ESR spectra of Etio and TPP V=O complexes dissolved in toluene measured at a temperature range from 20 oC to liquid nitrogen temperature (-196 oC). Both V=O complexes presented a typical isotropic spectrum at 20 oC, indicating that both they were in a high tumbling and free stage. Furthermore, both V=O complexes gave a typical anisotropic spectrum at -120 oC or below, suggesting that they were in a very low tumbling and constrained state under such conditions.33, 36, 37 The

7

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

spectra changed from isotropic to anisotropic through a series of transition spectra from 20 to -120 oC, indicating that the V=O tumbling rate decreased gradually, as reported in the literature.33 There is a significant change in the ESR spectra from isotropic to anisotropic when the temperature decreases from -90 to -100 oC for both complexes, because the melting point of toluene is -95 oC, at which the solidification of the solution constrains the mobility of the complexes, resulting in changes in the ESR spectrum (from isotropic to anisotropic in Figure 2). When comparing the ESR spectra of Etio and TPP V=O complexes in detail, it was further found that the change in the ESR spectra occurred at a slightly higher temperature for the TPP V=O complex than for Etio V=O complex to give similar spectrum at a temperature range from -40 to -70 oC. This may be due to the larger molecular size of the former compared with the latter, as shown in Figure 3, indicating that the ESR spectrum is sensitive to the structure and environment of the V=O complexes. Consequently, the series of transitional ESR spectra were used as references to estimate the mobility state of the V=O complexes in AR and its resin and asphaltene fractions. 3.2. ESR Spectra of V=O Complexes in LF-AR and its Resin and Asphaltene Fractions without Solvent As shown in Figure 5, the V=O complexes in LF-AR and its resin and asphaltene fractions at 20 oC presented similar anisotropic spectra with 16 lines, differing from the TPP V=O complex in toluene measured at the same temperature, indicating that these V=O complexes in LF-AR, resin and asphaltene could not freely mobile, possibly due to constraints imposed by the surrounding matrixes, as illustrated in the model of the metal complex in petroleum or asphaltene in Figure 4, which corresponds to the spectra of TPP or Etio V=O complex dissolved in toluene at -120 oC. Some solvent or higher 8

ACS Paragon Plus Environment

Page 8 of 35

Page 9 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

measurement temperature is required to moderate constraints on V=O complexes by the surrounding matrixes, and thus to liberate the V=O complexes from them to give an isotropic spectrum. 3.3. Comparing ESR Spectra of V=O Complexes in LF-AR Dissolved in THF and Toluene Figure 6 illustrates a series of ESR spectra of the V=O complexes in LF-AR dissolved in THF and toluene measured at 20 and 50 oC at concentrations from 40 to 80 wt %. The spectra profiles were modified slightly by the THF, when the concentration was 40 wt %, corresponding to the spectrum measured at -110 oC for TPP V=O dissolved in toluene. There is no obvious difference in the ESR spectra with and without solvent at a sample concentration of 80 wt %. The ESR spectra of the V=O complexes in LF-AR dissolved in THF were similar to those in toluene, indicating that the moderating effects of THF and toluene on the constraint on the V=O complexes by the surrounding molecules were similar under such conditions. ESR parameter values of the V=O complexes in LF-AR dissolved in THF and toluene measured at 20 and 50 oC were obtained and the results are summarized in Table 3. Both THF and toluene reduced the B parameter value in comparison with the values without any solvent. The B parameter value decreased at lower concentrations of AR and increasing measurement temperature, indicating that both addition of the solvents and an increase in the temperature increased the mobility of the V=O complexes in the matrix. However, no obvious difference in the mobility of the V=O complexes in LF-AR dissolved in THF and toluene was observed on the basis of the B parameter values. 3.4. ESR Spectra of V=O Complexes in LF-AR, Resin and Asphaltene Dissolved in Toluene 9

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 7 shows a series of ESR spectra of the V=O complexes in LF-AR dissolved in toluene at low concentrations, from 0.5 to 5.0 wt % measured at 20 and 100 oC, respectively. At a measurement temperature of 20 oC, the V=O complexes in LF-AR dissolved in toluene always displayed anisotropic spectra, even at a concentration as low as 2.0 wt %, indicating a slow tumbling rate (nearly rigid state), regardless of the concentrations. Such profiles were similar to that of the TPP V=O complex dissolved in toluene measured at -105 oC. The V=O complexes in LF-AR dissolved in toluene measured at 100 oC show isotropic spectra up to a concentration of 4.0 wt %, corresponding to that of the TPP V=O complex dissolved in toluene measured at 20 oC, indicating that the V=O tumbled rapidly, overcoming the constraints of the surrounding matrixes at the low concentration. High temperatures are favorable to liberating V=O complexes from the surrounding matrix. Figure 8 shows the ESR spectra of the V=O complexes in LF-AR resin dissolved in toluene at concentrations from 2 to 20 wt % measured at 20 oC and from 30 to 40 wt % measured at 100 oC. It is clear that the ESR spectra of V=O complexes in the resin dissolved in toluene at concentrations below 5 wt % at 20 oC were isotropic, whereas the spectra started to change to anisotropic spectra at a concentration above 8 wt %. The spectra at concentrations of 5, 8 and 20 wt % correspond to those of the TPP V=O complex dissolved in toluene at 20, -40 and -50 oC, respectively. The V=O complexes in LF-AR resin dissolved in toluene measured at 100 oC with concentrations up to 30 wt % still showed normal isotropic spectra. Increasing the concentration to 35 and 40 wt % reduced the V=O tumbling rate, with spectra corresponding to those of the TPP V=O complex dissolved in toluene at -20 and -50 oC, respectively. The concentration of resin in toluene influenced the V=O mobility, similar to AR. 10

ACS Paragon Plus Environment

Page 10 of 35

Page 11 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Figure 9 shows ESR spectra of the V=O complexes in LF-AR asphaltene dissolved in toluene at concentrations from 0.1 to 0.5 wt % measured at 20 and 100 oC, respectively. There was almost no ESR signal with the LF-AR asphaltene in toluene at a concentration of 0.1 wt % measured at 20 oC, possibly due to the too low V concentration to be detected. The ESR spectrum became definite as the concentration increased. The LF-AR asphaltene in toluene at 0.5 wt % showed a spectrum of V=O complexes with a very slow tumbling rate, corresponding to that at -100 oC of the TPP V=O complex dissolved in toluene. Meanwhile, the V=O complexes in the same sample also gave a spectrum showing a slow tumbling rate, even at 100 oC, corresponding to that of the TPP V=O complex dissolved in toluene at -40 oC, indicating strong constraints on the V=O complexes in LF-AR asphaltene. 3.5. ESR Spectra of V=O Complexes in KEC-AR, Its Resin and Asphaltene Dissolved in Toluene The ESR spectra of the V=O complexes in KEC-AR dissolved in toluene at the concentrations from 2 to 10 wt % measured at temperatures of 20 and 100 oC are shown in Figure 10. The V=O complexes in KEC-AR at 20 oC showed spectra with a low tumbling rate, regardless of the concentration from 2 to 10 wt %, corresponding to that of the TPP V=O complex dissolved in toluene at -60 to -90 oC. However, the KEC-AR at concentrations of 2 to 8 wt % at 100 oC showed isotropic spectra. The highest concentration of KEC-AR that gave a normal isotropic spectrum was higher than that of LF-AR at the same temperature of 100 oC, as summarized in Table 2. Figure 11 shows ESR spectra of the V=O complexes in KEC-AR resin dissolved in toluene at various concentrations measured at 20 and 100 oC, respectively. The V=O complexes in KEC-AR resin dissolved in toluene gave an isotropic spectrum at a 11

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

concentration of 5 wt % measured at 20 oC, the concentration giving an isotropic spectrum was similar to that of LF-AR resin dissolved in toluene. However, the V=O complexes in KEC-AR resin dissolved in toluene showed normal isotropic spectra at 100 oC at concentrations up to 40 wt % , obviously higher than that of LF-AR resin dissolved in toluene. The KEC-AR resin dissolved in toluene at 45 wt % started to transform the spectrum from isotropic to anisotropic, corresponding to that of the TPP V=O complex dissolved in toluene measured at -20 oC, and a concentration of 50 wt % gave a spectrum with a slow tumbling rate, corresponding to that of the TPP V=O complex dissolved in toluene measured at -50 oC. Figure 12 shows the ESR spectra of the V=O complexes in KEC-AR asphaltene dissolved in toluene measured at 20 and 100 oC at different concentrations from 0.3 to 1.0 wt %. KEC-AR asphaltene in toluene showed the ESR spectra with some differences at 20 and 100 oC. The spectra under these two temperatures at a concentration of 1.0 wt % correspond to those of the TPP V=O complex dissolved in toluene at -100 and -40 oC, respectively. Like the asphaltene from LF-AR, improvement in the mobility of the V=O complexes in asphaltene from KEC-AR by reducing the sample concentration in toluene to give an isotropic spectrum is also difficult, even at a measurement temperature of 100 oC. 4. DISCUSSION 4.1. V=O Complexes in Solvents The ESR spectra of V=O complexes in the petroleum fractions reflect strong constraints by the surrounding environment, where the rotational mobility of the V=O complexes is governed by their structure and peripheral organic components, including the solvent. For the ESR spectra of V=O complexes, the isotropic and anisotropic 12

ACS Paragon Plus Environment

Page 12 of 35

Page 13 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

natures of the spectra determined by the V=O tumbling rate reflect principally the V rotational mobility.33 Thus, the V=O complex in a rigid state shows an anisotropic spectrum, while the spectrum becomes more isotropic with increasing V=O rotational mobility, this moderates the directional anisotropy when the V=O complex is present in a liquid state at differing degrees of Brownian motion, as observed with the TPP V=O complex dissolved in toluene over a wide temperature range from -196 to 20 oC. The V=O complex in the liquid state and with low or medium molecular weight tends to give an isotropic spectrum, even at low measurement temperatures. In contrast, an V=O complex with high molecular weight, such as those with large substituents bounded to the tetradentate ligand or with the axial ligands22, 23 to the central V=O ion (e.g. V=O complexes in petroleum), shows an anisotropic spectrum, even in a solvent. Standard Etio and TPP V=O complexes dissolved in toluene show a series of spectra from free rotation to a rigid state with variable tumbling rates, where complete dissolution of the V=O complex in toluene provided an isotropic spectrum due to the free rotational mobility of the V=O complex. V=O complexes with different pendant ligands, such as Etio and TPP V=O complexes, can give different spectra when measured at -40 and -70 oC. However, the difference of the spectra of Etio and TPP V=O complexes is smaller when compared with those of V=O complexes in petroleum, so the surrounding molecules of V=O complexes in petroleum may exert more influence on V=O rotational mobility.

4.2. Effect of ARs and Their Fractions on V=O Complexes Mobility The V=O complexes in AR shows anisotropic spectra that are obviously different from those of TPP V=O complexes in toluene under the same measurement conditions. The differences can be ascribed to the greater variety in core ligands, their pendants, and surrounding molecules in the AR matrix. Regardless of the ligands, V=O complexes in petroleum appear molecularly dispersed in petroleum and in a liquid state. The 13

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

surrounding molecules of the matrix chemically or physically bound to the V=O complexes, as suggested in the Yen model,36 may govern the mobility and solubility of the V=O complexes. The different ESR spectra of TPP V=O complexes in toluene and V=O complexes in AR indicate that the V=O complexes in AR are strongly constrained by the surrounding matrix molecules via the non-covalent bonds or other interactions, resulting in their very slow rotational mobility. A recent FT-ICR-MS8 study showed several V=O porphyrins in crude oil, and different structures and derivatives of the V=O complexes are present in different crude oils. They include complexes with low or medium molecular weight, and which are smaller than TPP V=O complex, they should show an isotropic spectrum when the petroleum fraction becomes liquid or dissolved completely in a solvent, particularly at a low concentration through overcoming the constraints of the surrounding matrix. The matrixes of AR, its resin and asphaltene are known to govern the degree of mobility of V=O complexes. Thus, whether the V=O complexes in a matrix show an isotropic or anisotropic spectrum depends on the degree of liberation of V=O complexes from the matrix. When ARs or their fractions are dissolved in toluene, whether an isotropic or anisotropic spectrum is observed depends on the composition and properties of the matrixes (ARs, their resin or asphaltene), the sample concentration in the solvent and the measurement temperature. The maximum concentration in toluene to give an isotropic spectrum decreases in the order of resin >AR > asphaltene, and the asphaltene dissolved in toluene hardly gives an isotropic spectrum, even at a very low concentration and measurement temperature of 100 oC, suggesting that the constraining intensity of the surrounding molecules in asphaltene is significantly stronger than in the 14

ACS Paragon Plus Environment

Page 14 of 35

Page 15 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

AR and resin. Resin dissolved in toluene at 20 oC can give an isotropic spectrum, and higher temperatures more readily result in an isotropic spectrum. Solubility of the petroleum fraction can reflect their molecular composition. Additionally, compositional distribution should be considered in the presence or absence of toluene. Smaller molecules in the matrix may play the role of a solvent in the same fraction. KEC-AR and its resin give isotropic spectra at higher concentrations than those of LF-AR and its resin, as summarized in Table 2, probably due to weaker constraints on the V=O complexes in KEC-AR and its resin in comparison with that in LF-AR and its resin. The present findings indicate that higher concentrations of petroleum fractions dissolved in toluene tend to give anisotropic spectra due to the very slow tumbling rate of the V=O complexes. There may be a series of stages in which the V=O complexes in the matrix are liberated to varying extents. Complete dissolution of the matrix in toluene at a low concentration may allow V=O free mobility, whereas solvation of the matrix may moderate partially the constraint on the V=O complexes, which are still restricted to some extent by the surrounding matrix. The extent of this constraint defines the spectrum profile, as seen in the experimental results. The constraint extent should be influenced by the solvent to matrix ratio and the measurement temperature. It is important to further clarify the solvated state of the matrix and the V=O tumbling rate in the matrix. This discussion may suggest that V=O complexes in petroleum behave in concert with co-existing molecules in its fraction, even under conditions of a hydrotreatment reaction. The combined behavior of the V=O complexes and the matrix should be taken into consideration, particularly for coke formation and catalyst deactivation in the HDM 15

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

process. Liberation of V=O complexes from the matrix, as indicated by their molecular mobility, is important in improving metal removal. 5. CONCLUSIONS In this study, we have demonstrated that matrixes of V=O complexes in AR, resin and asphaltene restrict their rotational mobility by comparing their ESR spectra with those of the standard porphyrin complexes (Etio or TPP V=O complex) dissolved in toluene at temperatures ranging from -196 to 20 oC. V=O complexes in resin dissolved in toluene can give an isotropic spectrum at 20 oC. The highest concentration of resin dissolved in toluene that gives an isotropic spectrum at 100 oC is obviously much higher than those of the AR. The V=O complexes in both the AR and the resin dissolved in toluene show isotropic spectra at 100 oC, whereas the V=O complexes in the asphaltene dissolved in toluene hardly give an isotropic spectra. This is probably due to the chemical constraints on the V=O complexes in asphaltene by the intermolecular interactions. We conclude that the constraint on V=O complexes in different fractions increases in the order of resin < AR < asphaltene. The roles of aromatic and resin fractions should be recognized when dissolving V=O complexes in the asphaltene fraction present in AR. KEC-AR and its resin dissolved in toluene at a measurement temperature of 100 oC show isotropic spectra of V=O complexes up to concentration of 8 and 40 wt %, respectively, which are obviously higher than those of LF-AR and its resin, suggesting stronger constraints of LF-AR and its resin on V=O complexes than those in KEC-AR and its resin, indicating stronger aggregation in LF-AR.

AUTHOR INFORMATION Corresponding Author

16

ACS Paragon Plus Environment

Page 16 of 35

Page 17 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

*E-mail: [email protected]; Tel: (081) 092-662-0410.

ACKNOWLEDGEMENTS The authors acknowledge Japan Cooperation Center Petroleum (JCCP), the Kuwait Oil Company (KOC) and the Kuwait Institution for Scientific Research (KISR) for collaboration on this joint project. Acknowledgement is also extended to Kuwait National Petroleum Company (KNPC) for the in-kind contribution and technical support.

REFERENCES (1) Maity, S. K.; Perez, V. H.; Ancheyta, J.; Rana, M. S. Catalyst deactivation during hydrotreating of Maya crude in a batch reactor. Energy & Fuels. 2007, 21, 636-639. (2) Ancheyta-Juárez, J.; Maity, S. K.; Betancourt-Rivera, G. Comparison of different Ni-Mo/alumina catalysts on hydrodemetallization of Maya crude oil. Applied Catalysis A: General. 2001, 216, 195-208. (3) Rana, M. S.; Ancheyta, J.; Rayo, P.; Maity, S. K. Effect of alumina preparation on hydrodemetallization and hydrodesulfurization of Maya crude. Catalysis Today. 2004, 98, 151-160. (4) Rana, M. S.; Ancheyta, J.; Maity, S. K.; Rayo, P. Characteristics of Maya crude hydrodemetallization and hydrodesulfurization catalysts. Catalysis Today. 2005, 104, 86-93. (5) Ancheyta, J.; Betancourt, G.; Marroquin, G.; Centeno, G.; Castaneda, L. C.; Rayo, P. Hydroprocessing of Maya heavy crude oil in two reaction stages. Applied Catalysis A: General. 2002, 233, 159-170. (6) Ancheyta, J.; Centeno, G.; Trejo, F.; Marroquin, G.; Garcıa, J. A.; Tenorio, E.; Torres,

17

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

A. Extraction and characterization of asphaltenes from different crude oils and solvents. Energy & Fuels. 2002, 16, 1121-1127. (7) Tojima, M.; Suhara, S.; Imamura, M.; Furuta, A. Effect of heavy asphaltene on stability of residual oil. Catalysis Today. 1998, 43, 347-351. (8) Trejo, F.; Centeno, G.; Ancheyta, J. Precipitation, fractionation and characterization of asphaltenes from heavy and light crude oils. Fuel. 2004, 83, 2169-2175. (9) Shirokoff, J. W.; Siddiqui, M. N.; Ali, M. F. Characterization of the structure of saudi crude asphaltenes by X-ray diffraction. Energy & Fuels. 1997, 11, 561-565. (10) Ancheyta, J.; Centeno, G.; Trejo, F.; Marroquin, G. Changes in asphaltene properties during hydrotreating of heavy crudes. Energy & Fuels. 2003, 17, 1233-1238. (11) Camocho-Bragado, G. A.; Santiago, P.; Marin-Almazo, M.; Espinosa, M.; Romero, E. T. Fullerenic structures derived from oil asphaltenes. Carbon. 2002, 40, 2761-2766. (12) Trejo, F.; Ancheyta, J.; Morgan, T. J.; Herod, A. A.; Kandiyoti, R. Characterization of asphaltenes from hydrotreated products by SEC, LDMS, MALDI, NMR, and XRD. Energy & Fuels. 2007, 21, 2121-2128. (13) Trejo, F.; Ancheyta, J.; Rana, M. S. Structural characterization of asphaltenes obtained from hydroprocessed crude oils by SEM and TEM. Energy & Fuels. 2009, 23, 429-439. (14) Baker, E. W.; Yen, T. F.; Dickie, J. P.; Rhodes, R. E.; Clark, L. F. Mass spectrometry of porphyrins.Ⅱ. Characterization of petroporphyrins. J. Am. Chem. Soc. 1967, 89, 3631-3639. (15) Sugihara, J. M.; Bean, R. M. Direct determination of metalloporphyrins in Boscan crude oil. J. Chem. Eng. Data. 1962, 7, 269-271. (16) Zhao, X.; Liu, Y.; Xu, C. M.; Yan, Y. Y.; Zhang, Y. H. Separation and 18

ACS Paragon Plus Environment

Page 18 of 35

Page 19 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

characterization of vanadyl porphyrins in venezuela orinoco heavy crude oil. Energy & Fuels. 2013, 27, 2874-2882. (17) Rodgers, R. P.; Hendrickson, C. L.; Emmett, M. R.; Marshall, A. G.; Greaney, M. Molecular characterization of petroporphyrins in crude oil by electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry. Canadian Journal of Chemistry. 2001, 79, 546-551. (18) Liu, H.; Mu, J.; Wang, Z. X.; Ji, S. F.; Shi, Q. Characterization of vanadyl and nickel porphyrins enriched from heavy residues by positive-ion electrospray ionization FT-ICR mass spectrometry. Energy & Fuels. 2015, 29, 4803-4813. (19) McKenna, A.M.; Purcell, J. M.; Rodgers, R. P.; Marshall, A. G. Identification of vanadyl porphyrins in a heavy crude oil and raw asphaltene by atmospheric pressure photoionization fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry. Energy & Fuels. 2009, 23, 2122-2128. (20) Qian, K. N.; Edwards, K. E.; Mennito, A. S.; Walters, C. C.; Kushnerick, J. D. Enrichment, resolution, and identification of nickel porphyrins in petroleum asphaltene by cyclograph separation and atmospheric pressure photoionization fourier transform ion cyclotron resonance mass spectrometry. Anal. Chem. 2010, 82, 413-419. (21) Kim, T.; Ryu, J.; Kim, M. J.; Shul, Y. G.; Jeon, Y.; Park, Jl. Characterization and analysis of vanadium and nickel species in atmospheric residues. Fuel. 2014, 117, 783-791. (22) Bencosme, C. S.; Romero, C.; Simoni, S. Axial interaction of vanadyl tetraphenylporphyrin with Lewis bases. Inorg. Chem. 1985, 24, 1603-1604. (23) Dechaine, G. P.; Gray, M. R. Chemistry and association of vanadium compounds in heavy oil and bitumen, and implications for their selective removal. Energy & Fuels. 19

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

2010, 24, 2795-2808. (24) Quirke, J. M. E.; Eglinton, G.; Palmer, S. E.; Baker, E. W. High-performance liquid chromatographic and mass spectrometric analyses of porphyrins from deep-sea sediments. Chem. Geol. 1982, 35, 69-85. (25) Fish, R. H.; Komlenic, J. J.; Wines, B. K. Characterization and comparison of vanadyl and nickel compounds in heavy crude petroleums and asphaltenes by reverse-phase and size-exclusion liquid chromatography/graphite furnace atomic absorption spectrometry. Anal. Chem. 1984, 56, 2452-2460. (26) Duyck, C. B.; Saint'Pierre, T. D.; Miekeley, N.; Cristina, T. High performance liquid chromatography hyphenated to inductively coupled plasma mass spectrometry for V and Ni quantification as tetrapyrroles. Spectrochimica Acta Part B: Atomic Spectroscopy. 2011, 66, 362-367. (27) Desprez, A.; Bouyssiere, B.; Arnaudguilhem, C.; Krier, G.; Vernex-Loset, L. Study of the size distribution of sulfur, vanadium, and nickel compounds in four crude oils and their distillation cuts by gel permeation chromatography inductively coupled plasma high-resolution mass spectrometry. Energy & Fuels. 2014, 28, 3730-3737. (28) Acevedo, S.; Guzman, K.; Labrador, H.; Carrier, H.; Bouyssiere, B. Trapping of metallic porphyrins by asphaltene aggregates: A size exclusion microchromatography with high-resolution inductively coupled plasma mass spectrometric detection study. Energy & Fuels. 2012, 26, 4968-4977. (29) Barbier, J.; Marques, J.; Caumette, G.; Merdrignac, I.; Bouyssiere, B. Monitoring the behaviour and fate of nickel and vanadium complexes during vacuum residue hydrotreatment and fraction separation. Fuel Processing Technology. 2014, 119, 185-189. 20

ACS Paragon Plus Environment

Page 20 of 35

Page 21 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

(30) Barbier, J.; Lienemann, C. P.; Masle, A. L.; Chatron-Michaud, P.; Guichard, B. New insights into residue desulfurization processes: molecular size dependence of catalytic performances quantified by size exclusion chromatography-ICP/MS. Energy & Fuels. 2013, 27, 6567-6574. (31) Dickson, F. E.; Kunesh, C. J.; McGinnis, E. L.; Petrakis, L. Use of electron spin resonance to characterize the vanadium(1V)-sulfur species in petroleum. Analytical chemistry. 1972, 44, 978-981. (32) Premovic, P.I.; Allard, T.; Nikolic, N.D.; Tonsa, I.R.; Pavlovic, M.S. Estimation of vanadyl porphyrin concentration in sedimentary kerogens and asphaltenes. Fuel. 2000, 79, 813-819. (33) Campbell, R. F.; Freed, J. H. Slow-motional ESR spectra for vanadyl complexes and their model dependence. J. Phys. Chem. 1980, 84, 2668-2680. (34) Wong, G. K.; Yen, T. F. An electron spin resonance probe method for the understanding of petroleum asphaltene macrostructure. Journal of Petroleum Science and Engineering. 2000, 28, 55-64. (35) Biktagirov, T. B.; Gafurov, M. R.; Volodin, M. A.; Mamin, G.V.; Orlinskii, S. B. Electron paramagnetic resonance study of rotational mobility of vanadyl porphyrin complexes in crude oil asphaltenes: probing the effect of thermal treatment of heavy oils. Energy & Fuels. 2014, 28, 6683-6687. (36) Tynan, E. C.; Yen, T. F. Association of vanadium chelates in petroleum asphaltenes as studied by ESR. Fuel. 1969, 43, 191-208. (37) Yamada, Y.; Ouchi, K.; Sanada, Y.; Sohma, J. Study of carbonaceous mesophase through the e.s.r. spectra of vanadyl chelates. Fuel. 1978, 57, 79-84. (38) Espinosa, P. M.; Campero, A.; Salcedo, R. Electron spin resonance and electronic 21

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

structure of vanadyl-porphyrin in heavy crude oils. Inorg. Chem. 2001, 40, 4543-4549. (39) Malhotra, V. M.; Buckmaster, H. A. 34 GHz e.p.r. study of vanadyl complexes in various asphaltenes. Fuel. 1985, 64, 335-341.

22

ACS Paragon Plus Environment

Page 22 of 35

Page 23 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

TABLES Table 1 Properties of LF-AR and KEC-AR. KEC-AR

Boiling point ( C)

LF-AR >360

Density(g/ml)

1.0081

0.9745

C (wt %)

82.57

83.84

H (wt %)

10.06

10.95

S (wt %)

3.44

3.19

N (wt %)

0.33

0.29

V (ppm)

152.33

70.57

Ni (ppm)

23.79

13.83

Saturate (wt %)

16.66

25.69

Aromatics (wt %)

51.94

48.75

Resin (wt %)

19.46

18.38

Asphaltene (wt %)

11.94

7.18

o

>360

Table 2 Highest concentrations of AR and resin dissolved in toluene giving isotropic spectrum at 20 and 100 oC. Temp. 20oC o

100 C

Lower Fare

KEC

AR

Resin

Asphaltene

AR

Resin

Asphaltene

---

5%

---

---

5%

---

4%

30%

---

8%

40%

---

23

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 35

Table 3 ESR parameters of V=O complexes in LF-AR dissolved in THF and toluene measured at 20 and 50 oC. Samples o o

20 C 50 C o

20 C Toluene 50 oC

o

20 C THF 50 oC

A// mT

A┴ mT

g//

g┴

B

LF-AR

17.14

6.00

1.9615

1.9943

5.10

LF-AR

17.12

5.99

1.9615

1.9943

5.10

LF-AR-T 40 %

16.95

6.02

1.9621

1.9931

4.37

LF-AR-T 60 %

17.02

6.01

1.9620

1.9936

4.63

LF-AR-T 80 %

17.07

5.98

1.9618

1.9941

4.94

LF-AR-T 40 %

16.88

6.05

1.9622

1.9924

4.05

LF-AR-T 60 %

17.01

6.05

1.9618

1.9931

4.40

LF-AR-T 80 %

17.06

6.00

1.9616

1.9937

4.73

LF-AR-THF 40 %

17.04

6.05

1.9622

1.9932

4.41

LF-AR-THF 60 %

17.05

6.03

1.9622

1.9937

4.66

LF-AR-THF 80 %

17.14

6.01

1.9616

1.9941

4.96

LF-AR-THF 40 %

16.92

6.09

1.9624

1.9925

4.07

LF-AR-THF 60 %

17.01

6.05

1.9621

1.9931

4.37

LF-AR-THF 80 %

17.10

6.01

1.9620

1.9937

4.69

24

ACS Paragon Plus Environment

Page 25 of 35

FIGURES

200

Intensity (a.u.)

Intensity (a.u.)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

250

300

350

400

450

200

250

mT

(a)

300

mT

350

400

(b)

Figure 1. ESR spectra of V=O complex (a) isotropic spectrum and (b) anisotropic spectrum

25

ACS Paragon Plus Environment

450

Energy & Fuels

o

-60 C o

Intensity (a.u.)

20 C

o

-70 C o

o

-80 C

o

-90 C

o

-100 C

-20 C

o

-30 C

o

-40 C o

-50 C Etio

200

250

300

350

400

450 200

250

300

mT

mT

350

400

450

o

-105 C o

-110 C o

-120 C o

-140 C o

-196 C

200

250

300

mT

350

400

450

o

o

-60 C o -70 C

o

-80 C

o

-90 C

20 C Intensity (a.u.)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 35

o

-20 C

o

-30 C

o

-100 C

o

-40 C o

-50 C TPP

200

250

300

mT

350

400

450 200

250

300

350

400

450

mT

o

-105 C o -110 C o -120 C o -140 C o

-196 C

200

250

300

350

400

450

mT

Figure 2. ESR spectra of V=O complexes in Etio and TPP porphyrins dissolved in toluene.

26

ACS Paragon Plus Environment

Page 27 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Me

Et

N

N

Et

Me

O

O V

N

V

N

N

N

Et

Me

N

N

Et

Me

(a) Vanadyl Etioporphyrin (Etio V=O)

(b) Vanadyl Tetraphenylporphyrin (TPP V=O)

Figure 3. Structure of the Etio and TPP porphyrins23 M M

(a)

(b)

O

S

O

NH

H N H

N

Aromatic carbon (c) M

Metal chelate

Saturated carbon

Aromatic carbon

Figure 4. Model of V complexes in petroleum or asphaltene (a) associate at the edge, (b) intercalated and (c) at the centre among a group of micelles36

27

ACS Paragon Plus Environment

Energy & Fuels

LF-AR LF-Resin LF-Asphaltene Intensity(a.u.)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 35

200

250

300

mT

350

400

450

Figure 5. ESR spectra of V=O complexes in LF-AR, its resin and asphaltene measured at 20 oC.

28

ACS Paragon Plus Environment

Page 29 of 35

o

o

LF-AR-THF 40% 50 C o LF-AR-THF 60% 50 C o LF-AR-THF 80% 50 C

Intensity(a.u.)

LF-AR-THF 40% 20 C o LF-AR-THF 60% 20 C o LF-AR-THF 80% 20 C

200

250

300

350

400

450

mT

(a)

o

LF-AR-T 40% 20 C o LF-AR-T 60% 20 C o LF-AR-T 80% 20 C

o

LF-AR-T 40% 50 C o LF-AR-T 60% 50 C o LF-AR-T 80% 50 C

Intensity(a.u.)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

200

250

300

350

400

450

mT

(b) Figure 6. ESR spectra of V=O complexes in LF-AR dissolved in (a) THF and (b) toluene

29

ACS Paragon Plus Environment

Energy & Fuels

LF-AR-T 2.0 % LF-AR-T 5.0 %

Intensity(a.u.)

LF-AR-T 0.5 % LF-AR-T 4.0 %

200

250

300

mT

350

400

450

(a)

LF-AR-T 0.5 % LF-AR-T 4.0 %

LF-AR-T 2.0 % LF-AR-T 5.0 %

Intensity(a.u.)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 35

200

250

300

350 mT

400

450

(b) Figure 7. ESR spectra of V=O complexes in LF-AR dissolved in toluene measured at (a) 20 oC and (b) 100 oC.

30

ACS Paragon Plus Environment

Page 31 of 35

LF-R-T 5.0 % LF-R-T 20.0 %

Intensity(a.u.)

LF-R-T 2.0 % LF-R-T 8.0 %

200

250

300

mT

350

400

450

(a)

LF-R-T 30.0 % LF-R-T 40.0 %

LF-R-T 35.0 %

Intensity(a.u.)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

200

250

300

mT

350

400

450

(b) Figure 8. ESR spectra of V=O complexes in LF-AR resin dissolved in toluene measured at (a) 20 oC and (b) 100 oC.

31

ACS Paragon Plus Environment

Energy & Fuels

LF-As-T 0.3 %

Intensity(a.u.)

LF-As-T 0.1 % LF-As-T 0.5 %

200

250

300

mT

350

400

450

(a)

LF-As-T 0.1% LF-As-T 0.5%

LF-As-T 0.3%

Intensity(a.u.)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 35

200

250

300

mT

350

400

450

(b) Figure 9. ESR spectra of V=O complexes in LF-AR asphaltene dissolved in toluene measured at (a) 20 oC and (b) 100 oC.

32

ACS Paragon Plus Environment

Page 33 of 35

KEC-AR-T 2.0 % KEC-AR-T 8.0 % Intensity(a.u.)

KEC-AR-T 5.0 % KEC-AR-T 10.0 %

200

250

300

mT

350

400

450

(a)

KEC-AR-T 2.0 %

KEC-AR-T 5.0 %

KEC-AR-T 8.0 %

KEC-AR-T 10.0 %

Intensity(a.u.)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

200

250

300

mT

350

400

450

(b) Figure 10. ESR spectra of V=O complexes in KEC-AR dissolved in toluene measured at (a) 20 oC and (b) 100 oC.

33

ACS Paragon Plus Environment

Energy & Fuels

KEC-R-T 8.0 %

Intensity(a.u.)

KEC-R-T 5.0 % KEC-R-T 10.0 %

200

250

300

mT

350

400

450

(a)

KEC-R-T 30.0 % KEC-R-T 45.0 %

KEC-R-T 40.0 % KEC-R-T 50.0 %

Intensity(a.u.)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 35

200

250

300

mT

350

400

450

(b) Figure 11. ESR spectra of V=O complexes in KEC-AR resin dissolved in toluene measured at (a) 20 oC and (b) 100 oC.

34

ACS Paragon Plus Environment

Page 35 of 35

KEC-As-T 0.5 %

Intensity(a.u.)

KEC-As-T 0.3 % KEC-As-T 1.0 %

200

250

300

mT

350

400

450

(a)

KEC-As-T 0.3% KEC-As-T 1.0%

KEC-As-T 0.5%

Intensity(a.u.)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

200

250

300

mT

350

400

450

(b) Figure 12. ESR spectra of V=O complexes in KEC-AR asphaltene dissolved in toluene measured at (a) 20 oC and (b) 100 oC.

35

ACS Paragon Plus Environment