Enhanced Oil Spill Remediation Using Environmentally Asymmetric Dicationic Ionic Liquids: Synthesis, Characterization, and Evaluation

: The disastrous consequences for society—economically, environmentally, and socially— caused by oil spills encouraged us to treat this problem. The target of this work is to synthesize new amphiphilic dicationic ionic liquids (Ia, Ib, and Ic) and evaluate them spectroscopically and gravimetrically as potential oil spill dispersants at different temperatures to cover cold and warm areas. The synthesized ILs were well characterized by different tools for analysis of their surface activity and thermal stability. Ia, Ib, and Ic showed good dispersion effects, which were recorded to be 5.32, 20.45, and 33.61% for Ia, Ib, and Ic, respectively, at 10 ◦ C and 12.28, 52.55, and 66.80% for Ia, Ib, and Ic, respectively, at 30 ◦ C with a dispersant-to-oil ratio (DOR) of 0.8:10 (wt.%). Acute toxicity tests were elucidated against Nile tilapia and Oreochromis niloticus ﬁsh and conﬁrmed their slight toxicity by determining a LC50 value greater than 100 ppm after 96 h, which recorded 13.25, 17.75, and 37.5 mg/L for Ia, Ib, and Ic, respectively. Overall, the new synthesized ILs can be represented as sustainable materials for toxic chemicals to disperse oil spills.


Introduction
Oil spills are a type of pollution that is released from petroleum hydrocarbons into the surrounding environment [1], especially the marine eco-system, as a result of enormous human activities [2].Oil spills can cause severe negative effects economically [3,4], environmentally, and socially [5,6].Consequently, accidents spark a lot of political and media interest [7], which makes the government and the political community struggle together to respond to the oil spill and take the appropriate action to protect the environment from serious disaster [8,9].
The current challenge is to clean up oil spills caused by accidents.Accidents involving oil spills will continue to occur because society relies heavily on various petroleum products [10].Oil spill treatment methods can be classified as mechanical, biological, or chemical.In mechanical treatment, skimmers are fixed or mobile devices used for the removal of floating and/or emulsified oil from the surface of the water [11].Biological treatment includes the addition of microbes, nutrients, and/or oxygen and the subsequent biodegradation of spilled oil [12].A diverse community of microorganisms is being formed that can use them as their sole carbon source to degrade petroleum hydrocarbons.Chemical treatments include dispersants sprayed on oil spills to break them up into small droplets.
The oil dispersants split the oil slicks by minimizing the oil-water interface and accelerating the breakdown of the hydrophobic oil portion.
The use of IL-based surfactants for oil spill remediation is a relatively novel strategy [13,14].Due to their prospective environmental and technological advantages of low volatility, low toxicity, high surface activity [15], and thermal stability [16], ILs have garnered a lot of attention [17].A new class of dicationic ionic liquids (DILs) has been introduced.DILs are more adaptable than monocationic ionic liquids (MILs).That is referring to the enormous possible changes in the molecular structure by changing the cation or anion type; furthermore, the side chain linker between the two head groups could be a rigid or flexible spacer [16].
Usually, DILs exhibit greater thermal stability compared with MILs.They display decomposition temperatures between 330 and 400 • C and 145 and 185 • C for DILs and MILs, respectively.DILs are more viscous than MILs [17], so they have a wide range of applications, including high-temperature lubrication, as catalysts in esterification reactions, asphaltene, and oil spill dispersion [14].
The FT-IR bands of the prepared compounds were examined using a Nicolet Ia-10 at scan resolutions of 4000-400 cm −1 and 4 cm and a scan rate of 32 cm/min, in that order. 1 H-NMR spectra were screened using BURKER 1 H-NMR spectroscopy in DMSO-d6 solvent.(400.19MHz and a 5-mm broad-band inverse Z gradient probe).Elemental analysis was achieved via an elemental analyzer.Thermogravimetric analysis was studied using a thermal analyzer at a heating rate of 10 • C/min.Samples are heated from ambient temperature up to 600 • C under nitrogen flow.The thermal degradation discussion was revealed at the point of 95% weight loss from the original weight.Surface tension parameters were obtained by a tensiometer (Du Nouty) with a platinum ring.Fresh aqueous solutions of the examined ILs (as surfactants) were prepared and measured within 0.01-0.00001M/L at the 25 • C concentration range.Compound I was previously prepared in our previous work [20,21] by stirring 0.1 mL of 2-methyl imidazole with potassium hydroxide dissolved in acetonitrile.Following the addition of 0.1 of bromododecane, which was completely miscible, drop by drop, the stirring stopped when the white precipitate was noticed (after about 3 h).By filtration, the white precipitate of KBr was eliminated, and the filtrate was vaporized in an oven under vacuum.ii.
The completion of the reaction was monitored and confirmed by the TLC technique in Scheme 1.
h).By filtration, the white precipitate of KBr was eliminated, and the filtrate was vaporized in an oven under vacuum.ii.
[22].Ia, Ib, and Ic (dicationic ILs) were synthesized by the compounding of compound I (1:1) with a, b, and c under refluxing for 3 h.The products were collected and dried via evaporation under vacuum (89.2, 87.5, and 89.4% yield).The completion of the reaction was monitored and confirmed by the TLC technique in Scheme 1.

Quantum Chemical Parameters
Synthesized ionic liquids (Ia, Ib, and Ic) differ from each other by the type of cation.Quantum chemical calculations and the geometric shape of prepared molecules were performed using density functional theory (DFT).The quantum chemical calculations were performed by utilizing the basis set 6-31G'(d,p), where Beck's three-parameter exchange functional was used along with the Lee-Yang-Parr nonlocal correlation function (B3LYP).The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) were studied for all synthesized compounds.The calculated quantum parameters, including energy gap (ΔE), chemical hardness (η), chemical softness (σ), dipole moment (µ), and electron affinity (A), were recorded.Chemical hardness (η) is defined as the resistance of an atom to the transfer of charge or change in its electron configuration and is estimated according to Equation (1).
The softness (σ), which is described as the capacity of the atom or atom group to receive electrons, can be estimated according to Equation (2).
The ionization potential (I) can be derived from EHOMO, while the electron affinity (A) related to ELUMO.

Quantum Chemical Parameters
Synthesized ionic liquids (Ia, Ib, and Ic) differ from each other by the type of cation.Quantum chemical calculations and the geometric shape of prepared molecules were performed using density functional theory (DFT).The quantum chemical calculations were performed by utilizing the basis set 6-31G'(d,p), where Beck's three-parameter exchange functional was used along with the Lee-Yang-Parr nonlocal correlation function (B3LYP).The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) were studied for all synthesized compounds.The calculated quantum parameters, including energy gap (∆E), chemical hardness (η), chemical softness (σ), dipole moment (µ), and electron affinity (A), were recorded.Chemical hardness (η) is defined as the resistance of an atom to the transfer of charge or change in its electron configuration and is estimated according to Equation (1).
The softness (σ), which is described as the capacity of the atom or atom group to receive electrons, can be estimated according to Equation (2).
The ionization potential (I) can be derived from E HOMO, while the electron affinity (A) related to E LUMO.

Evaluation ot the Prepared ILs for Oil Spill Remidiation Efficiency Test
The competence of Ia, Ib, and Ic as oil dispersants at different temperatures and concentrations was illustrated using the same method as discussed in our previous work [21].Figure 1 shows how 250 mL of seawater was poured into a separate funnel that was kept at a constant temperature, and then 5 mL of crude oil was added to the water surface for 1 min.The amounts of different concentrations of the studied dispersants were determined, the water with oil and dispersant was shaken, 50 mL of oily water emulsion was taken in a measuring cylinder, and finally the oil was extracted using chloroform.The extracted oil in chloroform was transferred into a 100 mL volumetric flask, completed with chloroform to the mark, and mixed well.By using UV spectroscopy, determine the absorption of the extract against the blank.The efficiency index was calculated according to the following Equation ( 5

Acute Toxicity Test
The compounds (Ia, Ib, and Ic) were tested against healthy Nile tilapia, Oreochromis niloticus, with an average body weight of 40.0 ± 2.0 g.Each treatment had three replicates for each compound.The lab conditions have to stay constant.Abnormal behavior and dead fish were noticed and recorded every day for four days.Then the value of the 96 h LC50 was calculated [24].

Characterization of the Synthesized ILs
FT-IR spectrum Figure 2 revealed that N-H broad stretching peaks appeared around 3416-3439 cm −1 .The stretching peaks in the range 3017-3050 cm −1 were resorted to by C-H aromatic stretching, which belongs to the vibrational motion in the imidazole ring [25], while the aliphatic stretching asymmetric C-H appeared between 2853 and 2922 cm −1 .C-C vibrations related to the imidazole ring were recorded between 1631 and 1642 cm −1 .The peaks between 1375 and 1380 cm −1 resorted to the C-N bending symmetrical mode of vibrations in imidazole rings.The vibrational bending (in plane) peaks of imidazole appeared around 1152-1183 cm −1 .While the values 725-773 cm −1 are attributed to the out-ofplane bending vibrations.The treated crude oil was studied by SARA analysis (saturates, aromatics, resin, and asphaltene), where asphaltene content must be determined and isolated.The asphaltene utilized in our study was extracted from West Bakr Petroleum Company's heavy crude oil according to IP-143 (the standard method).A specific amount (10 g) of the crude oil was refluxed for 1 h with 300 mL of n-heptane.Allow the mixture to cool for 2 h in a dark place before filtering.Inorganic materials, waxy substances, and asphaltene were still on the filter paper, which was folded and placed in soxhelt to be washed with heated n-heptane to get rid of any waxy substances.The asphaltene was dissolved in hot toluene.Finally, evaporate the toluene and calculate the asphaltene percent according to the following Equation ( 6): where M and G refer to the mass of asphaltene and crude oil in grams, respectively.Maltene, or deasphalted oil, is then separated into saturates, aromatics, and resins according to ASTM D 2007.Saturates and aromatics were analyzed by GC and HPLC spectroscopy.The saturates were analyzed by a gas chromatography mass spectrometer (GC-MS, QP 2010) using a DB-5 capillary column (containing 5% phenyl and 95% polydimethyl siloxane) and measured at 30 m × 0.32 mm × 0.23 µm film thickness.The MS scan parameters were an m/z range of 40-500, a scan interval of 0.5 s, a scan speed of 1000 amu s −1 and a detector voltage of 1.1 kV.It is the used carrier gas, and with a constant linear velocity of 1.1 mL min −1 and a programmed temperature (of 40-280 • C at 2 • C min −1 ) the samples were held for 30 min at 280 • C.
Aromatics were analyzed to detect PAHs by HPLC (Agilent Technologies 1200 Series HPLC) with a UV detector using a ZORBAX 5 µm, 4.6 × 250 mm NH 2 -column, which separates the aromatic fractions into sub-fractions.N-heptane was used as the mobile Separations 2023, 10, 397 6 of 16 phase in a column with a flow rate of 3 mL/min.The delay and elution time were detected, and the aromatic sub-fractions were noticed in the HPLC figures.All the HPLC systems were kept at room temperature, and the data handling was justified by Agilent software (Santa Clara, CA, USA) [23].

Acute Toxicity Test
The compounds (Ia, Ib, and Ic) were tested against healthy Nile tilapia, Oreochromis niloticus, with an average body weight of 40.0 ± 2.0 g.Each treatment had three replicates for each compound.The lab conditions have to stay constant.Abnormal behavior and dead fish were noticed and recorded every day for four days.Then the value of the 96 h LC50 was calculated [24].

Characterization of the Synthesized ILs
FT-IR spectrum Figure 2 revealed that N-H broad stretching peaks appeared around 3416-3439 cm −1 .The stretching peaks in the range 3017-3050 cm −1 were resorted to by C-H aromatic stretching, which belongs to the vibrational motion in the imidazole ring [25], while the aliphatic stretching asymmetric C-H appeared between 2853 and 2922 cm −1 .C-C vibrations related to the imidazole ring were recorded between 1631 and 1642 cm −1 .The peaks between 1375 and 1380 cm −1 resorted to the C-N bending symmetrical mode of vibrations in imidazole rings.The vibrational bending (in plane) peaks of imidazole appeared around 1152-1183 cm −1 .While the values 725-773 cm −1 are attributed to the out-of-plane bending vibrations.

Acute Toxicity Test
The compounds (Ia, Ib, and Ic) were tested against healthy Nile tilapia, Oreochromis niloticus, with an average body weight of 40.0 ± 2.0 g.Each treatment had three replicates for each compound.The lab conditions have to stay constant.Abnormal behavior and dead fish were noticed and recorded every day for four days.Then the value of the 96 h LC50 was calculated [24].

Characterization of the Synthesized ILs
FT-IR spectrum Figure 2 revealed that N-H broad stretching peaks appeared around 3416-3439 cm −1 .The stretching peaks in the range 3017-3050 cm −1 were resorted to by C-H aromatic stretching, which belongs to the vibrational motion in the imidazole ring [25], while the aliphatic stretching asymmetric C-H appeared between 2853 and 2922 cm −1 .C-C vibrations related to the imidazole ring were recorded between 1631 and 1642 cm −1 .The peaks between 1375 and 1380 cm −1 resorted to the C-N bending symmetrical mode of vibrations in imidazole rings.The vibrational bending (in plane) peaks of imidazole appeared around 1152-1183 cm −1 .While the values 725-773 cm −1 are attributed to the out-ofplane bending vibrations.Elemental Analysis: Table 7 indicated that the experimental values of the prepared ILs are in acceptable agreement with the theoretically calculated values.Thermo-gravimetric Analysis (TGA) and DTG composite pictures of Ia, Ib, and Ic recorded high onset temperatures, which related to the resistance of these materials toward thermal degradation.Figure 6 confirmed that Ia, Ib, and Ic are thermally stable, and the actual first steps of decomposition were recorded as 247, 223, and 239 • C, respectively, while the final decomposition was up to 600 • C, which presented endothermal phenomena until reaching the maximum mass loss percent (100% wt.loss of the sample) [26].The position of the methyl group affects the thermal stability; in the meta position (Ic), the compound showed higher thermal stability than in the para position (Ib) because of the effect of mesomeric resonance, which prevents C-bond leakage [27].Thermo-gravimetric (TGA) and DTG composite pictures of Ia, Ib, and Ic recorded high onset temperatures, which related to the resistance of these materials toward thermal degradation.Figure 6 confirmed that Ia, Ib, and Ic are thermally stable, and the actual first steps of decomposition were recorded as 247, 223, and 239 °C, respectively, while the final decomposition was up to 600 °C, which presented endothermal phenomena until reaching the maximum mass loss percent (100% wt.loss of the sample) [26].The position of the methyl group affects the thermal stability; in the meta position (Ic), the compound showed higher thermal stability than in the para position (Ib) because of the effect of mesomeric resonance, which prevents C-bond leakage [27].

Surface Tension and Critical Micelle Concentration (CMC)
The critical micelle concentrations (CMC) of Ia, Ib, and Ic were determined at 25 °C, as shown in Figure 7.The synthesized compounds' CMC has a low value and significant surface activity, as can be observed.Table 8 displays their CMC and surface tensions at CMCs.The adsorption of the synthesized compounds was studied using surface tension measurements.These results demonstrated a decrease in surface tension, indicating a high (Ia) (Ib) (Ic)

Surface and Physical Characteristics of ILs 3.2.1. Surface Tension and Critical Micelle Concentration (CMC)
The critical micelle concentrations (CMC) of Ia, Ib, and Ic were determined at 25 • C, as shown in Figure 7.The synthesized compounds' CMC has a low value and significant surface activity, as can be observed.Table 8 displays their CMC and surface tensions at CMCs.The adsorption of the synthesized compounds was studied using surface tension measurements.These results demonstrated a decrease in surface tension, indicating a high effect on surface activity.The results revealed that Ia, Ib, and Ic have significant power to lower aqueous system surface tension [21,28].With the addition of the surfactant, the surface tension below the CMC was reduced.The surface tension was in the order Ic > Ia > Ib, while the CMC was in the order Ia > Ib > Ic.CMC is useful since it displays surfactant assembly proclivity in water.The micellization is driven by hydrophobic forces opposing electrostatic repulsion among the ionic head groups at the micelle surface.Micelles are known to be disorganized assemblies with a mobile hydrocarbon chain within.The Ic has a methyl group at position 3 of the pyridine ring, making it the most reactive.Micelles are generated by zwitter-ionic species with a characteristic nucleophilic group covalently connected to their portions [30].It is well known that the intermolecular interactions of monomeric surfactants and their interactions with solvents primarily regulate the aggregation behavior of typical monocationic surfactants in solution.In contrast to conventional monocationic surfactants, dicationic surfactants' aggregation behavior in solution is governed by the cooperativeness of both intra-and inter-molecular contacts, in addition to their interactions with solvents.
The dicationic IL molecules must arrange the conformation in an aggregate state when the conformation is found in the free state.In this case, the two alkyl chains are as close as possible, and the spacer is gradually incorporated into the aggregate hydrophobic core as the spacer's hydrophobicity increases.
Additionally, in the aqueous solutions of dicationic ILs, the spacer causes a strong inclination to form a wide range of aggregate forms [21], such as worm-like micelles and tubule aggregates.An acceptable explanation for the development of worm-like micelles could be found in the surfactant packing parameter theory (SPP).By comparing monocationic with dicationic ILs as surfactants, it was found that they typically have higher SPP values, which causes the spontaneous formation of rod-or worm-shaped micellar aggregates even just above the CMC [29].
With the addition of the surfactant, the surface tension below the CMC was reduced.The surface tension was in the order Ic > Ia > Ib, while the CMC was in the order Ia > Ib > Ic.CMC is useful since it displays surfactant assembly proclivity in water.The micellization is driven by hydrophobic forces opposing electrostatic repulsion among the ionic head groups at the micelle surface.Micelles are known to be disorganized assemblies with a mobile hydrocarbon chain within.The Ic has a methyl group at position 3 of the pyridine ring, making it the most reactive.Micelles are generated by zwitter-ionic species with a characteristic nucleophilic group covalently connected to their portions [30].
The working action mechanism of cationic surfactants mainly depends on the adsorption of hydrophilic groups, which are represented by positively charged nitrogen atoms in Ia, Ib, and Ic, onto a polar hydrophobic group on a nonpolar phase.At lower concentrations, a higher Pc20 represents a 20 mN/m decrease in surface tension.The low value of πCMC is directly related to the low surface activity.As shown in Table 8, πCMC varies between 40 and 38 dyne/cm.According to Gibb's adsorption isotherm, the relation between the surface tension and logarithmic concentration below the CMC region was used to calculate the surfactant molecules that adsorbed at the surface (max).
The IL solution concentration affects the number of ionic species (n).If the surfactant is present near or on the surface, it will decrease the surface energy.By increasing the ILs' activity, the surface tension will decrease, and the maximum will be positive, as shown in Table 8.Amin was recorded for Ic and could be calculated according to the following Equation ( 8): A min = 10 16 /N.Γ max (8) where N, represent an Avogadro number.

Thermodynamics of Micellization and Adsorption
∆G o mic provides details about the aggregation process nature gives information about the nature of the aggregation process.The following equation illustrates the gemini surfactant's energy as shown: Negative values of ∆G o mic indicate that the process of micellization is spontaneous.Thermodynamic quantities relevant to micellization are the changes of free energy.
At ambient temperature and atmospheric pressure, the values of ∆G o mic are always negative.The thermodynamic behavior of micellization was found, as Ic is less negative than Ib and Ia.The free energy of adsorption at the air-liquid interface (a measure of the free energy of transfer per mole of surfactant at a unit concentration from bulk to surface at unit pressure) was calculated using Equation ( 10): A negative value of ∆G o ads indicates that the adsorption process is spontaneous.Furthermore, the higher values of ∆G o ads compared to ∆G o mic , Table 8 show more spontaneous adsorption at the solution-to-air interface than micelization does in the bulk.

Quantum Chemical Calculations
The optimized geometry and the frontier molecular orbital density of the investigated compounds are shown in S1.The chemical reactivity of the prepared compounds was explained by the value of ∆E.Molecules that have a small (∆E) are less stable and more reactive.(σ) and (η) are properties that are related to the reactivity of the compounds.Small (σ) values correspond to soft molecules that can easily donate electrons to the acceptor system.
Consequently, adsorption on asphaltene surfaces could take place easily for the molecule that has the highest (σ) value, postponing the asphaltene dispersion.(µ) is another index factor that can illustrate the effect of prepared compounds on asphaltenes, which can be defined as the measurements of the bond polarity in molecule S2.According to (σ) and (η) values, the studied compounds were considered soft molecules.
The results in Table 9 indicate that Ia, Ib, a Ic have an appreciable (µ) that facilitates the interaction with asphaltene [31].Data showed that Ic is the most effective dispersant that is compatible with (A) (1.119 eV) due to its high electron density and high aromaticity, which accelerate the penetration of crude oil and accelerate the adsorption on asphaltene surfaces and their aggregation.

Evaluation of the Synthesized ILs as Oil Dispersants:
The evaluation of the synthesized compounds as dispersants was discussed using different ratios of dispersants to crude oil, wt.% (DOR), at specific temperatures (10 and 30 • C).

Effect of Temperatures
It was noticed that the dispersion efficiency of three newly synthesized compounds increased with increasing temperature for the same dispersant.In the case of cold sea water, Table 10 (at 10 • C) shows that Ia, Ib, and Ic recorded 5.32, 20.45, and 33.61% as dispersion efficiency, respectively.Whereas, the dispersion efficiency increased in warm water (30 • C), which recorded 12.28, 52.55, and 66.80% for Ia, Ib, and Ic in order.Whereas, the dispersion efficiency increased in warm water (30 • C), which recorded 12.28, 52.55, and 66.80% for Ia, Ib, and Ic in order in Table 11.  10 and 11 indicated that the efficiency of the synthesized dispersants increased with increasing the DOR, wt.%, from (0.4:10) to (0.8:10), started slightly by using (1.2:10), and decreased again when (1.6:10) was used.So, the optimum ratio in this study was 0.8:10.

Effect of Cation Structure
Three derivatives of pyridine were used: pyridine, 3-methylpyridine, and 4-methylpyridine, and the results showed that the dispersion efficiency ranked as Ia < Ib < Ic which refers to the presence of CH 3 as an electron-releasing group at position 4, which activates the pyridinum ring more than in position 3 Figure 8.In para-position, the hyper conjugation of the methyl group adjacent to the pyridine ring increases the electron density of N, so the effects of the N atom increase.
The amphiphilic ILs might be completely absorbed on the surface of the oil droplet, and the dispersion efficiency increases due to the decrease in oil-water interfacial tension, which facilitates the formation of small droplets in large quantities to be entered through the water column, resulting in an emulsion that seems to have a sun-like structure.The presence of non-polar ends in the synthesized ILs would help to form a stable emulsion via van der Waals forces [32,33].Increasing the IL concentrations results in an increase in the IL molecules adsorbed at the oil-water interface.According to the nature of IL molecules, they provide electrostatic and steric barriers to the coalescence of the dispersed oil drops.

Effect of Cation Structure
Three derivatives of pyridine were used: pyridine, 3-methylpyridine, and 4methylpyridine, and the results showed that the dispersion efficiency ranked as Ia < Ib < Ic which refers to the presence of CH3 as an electron-releasing group at position 4, which activates the pyridinum ring more than in position 3 Figure 8.In para-position, the hyper conjugation of the methyl group adjacent to the pyridine ring increases the electron density of N, so the effects of the N atom increase.The amphiphilic ILs might be completely absorbed on the surface of the oil droplet, and the dispersion efficiency increases due to the decrease in oil-water interfacial tension, which facilitates the formation of small droplets in large quantities to be entered through the water column, resulting in an emulsion that seems to have a sun-like structure.The presence of non-polar ends in the synthesized ILs would help to form a stable emulsion via van der Waals forces [32,33].Increasing the IL concentrations results in an increase in the IL molecules adsorbed at the oil-water interface.According to the nature of IL molecules, they provide electrostatic and steric barriers to the coalescence of the dispersed oil drops.

Stability of Crude Oil & Water Emulsions
In order to evaluate the performance of the prepared samples as oil spill dispersants regarding benchmark industrial standards, samples of crude oil and water were emulsified by Tween-80 (as a commercial oil spill dispersant) [34], Ia, Ib, and Ic, to evaluate the emulsion stability at 30 °C with (0.8:10 wt.%) DOR.Sample photos were taken immediately after sample preparation (at 0 min) and after 30 min.The results are shown in Figure 9. Firstly, all samples formed a completely homogenous emulsion.Consequently, the prepared emulsion with Tween-80 separated into two clear layers after 30 min.The upper, darker layer is oil-rich, while the lower layer is water-rich and lighter in color.The formed emulsions with Ia, Ib, and Ic were gradually stable, and Ic showed the most stability due to the highest dispersion and emulsification capacity compared to others.

Stability of Crude Oil & Water Emulsions
In order to evaluate the performance of the prepared samples as oil spill dispersants regarding benchmark industrial standards, samples of crude oil and water were emulsified by Tween-80 (as a commercial oil spill dispersant) [34], Ia, Ib, and Ic, to evaluate the emulsion stability at 30 • C with (0.8:10 wt.%) DOR.Sample photos were taken immediately after sample preparation (at 0 min) and after 30 min.The results are shown in Figure 9. Firstly, all samples formed a completely homogenous emulsion.Consequently, the prepared emulsion with Tween-80 separated into two clear layers after 30 min.The upper, darker layer is oil-rich, while the lower layer is water-rich and lighter in color.The formed emulsions with Ia, Ib, and Ic were gradually stable, and Ic showed the most stability due to the highest dispersion and emulsification capacity compared to others.

Group Composition of Undispersed Crude Oil
The treated crude oil with the different synthesized dispersants at 30 °C with the optimum ratio (0.8-10) wt.% was recovered, weighted accurately, and separated into its components via SARA analysis.Gravimetrically, the dispersion efficiency was determined as shown in Table 12.All results obtained through gravimetric and UV spectroscopy for efficiency determination led to the same trend with accepted values.
Results revealed that the higher asphaltene and resin content is directly proportional to the dispersion [35].The highest percentage of asphaltene and resin was recorded in dispersant Ic Table 13.

Group Composition of Undispersed Crude Oil
The treated crude oil with the different synthesized dispersants at 30 • C with the optimum ratio (0.8-10) wt.% was recovered, weighted accurately, and separated into its components via SARA analysis.Gravimetrically, the dispersion efficiency was determined as shown in Table 12.All results obtained through gravimetric and UV spectroscopy for efficiency determination led to the same trend with accepted values.
Results revealed that the higher asphaltene and resin content is directly proportional to the dispersion [35].The highest percentage of asphaltene and resin was recorded in dispersant Ic Table 13.The dispersion efficiency was also understood by studying the alkane's distributions (Figure 10).It was noticed that the gas chromatograph of crude oil is considered a fingerprint to be compared with the treated crude oil using Ia, Ib, and Ic.Clearly, it was found that there is a change in GC profiles; otherwise, there is a reduction in peak heights or disappearance of the number of paraffin peaks.Undispersed, treated crude oil with Ic dispersant recorded the highest reduction.3.6.2.HPLC for Polyaromatic Hydrocarbons (PAHs) Figure 11 illustrates the PAH fingerprints of the studied (blank) and treated crude oils via Ia, Ib, and Ic [36].PAHs were characterized by the formation of high-molecularweight aromatic rings (di-, tri-, tetra-, penta-, and hexa-).The presence of aromatics with high molecular weight aromatic rings (tetra-, penta-, and hexa-) is more carcinogenic than that of those of smaller molecular weight (di and tri).The disappearance of higher PAHs in undispersed, treated crude oil is referred to as the degradation of them as an effect of dispersion.HPLC results showed that the efficiency of dispersants' effect on PAH distribution increased, whereas penta and hexa aromatics (carcinogenic parts) decreased.

HPLC for Polyaromatic Hydrocarbons (PAHs)
Figure 11 illustrates the PAH fingerprints of the studied (blank) and treated crude oils via Ia, Ib, and Ic [36].PAHs were characterized by the formation of high-molecular-weight aromatic rings (di-, tri-, tetra-, penta-, and hexa-).The presence of aromatics with high molecular weight aromatic rings (tetra-, penta-, and hexa-) is more carcinogenic than that of those of smaller molecular weight (di and tri).The disappearance of higher PAHs in undispersed, treated crude oil is referred to as the degradation of them as an effect of dispersion.HPLC results showed that the efficiency of dispersants' effect on PAH distribution increased, whereas penta and hexa aromatics (carcinogenic parts) decreased.
weight aromatic rings (di-, tri-, tetra-, penta-, and hexa-).The presence of aromatics with high molecular weight aromatic rings (tetra-, penta-, and hexa-) is more carcinogenic than that of those of smaller molecular weight (di and tri).The disappearance of higher PAHs in undispersed, treated crude oil is referred to as the degradation of them as an effect of dispersion.HPLC results showed that the efficiency of dispersants' effect on PAH distribution increased, whereas penta and hexa aromatics (carcinogenic parts) decreased.

Acute Toxicity of the Prepared Dispersants
The toxicity of Ia, Ib, and Ic was tested to evaluate their toxicity against fish (Nile tilapia and Oreochromis niloticus).Results indicated an acute fish toxicity test (LC50) of concentrations greater than 100 ppm after 96 h of observation.These values indicate that the LC50s of Ia, Ib, and Ic record 13.25, 17.75, and 37.5 mg/L, respectively, which are located between 10 and 100 mg/L [32], providing that the synthesized dispersants are "slightly toxic" under the acute toxicity test.By comparing the LC50 of the prepared ILs with that of commercial Corexit 9500 (96 h LC50 = 35.9mg/L), it was found that the synthesized dispersants were less toxic.The results of the toxicity test and assessment of ILs against marine toxicity confirmed the ability to use the synthesized ILs as potential dispersants and safety applications for oil remediation with a slight toxicity effect on the marine ecosystem [37].

Conclusions
In the presented work, a series of dicationic ionic liquids (Ia, Ib, and Ic) were synthesized with different pyridine derivatives and well characterized via FT-IR, 1 H-NMR, and elemental analysis.They recorded good physico-chemical properties (surface activity and thermal stability).Density functional theory (DFT) parameters were studied to investigate the geometry optimization of electronic structures, energy gap (∆E), reactivity, hardness, and softness.The calculated parameters showed that Ic is the most effective dispersant; due to its high electron density and high aromaticity, it accelerates the penetration of crude oil, the adsorption on asphalt surfaces, and its aggregation.
The evaluation of Ia, Ib, and Ic as oil spill dispersants occurred by spectroscopic (U.V. spectroscopy) and gravimetric techniques at different temperatures (10 and 30 • C).The experimental data showed that the synthesized ILs ranked as Ia < Ib < Ic which refer to the presence of the electron-releasing group (CH 3 ) at position 4, activating the pyridinium ring more than in position 3. Dispersion efficiency increased with increasing the ratio of dispersant to crude oil; the optimum ratio of dispersants to crude oil (DOR) in this study was (0.8:10) wt.% and recorded 33.61% at 10 • C and 66.0% at 30 • C. Undispersed, treated crude oil with dispersants was studied via GC and HPLC techniques.The results showed that the dispersants affected saturated and PAH distributions.
The toxicity of Ia, Ib, and Ic was tested to evaluate their toxicity against fish, and the studies indicated that they are "slightly toxic" under the acute toxicity test (LC50 between 10 and 100 mg/L).Generally, the excellent oil spill dispersion and slight toxicity of the newly synthesized ILs could provide a real chance to be applied in industrial oil spill remediation.
1 H-NMR spectrum for Ia, Ib, and Ic were confirmed and detected, as shown in Figures 3-5 and Tables 4-6.Separations 2023, 10, x FOR PEER REVIEW 7 of 17 1 H-NMR spectrum for Ia, Ib, and Ic were confirmed and detected, as shown in Figures 3-5 and Tables 4-6.

Figure 5 .
Figure 5.Chemical structure of Ic.Elemental Analysis: Table7indicated that the experimental values of the prepared ILs are in acceptable agreement with the theoretically calculated values.

Separations 2023 , 17 Figure 10 .
Figure 10.C number of blank and treated crude oil with Ia, Ib & Ic.

Figure 10 .
Figure 10.C number of blank and treated crude oil with Ia, Ib & Ic.

Figure 11 .
Figure 11.PAHs of blank and treated crude oil with Ia, Ib & Ic.

Figure 11 .
Figure 11.PAHs of blank and treated crude oil with Ia, Ib & Ic.

Table 2 .
Physical characteristics of sea water.

Table 3 .
Inorganic Chemical Constituents of sea water.

Table 6 .
Table 7 indicated that the experimental values of the prepared ILs are in acceptable agreement with the theoretically calculated values.

Table 8 .
Surface parameters of the synthesized ILs.

Table 8 .
Surface parameters of the synthesized ILs.

Table 12 .
Gravimetric determination of efficiency.

Table 12 .
Gravimetric determination of efficiency.

Table 13 .
SARA analysis of treated crude oil.