Next Article in Journal
Identification of Phytoconstituents in Leea indica (Burm. F.) Merr. Leaves by High Performance Liquid Chromatography Micro Time-of-Flight Mass Spectrometry
Previous Article in Journal
All Solid-State Poly (Vinyl Chloride) Membrane Potentiometric Sensor Integrated with Nano-Beads Imprinted Polymers for Sensitive and Rapid Detection of Bispyribac Herbicide as Organic Pollutant
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Screening of Suitable Ionic Liquids as Green Solvents for Extraction of Eicosapentaenoic Acid (EPA) from Microalgae Biomass Using COSMO-RS Model

by
Shiva Rezaei Motlagh
1,
Razif Harun
1,*,
Dayang Radiah Awang Biak
1,
Siti Aslina Hussain
1,
Wan Azlina Wan Ab Karim Ghani
1,
Ramin Khezri
1,
Cecilia Devi Wilfred
2 and
Amal A. M. Elgharbawy
3
1
Department of Chemical and Environmental Engineering, Faculty of Engineering, University Putra Malaysia, UPM, Serdang, Selangor 43400, Malaysia
2
Department of Fundamental and Applied Sciences, Centre of Research in Ionic Liquids (CORIL), Universiti Teknologi Petronas, UTP, Bandar Seri Iskandar, Perak 32610, Malaysia
3
International Institute for Halal Research and Training (INHART), International Islamic University Malaysia, Gombak, Kuala Lumpur 50728, Malaysia
*
Author to whom correspondence should be addressed.
Molecules 2019, 24(4), 713; https://doi.org/10.3390/molecules24040713
Submission received: 22 December 2018 / Revised: 18 January 2019 / Accepted: 20 January 2019 / Published: 16 February 2019
(This article belongs to the Section Computational and Theoretical Chemistry)

Abstract

:
Omega-3 poly unsaturated fatty acids (PUFA) particularly eicosapentaenoic acid (EPA), and docosahexaenoic acid (DHA), have many health benefits including reducing the risk of cancer and cardiovascular disease. Recently, the use of ionic liquids (ILs) in lipid extraction from microalgae provides the potential to overcome common drawbacks, offers several other benefits. To date, very limited researches are available to focus on extracting microalgae lipid and PUFA in particular by using ILs. The objective of current work is to screen the potential ILs that can be applied in EPA extraction. In this study, fast ILs screening was performed with the help of a conductor like screening model for real solvents (COSMO-RS) and the ILs with higher capacity values for use in extraction of EPA were compared. According to the results, the highest capacity for EPA extraction among 352 screened cation/anion combinations belongs to [TMAm][SO4]. It is expected to achieve a higher yield of EPA once applying this combination as the solvent in the process of extraction. ILs with small anions were observed to have higher capacities, as well possessing higher charge density compared to larger ones, and therefore, they are more preferable for extraction purposes. Moreover, shorter alkyl chain cations are preferred when using imidazolium-based IL, which agrees with experimental data.

1. Introduction

Health benefits provided by omega-3 polyunsaturated fatty acids (PUFAs), including eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) in particular are noteworthy to the society as it reduces the risk of cancer, autoimmune, cardiovascular disease, inflammatory disorders, cystic fibrosis, disrupted neurological function, bowel disease, and mental illness [1,2,3,4,5,6]. The use of PUFA as health prospects has gained great attention in last years. Current dietary health organization recommends a daily intake of 250–500 mg of EPA + DHA for primary and secondary prevention of coronary heart disease [7,8]. Additionally, according to Shahidi and Ambigaipalan human bodies requires the supplementary sources of EPA and DHA as they are vulnerable to shortages of necessary enzymes to place a double bond at the omega-3 position and hence are unable to synthesize them to a sufficient extent [9]. Those are one of the two classes of essential fatty acids that have such characteristics and therefore must be taken from diet and supplements. Thus, the external source of omega-3 PUFA is necessary.
Lately, the main supply of both EPA and DHA are provided from fish, however, the intake of fish oil for long period may cause a deficiency of vitamin E due to the high level of vitamin A and D involved in fish lipids [10,11]. Many harmful contaminants such as methylmercury, copper, and organic pollutants as polychlorinated biphenyls (PCBs) or dioxins are also found in some species of fish especially in salmon, sardine, anchovy, and tuna, which may impose toxic effects on human health [11,12], and therefore, alternative sources of omega-3 PUFA are required. Microalgae oil and fish oil have the same amount of omega-3 but microalgae oil has more advantages in terms of health and less toxicity compared to fish oil [13]. The oil derived from microalgae contains squalene and phytosterols, which indicates several health benefits to the human body and it not having an unpleasant smell and cholesterol. The process of converting microalgae into omega-3 PUFA mainly consists of four steps, including microalgae cultivation, harvesting, cell disruption (lipid extraction), and transesterification (GC/FID analysis). There are numerous methods that have been used for lipid extraction, but most of them using conventional solvents, which are mostly volatile, aromatic, thermally unstable, and environmentally unfriendly, in addition to causing health and safety problems [14]. For instance, the use of solvents such as hexane or isopropanol in Soxhlet extraction likewise chloroform, methanol, and water in the Bligh and Dyer method (both methods are considered to be the most widely used techniques for lipid extraction) introduced a few concerns due to the toxicity associated with the solvents that brings biosafety issues and is required to be replaced by environmentally friendly, biocompatible, and less toxic solvents.
Therefore, solvent extraction using ionic liquids (ILs) are preferable and has potential in the extraction of the microalgae lipids as an alternative to traditional solvents that use toxic solvent, long duration time, and demand a huge amount of energy [15]. In solvent extraction, it is crucial to have high capacity to extract the targeted solute. ILs are salts with low melting point and consist of large asymmetric organic cations and smaller organic or inorganic anions [16]. The cations of ILs are generally composed of nitrogen contacting ring like imidazolium, pyrrolidinium to which different functional groups can be added [17]. Additionally, ILs are often considered to be “green” solvents because of their very low vapor pressure, great thermal stability with a broad electrochemical window, and high ionic conductivity compared to conventional solvents [18]. In a study by Pan et al. [19], the microwave-assisted extraction with [BMIM][HSO4] as IL was observed to have higher yield rather than the conventional method of oil bath using an organic solvent for the extraction of lipid from Nannochloropsis salina, Chlorella sorokiniana, and Galdieria sulphuraria. Moreover, Kim et al. [20] used the mixture of IL like [BMIM][CF3SO3] with methanol to extract lipid (biodiesel purpose) from microalgae C.vulgaris, which obtained 19%. According to the literature, ILs in general have prominent advantages over organic solvent like hexane, chloroform, however, there are only a few studies available to use ILs for extraction of long chain fatty acids such as EPA and DHA. It is noticeable that using ILs for processing is costly and time-consuming, and thus a reasonable screening of different cations and anions combination is obligatory prior to selecting the appropriate ILs. Therefore, it seems crucial to predict ILs behavior before the experimental application. The conductor like screening model for real solvent (COSMO-RS) proposed by Klamt is a model that can be used for such fast prior prediction [21,22]. COSMO-RS is a new method to predict the thermodynamic properties of fluid and liquid mixture based on quantum chemistry concept of density functional theory (DFT) [23]. Chosen the suitable IL for the extraction efficiency, some thermodynamic properties are important. The solvent capacity, selectivity at infinite dilution and performance index was used many times by previous studies for solvent liquid-liquid extraction prediction, which was calculated under the activity coefficient at the infinite dilution. For instance, Mansoure et al. [24] used the COSMO-RS tool for the prediction of ammonium ILs capacity and then selectivity values for β-carotene separation from n-hexane. Additionally, the study reported by Zeeshan et al. [25] COSMO-RS software used for the calculation of activity coefficient at infinite dilution and then selectivity, capacity, and performance index to predict different types of ILs potential for asphaltene extraction as a part of liquid-liquid extraction. Moreover, in a recent study by Xing et al. [26], it has been used COSMO-RS for prediction and exploring the potential separation mechanism of long chain fatty acids (EPA, DHA, CLAs, OA, and SA) through the calculation of different solvents capacity, selectivity, and performance index values. In another study, adjustable parameters were re-optimized in COSMO-RS to fit for the systems containing ionic liquids (ILs), whereby a vast numbers of activity coefficients and at infinite dilution and CO2 solubility were collected from references and used as a training set. The results showed that the predicted results by COSMO-RS model with the new optimized parameters are in agreement with experimental data [27]. As pointed out by previous literatures, COSMO-RS only require the structural information of interacting species in order to calculate their related thermodynamic properties through a set of defined mathematical equations and hence one does not need to perform any experimentations in advance [28,29]. It is noteworthy to mention that, although COSMO-RS is a powerful screening tool, this software has been addressed by a number of authors as a mean field model that does not account for solvent molecular structure and correlations. Therefore, it is suggested that when using COSMO-RS, the important effects of explicit interactions between organic molecules and water, which can lead to charge transfer between the solute and water, must also be considered [30,31,32,33,34].
While the ILs screening potential seems crucial for extraction of omega-3 (EPA, DHA and ALA) compounds from biomass, however to the best of the author’s knowledge, no significant report has been published yet to concern the application of COSMO-RS in that particular regard. Therefore, this study strictly help and open a new horizon about ILs interaction toward valuable omega-3 compound for researchers who seek to replace the conventional solvent with proper green solvent (ILs) without any requirement to perform expensive and time consuming lab experimentations. Moreover, what makes the COSMO-RS screening even more valuable is that it provides the facility to study over the variety of 352 combination of cations and anions as ILs to act as solvents in interaction with EPA compounds.
In this study, COSMO-RS has been used as the screening methodology to select the most suitable ILs as a part of solid-liquid extraction of EPA from microalgae. This work aims to predict infinite capacity values of ILs for EPA extraction and without the use of any experimental data.

2. Results and Discussion

The principle of COSMO-RS as a physical model considers the interaction energy of polarization charge densities between the sigma bonds. Net surface charge of EPA is calculated in terms of Sigma profile (σ-profile), which is related to the charge density profiles of molecular surface and Sigma potential (σ-potential) in order to clear its chemical nature. The polarity of the molecule is typically determined by the help of molecular electrostatic potentials. The best approach to investigate the molecular similarity is through comparing the electron density or electrostatic potential that is created by molecules, considering all space within the van der Waals surfaces of the molecules or within all regions outside the molecular cores [35]. The interactions of the electrostatic (polarity) among IL systems will integrate the hydrogen-bonding and hydrophobicity of a structure in the method of COSMO-RS. Therefore, it demonstrates multiple solvation interactions of ILs and presents the quantity and clear visual of polarity distribution on molecular surfaces with the help of sigma (σ) profiles [36]. According to the σ-profile, higher absolute values of σ indicates the stronger compound as a hydrogen bond donor or a hydrogen bond acceptor [37].
The chemical structure of EPA and its sigma surface are illustrated in Figure 1. The results of σ-profile and σ-potential are indicated in Figure 2 and Figure 3, respectively. Significant terms related to qualitative description of EPA molecules including hydrogen bonding, polarity, and lipophilicity/hydrophilicity can be visualized with the help of COSMO-RS 3D screening charge distribution (sigma surface). Charge distribution is coded in colors such as red color, which represents the negative (−ve) charge, whereas blue and green represents the positive (+ve) and neutral (0) charge on the molecular surface, respectively. The screening charge density surface (s-surface) contains all relevant information for COSMO-RS to calculate the chemical potential.
The capacity values related to five types of cation-base ILs including imidazolium, pyridinium, pyrrolidinium, piperidinium, and tetra-methyl ammonium with 22 types of anions were screened at 298.15 K and the results are summarized in Figure 4a–e. The results are followed by recommendations for the selection of suitable ILs to be applied in EPA extraction.
Generally, sigma profile divided into three important parts. The first part, where σ < −0.01 e/nm2, shows the hydrogen bond donor part. The second part extending between −0.01 < σ < +0.01 e/nm2, which belongs to the non-polarized area and the third part, where σ > +0.01 e/nm2, shows the hydrogen bond acceptor part. It is remarkable that there are a series of sharp a peak in the non-polarized area and there is also a slight peak in hydrogen bond acceptor and donor. Consequently, EPA is highly non-polar and is due to all of the peaks merely crowded in the non-polar region.
Chemical potential or sigma potential (μ(σ)) demonstrate the deep understanding about which ILs structures have a high impact on interaction with EPA molecules. Sigma potential (μ(σ)) like sigma profile is divided into three important areas. The first zone where σ < −0.01 e/nm−2 belongs to the interaction with hydrogen bond donors. The second zone where −0.01< σ < 0.01 e/nm−2 is related to non-polar area and the third zone, where σ > 0.01 e/nm−2 belongs to the hydrogen bond acceptor part. Sigma potential shows that EPA can interact with a polar compound due to the peaks located in negative values on the y axis from 0 < (μ(σ)) < −0.7. Additionally it has high attraction toward both hydrogen-bond donor and acceptor molecules.

Investigation of the Effect of Cations Alkyl Chain Length on EPA Extraction

In this study, the capacity values represent the interaction between the chemical structure of ILs and EPA molecules. Further, the COSMO-RS method is used to predict the behaviour of 22 types of anions through the investigation of the effect of increasing cation alkyl chain length base ILs towards EPA in the case of imidazolium, pyridinium, pyrrolidinium, piperidinium, and tetraalkyl ammonium. From the screening, it can be seen that increasing the alkyl chain length in the case of imidazolium cations, had a negative impact on the extraction capacity of the IL (Figure 4 and Supplementary Materials).
This seems true as the capacity could be arranged in the order [EMIM] > [BMIM] > [HMIM]> [OMIM]. However, the IL capacity is remarkably high when combined with the anions, SO42− > Cl > Br. It is speculated that the interaction between the cation and the anion in these four exceptions is due to the absence of any alky chain in the anion, which makes the interaction (cation-anion) stronger. In addition, the charge density on the cation decreases and the molar volume increases according to the increase in alkyl chain length. The flexible alkyl chain may avoid the ions to pack effectively and cause a decrease in density as the result. It is expected that the increase in alkyl chain length and higher flexibility as the result, may interrupt the movement of IL component and avoid them to pass by each other. Furthermore it may increase the dispersion forces due to the enlarged molecular volume [38].
The inorganic ion’s influence seems to be well described by the Hofmeister series. The Hofmeister effect has been studied for decades and in several aspects related to interfacial phenomena such as the behaviour of colloids, the denaturation of proteins, and the structure of microemulsions. The mechanism of the Hofmeister effect is not yet clear, however the notable fact has been proven that the interface within the hydrophilic and hydrophobic environments can be stabilized with presence of ions in the solution and their connections with the charged or polar head groups of amphiphilic molecules [39]. It is also noteworthy that the small anions can easily be replaced by the EPA hydroxyl group.
On the other hand, increasing the alkyl chain length attached to the cation, 3-methyl-pyridinium, contributed to the reduction of the extraction capacity in general, regardless of the anion used. Together, chloride followed by sulfate was the best among other anions tested in the screening. While the capacity value for the pair 1-ethyl-3-methyl-pyridinium and SO4 was the highest, a general conclusion cannot be drawn based on the simulation. However, it can be accepted that according to properties of S=O as a good proton acceptor in forming H-bonds with OH from EPA (Figure 5), therefore the hydrogen bond (S=OH) between the cation and anion pair is strong. The regarded fact is even more emphasized while dealing with imidazolium-based IL with SO42− as the anion. The different changing patterns in hydrogen bond strength suggest the possible existence of selective interactions among water, anion, and acidic protons on the imidazolium ring [40].
EPA has two hydrogen bond acceptors and one hydrogen bond donor, and therefore, possibilities to form H-bonding with another EPA molecule is likely to occur during extraction, thus one molecule could attract another molecule, as a chain reaction. The cation, in this case, has no protonated hydrogen, so no hydrogen bond could be from between the alkyl group and the fatty acid. However, it is possible to form H-bond with the nitrogen due to the lone pair. This makes the proposed scheme a likelihood during extraction.
While poor capacity was observed in case of Br with [EMPyr], higher values were recorded when paired with [EMPyrro] in the order: SO42− > Cl > Br > propanoate. Likewise, [CnMPIP] has a similar trend when paired with the anions: SO42−, Cl, Br, and propanoate, although the cyclic nature of the cation might slightly reduce the capacity.
On the other hand, one exception was observed with tetramethylammonium-IL. With the large values of capacity when the anions are paired with SO42−, Cl, Br and propanoate, nitrate [NO3] has a different trend from previous cases. The nitrate anion has an electron to donate as well as the lone pair on the nitrogen, which makes it into a potential electron donor (able to bind). Moreover, looking into the tetra-methyl ammonium cation, it has positively charged nitrogen as all electrons are occupied with methyl groups. This may explain the low capacity of the pair (tetra-methyl ammonium nitrate) in comparison to Cl, Br which possess higher electronegativity.
It may be predicted that in case larger alkyl groups, the anion can be at a more stable position, which could hinder its capacity for fatty acids extraction. Another suggestion is that ILs with hydrophobic anions (such as NTF2) and long alkyl chain are a better choice for the extraction of non-polar substances, such as EPA, as the long carbon chain reduced the polarity of the FFA. On the other hand, ILs such as [EMIM]Cl are highly hydrophilic.
In this situation, the correlation between n-3 PUFA extraction and ionic liquids structures (aromatic or delocalized cation) for more enhanced PUFA extraction was suggested by Cheong et al. (2011) [41]. The hydrophilic or hydrophobic characteristics of ILs is significantly dependent on the structures of cations and anions, however the water miscibility of ILs are highly associated with the anion portion. The IL based on [PF6] and [Tf2N] are typically water-immiscible and are therefore desirable for forming biphasic systems. The regarded ILs are widely used for the target of extraction [42]. Li et al., (2009) was the first to report a study on the recovery of n-3 PUFA methyl esters from fish oil in which the ILs were combined with silver salts, and as the result, the hydrophobic ILs with large anions of low lattice energy were selected as the most desirable output [43,44]. Viscosity plays a crucial role, in the extraction of which, lower viscosity extensively improves the process. Moreover from the results, it was observed that the extraction capability of pyrrolidinium and tetra ammonium-based ILs were considerably lower than imidazolium and pyridinium-based ILs as it was below 20% for n-3 PUFAs and lower than 6% for n-3 PUFAEEs.
It can be noticed from these findings that the enrichment efficiency of n-3 PUFAs are strongly correlated to the aromatic structure of ILs. The correlation can be explained as the IL with aromatic structure is able to form temporary weak π-π interactions with polyunsaturated bonds in n-3 PUFA and enriches the component in the IL phase. It is found that ILs with cations containing aromatic rings, namely, the N,N′-dialkylimidazolium ([BMIM]PF6, [OMIM]PF6, [BMIM]BF4, [OMIM]BF4, and [BMIM]OTF) and N-alkylpyridinium ([BuMePyr]DCN, [BuMePyr]OTF, and [MeOcPyr]BF4) types have higher selectivity and extraction capability for n-3 PUFA and n-3 PUFAEE, as compared to those without aromatic structure. It can be concluded that the selectivity increases and hence the extraction capability improves with alkyl chain length from C4 to C8. However, it is not yet accurate to select the IL with higher selectivity merely by looking through cations. For example, in a study by Reference [41], same cation was used with different anions in screening ILs and it was observed that smaller anions have better extraction capability for PUFA. Fluoride (F) has the highest electronegative value and has the smallest anion size among the other anions in the comparison, [F = Ne = 1s22s22p6], which might be one of the reasons of its extraction capability, besides the stable electron configuration (the corresponding noble gas), although the COMSO-RS analysis demonstrated the opposite.
Table 1 summarize Figure 4a–e for the ILs with high potential of EPA extraction in terms of infinite dilution activity coefficient and capacity values. In the following table, there is shortlist of the ILs with high capacity and lower activity coefficient at infinite dilution. Activity coefficient at infinite dilution is a term firstly introduced and discussed by Reference [45]. In this study, the activity coefficient is a term to account for the effective concentration of a solute in a solution of which its chemical potential can be expressed by the activity of the substances present in the solution. When the concentration of the solute tends to zero, the terms are hence called the activity coefficient at infinite dilution [46,47]. Activity coefficient at infinite dilution provides a deep vision on the degree of non-ideal behaviour in ILs–EPA compounds mixtures. Low activity coefficient values show high solubility of the ILs with weak solute-ILs interaction and vice versa, as described by References [48,49].
According to the listed ILs in Table 1, mono-cations were the selected cations since the anion (SO42−) with the charge of (−2) is outstanding and hence the overall charge for the IL molecule will be (−1). The presence of free charge is expected to increase the polarity and therefore increase the solvent extractive particularly.
The final conclusion about the rate of mass transfer, the solubility of all components of the model solution, the distribution coefficients, and selectivity must be determined in order to decide the suitable IL for extraction [50].

3. Methodology

3.1. Computational Details of COSMO-RS

The molecules geometries were optimized by turbo mole 6.3 program package (TURBOMOLE GmbH (ltd), Univerisity of Karlsruhe Karlsruhe, Germany) that follows the density functional theory (DFT), to generate COSMO files using Becke and Perdew (BP) functional B88-86 with a triple-zeta valence with polarization (TZVP), and the resolution of identity standard (RI) approximation as a basis set [51].
COSMO-RS is a theory with the purpose of quantum theory, surface interactions, dielectric continuum models, and statistical thermodynamic properties that has been described for the first time in a work by Reference [23]. At first, in 1998 the application of COSMO-RS was limited by calculation of the activity coefficient in infinite dilution for vapour-liquid equilibria of binary mixtures. Then, the usage of COSMO-RS go further to calculate all kinds of phase equilibrium such as vapour-liquid, liquid-liquid, and solid-liquid predictions [52,53].
In this study, the COSMO-RS approach has been used to screen 352 ILs (22 anions-16 cations combination) and to calculate their related capacity values. Selected ILs (with higher capacity values) are intended to be applied as solvents in EPA extraction. Figure 6 illustrates the following steps for screening, and therefore, the selection of an appropriate IL to be used in the extraction of the total EPA compound.

3.2. Data Collection and Screening Process

The chemical structure of EPA which is illustrated in Table 2 is first optimized by TZVP according to thermodynamic properties and inserted to the data bank of COSMO-RS. Next, cations were selected from imidazolium, pyridinium, pyrrolidinium, piperidinium, and tetra-methyl ammonium based-ILs and hydrophilic and hydrophobic anions. The screened anions and cations based ILs were referred from the lipid extraction literature as a basis, as shown in Table 3 and Table 4. Then, the software estimates the sigma profile and sigma potential of EPA and calculates the related capacity values for each IL toward EPA compound.

3.3. COSMO-RS Assumptions

Solvent extraction process is a robust and economical methodology to be applied for extraction of EPA. There are three assumptions that were made during the screening process including: (a) The liquid state is incompressible, (b) all parts of the molecular surfaces can be in contact with each other, and (c) only pairwise interactions of molecular surface patches are allowed.
With respect to the stated assumptions, screening of charge density and capacity value at infinite dilution were performed by COSMO-RS. The significance of capacity value at infinite dilution is where it defined whether the targeted solvent has adequate strength in extraction. The screening charge density is defined as sigma profile of mixture ( p s ( σ ) ) and it can be calculated by Equations (1) and (2) [23,54,55]:
p s ( σ ) = i = s X i   p X i ( σ )
p X i ( σ ) = n i ( σ ) n i = A i ( σ ) A i
In Equation (1), X i is the mole fraction of compound X, and p X i ( σ ) is the sigma profile of compound X. Additionally, in Equation (2), n i ( σ ) is the number of the segment with the surface charge density σ , A i ( σ ) ) is the total surface area of all segment with particular charge density σ .
The sigma potential ( μ s ( σ ) ) rely on composition, temperature (T), sigma profile ( P s ) ,   and electrostatic interaction energies ( E m i s f i t ) , which is determined by Equation (3):
μ s ( σ ) = R T a e f f l n [ P s ( σ ) e x p ( a e f f R T ( μ s ( σ ) E m i s f i t ( σ , σ ) E H B ( σ , σ ) ) ) d σ ]
Equations (4) represent the electrostatic interaction energies [56].
E m i s f i t ( σ , σ ) = a e f f α 2 ( σ + σ ) 2
where R is ideal gas constant, T is also absolute temperature, a e f f is the effective area of contact, E H B is hydrogen bonding energy, E m i s f i t   is electrostatic interaction energy of two segments per unit area, α′ is an interaction parameter, σ and σ′ are two different screening charge densities for solute molecules that are in contact with each other, and E m i s f i t   is the electrostatic interaction energy of two segments per unit area.
The solvent capacity accounts for the quantity of solute that has been removed from the mixture by the extracting solvent [57]. The solvent capacity at infinite dilution, is numerically the maximum amount of solute (EPA) that can be dissolved in the solvent (IL). The capacity values of ILs for the extraction of EPA calculated by Equation (5).
C = ( 1 γ ) I o n i c   l i q u i d   p h a s e
where, γ is the activity coefficient of solvent at infinite dilution. The value of capacity indicates the amount of solvent that is used during the separation process. The higher solvent capacity results in a lower amount of solvent used. Therefore, the solvent with highest capacity value tends to be more desirable to choose as a separating agent.

4. Conclusions

Overall, the COSMO-RS method was used to screen the suitable ILs for microalgae lipid extraction based on the findings from this study. The highest capacity for EPA extraction among 352 cation/anion combinations was belong to [TMAm][SO4] followed by [BMPIP][SO4]   <   [EMPyr][SO4]   <   [EMIM][SO4] < [EMPyrro][SO4]   <   [TMAm][SO4] by the order of low to high. It is expected to apply this combination ([TMAm][SO4]) to produce a high yield of EPA from microalgae. It is expected to have a better extraction with the small anions that possess high charge density compared to larger ones. Moreover, shorter alkyl chain cations are preferred when using imidazolium-based IL, which agrees with experimental data. Overall, we may conclude that other factors should be addressed to select the IL for extraction such as the rate of mass transfer, the solubility, selectivity, and viscosity. Indeed, screening for the most suitable ILs for extraction reduces the time required and therefore would be more efficient in terms of resources.

Supplementary Materials

The following are available online at https://www.mdpi.com/1420-3049/24/4/713/s1.

Author Contributions

Writing‒original draft preparation, methodology, software, analysis, S.R.M.; supervision, R.H., D.R.A.B., S.A.H., W.A.W.A.K.G. and C.D.W.; writing‒review and editing, R.K. and A.A.M.E.

Funding

This research was funded by “Matching grant” with grant number 9300432.

Acknowledgments

This work has been supported by the Department of Chemical and Environmental Engineering, Universiti Putra Malaysia (UPM) and Centre of Ionic Liquids, Universiti Teknologi Petronas (UTP). The author would like to acknowledge the ionic liquid research centre in Universiti Teknologi Petronas (UTP) to provide the facility to access COSMO-RS software as the main methodology for this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chowdhury, R.; Stevens, S.; Gorman, D.; Pan, A.; Warnakula, S.; Chowdhury, S.; Ward, H.; Johnson, L.; Crowe, F.; Hu, F.B.; et al. Association between fish consumption, long chain omega 3 fatty acids, and risk of cerebrovascular disease: Systematic review and meta-analysis. BMJ 2012, 345, e6698. [Google Scholar] [CrossRef] [PubMed]
  2. Fabian, C.J.; Kimler, B.F.; Hursting, S.D. Omega-3 fatty acids for breast cancer prevention and survivorship. Breast Cancer Res. 2015. [Google Scholar] [CrossRef] [PubMed]
  3. Lazo, M.; Manzano-patin, A.; Cortez-lugo, M.; Romieu, I.; Te, M.M.; Hernandez-Avila, M.; Holguin, F.; Julien, P.; Be, M.C. Omega-3 Fatty Acid Prevents Heart Rate Variability Reductions Associated with Particulate Matter. Am. J. Respir. Crit. Care Med. 2005, 132, 1534–1540. [Google Scholar]
  4. Marik, P.E.; Varon, J. Omega-3 dietary supplements and the risk of cardiovascular events: A systematic review. Clin. Cardiol. 2009, 32, 365–372. [Google Scholar] [CrossRef] [PubMed]
  5. Rubio-Rodríguez, N.; Beltrán, S.; Jaime, I.; de Diego, S.M.; Sanz, M.T.; Carballido, J.R. Production of omega-3 polyunsaturated fatty acid concentrates: A review. Innov. Food Sci. Emerg. Technol. 2010, 11, 1–12. [Google Scholar] [CrossRef]
  6. Ruxton, C.H.S.; Reed, S.C.; Simpson, M.J.A.; Millington, K.J. The health benefits of omega-3 polyunsaturated fatty acids: A review of the evidence. J. Hum. Nutr. Diet. 2004, 17, 449–459. [Google Scholar] [CrossRef]
  7. Nichols, P.D.; McManus, A.; Krail, K.; Sinclair, A.J.; Miller, M. Recent advances in omega-3: Health benefits, Sources, Products and bioavailability. Nutrients 2014, 6, 3727–3733. [Google Scholar] [CrossRef]
  8. Ozogul, Y.; Şimşek, A.; Balikçi, E.; Kenar, M. The effects of extraction methods on the contents of fatty acids, especially EPA and DHA in marine lipids. Int. J. Food Sci. Nutr. 2012, 63, 326–331. [Google Scholar] [CrossRef]
  9. Shahidi, F.; Ambigaipalan, P. Omega-3 Polyunsaturated Fatty Acids and Their Health Benefits. Annu. Rev. Food Sci. Technol. 2018, 9, 345–381. [Google Scholar] [CrossRef]
  10. Adarme-vega, T.C.; Lim, D.K.Y.; Timmins, M.; Vernen, F.; Li, Y.; Schenk, P.M. Microalgal biofactories: A promising approach towards sustainable omega-3 fatty acid production. Microb. Cell Fact. 2012, 11, 96. [Google Scholar] [CrossRef]
  11. Davidson, M.H.; Burns, J.H.; Subbaiah, P.V.; Conn, M.E.; Drennan, K.B. Marine oil capsule therapy for the treatment of hyperlipidemia. Arch. Intern. Med. 1991, 151, 1732–1740. [Google Scholar] [CrossRef] [PubMed]
  12. Swanson, D.; Block, R.; Mousa, S.A. Omega-3 fatty acids EPA and DHA: Health benefits throughout life. Adv. Nutr. 2012, 3, 1–7. [Google Scholar] [CrossRef]
  13. Conchillo, A.; Valencia, I.; Puente, A.; Astiasarán, D.A.I. Componentes Funcionales en Aceites de pescado y de Alga. 2006. Available online: http://www.redalyc.org/articulo.oa?id=309226444013 (accessed on 23 November 2005).
  14. Nelson, W.M. Chapter 3 Are Ionic Liquids Green Solvents? Ionic Liquids as Greener Solvents. ACS Symp. Ser. 2002, 9, 30–41. [Google Scholar]
  15. Berthod, A.; Ruiz-Ángel, M.J.; Carda-Broch, S. ARTICLE IN PRESS G Model Recent advances on ionic liquid uses in separation techniques. J. Chromatogr. A 2017, 1559, 1–15. [Google Scholar]
  16. Shiddiky, M.J.A.; Torriero, A.A.J. Application of ionic liquids in electrochemical sensing systems. Biosens. Bioelectron. 2011, 26, 1775–1787. [Google Scholar] [CrossRef]
  17. Keskin, S.; Kayrak-Talay, D.; Akman, U.; Hortaçsu, Ö. A review of ionic liquids towards supercritical fluid applications. J. Supercrit. Fluids 2007, 43, 150–180. [Google Scholar] [CrossRef]
  18. Biczak, R.; Pawłowska, B.; Bałczewski, P.; Rychter, P. The role of the anion in the toxicity of imidazolium ionic liquids. J. Hazard. Mater. 2014, 274, 181–190. [Google Scholar] [CrossRef] [PubMed]
  19. Pan, J.; Muppaneni, T.; Sun, Y.; Reddy, H.K.; Fu, J.; Lu, X.; Deng, S. Microwave-assisted extraction of lipids from microalgae using an ionic liquid solvent [BMIM][HSO4]. Fuel 2016, 178, 49–55. [Google Scholar] [CrossRef]
  20. Kim, Y.; Choi, Y.; Park, J.; Lee, S.; Yang, Y.; Joo, H.; Park, T.; Hwan, Y.; Hyun, S. Bioresource Technology Ionic liquid-mediated extraction of lipids from algal biomass. Bioresour. Technol. 2012, 109, 312–315. [Google Scholar] [CrossRef]
  21. Eckert, F.; Klamt, A.; Gmbh, C.; Kg, C.; Leverkusen, D. Fast Solvent Screening via Quantum Chemistry: COSMO-RS Approach. AIChE J. 2002, 48, 369–385. [Google Scholar] [CrossRef]
  22. Klamt, A. The COSMO and COSMO-RS solvation models. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2018, 8, e1338. [Google Scholar] [CrossRef]
  23. Klamt, A. Conductor-like Screening Model for Real Solvents: A New Approach to the Quantitative Calculation of Solvation Phenomena. J. Phys. Chem. 1995, 99, 2224–2235. [Google Scholar] [CrossRef]
  24. Rajabi, M.S.; Moniruzzaman, M.; Mahmood, H.; Sivapragasam, M.; Bustam, M.A. Extraction of β-carotene from organic phase using ammonium based ionic liquids. J. Mol. Liq. 2017, 227, 15–20. [Google Scholar] [CrossRef]
  25. Rashid, Z.; Wilfred, C.D.; Gnanasundaram, N.; Arunagiri, A.; Murugesan, T. Screening of ionic liquids as green oilfield solvents for the potential removal of asphaltene from simulated oil: COSMO-RS model approach. J. Mol. Liq. 2018, 255, 492–503. [Google Scholar] [CrossRef]
  26. Xing, H.; Zhang, X.; Yang, Q.; Liu, R.; Bao, Z.; Su, B.; Yang, Y.; Ren, Q. Separation of long chain fatty acids with different number of unsaturated bonds by fractional extraction: Experimental and COSMO-RS study. Food Chem. 2014, 143, 411–417. [Google Scholar] [CrossRef]
  27. Han, J.; Dai, C.; Yu, G.; Lei, Z. Parameterization of COSMO-RS model for ionic liquids. Green Energy Environ. 2018, 3, 247–265. [Google Scholar] [CrossRef]
  28. Durand, M.; Molinier, V.; Kunz, W.; Aubry, J. Classification of Organic Solvents Revisited by Using the COSMO-RS Approach. Chem. Eur. J. 2011, 5155–5164. [Google Scholar] [CrossRef]
  29. Zhou, T.; Qi, Z.; Sundmacher, K. Model-based method for the screening of solvents for chemical reactions. Chem. Eng. Sci. 2014, 115, 177–185. [Google Scholar] [CrossRef]
  30. Hassanali, A.A.; Zhong, D.; Singer, S.J. An AIMD study of CPD repair mechanism in water: Role of solvent in ring splitting. J. Phys. Chem. B 2011, 115, 3860–3871. [Google Scholar] [CrossRef]
  31. Jong, K.; Grisanti, L.; Hassanali, A. Hydrogen bond networks and hydrophobic effects in the amyloid β30–35 chain in water: A molecular dynamics study. J. Chem. Inf. Model. 2017, 57, 1548–1562. [Google Scholar] [CrossRef]
  32. Ufimtsev, I.S.; Luehr, N.; Martinez, T.J. Charge transfer and polarization in solvated proteins from ab initio molecular dynamics. J. Phys. Chem. Lett. 2011, 2, 1789–1793. [Google Scholar] [CrossRef]
  33. van der Vaart, A.; Merz, K.M. The role of polarization and charge transfer in the solvation of biomolecules. J. Am. Chem. Soc. 1999, 121, 9182–9190. [Google Scholar] [CrossRef]
  34. Ansari, N.; Dandekar, R.; Caravati, S.; Sosso, G.C.; Hassanali, A. High and low density patches in simulated liquid water. J. Chem. Phys. 2018, 149, 204507. [Google Scholar] [CrossRef] [PubMed]
  35. Mohan, N.; Suresh, C.H. A Molecular Electrostatic Potential Analysis of Hydrogen, Halogen, and Dihydrogen Bonds. J. Phys. Chem. A 2014, 118, 1697–1705. [Google Scholar] [CrossRef] [PubMed]
  36. Guo, Z.; Lue, B.; Thomasen, K.; Meyer, A.S.; Xu, X. Predictions of flavonoid solubility in ionic liquids by COSMO-RS: Experimental verification, structural elucidation, and solvation characterization. Green Chem. 2007, 9, 1362–1373. [Google Scholar] [CrossRef]
  37. Man, M.S.; Abdullah, M.A.M.; Abdullah, S.B.; Yaacob, Z. Screening Cation and Anion of Ionic Liquid for Dissolution of Silicon Dioxide using COSMO-RS. Indian J. Sci. Technol. 2017, 10. [Google Scholar] [CrossRef]
  38. Hawker, R.R.; Haines, R.S.; Harper, J.B. Variation of the Cation of Ionic Liquids: The Effects on Their Physicochemical Properties and Reaction Outcome. Targets Heterocycl. Syst. Prop. 2015, 18, 141–213. [Google Scholar] [CrossRef]
  39. Shahriari, S.; Neves, C.M.S.S.; Freire, M.G.; Coutinho, J.A.P. Role of the Hofmeister series in the formation of ionic-liquid-based aqueous biphasic systems. J. Phys. Chem. B 2012, 116, 7252–7258. [Google Scholar] [CrossRef]
  40. Zhang, Q.G.; Wang, N.N.; Yu, Z.W. The hydrogen bonding interactions between the ionic liquid 1-Ethyl-3-methylimidazolium ethyl sulfate and water. J. Phys. Chem. B 2010, 114, 4747–4754. [Google Scholar] [CrossRef]
  41. Cheong, L.; Guo, Z.; Yang, Z.; Chua, S.; Xu, X. Extraction and Enrichment of n-3 Polyunsaturated Fatty Acids and Ethyl Esters through Reversible π À π Complexation with Aromatic Rings Containing Ionic Liquids. J. Agric. Food chem. 2011, 59, 8961–8967. [Google Scholar] [CrossRef]
  42. Fu, H.; Li, M.; Ni, R.; Lo, Y.M. Enzymatic catalysis for sustainable production of high omega-3 triglyceride oil using imidazolium-based ionic liquids. Food Sci. Nutr. 2018, 6, 2020–2027. [Google Scholar] [CrossRef] [PubMed]
  43. Li, M.; Pham, P.J.; Wang, T.; Pittman, C.U., Jr.; Li, T. Selective extraction and enrichment of polyunsaturated fatty acid methyl esters from fish oil by novel π-complexing sorbents. Sep. Purif. Technol. 2009, 66, 1–8. [Google Scholar] [CrossRef]
  44. Li, M.; Pittman, C.U., Jr.; Li, T. Extraction of polyunsaturated fatty acid methyl esters by imidazolium-based ionic liquids containing silver tetrafluoroborate—Extraction equilibrium studies. Talanta 2009, 78, 1364–1370. [Google Scholar] [CrossRef] [PubMed]
  45. Orawiec, M. Activity Coe ffi cients at In fi nite Dilution of Organic Solutes and Water in Tributylethylphosphonium Diethylphosphate Using Gas—Liquid Chromatography: Thermodynamic Properties of Mixtures Containing Ionic Liquids. J. Chem. Eng. Data 2016, 61, 1793–1802. [Google Scholar]
  46. Heintz, A.; Sergey, P.; Verevkin, D.O. Thermodynamic Properties of Mixtures Containing Ionic Liquids. 8. Activity Coefficients at Infinite Dilution of Hydrocarbons, Alcohols, Esters, and Aldehydes in 1-Hexyl-3-methylimidazolium Bis(trifluoromethylsulfonyl) Imide Using Gas—Liquid Chromatograp. Chem. Eng 2006, 51, 434–437. [Google Scholar] [CrossRef]
  47. Joshua, J.; Magaret, S.; Pranesh, M.; Chellappan, K.; Cecilia, W.; Mutalib, M.I.A. Separation and Puri fi cation Technology Extractive desulfurization of model fuel oil using ester functionalized imidazolium ionic liquids. Sep. Purif. Technol. 2018, 196, 115–123. [Google Scholar]
  48. Sobota, M.; Vrbka, P. 1-Ethyl-3-methyl-imidazolium Nitrate. J. Phys. Chem. B 2009, 113, 4323–4332. [Google Scholar] [CrossRef]
  49. Marciniak, A.; Letcher, T.M. Activity Coefficients at Infinite Dilution and Physicochemical Properties for Organic Solutes and Water in the Ionic Liquid 1-Ethyl-3-methylimidazolium. J. Solut. Chem. 2014, 44, 413–430. [Google Scholar]
  50. Gabrić, B.; Sander, A.; Bubalo, M.C.; Macut, D. Extraction of S- and N-Compounds from the Mixture of Hydrocarbons by Ionic Liquids as Selective Solvents. Sci. World J. 2013, 2013, 512953. [Google Scholar]
  51. Schafer, A.; Klamt, A.; Sattel, D.; Lohrenz, J.C.W.; Eckert, F. COSMO Implementation in TURBOMOLE: Extension of an efficient quantum chemical code towards liquid systems COSMO-Implementation in TURBOMOLE: Extension of an efficient quantum chemical code towards liquid systems. Phys. Chem. Chem. Phys. 2000, 2, 2187–2193. [Google Scholar] [CrossRef]
  52. Klamt, A.; Jonas, V.; Bu, T.; Lohrenz, J.C.W. Le Refinement and Parametrization of COSMO-RS. J. Phys. Chem. A 1998, 5639, 5074–5085. [Google Scholar] [CrossRef]
  53. Lin, S.; Sandler, S.I. A Priori Phase Equilibrium Prediction from a Segment Contribution Solvation Model. Ind. Eng. Chem. Res. 2002, 41, 899–913. [Google Scholar] [CrossRef]
  54. Mullins, E.; Oldland, R.; Liu, Y.A.; Wang, S.; Sandler, S.I.; Chen, C.C.; Zwolak, M.; Seavey, K.C. Sigma-profile database for using COSMO-based thermodynamic methods. Ind. Eng. Chem. Res. 2006, 45, 4389–4415. [Google Scholar] [CrossRef]
  55. Anantharaj, R.; Banerjee, T. COSMO-RS based predictions for the desulphurization of diesel oil using ionic liquids: Effect of cation and anion combination. Fuel Process. Technol. 2011, 92, 39–52. [Google Scholar] [CrossRef]
  56. Zhou, T.; Chen, L.; Ye, Y.; Chen, L.; Qi, Z.; Freund, H.; Sundmacher, K. An overview of mutual solubility of ionic liquids and water predicted by COSMO-RS. Ind. Eng. Chem. Res. 2012, 51, 6256–6264. [Google Scholar] [CrossRef]
  57. Tumba, K.A. Infinite Dilution Activity Coefficient Measurements of Organic Solutes in Fluorinated Ionic Liquids by Gas-Liquid Chromatography and the Inert Gas Kaniki Armel Tumba. Ph.D. Thesis, Mangosuthu University, Umlazi, South Africa, 2010. [Google Scholar]
Sample Availability: Samples of the compounds are not available from the authors.
Figure 1. (A) EPA chemical structure and (B) sigma surface.
Figure 1. (A) EPA chemical structure and (B) sigma surface.
Molecules 24 00713 g001
Figure 2. The σ-profile of EPA was obtained by COSMO-RS.
Figure 2. The σ-profile of EPA was obtained by COSMO-RS.
Molecules 24 00713 g002
Figure 3. The σ-potential of EPA was obtained by COSMO-RS.
Figure 3. The σ-potential of EPA was obtained by COSMO-RS.
Molecules 24 00713 g003
Figure 4. The capacity values at infinite dilution of (a) Imidazolium, (b) pyridinium, (c) Pyrrololidinium, (d) piperidinium, and (e) tetra-methyl ammonium cations alkyl chain length with 22 anion types for EPA extraction at T = 298.15 K.
Figure 4. The capacity values at infinite dilution of (a) Imidazolium, (b) pyridinium, (c) Pyrrololidinium, (d) piperidinium, and (e) tetra-methyl ammonium cations alkyl chain length with 22 anion types for EPA extraction at T = 298.15 K.
Molecules 24 00713 g004aMolecules 24 00713 g004bMolecules 24 00713 g004c
Figure 5. Proposed hydrogen bonding between the IL, [EMIM]+SO42−, EPA, and water available in the system during extraction.
Figure 5. Proposed hydrogen bonding between the IL, [EMIM]+SO42−, EPA, and water available in the system during extraction.
Molecules 24 00713 g005
Figure 6. COSMO-RS calculation steps.
Figure 6. COSMO-RS calculation steps.
Molecules 24 00713 g006
Table 1. Summarize the selective ILs for EPA extraction.
Table 1. Summarize the selective ILs for EPA extraction.
#Shortlisted ILs Infinite Dilution Activity Coefficient at 298.15 KInfinite Dilution Capacity at 298.15 K
1[TMAm]SO41.00 × 10−119.97 × 1010
2[TMAm]Cl9.52 × 10−101.05 × 109
3[TMAm]Br3.61 × 10−827677690
4[EMPyrro]SO44.40 × 10−72273253
5[EMPyrro]Cl1.68 × 10−559473.56
6[EMIM]SO41.87 × 10−553360.83
7[MPPIP]SO44.17 × 10−523995.28
8[EmPyr]SO47.10 × 10−514084.39
9[BMPyrro]SO49.03 × 10−511073.19
10[BPPIP]SO42.13 × 10−44689.32
11[EMPyrro]Br2.41 × 10−44147.65
12[EMIM]Cl6.61 × 10−41513.22
13[TMAm]propanoate6.97 × 10−41434.47
14[TMAm]NO38.39 × 10−41191.69
15[EMPyrro]propanoate8.67 × 10−41153.16
16[HMPyrro]SO49.62 × 10−41039.60
17[MPPIP]Cl1.01 × 10−3991.30
18[BMIM]SO41.28 × 10−3781.77
19[HPPIP]SO41.28 × 10−3781.45
20[BmPyr]SO41.37 × 10−3727.47
21[EmPyr]Cl1.74 × 10−3575.71
Table 2. Structure of EPA.
Table 2. Structure of EPA.
Shorthand SignSynthetic NameTrivial NameFormulaChemical Structure
20:5Eicosapentanoic acidEPAC20H30O2 Molecules 24 00713 i001
Table 3. The screened anions used in this study.
Table 3. The screened anions used in this study.
No.Anions/AbbreviationsChemical Structure
1Chloride [Cl]Cl
2Bromide [Br]Br
3Tetrafluoroborate [BF4] Molecules 24 00713 i002
4Hexafluorophosphate [PF6] Molecules 24 00713 i003
5Nitrate [NO3] Molecules 24 00713 i004
6Dicyanamide [DCN] Molecules 24 00713 i005
7Diethylphosphate [C4H10O4P] Molecules 24 00713 i006
8Tetrachloro aluminate [AlCl4] Molecules 24 00713 i007
9Methyl sulfate [CH4O4S] Molecules 24 00713 i008
10Thiocyanate [SCN] Molecules 24 00713 i009
11Methane sulfonate [C2H6O3S] Molecules 24 00713 i010
12Ethyl sulfate [C2H6O4S] Molecules 24 00713 i011
13Benzoate [C7H5O2] Molecules 24 00713 i012
14Sulfate [SO4] Molecules 24 00713 i013
15Hydrogen sulphate [HSO4] Molecules 24 00713 i014
16Dimethyl phosphate [C2H7O4P] Molecules 24 00713 i015
17Propanoate [C3H5O2] Molecules 24 00713 i016
18Toluene-4-Sulfonate [C7H7O3S] Molecules 24 00713 i017
19Tri fluoro methane-Sulfonate [CF3SO3] Molecules 24 00713 i018
20Bis(Trifluoromethyl)Imide [NHC2F6] Molecules 24 00713 i019
21Trifluoroacetate [CF3CO2] Molecules 24 00713 i020
22Bis(trifluoromethylsulfonyl) imide [TF2N] Molecules 24 00713 i021
Table 4. The screened cations used in this study.
Table 4. The screened cations used in this study.
No.Cations/AbbreviationsChemical Structures
11-ethyl-3-methyl imidazolium [EMIM] Molecules 24 00713 i022
21-butyl-3-methyl imidazolium [BMIM] Molecules 24 00713 i023
31-hexyl-3-methyl imidazolium [HMIM] Molecules 24 00713 i024
41-octyl-3-methyl imidazolium [OMIM] Molecules 24 00713 i025
51-ethyl-3-methyl pyridinium [EMPyr] Molecules 24 00713 i026
61-butyl-3-methyl pyridinium [BMPyr] Molecules 24 00713 i027
71-hexyl-3-methyl pyridinium [HMPyr] Molecules 24 00713 i028
81-octhyl-3-methyl pyridinium [OMPyr] Molecules 24 00713 i029
91-butyl-1-methyl pyrrolidinium [BMPyrro] Molecules 24 00713 i030
101-ethyl-1-methyl pyrrolidinium [EMPyrro] Molecules 24 00713 i031
111-hexyl-1-methyl pyrrolidinium [HMPyrro] Molecules 24 00713 i032
121-methyl-1-octyl pyrrolidinium [MOPyrro] Molecules 24 00713 i033
131-butyl-1-methyl piperidinium [BMPIP] Molecules 24 00713 i034
141-hexyl-1-methyl piperidinium [HMPIP] Molecules 24 00713 i035
151-methyl-1-propyl piperidinium [MPPIP] Molecules 24 00713 i036
16Tetra-methyl ammonium [TMAm] Molecules 24 00713 i037

Share and Cite

MDPI and ACS Style

Rezaei Motlagh, S.; Harun, R.; Awang Biak, D.R.; Hussain, S.A.; Wan Ab Karim Ghani, W.A.; Khezri, R.; Wilfred, C.D.; Elgharbawy, A.A.M. Screening of Suitable Ionic Liquids as Green Solvents for Extraction of Eicosapentaenoic Acid (EPA) from Microalgae Biomass Using COSMO-RS Model. Molecules 2019, 24, 713. https://doi.org/10.3390/molecules24040713

AMA Style

Rezaei Motlagh S, Harun R, Awang Biak DR, Hussain SA, Wan Ab Karim Ghani WA, Khezri R, Wilfred CD, Elgharbawy AAM. Screening of Suitable Ionic Liquids as Green Solvents for Extraction of Eicosapentaenoic Acid (EPA) from Microalgae Biomass Using COSMO-RS Model. Molecules. 2019; 24(4):713. https://doi.org/10.3390/molecules24040713

Chicago/Turabian Style

Rezaei Motlagh, Shiva, Razif Harun, Dayang Radiah Awang Biak, Siti Aslina Hussain, Wan Azlina Wan Ab Karim Ghani, Ramin Khezri, Cecilia Devi Wilfred, and Amal A. M. Elgharbawy. 2019. "Screening of Suitable Ionic Liquids as Green Solvents for Extraction of Eicosapentaenoic Acid (EPA) from Microalgae Biomass Using COSMO-RS Model" Molecules 24, no. 4: 713. https://doi.org/10.3390/molecules24040713

APA Style

Rezaei Motlagh, S., Harun, R., Awang Biak, D. R., Hussain, S. A., Wan Ab Karim Ghani, W. A., Khezri, R., Wilfred, C. D., & Elgharbawy, A. A. M. (2019). Screening of Suitable Ionic Liquids as Green Solvents for Extraction of Eicosapentaenoic Acid (EPA) from Microalgae Biomass Using COSMO-RS Model. Molecules, 24(4), 713. https://doi.org/10.3390/molecules24040713

Article Metrics

Back to TopTop