skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Uncovering the Critical Factors that Enable Extractive Desulfurization of Fuels in Ionic Liquids and Deep Eutectic Solvents from Simulations
Environmental regulatory agencies have implemented stringent restrictions on the permissible levels of sulfur compounds in fuel to reduce harmful emissions and improve air quality. Problematically, traditional desulfurization methods have shown low effectiveness in the removal of refractory sulfur compounds, e.g., thiophene (TS), dibenzothiophene (DBT), and 4-methyldibenzothiophene (MDBT). In this work, molecular dynamics (MD) simulations and free energy perturbation (FEP) have been applied to investigate the use of ionic liquids (ILs) and deep eutectic solvents (DESs) as efficient TS/DBT/MDBT extractants. For the IL simulations, the selected cation was 1-butyl-3-methylimidazolium [BMIM], and the anions included chloride [Cl], thiocyanate [SCN], tetrafluoroborate [BF4], hexafluorophosphate [PF6], and bis(trifluoromethylsulfonyl)amide [NTf2]. The DESs were composed of choline chloride with ethylene glycol (CCEtg) or with glycerol (CCGly). Calculation of excess chemical potentials predicted the ILs to be more promising extractants with energies lower by 1-3 kcal/mol compared to DESs. Increasing IL anion size was positively correlated to enhanced solvation of S-compounds, which was influenced by energetically dominant solute-anion interactions and favorable solute-[BMIM] pi-pi stacking. For the DESs, the solvent components offered a range of synergistic, yet comparatively weaker electrostatic interactions that included hydrogen bonding and cation-pi interactions. An in-depth analysis of the structure of IL and DES systems is presented, along with a discussion of the critical factors behind experimental trends of S-compound extraction efficiency.  more » « less
Award ID(s):
2102038
PAR ID:
10512887
Author(s) / Creator(s):
;
Publisher / Repository:
American Chemical Society
Date Published:
Journal Name:
The Journal of Physical Chemistry B
Volume:
127
Issue:
28
ISSN:
1520-6106
Page Range / eLocation ID:
6354 to 6373
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The conversion of biomass to 5-hydroxymethylfurfural (HMF) holds substantial promise as a renewable energy source. Notably, HMF can be transformed into 2,5-bis(hydroxymethyl)furan (BHMF), a crucial reactant in biofuel production, but requires harsh operating conditions, H2, and precious metal catalysts. A recently reported Cannizzaro reaction of HMF to BHMF, characterized by its efficiency, mild conditions, and eco-friendliness, instead employed ionic liquids (ILs) to achieve high yields. In this study, combined quantum mechanical and molecular mechanical (QM/MM) simulations in conjunction with Metropolis Monte Carlo statistical mechanics and free-energy perturbation theory utilized M06-2X/6-31+G(d), PDDG/PM3, and the OPLS-VSIL force field to uncover important solute–solvent interactions present in the HMF to BHMF reaction pathway. The Cannizzaro reaction was examined in water and in five ILs composed of the 1-butyl-3-methylimidazolium [BMIM] cation coupled to hexafluorophosphate, tetrafluoroborate, thiocyanate, chloride, and bromide. Energetic and structural analysis of the rate-determining hydride transfer between HMF and the hydride-donor anion HMFOH− attributed the enhanced reactivity to highly organized solvent interactions featuring (1) hydrogen bonding between the ring protons of [BMIM] and the negatively charged carbonyl oxygen atoms on the transition structure, (2) favorable electrostatic interactions between the IL anions and solute hydroxyl groups, and (3) beneficial π–π stacking interactions between [BMIM] and the two aromatic rings present in HMF and HMFOH−. 
    more » « less
  2. Removing sulfur compounds from petroleum is essential in mitigating harmful emissions linked to fuel combustion. Problematically, thiophenic compounds that are readily present in fossil fuels resist conventional desulfurization methods. Extractive Desulfurization (EDS) provides an attractive alternative that can be selective for heterocyclic sulfur compounds, i.e., thiophene and its derivatives, through the choice of solvent employed. To this end, a burgeoning field has developed around the use of ionic liquids (ILs) given their ability to be fine-tuned with varying levels of polarity and solubility to suit the specific requirements of the desulfurization process. Hundreds of experimental studies featuring the use of IL technologies have provided encouraging trends for the design of more efficient extractants; however, conflicting data and a lack of a definitive understanding of important solute-ion interactions has presented challenges in advancing the field. More recently, computational investigations have been employed to unravel these key interactions and to inform design principles for future high-performance IL extractants. The myriad of intermolecular forces, e.g., coulombic, dispersive, and steric, and their subtle interplay present in IL-mediated EDS processes are prime for study using computational methodologies that include quantum mechanics (QM) at the ion-solute interaction level, molecular dynamics (MD) for the simulation of bulk-phase solvent properties, and the Conductor-Like Screening Model (COSMO) for high-throughput screening. This minireview summarizes computational advances and findings in the field of IL-mediated EDS in a format suitable for theoreticians and experimental chemists alike with discussions provided of future directions for the field. 
    more » « less
  3. Abstract This study examines the activity of chemisorbed CO2 species in the microenvironment formed by bifunctional ionic liquids (ILs) in the reactive capture and conversion (RCC) of CO2 to CO on silver. Comparative electroanalytical measurements with imidazolium based ILs were performed to probe the impact of electrostatic interactions, anion and cation basicity, and hydrogen bonding on RCC. Particularly, ILs with 1-ethyl,3-methylimidazolium ([EMIM]+) and 1-ethyl, 2,3-methylimidazolium ([EMMIM]+) cations and aprotic heterocyclic anions of 2-cyanopyrrolide ([2-CNpyr]) and 1,2,4-triazolide ([1,2,4-Triz]) were examined for RCC. It was found that anion–CO2 carbamate complexes facilitate RCC at significantly lower overpotentials compared to cation–CO2 carboxylate complexes. Additionally, [EMIM]+ was found to better stabilize anion–CO2 complexes than [EMMIM]+. Furthermore, it was found that 2-CNpyrH that naturally forms in CO2 absorption competes for electrode surface adsorption with the anion–CO2 carbamate complex, thereby reducing the electrochemical activity of the anion–CO2 complex. These results highlight the importance of IL structure in tuning the interfacial interactions and suggest that ILs with anion-dominated CO2 chemisorption enhances CO2 utilization in RCC applications. 
    more » « less
  4. null (Ed.)
    Ionic liquids (ILs) exhibit unique properties that have led to their development and widespread use for a variety of applications. Development efforts have generally focused on achieving desired macroscopic properties via tuning of the IL through variation of the cations and anions. Both the macroscopic and microscopic properties of an IL influence its tunability and thus feasibility of use for selected applications. Works geared toward a microscopic understanding of the nature and strength of the intrinsic cation-anion interactions of ILs have been limited to date. Specifically, the intrinsic strength of the cation-anion interactions in ILs is largely unknown. In previous work, we employed threshold collision-induced dissociation (TCID) approaches supported and enhanced by electronic structure calculations to determine the bond dissociation energies (BDEs) and characterize the nature of the cation-anion interactions in a series of four 2:1 clusters of 1-alkyl-3-methylimidazolium cations with the hexafluorophosphate anion, [2C n mim:PF 6 ] + . To examine the effects of the 1-alkyl chain on the structure and energetics of binding, the cation was varied over the series: 1-ethyl-3-methylimidazolium, [C 2 mim] + , 1-butyl-3-methylimidazolium, [C 4 mim] + , 1-hexyl-3-methylimidazolium, [C 6 mim] + , and 1-octyl-3-methylimidazolium, [C 8 mim] + . The variation in the strength of binding among these [2C n mim:PF 6 ] + clusters was found to be similar in magnitude to the average experimental uncertainty in the measurements. To definitively establish an absolute order of binding among these [2C n mim:PF 6 ] + clusters, we extend this work again using TCID and electronic structure theory approaches to include competitive binding studies of three mixed 2:1 clusters of 1-alkyl-3-methylimidazolium cations and the hexafluorophosphate anion, [C n-2 mim:PF 6 :C n mim] + for n = 4, 6, and 8. The absolute BDEs of these mixed [C n-2 mim:PF 6 :C n mim] + clusters as well as the absolute difference in the strength of the intrinsic binding interactions as a function of the cation are determined with significantly improved precision. By combining the thermochemical results of the previous independent and present competitive measurements, the BDEs of the [2C n mim:PF 6 ] + clusters are both more accurately and more precisely determined. Comparisons are made to results for the analogous [2C n mim:BF 4 ] + and [C n-2 mim:PF 6 :C n mim] + clusters previously examined to elucidate the effects of the [PF 6 ] - and [BF 4 ] - anions on the binding. 
    more » « less
  5. Ionic liquids (ILs) are gaining attention as protein stabilizers and refolding additives. However, varying degrees of success with this approach motivates the need to better understand fundamental IL-protein interactions. A combination of experiment and simulation is used to investigate the thermal unfolding of lysozyme in the presence of two imidazolium-based ILs (1-ethyl-3-methylimidazolium ethylsulfate, [EMIM][EtSO 4 ] and 1-ethyl-3-methylimidazolium diethylphosphate, [EMIM][Et 2 PO 4 ]). Both ILs reduce lysozyme melting temperature Tm , but more gradually than strong denaturants. [EMIM][Et 2 PO 4 ] lowers lysozyme Tm more readily than [EMIM][EtSO 4 ], as well as requiring less energy to unfold the protein, as determined by the calorimetric enthalpy ΔH. Intrinsic fluorescence measurements indicate that both ILs bind to tryptophan residues in a dynamic mode, and furthermore, molecular dynamics simulations show a high density of [EMIM] + near lysozyme’s Trp62 residue. For both ILs approximately half of the [EMIM] + cations near Trp62 show perfect alignment of their respective rings. The [EMIM] + cations, having a "local" effect in binding to tryptophan,likely perturb a critically important Arg-Trp-Arg bridge through favorable π − π and cation-π interactions. Simulations show that the anions, [EtSO 4 ] - and [Et 2 PO 4 ] - , interact in a "global" manner with lysozyme, due to this protein’s strong net positive charge. The anions also determine the local distribution of ions surrounding the protein. [Et 2 PO 4 ] - is found to have a closer first coordination shell around the protein and stronger Coulomb interactions with lysozyme than [EtSO 4 ] - , which could explain why the former anion is more destabilizing. Patching of ILs to the protein surface is also observed, suggesting there is no universal IL solvent for proteins, and highlighting the complexity of the IL-protein environment. 
    more » « less