Quantum chemical calculation, molecular dynamics simulation and process design for separation of heptane - butanol using ionic liquids extraction

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.

Process Simulation and Optimization on Ionic Liquids

Ionic liquids (ILs) are promising alternative compounds that enable the development of technologies based on their unique properties as solvents or catalysts. These technologies require integrated product and process designs to select ILs with optimal process performances at an industrial scale to promote cost-effective and sustainable technologies. The digital era and multiscale research methodologies have changed the paradigm from experiment-oriented to hybrid experimental-computational developments guided by process engineering. This Review summarizes the relevant contributions (>300 research papers) of process simulations to advance IL-based technology developments by guiding experimental research efforts and enhancing industrial transferability. Robust simulation methodologies, mostly based on predictive COSMO-SAC/RS and UNIFAC models in Aspen Plus software, were applied to analyze key IL applications: physical and chemical CO2 capture, CO2 conversion, gas separation, liquid-liquid extraction, extractive distillation, refrigeration cycles, and biorefinery. The contributions concern the IL selection criteria, operational unit design, equipment sizing, technoeconomic and environmental analyses, and process optimization to promote the competitiveness of the proposed IL-based technologies. Process simulation revealed that multiscale research strategies enable advancement in the technological development of IL applications by focusing research efforts to overcome the limitations and exploit the excellent properties of ILs.

Ionic liquids as a useful tool for tailoring active pharmaceutical ingredients

Ionic liquids (ILs) have been widely used in biomedical and pharmaceutical fields as solvents or permeation enhancers. Recently, more and more researchers focused on optimizing the physicochemical properties of active pharmaceutical ingredient (API) by ILs technology. Converting APIs into ILs (API-ILs) has shown great potential for drug delivery by eliminating polymorphism, tailoring solubility, improving thermal stability, increasing dissolution, controlling drug release, modulating the surfactant properties, enhancing permeability of APIs and modulating cytotoxicity on tumor cells. In addition, API-ILs are also used in various formulations as active ingredients, such as solutions, emulsions, even tablets or nanoparticles. This paper aims to review current status of API-ILs, including the rational and design, preparation and characterization, the improvement on the physicochemical characteristics of APIs, the compatibility of API-ILs with various formulations, and the future prospects of API-ILs in biomedical and pharmaceutical fields.

Relating the structure and dynamics of ionic liquids under shear by means of reverse non-equilibrium molecular dynamics simulations

The effect of the shear rate on the viscosity and the structure of 1-ethyl-3-methylimidazolium based ionic liquids with three different anions (tetrafluoroborate, dicyanamide, and bis(trifluoromethylsulfonyl)imide) was studied by means of reverse non-equilibrium molecular dynamics (RNEMD) simulations using a polarizable force field. The three liquids display a Newtonian plateau followed by a shear thinning regime at shear rates of the order of GHz. Even though the main features of the liquid structure remains under shear, systematic changes were noticed at the GHz rates, with coordination shells becoming more diffuse as noticed by the reduction in the difference between consecutive maxima and minima in the radial distribution function. Interestingly, these structural changes with the shear rate can be precisely fitted using the Carreau equation, which is a well-known expression for the shear rate dependence of the viscosity. The fitting parameters for different distributions can be used to explain qualitatively the shear thinning behavior of these liquids. In the GHz range, the cations and, in a minor extension, some anions, tend to assume preferentially a parallel orientation with the flux, which contributes to the shear thinning behavior and may have consequences for adhesion in applications as lubricants.