Molecular modeling insights in the extraction of benzene from hydrocarbon stream using deep eutectic solvent

Atomistic molecular dynamics simulation and COSMO-SAC approach for enhanced 1,3-propanediol extraction with imidazolium-based ionic liquids

CONTEXT

1,3-Propanediol (1,3-PDO) is a key chemical in various industries, including pharmaceuticals and material sciences, and is projected to see significant market growth. However, the current challenges in its downstream processing, particularly in terms of cost and efficiency, highlight the need for innovative solutions. Our study delves into using ionic liquids (ILs) as a potential alternative, aiming to address these critical separation challenges more sustainably and efficiently. In this study, we utilized molecular dynamics (MD) simulations and the COSMO-SAC to examine 1,3-propanediol (1,3-PDO) extraction using four imidazolium-based ionic liquids with 1-butyl-3-methylimidazolium [Bmim] cation and with different anions bis(pentafluoroethanesulfonyl)imide [NPF2]-, bis(trifluoromethylsulfonyl)imide [NTF2]-, thiocyanate [SCN]-, and trifluoromethanesulfonate [TFO]-. Molecular dynamics simulations, incorporating analysis of radial distribution functions (RDF) and spatial distribution functions (SDF), revealed that [Bmim][SCN] and [Bmim][TFO] exhibit enhanced interactions with 1,3-PDO. Notably, [Bmim][SCN] formed the most hydrogen bonds, averaging 1.639 per molecule, due to its coordinating [SCN]- anion. This was in contrast to the fewer hydrogen bonds formed by non-coordinating anions in [Bmim][NPF2] and [Bmim][NTF2]. In ternary systems, [Bmim][SCN] and [Bmim][TFO] demonstrated superior selectivity for 1,3-PDO extraction compared to the other ionic liquids, with selectivity values around 29. These findings, supported by COSMO-SAC predictive modeling, highlight the potential of [Bmim][SCN] as a promising candidate for 1,3-PDO extraction, emphasizing the importance of anion selection in optimizing ionic liquid properties for this application.

METHODS

In our study, we employed MD simulations, incorporating the OPLS-AA force field, and COSMO-SAC to investigate the extraction of 1,3-PDO using imidazolium-based ionic liquids: [Bmim][NTF2], [Bmim][NPF2], [Bmim][SCN], and [Bmim][TFO]. The MD simulations were conducted using LAMMPS software, focusing on elucidating the RDF, SDF, and hydrogen bonding. Analysis of the distribution coefficient (β) and selectivity (S) for the ternary mixture was also conducted. These aspects of the simulation were analyzed using TRAVIS and VMD software. Additionally, the COSMO-SAC model was employed to determine the activity coefficients of 1,3-PDO in the ionic liquids, with molecular optimization conducted using Gaussian16 and sigma profile calculations performed using COSMO-SAC.

Molecular Dynamic Insights into the Distinct Solvation Structures of Aromatic and Aliphatic Compounds in Monoethanolamine-Based Deep Eutectic Solvents

Deep eutectic solvents (DESs) are developing as an alternate medium for aromatic extraction, especially benzene and thiophene from aliphatic hydrocarbon mixtures. In this work, molecular dynamics (MD) simulations were first used to investigate the solvation structure of benzene, thiophene, and n-hexane in monoethanolamine-based DESs. It reveals the liquid structures in the adjacent neighbor shells, which is a function of electron-withdrawing sulfur attached to thiophene and the π-electron cloud of benzene. The intermolecular forces between aromatic, aliphatic, and DES components are analyzed in van der Waals and hydrogen bond interactions. The chloride ions serve as a charge carrier bridge between choline and monoethanolamine precursors. The solvation of benzene, thiophene, and n-hexane in the DESs depends on volume expansion and minor solvent structural changes. Density functional theory results provided information on the mechanism of short-range interactions between organic solutes and studied DES. It aids in understanding the structural orientations of a DES with the addition of solutes, essential to the formation of DES. The solvation shell structure and characteristics were investigated in tandem with the possibility of benzene and thiophene clustering. The ¹H NMR and 2D ¹H-¹H-NOESY were used to investigate the intermolecular interactions between benzene, thiophene, and n-hexane with monoethanolamine-based solvents. It concludes that high-ordered DES1 is more inclined to higher solubility than lower-ordered ones with a higher molar ratio of monoethanolamine. The solvation was reduced because the entropy gain was not maximized in the lower ordered DESs.

Computer Simulations of Deep Eutectic Solvents: Challenges, Solutions, and Perspectives

Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications, such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can predict and reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.

Evaluation of hybrid solvents featuring choline chloride-based deep eutectic solvents and ethanol as extractants for the liquid-liquid extraction of benzene from n-hexane: towards a green and sustainable paradigm

Deep eutectic solvents (DESs) have high viscosities, but known to be mitigated by addition of suitable co-solvent. The effect of such co-solvent on the extraction efficiency of the hybrid solvent is hardly known. This study examined the effect of ethanol on three choline chloride-based DESs (glyceline, reline, and ethaline) by mixing each in turn with ethanol in various volume proportions. The hybrid solvents were evaluated for the extraction of benzene from n-hexane. Pseudo-ternary liquid–liquid equilibrium data were obtained using the refractive index method at 303 K and 1 atm for the systems, n-hexane (1) + benzene (2) + hybrid solvent (glyceline/ethanol, ethaline/ethanol, reline/ethanol) (3), and used to evaluate distribution coefficient (D) and selectivity (S). Furthermore, the physicochemical properties of the hybrid solvents were also determined. The results indicate increase in selectivity with increasing ethanol addition up to 50% and decrease with further addition. All hybrid solvents with 50% ethanol outperform sulfolane and are suitable replacement for same as green and sustainable extractant for aromatics from aliphatics. The glyceline + 50% ethanol emerged the overall best with 49.73% elevation in selectivity and 41.15% reduction in viscosity relative to the neat glyceline. The finding of this study is expected to fillip the drive for paradigm shift in petrochemical industries.

Liquid-Liquid Equilibrium of Deep Eutectic Solvent-Aromatic-Aliphatic Ternary Systems: Experimental Study with COSMO Model Predictions

Common solvents used for aromatic extraction from aliphatics typically degrade into toxic compounds, while green alternatives perform poorly compared to the state-of-the-art solvents. Deep eutectic solvents (DES) are a novel solvent type made of hydrogen bond donors (HBD) and hydrogen bond acceptors (HBA). DES have been applied in various applications, including advanced separations. In this study, DES were studied experimentally and using the Conductor-like Screening Model (COSMO) to separate benzene from cyclohexane as model compounds for an aromatic:aliphatic system. Both equilibrium and kinetic studies were performed to determine the liquid liquid equilibrium (LLE) and mass transfer rate for the DES-based separation. Selected HBAs including tetrabutylammonium bromide (N4444Br), tetrahexylammonium bromide (N6666Br), choline chloride (ChCl), and methyltriphenylphosphonium bromide (METPB) were paired with HBDs including ethylene glycol (EG) and glycerol (Gly). COSMO was used, with adjustments to reflect DES specific interactions, to predict the liquid-liquid equilibrium (LLE). COSMO results showed that ChCl and N6666Br-based DES extracted too little benzene or too much cyclohexane, respectively, to be considered for experimental evaluation. Overall, the COSMO model predictions for LLE of EG-based DES were very accurate, with root-mean-square deviations (RMSD) below 1% for both N4444Br:EG and METPB:EG. The glycerol systems were less accurately modeled, with RMSD’s of 4% for N4444Br:Gly and 6% for METPB:Gly. The lower accuracy of glycerol system predictions fmay be due to limitations in COSMO for handling glycerol’s influence on polarizability in the DES that is not seen in EG-based DES. Mass transfer kinetics were determined experimentally for DES and the results were fit to a first order kinetics model. METPB:Gly had the highest mass transfer coefficient at 0.180 min−1, followed by N4444Br:EG at 0.143 min−1. N4444Br:Gly and METPB:EG had the lowest mass transfer coefficients at 0.096 min−1 and 0.084 min−1, respectively. It was found that mass transfer rate was not directly related to maximum benzene solubility, as N4444Br:EG and METPB:Gly had the highest and lowest benzene removal, respectively, but had similar mass transfer coefficients.