- 1Technical Chemistry, Department of Chemistry, Chemical-Biological Centre, Umeå University, Umeå, Sweden
- 2Physical Chemistry Department, University of Tabriz, Tabriz, Iran
- 3Division of Material Science, Luleå University of Technology, Luleå, Sweden
- 4Industrial Chemistry & Reaction Engineering, John Gadolin Process Chemistry Centre, Åbo Akademi University, Turku, Finland
The density, sound velocity, and viscosity of 1-ethyl-3-methylimidazolium chloride [C2mim]Cl in pure water and aqueous solutions of some electrolytes such as potassium chloride, potassium carbonate, and potassium phosphate (weight fraction of salt fixed at ws = 0. 11) have been measured over a wide range of temperatures from 298.15 to 318.15 K. The obtained experimental data have been used to compute various volumetric, compressibility, and viscometric parameters, e.g., apparent molar properties, limiting apparent molar and transfer properties. The co-sphere overlap model was employed to describe the dominant intermolecular interactions in the ternary solutions. Additionally, the structure making/breaking nature of the [C2mim]Cl in the ternary solutions has been discussed in terms of Hepler's constant and the temperature derivative of viscosity B-coefficient (dB/dT). The activation free energy of solvent and solute, activation enthalpy, and activation entropy have been calculated by the application of transition state theory. The calculated parameters have been interpreted in the sense of solvent–solute and solute–solute interactions. The Fourier transform infrared (FTIR) studies also have been done for the studied systems. Volumetric, acoustic, viscometric, and spectroscopic studies can render some evidence and help to understand the aqueous solution behavior of ionic liquids.
Introduction
Ionic liquids (ILs) have been utilized and studied in the last decades in different fields like organic catalysis, biomass treatment/conversion, gas separation, electrolytes in the extraction of diverse substances, and reaction media, due to their inherent advantages and unique properties such as negligible vapor pressure, high thermal and chemical stability, designable ability, and excellent solubility (Wang et al., 2014; Sarmad et al., 2017; Cao et al., 2018; Xu et al., 2018; Hui et al., 2019; An et al., 2020; Wu et al., 2020).
They are considered as potential environmentally friendly solvents to replace the common volatile organic solvents and thus reduce environmental footprints (Gomes et al., 2019; Kaur et al., 2020). In industrial applications, instead of pure compounds, frequently, mixtures of ILs with different solvents are used, and therefore, investigation of the behavior of binary or ternary mixtures of ILs is of great importance to unravel their structural characteristics which can provide valuable information regarding solute–solute and solute–solvent interactions (Zafarani-Moattar and Sarmad, 2012).
Imidazolium-based ILs have been gaining more attention in the field of materials, physics, and chemistry as pure solvents or cosolvents in aqueous/non-aqueous systems and biphasic mixtures (Wang et al., 2007; Zafarani-Moattar and Sarmad, 2011). Most of the studies have been related to the synthesis of ILs and their chemical reactions. Indeed, investigation of the physicochemical properties of the systems containing ILs in a wide range of pressures and temperatures is in demand. These properties are essential for optimizing and designing the systems containing IL and solvent in both laboratory and industrial scales. Additionally, for the evaluation of molecular interactions, thermophysical properties and their deviations from ideality are important (Pal and Gaba, 2008).
Kaur et al. studied volumetric and acoustic properties of choline acetate IL in α, ω-alkane diols, at different temperatures (Kaur et al., 2020). Volumetric, acoustic, and spectroscopic properties of an aqueous ternary mixture of IL 1-hexyl-3-methylimidazolium bromide and amino acids (L-alanine and L-phenylalanine) were investigated by Kumar and Sharma (2020). Volumetric, acoustic, and infrared spectroscopic study of amino acids in aqueous solutions of pyrrolidinium-based IL was also carried out by Kumar et al. (2020). Sharma et al. (2020) studied molecular interactions of L-histidine in aqueous imidazolium-based IL solutions using volumetric, acoustic, and viscometric approaches, at different temperatures. On the other hand, Zafarani-Moattar et al. (2020) investigated the volumetric, acoustic, and viscometric properties of some choline amino acid ILs in aqueous polypropylene glycol 400 and polyethylene glycol 400 solutions.
Marcinkowski et al. (2019) studied interactions of N-alkyl-N-methylmorpholinium-based ILs with acetonitrile using molecular dynamic simulations along with density and velocity of sound measurements.
The volumetric and acoustic behaviors of a protic IL, 2-hydroxy ethylammonium lactate, was investigated in different hydroxylic media, e.g., water, methanol, and ethanol, at various temperatures (288.15–323.15 K) and under atmospheric pressure by Barros et al. (2018). The influence of temperature, alkyl chain length, and nature of the anion on the volumetric and acoustic properties of imidazolium-based ILs with alkyl nitriles was studied as well (Bhanuprakash et al., 2020).
Interactions of 1-butyl-1-methylpyrrolidinium tetrafluoroborate and 1-butyl-1-methyl piperidinium tetrafluoroborate with acetonitrile and water were investigated by density and speed of sound measurements comparatively and at a different temperature, by Sahin and Ayranci (2019). Gaba et al. (2019) studied the molecular interactions of L-glutamic acid and L-aspartic acid in aqueous imidazolium-based IL solutions, at different temperatures, by measuring their volumetric and acoustic properties. Furthermore, the apparent molar volume, apparent molar isentropic compressibility, and corresponding limiting values were computed using the experimental density and speed of sound data. These quantities are helpful when investigating the solvation behavior and various interactions present in the ternary amino acid + IL + water solutions.
To understand the intermolecular interaction existing in the ternary systems of citrate salts + water + 1-butyl-3-methyl imidazolium tetrafluoroborate, Kumar and Katal (2018b) studied the volumetric and compressibility properties at 0.1 MPa and different temperatures. The structure-making ability of citrate salts in IL solution is observed.
Zafarani and Sarmad investigated the effect of tripotassium phosphate, dipotassium hydrogen phosphate, and potassium dihydrogen phosphate on aqueous solution behavior of 1-ethyl-3-methylimidazolium bromide at different temperatures. The negative dB/dT values of IL in water and aqueous electrolyte solutions revealed that IL acts as a structure maker in water and in aqueous electrolyte solutions (Zafarani-Moattar and Sarmad, 2010).
Volumetric, acoustic, and transport properties of IL, 1-butyl-3-methyl imidazolium chloride in aqueous lithium bromide solutions at different temperatures (298.15–318.15 K) were studied as well. The calculated volumetric interaction parameters indicated the existence of hydrophilic–ionic or ionic–ionic interactions in the studied systems and suggested that the interactions between the IL and lithium bromide are mainly pairwise (Zafarani-Moattar et al., 2014).
The interactions of some phosphate and citrate salts with 1-hexyl-3-methylimidazolium chloride as a function of temperature have been investigated by a combination of volumetric and acoustic measurements by Kumar and coworkers. The absorption values revealed the existence of interactions between citrate/phosphate salt, IL, and water molecules, which supported results obtained from volumetric and acoustic data (Kumar and Chadha, 2014, 2015).
Volumetric and compressibility behaviors of binary and ternary aqueous solutions of 1-hexyl-3-methylimidazolium chloride in methyl potassium malonate and methyl potassium malonate were studied by Sadeghi and Mahdavi (2012); the obtained results have been interpreted in terms of the solute–solute and solute–water interactions.
Ternary systems composed of phosphate salts + 1-butyl-3-methyl imidazolium tetrafluoroborate + water and citrate salts + 1-butyl-3-methyl imidazolium tetrafluoroborate + water at 0.1 MPa and different temperature have been investigated by Kumar and Katal (2018a,c). The pair and triplet interaction coefficients have been computed from transfer parameters and have demonstrated the supremacy of pairwise interactions in the system. The possible intermolecular interactions such as hydrophilic–hydrophobic, hydrophilic–hydrophilic, or hydrophobic–hydrophobic were investigated as well.
In the present study, density, speed of sound, and viscosity of 1-ethyl-3-methyl imidazolium chloride [C2mim]Cl in pure water as well as in an aqueous solution of potassium chloride, potassium carbonate, and potassium phosphate (ws = 0.11) at T = 298.15–318.15 K were measured. The apparent molar volume, isentropic compressibility, apparent molar isentropic compressibility, and relative viscosity were computed using the experimental data. Transfer apparent molar volume, apparent isentropic compressibility, and viscosity B-coefficient of IL from pure water to aqueous electrolyte solutions were also calculated which provides valuable information regarding the solute–solute and solute–solvent interactions that exist in aqueous IL + electrolyte solutions. The structure-making/breaking nature of the [C2mim]Cl in the ternary solutions has been discussed in terms of Hepler's constant and temperature derivative of viscosity B-coefficient (dB/dT). The activation free energy, enthalpy, and entropy were also calculated by the application of transition state theory. The calculated parameters were interpreted in the sense of solvent–solute and solute–solute interactions. The Fourier transform infrared (FTIR) studies and thermal analysis were also performed for the studied systems. The temperature dependency of viscosity was successfully fitted to the Vogel–Fulcher–Tammann (VFT) equation.
Materials and Methods
All chemicals (electrolytes and IL) were obtained from Merck with a purity >98 (wt%). The electrolytes were used as received without further purification, but they were dried at 100–110°C overnight. The water content in the IL was determined using a microprocessor-based automatic Karl–Fischer titrator. The water content found for [C2mim]Cl was 1.16% by mass. This water content in the IL was taken into account during the preparation of the aqueous solutions of IL. For the preparation of solutions, double-distilled and deionized water was used.
All the solutions were prepared afresh and scaled using an analytical balance (Shimadzu, 321-34553, Shimadzu Co., Japan) with precision (10−7 kg). All the solutions were kept tightly sealed to minimize the absorption of atmospheric moisture. Measurements were performed immediately after the preparation of solutions.
Density Measurements
The density and sound velocity of mixtures were measured using a commercial density and sound velocity measurement apparatus (Anton Paar DSA 5000 densimeter and sound velocity analyzer). Both of the speed of sound and density are extremely sensitive to temperature, and therefore, the temperature was kept constant within ±10−3 K. The experimental uncertainties of density and ultrasonic velocity measurements were estimated to be in the range of ±3.0 ± 10−6 g·cm−3 and ±0.1 m·s−1, respectively. The apparatus was calibrated with double-distilled deionized and degassed water and dry air at atmospheric pressure.
Viscosity Measurements
The viscosity measurements were carried out with a suspended level Ubbelohde-type capillary viscometer, mounted in a water thermostat (Julabo, MD-18V, Germany) which was calibrated with water at five different temperatures (298.15–318.15 K). The flow time of a constant volume of liquid through a capillary was measured with an electronic stopwatch with a resolution of 0.01 s. An average of at least four readings of flow time with a variation not exceeding ±0.1 s was calculated for each solution. The temperature of the thermostat bath was controlled within ±0.01 K. The measured viscosities were found to be accurate within ±0.001 mPa. Also, the dynamic viscosity, η, was calculated by applying the following relation:
here, t, K, θ, and d are referred to as the flow time, the viscometer constant, the Hagenbach correction factor, and the density, respectively. The viscometer constant was determined by calibrating with distilled water at working temperatures (Marsh, 1987). The uncertainty for the dynamic viscosity determination was estimated to be around ±0.5%.
Thermogravimetric Analysis (TGA)
STA 449C Jupiter (NETZSCH, Germany) was employed to perform the TGA. Five to 10 mg of samples was heated under argon atmosphere from 20 to 500°C, at a heating rate of 10°C·min−1. NETZSCH TA software was used for the processing of the obtained data.
Fourier Transform Infrared (FTIR) Analysis
The attenuated total reflectance Fourier transforms infrared spectroscopy (ATR-FTIR) was utilized to characterize the [C2mim]Cl in pure water and aqueous electrolyte solutions. The Bruker Vertex 80v FTIR spectrometer with a DTGS detector was applied. The FTIR analysis was performed at room temperature for all the studied systems.
Results and Discussions
Volumetric Properties
The experimental density and sound velocity for binary aqueous [C2mim]Cl solution along with ternary [C2mim]Cl + aqueous electrolyte solution at T = 298.15–318.15 K are reported in Supplementary Table 1. Apparent molar volumes, Vϕ, of [C2mim]Cl in pure water and aqueous electrolyte solutions (ws = 0.11), at different temperatures, were determined using the experimental density data via the following equation:
where M and m are the molar mass (kg·mol−1) and the molality (mol·kg−1) of [C2mim]Cl, respectively; ρ and ρ0 are the densities (kg·m−3) of ternary solutions ([C2mim]Cl + H2O + electrolyte) and the solvent (H2O + electrolyte), respectively. The apparent molar volume in the aqueous electrolyte solutions is larger than that of pure water and increases by increasing the temperature. This means that the size of the solute [C2mim]Cl in aqueous electrolyte solutions is larger than in pure water. This can be attributed to the phenomena that in aqueous solutions, [C2mim]Cl dissociated to anion and cation, whereas in the presence of potassium salts, they formed ion pairs with the ions of the salts which are physically bonded together and, therefore, the electrostatic interactions between water molecules and the ions decrease. In the studied electrolyte solutions, the apparent molar volume decreased by increasing the charge of anions as follows: KCl > K2CO3 > K3PO4.
The limiting apparent molar volumes, at infinite dilution (), were estimated by applying the Redlich–Mayer-type equation to the corresponding data via Equation (3) (Redlich and Meyer, 1964).
The values along with adjustable parameters (SV, bV) and standard deviations are summarized in Table 1. The magnitude of values increases with an increase in temperature as can be depicted in Figure 1, indicating that the hydration effects in solutions are strongly sensitive to temperature. The limiting apparent molar volumes render beneficial information regarding solvent–solute interactions. The reason is that at infinite dilution, each ion is surrounded just by solvent molecules and is distant with other ions. Consequently, the limiting apparent molar volume is a measure of ion–solvent interaction and uninfluenced by the ion–ion interactions. The calculated limiting apparent molar volumes were used to determine the limiting molar volumes of transfer of IL from water to aqueous electrolyte solutions (Rajagopal and Jayabalakrishnan, 2009; Zafarani-Moattar and Sarmad, 2010; Kaczkowska et al., 2020).
Table 1. Limiting apparent molar volumes along with the coefficients of Equation (3) for the studied systems.
The obtained values are reported in Table 2 and depicted in Figure 2. Evidently, as can be seen for all studied systems, the values are negative and generally decrease with increasing temperature. The negative values suggest that the hydrophobic–ionic interactions predominate over the hydrophilic–ionic interactions.
Table 2. Limiting molar volume, limiting molar isentropic compressibility, and B-coefficient of transfer of [C2mim]Cl from pure water to aqueous electrolyte solutions (ws = 0.11).
Figure 2. Limiting apparent molar volumes of transfer of [C2mim]Cl from pure water to aqueous electrolyte solutions.
In the studied ternary system ([C2mim]Cl + H2O + electrolyte), three types of interactions possibly exist between the solvent and solute: (1) ion–ion interactions among potassium ions of electrolytes and chloride ions of IL, (2) ion–ion interactions among imidazolium ion of IL and phosphate, carbonate, and chloride ions of electrolytes, (3) ion–apolar group interactions. In the systems in which type 3 interactions are dominant, the limiting molar volumes of transfer are negative. According to the co-sphere overlap model, the overlap of two co-spheres of hydrophobic and ionic hydration releases some water molecule from the solvation sphere to the bulk that leads to a negative volume contribution (Scheme 1) (Friedman and Krishnan, 1973; Lin et al., 2006).
Scheme 1. The interaction of two co-spheres (Lin et al., 2006).
To identify the hydrophilic or hydrophobic features of [C2mim]Cl in electrolyte solutions, the temperature dependency of the obtained values was correlated via the following quadratic equation (Rajagopal and Jayabalakrishnan, 2009):
here a, b, and c are empirical coefficients that have been determined so Equation (5) has the following form for the studied systems.
According to Hepler's theory, Hepler's constant is related to the hydrophilic or hydrophobic feature of the solute and is a better proof in describing whether a solute acts as a breaker or maker of a certain structure (Helper, 1969). If the sign of Hepler's constant is positive, it means that the solute is hydrophobic and is behaving as making (developing) the structure. Vice versa, the negative sign shows that the solute is hydrophilic and acts as a breaker of the structure. The obtained values for all the studied systems are positive (Equation 6) and revealed that [C2mim]Cl in pure water and electrolyte solutions is a structure maker.
Differentiating Equation 5 to temperature results in the limiting apparent molar expansibility ().
The obtained values of for studied systems are reported in Table 3. At all temperatures, values of IL in water and electrolyte solutions are positive and increased with temperature due to the increase of volume. This behavior can be attributed to the solvation and shrinkage of the solvent around the ions in aqueous electrolyte solutions. Some of the water molecules upon heating the solution may be released from the hydration layer, which can increase the volume of solution slightly faster than that of pure water, so will be positive. Under ambient conditions, the expansivity of water is around 0.25 × 10−3 K−1, which is less than that in aqueous IL solutions (Gu and Brennecke, 2002).
Table 3. Limiting apparent molar expansibility and isobaric thermal expansion coefficients for studied systems ([C2mim]Cl + H2O + electrolyte; ws = 0.11) at different temperatures.
The limiting apparent molar expansibility along with limiting apparent molar volume was used to determine the isobaric thermal expansion coefficient (αP) through Equation (8).
The obtained values of αP are listed in Table 3. The value of αP in water was higher than that in electrolyte solutions.
Acoustic Properties
The isentropic compressibility was determined using the experimental density and speed of sound by applying Laplace–Newton's equation.
For all the studied systems, the isentropic compressibility decreased by the molality of IL, at all studied temperatures, and the reason may be hydration of ions and breaking up of the three-dimensional network of water molecules. The compressibility isotherms for the system [C2mim]Cl + H2O + K3PO4, as an example, is demonstrated in Figure 3. As can be seen, the isotherms intersect at mIL = 0.921 mol·kg−1; at this point, the isentropic compressibility was equal and temperature independent. The same behavior has been observed in the electrolyte and non-electrolyte solutions (Schmelzer et al., 2004; Wahab and Mahiuddin, 2004).
Figure 3. Isentropic compressibility of [C2mim]Cl + H2O + K3PO4 as a function of [C2mim]Cl molality at different temperatures.
According to Equation (10), the isentropic compressibility of electrolyte solutions is the sum of two contributions (Wahab and Mahiuddin, 2004):
We can see (Figure 3) that up to the converging point, κs (solvent intrinsic) is the main contributor to the total value of isentropic compressibility, and after that, the solute intrinsic part is the major contributor. At the coverage point, dκs/dT = 0, and it seems that at this point the inherent structure of water is totally broken down, and therefore, the water molecules are assimilated in the ions' primary hydration shell. Above the coverage point, dκs/dT > 0, so the contribution of the solvent intrinsic properties is almost naught, in contrast to the case when we are below the coverage point, dκs/dT < 0. Indeed, the presence of an intersection point can be attributed to the formation of a clathrate-like structure in an aqueous IL solution similar to tetraethyl ammonium salts, which exhibited a similar behavior (Jerie et al., 1999).
The apparent molar isentropic compressibility of the studied systems was computed by employing Equation (11).
Here, κs0 and κs, respectively, represent the isentropic compressibility of a pure solvent and a mixture. Since the concentration dependency of the apparent molar isentropic compressibility was non-linear, the obtained values were fitted to the Redlich–Mayer-type equation as below:
The adjustable parameters of the above equation along with limiting apparent isentropic compressibility, at infinite dilution (), are tabulated in Table 4. The obtained values for [C2mim]Cl in aqueous solutions, at all temperatures, were negative and smaller than those obtained for the electrolyte solutions. In the studied electrolyte solutions, the values at lower temperatures are negative, which means that the water molecules around IL presented higher resistances to compression when compared with bulk. Also, the values of increased by increasing the temperature. The negative values at lower temperatures can be attributed to the attractive interactions between the IL and solvent due to the hydration of the ions. The values of are affected by the size of the molecules and the degree of penetration of solvent molecules. The size of the molecules has a positive contribution to , and thereby, large molecules due to the free intermolecular spaces possess inherent compressibility more so than in the case of a mixture. The penetration of the solvent molecules has a negative contribution; this can result in the formation of intramolecular free spaces because the IL molecules interacted with the solvent molecules. Upon increasing temperature, the electrostriction was reduced, and the interactions of IL with solvent were weakened. Thus, the solvent molecules around the IL are released into the bulk; and therefore, the medium becomes more compressible, and the value of is increased.
Table 4. Limiting apparent molar isentropic compressibility along with the coefficients of Equation (6) for the studied systems at different temperatures.
The limiting molar compressions of transfer () from water to aqueous electrolyte solutions have been determined using the following relation:
The obtained values are listed in Table 2, and for all studied systems, the values of are positive and decrease with increasing temperature.
Viscosity
The relative viscosities of [C2mim]Cl in aqueous and electrolyte solutions were determined using the experimental viscosity data (Supplementary Table 2) and employing the following equation.
Herein, η and η0 are the viscosities of the solution and solvent, respectively. The viscosity B-coefficients were obtained by fitting the relative viscosity (in the linear region) to the Jones–Dole equation by a least-squares method (Jones and Dole, 1929).
In the above equation, C is the molarity of solutions (mol·dm−3). The obtained B-coefficients are presented in Table 5. For all the studied systems, B-coefficients are positive and decreased by increasing temperature as well as decreasing the charge of the anion (KCl < K2CO3 < K3PO4).
The viscosity B-coefficients are of great importance in the investigation of the transport properties of solutions. They render valuable information regarding the solvation of the solute and their influence in the structure of the solvent surrounding the solute molecules. Additionally, some activation parameters can be determined using B-coefficients. The viscosity B-coefficients depend on the size of the solvent and solute molecules as well as the interactions between the solvent–solute (Jenkins and Marcus, 1995; Rajagopal and Jayabalakrishnan, 2009). Positive and large B-coefficient values demonstrate a structure-making phenomenon, thus rendering the hydrophobicity of the solute molecules and giving rise to hydrogen-bonded action of solute on solvent molecules. Table 5 illustrates the structure-making ability of [C2mim]Cl in both aqueous and electrolyte solutions, and the values reflect the presence of strong ion–solvent interactions. The temperature derivatives of the B-coefficients (dB/dT) of the studied systems are negative since their values decrease by increasing temperature. The sign of dB/dT also provides some information regarding the structure-making/breaking ability of a solute in a solvent media (Sarma and Ahluwalia, 1973). For [C2mim]Cl in aqueous and electrolyte solutions, the values of dB/dT are negative and indicate the structure-making feature of [C2mim]Cl. These results are in excellent agreement with those concluded from Hepler's constant.
The viscosity B-coefficients were used to calculate ΔtB transfer as follows:
The reported values of ΔtB in Table 2 and Figure 4, for all cases, are positive and increase by increasing the temperature as well as decreasing the charge of electrolytes' anion. The positive values of ΔtB represented the more structured medium in the presence of aqueous electrolyte solutions.
Figure 4. Viscosity B-coefficient of transfer of [C2mim]Cl from pure water to aqueous electrolyte solutions.
The viscosity B-coefficients were also used to evaluate the free energy of activation per moles of the solvent and solute by employing the transition state theory suggested by Feakins et al. (1993) and Glasstone et al. (1941):
Here, and are the limiting apparent molar volumes of the solute and solvent, respectively. and denote the activation free energies of the solute and solvent, respectively. h is Planck's constant, NA is Avogadro's number, η0 is the viscosity of the solvent, and R is the gas constant. The activation entropy and enthalpy for different systems could be estimated using the following relations:
The calculated values of activation free energies, entropies, and enthalpies are listed in Table 6. According to the reported values in Table 6, values are positive and also larger than for all the studied systems, which signify the structure-maker ability of [C2mim]Cl in aqueous electrolyte solutions and are therefore supporting our conclusions through Hepler's constant and dB/dT. Larger values demonstrated the presence of stronger ion–solvent interactions as well. Indeed, the formation of a transition state occurred, followed by the deformation of the intermolecular forces in the solvent structure. Similar results were reported for the amino acids in aqueous salbutamol sulfate solution and glycine in aqueous transition metal chloride solutions (Mishra and Gautam, 2001; Rajagopal and Jayabalakrishnan, 2009).
Table 6. The activation free energy of solvent, solute, entropy, and enthalpy for the studied systems at different temperatures.
Based on the transition state theory, each solvent molecule in 1 mol of the solution must pass across the transition state and interact with solute molecules (more/less strongly) (Feakins et al., 1993). The activation free energy () incorporates the free energy of transfer of the solute from base state to the transition state along with the free energy of solute across its viscous transition state . is calculated using Equation (22) (Feakins et al., 1993; Yan et al., 2002).
The obtained values for can be found in Table 7, and they are for all studied systems positive and larger than . This realization revealed that for the formation of a transition state, more solute–solvent bonds should be broken; in other words, the formation of the transition state is less favored (Mishra and Gautam, 2001; Yan et al., 2002; Rajagopal and Jayabalakrishnan, 2009).
Table 7. Thermodynamic activation parameter transfer of [C2mim]Cl from water to aqueous electrolyte solutions (ws = 0.11).
According to Table 6, values are positive for all systems, which means that the achievement of the transition state for viscous flow has occurred along with bond breaking and order decrease. These results are in good agreement with conclusions drawn from apparent molar volume and viscosity, so the results confirm each other.
The temperature dependency of the viscosity was fitted with the Vogel–Fulcher–Tammann (VFT) equation, which is capable of describing the relation between temperature and dynamic viscosity, at constant molar composition, and can be expressed as follows (Verdía et al., 2014; González et al., 2015; Yang et al., 2017; Cai et al., 2020):
A, B, and T0, which are adjustable parameters along with the correlation coefficients (R2), are given in Supplementary Table 3. As depicted in Figure 5, the viscosity of a solution, at a fixed molality, decreased with increasing temperature. Since the viscosity of a liquid is a measure of the cohesion of the molecules, it is sensitive to temperature. Consequently, the higher the temperature, the stronger the molecules' vibrations, the smaller the molecules' cohesion, and consequently, the lower the viscosity. The viscosity of all solutions increased by increasing molality, due to the smaller distances between anions and cations, so electrostatic interactions became stronger and resulted in increased viscosity.
Figure 5. The viscosity of [C2mim]Cl + H2O + K3PO4 as a function of [C2mim]Cl molality at different temperatures.
FTIR Analysis
For the structure determination of all kinds of compounds including inorganic or organic, FTIR spectroscopy is the most important characterization technique and can be employed. This technique can be applied for both a single compound and a mixture of compounds in a powdered, solid, liquid, or paste state. The FTIR method supplies insight into the molecular interactions that take place in the systems, based on the vibrations of the atoms of a molecule. The FTIR apparatus sends the infrared radiations toward the sample, and some of them pass through the sample whereas the other part is absorbed by the molecules of the sample. The absorbed radiations are transformed into vibrational or rotational energy by the molecules of the sample. The resulting spectra represent the fingerprint of the molecule. The FTIR spectra of [C2mim]Cl + H2O and [C2mim]Cl + H2O + electrolytes were recorded in the wavenumber (4.000–400 cm−1), and the values are depicted in Figure 6.
The aqueous solution of [C2mim]Cl gives rise to major peaks at 1.165 cm−1, corresponding to the stretching vibration of the amine group (C–N), 1.350–1.454 cm−1 bending of the C–H, 1.575 cm−1 stretching vibration of C=C, and 1.633 cm−1 stretching vibration of N=C. The broad peak corresponding to the O–H stretching band for water appears at 3.356 cm−1 (Seethalakshmi et al., 2012).
As shown in Figure 6, the spectrum of [C2mim]Cl + H2O + K3PO4 shows a band around 991 cm−1, which is related to the asymmetric stretching P–O in aqueous ions besides the earlier-mentioned peaks for [C2mim]Cl (Jastrzebski et al., 2011; Zhang et al., 2018). Furthermore, the spectrum of [C2mim]Cl + H2O + K2CO3 comprises characteristic peaks for [C2mim]Cl as well as a peak at 1,381 cm−1 that belongs to the asymmetric stretching of CO2 in a carbonate structure (Toops et al., 2005; Liu et al., 2017).
According to Figure 6, there is no significant difference between the spectra of aqueous [C2mim]Cl and [C2mim]Cl + H2O + KCl due to the dissociation of KCl in the aqueous medium, which does not show any absorbance to the infrared radiations (Max and Chapados, 1999).
The broad O–H peaks in all the studied samples arise from the water molecules in the solution around [C2mim]Cl and also K3PO4, K2CO3, and KCl molecules. However, a small change in FTIR wavenumber of O–H peak in different systems may be interpreted by some structural change or possible interactions between [C2mim]Cl and electrolyte molecules. A similar behavior was reported recently for the IL + water + amino acid systems (Kumar and Sharma, 2020; Kumar et al., 2020).
Thermal Properties
The thermal behavior of [C2mim]Cl in water and aqueous electrolyte solutions was also investigated, and the results are presented in Figure 7. According to TGA–DTA diagrams, the studied systems undergo a two-step thermal decomposition process (Efimova et al., 2015). The first step occurred around 70–80°C, which can be attributed to the dehydration of the samples. The decomposition of [C2mim]Cl takes place at 286.7°C, and the addition of different electrolytes shifted the decomposition of [C2mim]Cl to lower temperatures 271.7, 274.2, and 276.7, respectively, for K3PO4, K2CO3, and KCl. Therefore, increasing the charge of anions resulted in the decomposition of IL at lower temperatures and decreased the thermal stability of IL more.
Conclusions
The density, speed of sound, and viscosity of [C2mim]Cl in water and aqueous KCl, K2CO3, and K3PO4 (ws = 0.11) solutions were measured at T = 298.15–318.15 K. From the experimental data, the apparent molar volume, apparent isentropic compressibility, and relative viscosity have been determined. The limiting apparent molar volume and isentropic compressibility along with viscosity B-coefficient were calculated as well, which were utilized to evaluate the corresponding transfer parameters for [C2mim]Cl from water to aqueous electrolyte solution. The obtained result exhibited a negative transfer volume of [C2mim]Cl from water to the aqueous electrolyte solutions, which decreased by increasing temperature. The isentropic compressibility of all studied systems decreased by increasing the concentration of [C2mim]Cl and temperature. The intersection point in isentropic compressibility isotherms may indicate the formation of a clathrate-like structure. Hepler's constant and dB/dT revealed that [C2mim]Cl in pure water and electrolyte solutions acts as a structure maker. The activation free energy and activation enthalpy values were also calculated, which confirm the structure maker behavior of [C2mim]Cl in aqueous electrolyte solutions. The FTIR studies revealed that a change in the wavenumber of the O–H peak in different systems might indicate some structural changes or possible interactions between [C2mim]Cl and electrolyte molecules. The temperature dependency of viscosity was satisfactorily fitted to the VFT equation. The viscosity of [C2mim]Cl in aqueous electrolyte solutions is larger than that in pure water and increases with increasing IL concentration and decreases with temperature. The TGA revealed that all the studied systems undergo a two-step thermal decomposition process. By the addition of different electrolytes to the aqueous IL solution, the thermal decomposition temperature of IL was shifted to lower temperatures. The higher the anion charge, the lower the decomposition temperature.
Data Availability Statement
All datasets generated for this study are included in the article/Supplementary Material.
Author Contributions
SS conducted and organized all the experiments, calculated and analyzed the data, and wrote the manuscript. MZ-M and J-PM supervised the experiments and analysis and edited and reviewed the manuscript. DN performed thermal analysis, analyzed data, and edited and reviewed the manuscript. All the authors contributed to the manuscript and approved the submitted manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The University of Tabriz Research Council, Bio4Energy program, and ÅForsk are gratefully acknowledged. This work is a part of the activities of the Technical Chemistry, Department of Chemistry, Chemical Biological Centre, Umeå University, Sweden, and the Johan Gadolin Process Chemistry Centre at Åbo Akademi University in Finland.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2020.593786/full#supplementary-material
References
An, X., Li, Z., Zhou, Y., Zhu, W., and Tao, D. (2020). Rapid capture and efficient removal of low-concentration SO2 in simulated flue gas by hypercrosslinked hollow nanotube ionic polymers. Chem. Eng. J. 394:124859. doi: 10.1016/j.cej.2020.124859
Barros, S., Andrade, R. S., and Iglesias, M. (2018). Effect of temperature on thermodynamic properties of protic ionic liquids: 2-hydroxy ethylammonium lactate (2-HEAL) short hydroxylic solvent. Int. J. Thermodyn. 21, 70–80. doi: 10.5541/ijot.303607
Bhanuprakash, P., Prathibha, R., Gardas, R. L., Sivakumar, K., and Jyothi, N. (2020). Effect of temperature, nature of anion and alkyl chain length on the volumetric and acoustic properties of ionic liquid [C4C1im][MeSO4] with alkyl nitriles. J. Mol. Liq. 302:112507. doi: 10.1016/j.molliq.2020.112507
Cai, G., Yang, S., Wang, X., Zhou, Q., Xu, J., and Lu, X. (2020). Densities and viscosities of binary mixtures containing the polyhydric protic ionic liquid (2-hydroxy-N-(2-hydroxyethyl)-N-methylethanaminium methanesulfonate) and water or alcohols. J. Sol. Chem. 49, 423–457. doi: 10.1007/s10953-020-00968-2
Cao, Z., Zhao, X., He, F., Zhou, Y., Huang, K., Zheng, A., et al. (2018). Highly efficient indirect hydration of olefins to alcohols using superacidic polyoxometalate-based ionic hybrids catalysts. Ind. Eng. Chem. Res. 57, 6654–6663. doi: 10.1021/acs.iecr.8b00535
Efimova, A., Pfützner, L., and Schmidt, P. (2015). Thermal stability and decomposition mechanism of 1-ethyl-3-methylimidazolium halides. Thermochim. Acta 604, 129–136. doi: 10.1016/j.tca.2015.02.001
Feakins, D., Bates, F. M., Waghorne, W. E., and Lawrence, K. G. (1993). Relative viscosities and quasi-thermodynamics of solutions of tert-butyl alcohol in the methanol–water system: a different view of the alkyl–water interaction. J. Chem. Soc. Faraday Trans. 89, 3381–3388. doi: 10.1039/FT9938903381
Friedman, H., and Krishnan, C. (1973). “Thermodynamics of ionic hydration,” in Aqueous Solutions of Simple ElectrolytesAnonymous, ed F. Franks (Boston, MA: Springer), 1–118.
Gaba, R., Pal, A., Sharma, D., Kumar, H., and Kumar, A. (2019). Molecular interactions of some non-essential amino acids in aqueous solutions of 1-methylimidazolium chloride at different temperatures. J. Mol. Liq. 279, 711–718. doi: 10.1016/j.molliq.2019.01.094
Glasstone, S., Laidler, K. J., and Eyring, H. (1941). The theory of rate processes. Nature 149, 509–510. doi: 10.1038/149509a0
Gomes, J. M., Silva, S. S., and Reis, R. L. (2019). Biocompatible ionic liquids: Fundamental behaviours and applications. Chem. Soc. Rev. 48, 4317–4335. doi: 10.1039/C9CS00016J
González, E. J., González, B., and Macedo, E. A. (2015). Effect of the relative humidity and isomeric structure on the physical properties of pyridinium based-ionic liquids. J. Chem. Thermodyn. 86, 96–105. doi: 10.1016/j.jct.2015.02.020
Gu, Z., and Brennecke, J. F. (2002). Volume expansivities and isothermal compressibilities of imidazolium and pyridinium-based ionic liquids. J. Chem. Eng. Data 47, 339–345. doi: 10.1021/je010242u
Helper, L. (1969). Studies on viscosities and densities of R4NX in ME water mixtures of different temperatures. Can. J. Chem. 47, 4613–4616.
Hui, W., Zhou, Y., Dong, Y., Cao, Z., He, F., Cai, M., et al. (2019). Efficient hydrolysis of hemicellulose to furfural by novel superacid SO4H-functionalized ionic liquids. Green Energy Environ. 4, 49–55. doi: 10.1016/j.gee.2018.06.002
Jastrzebski, W., Sitarz, M., Rokita, M., and Bulat, K. (2011). Infrared spectroscopy of different phosphates structures. Spectrochim. Acta A Mol. Biomol. Spectrosc. 79, 722–727. doi: 10.1016/j.saa.2010.08.044
Jenkins, H. D. B., and Marcus, Y. (1995). Viscosity B-coefficients of ions in solution. Chem. Rev. 95, 2695–2724. doi: 10.1021/cr00040a004
Jerie, K., Baranowski, A., Jákli, G., and Gliński, J. (1999). Structure of aqueous solutions of tetraethylammonium chloride investigated by positron annihilation and ultrasonic methods. J. Radioanal. Nucl. 240, 223–229. doi: 10.1007/BF02349158
Jones, G., and Dole, M. (1929). Density and viscosity studies on some R4NX in DME water mixtures at different temperatures. J. Am. Chem. Soc. 51, 2950–2956. doi: 10.1021/ja01385a012
Kaczkowska, E., Wawer, J., Tyczynska, M., Jozwiak, M., Borun, A., and Krakowiak, J. (2020). The hydration properties of protein stabilizer, trimethylamine-N-oxide in aqueous solutions of N-methylacetamide—The volumetric and compressibility studies between 288.15 and 308.15 K. Thermochim. Acta 685:178535. doi: 10.1016/j.tca.2020.178535
Kaur, A., Kaur, N., and Chaudhary, G. R. (2020). Volumetric and acoustic approach for investigating molecular interactions of choline acetate ionic liquid in α, ω-alkanediols at different temperatures. J. Mol. Liq. 315:113330. doi: 10.1016/j.molliq.2020.113330
Kumar, H., and Chadha, C. (2014). Densities, sound speed, and UV absorption studies of trisodium citrate and tripotassium citrate in aqueous solutions of 1-hexyl-3-methylimidazolium chloride [C6mim][cl]. J. Chem. Eng. Data 59, 4049–4061. doi: 10.1021/je5006906
Kumar, H., and Chadha, C. (2015). Interionic interactions of aqueous mixtures of ionic liquid 1-hexyl-3-methylimidazolium chloride with phosphate salts at T=(288.15 to 308.15) K: volumetric, acoustic and UV absorption studies. J. Chem. Eng. Data 60, 2937–2950. doi: 10.1021/acs.jced.5b00326
Kumar, H., and Katal, A. (2018a). Spectroscopic and temperature dependent physicochemical studies on interactional behavior of phosphate salts in aqueous ionic liquid (1-butyl-3-methyl imidazolium tetrafluoroborate) solutions. J. Mol. Liquids 251, 273–285. doi: 10.1016/j.molliq.2017.12.092
Kumar, H., and Katal, A. (2018b). Temperature dependent physicochemical and spectroscopic (FT-IR) studies of citrate salts (trilithium citrate and triammonium citrate) in aqueous ionic liquid [C4mim][BF4](1-butyl-3-methyl imidazolium tetrafluroborate) solutions. J. Mol. Liquids 256, 148–162. doi: 10.1016/j.molliq.2018.02.011
Kumar, H., and Katal, A. (2018c). Volumetric, acoustic and spectroscopic studies (FT-IR) of trisodium (TSC) and tripotassium citrate (TPC) in aqueous solution of ionic liquid 1-butyl-3-methylimidazolium tetrafluoro borate [C4mim][BF4] at different temperatures. J. Chem. Eng. Data 116, 85–96. doi: 10.1016/j.jct.2017.08.025
Kumar, H., and Sharma, R. (2020). Influence of 1-hexyl-3-methylimidazolium bromide ionic liquid on the volumetric and acoustic properties of amino acids (L-alanine and L-phenylalanine) at different temperatures. J. Mol. Liq. 304:112666. doi: 10.1016/j.molliq.2020.112666
Kumar, H., Singh, G., Kataria, R., and Sharma, S. K. (2020). Volumetric, acoustic and infrared spectroscopic study of amino acids in aqueous solutions of pyrrolidinium based ionic liquid, 1–butyl−1–methyl pyrrolidinium bromide. J. Mol. Liq. 303:112592. doi: 10.1016/j.molliq.2020.112592
Lin, G., Bian, P., and Lin, R. (2006). The limiting partial molar volume and transfer partial molar volume of glycylglycine in aqueous sodium halide solutions at 298.15 K and 308.15 K. J. Chem. Thermodyn. 38, 144–151. doi: 10.1016/j.jct.2005.04.012
Liu, Y., Guo, L., Zhao, D., Li, X., Gao, Z., Ding, T., et al. (2017). Enhanced activity of CuO/K2CO3/MgAl2O4 catalyst for lean NOx storage and reduction at high temperatures. RSC Adv. 7, 27405–27414. doi: 10.1039/C7RA03200E
Marcinkowski, Ł., Smiechowski, M., Szepiński, E., Kloskowski, A., and Warmińska, D. (2019). Interactions of N-alkyl-N-methylmorpholinium based ionic liquids with acetonitrile studied by density and velocity of sound measurements and molecular dynamics simulations. J. Mol. Liq. 286:110875. doi: 10.1016/j.molliq.2019.04.152
Marsh, K. N. (1987) Recommended Reference Materials for the Realization of Physicochemical Properties. Oxford: Blackwell.
Max, J., and Chapados, C. (1999). Influence of anomalous dispersion on the ATR spectra of aqueous solutions. Appl. Spectrosc. 53, 1045–1053. doi: 10.1366/0003702991948044
Mishra, A., and Gautam, S. (2001). Viscometric and volumetric studies of some transition metal chlorides in glycine water solution. Indian J. Chem. 40, 100–104.
Pal, A., and Gaba, R. (2008). Densities, excess molar volumes, speeds of sound and isothermal compressibilities for {2-(2-hexyloxyethoxy) ethanol n-alkanol} systems at temperatures between (288.15 and 308.15) K. J. Chem. Thermodyn. 40, 750–758. doi: 10.1016/j.jct.2008.01.015
Rajagopal, K., and Jayabalakrishnan, S. (2009). Volumetric and viscometric studies of 4-aminobutyric acid in aqueous solutions of salbutamol sulphate at 308.15, 313.15 and 318.15 K. Chin. J. Chem. Eng. 17, 796–804. doi: 10.1016/S1004-9541(08)60279-7
Redlich, O., and Meyer, D. M. (1964). The molal volumes of electrolytes. Chem. Rev. 64, 221–227. doi: 10.1021/cr60229a001
Sadeghi, R., and Mahdavi, A. (2012). (Vapour liquid) equilibria, volumetric and compressibility behaviour of binary and ternary aqueous solutions of 1-hexyl-3-methylimidazolium chloride, methyl potassium malonate, and ethyl potassium malonate. J. Chem. Thermodyn. 47, 347–357. doi: 10.1016/j.jct.2011.11.014
Sahin, E. M., and Ayranci, E. (2019). Comparison of solution behaviors of two structurally related ionic liquids in water and in acetonitrile as studied by volumetric and acoustic properties. J. Sol. Chem. 48, 1503–1518. doi: 10.1007/s10953-019-00931-w
Sarma, T., and Ahluwalia, J. (1973). Experimental studies on the structures of aqueous solutions of hydrophobic solutes. Chem. Soc. Rev. 2, 203–232. doi: 10.1039/cs9730200203
Sarmad, S., Mikkola, J., and Ji, X. (2017). Carbon dioxide capture with ionic liquids and deep eutectic solvents: a new generation of sorbents. ChemSusChem 10, 324–352. doi: 10.1002/cssc.201600987
Schmelzer, C., Zwirbla, W., Rosenfeld, E., and Linde, B. (2004). Acoustic investigations of pseudo-stable structures in aqueous solutions of polyethylene glycols. J. Mol. Struct. 699, 47–51. doi: 10.1016/j.molstruc.2004.04.027
Seethalakshmi, K., Rani, E. J. V., Padmavathy, R., and Radha, N. (2012). FT-IR spectral analysis of imidazolium chloride. Int. J. Curr. Res. Rev. 4, 31–36.
Sharma, S., Singh, M., Sharma, S., Singh, J., Sharma, A. K., and Sharma, M. (2020). Molecular interactions of L-histidine in aqueous ionic liquid [C4mim][BF4] solution at different temperatures: volumetric, acoustic and viscometric approach. J. Mol. Liq. 303:112596. doi: 10.1016/j.molliq.2020.112596
Toops, T. J., Smith, D. B., Epling, W. S., Parks, J. E., and Partridge, W. P. (2005). Quantified NOx adsorption on pt/K/gamma-al2O3 and the effects of CO2 and H2O. Appl. Cat. B Environ. 58, 255–264. doi: 10.1016/j.apcatb.2004.10.022
Verdía, P., Hernaiz, M., González, E. J., Macedo, E. A., Salgado, J., and Tojo, E. (2014). Effect of the number, position and length of alkyl chains on the physical properties of polysubstituted pyridinium ionic liquids. J. Chem. Thermodyn. 69, 19–26. doi: 10.1016/j.jct.2013.09.028
Wahab, A., and Mahiuddin, S. (2004). Isentropic compressibility, electrical conductivity, shear relaxation time, surface tension, and raman spectra of aqueous zinc nitrate solutions. J. Chem. Eng. Data 49, 126–132. doi: 10.1021/je0302001
Wang, J., Wang, H., Zhang, S., Zhang, H., and Zhao, Y. (2007). Conductivities, volumes, fluorescence, and aggregation behavior of ionic liquids [C4mim][BF4] and [C n mim] br (n= 4, 6, 8, 10, 12) in aqueous solutions. J. Chem. Phys. B 111, 6181–6188. doi: 10.1021/jp068798h
Wang, X., Chi, Y., and Mu, T. (2014). A review on the transport properties of ionic liquids. J. Mol. Liq. 193, 262–266. doi: 10.1016/j.molliq.2014.03.011
Wu, X., Sang, X., Li, Z., and Tao, D. (2020). Study on physicochemical properties and basicity of carbanion-functionalized ionic liquids. J. Mol. Liq. 312:113405. doi: 10.1016/j.molliq.2020.113405
Xu, A., Chen, L., and Wang, J. (2018). Functionalized imidazalium carboxylates for enhancing practical applicability in cellulose processing. Macromolecules 51, 4158–4166. doi: 10.1021/acs.macromol.8b00724
Yan, Z., Wang, J., and Lu, J. (2002). Viscosity behavior of some α-amino acids and their groups in water–sodium acetate mixtures. Biophys. Chem. 99, 199–207. doi: 10.1016/S0301-4622(02)00200-4
Yang, F., Ma, Q., Wang, X., and Liu, Z. (2017). Influence of aprotic cosolvents on the thermophysical properties of imidazolium-based ionic liquid. J. Chem. Eng. Data 62, 1628–1638. doi: 10.1021/acs.jced.7b00002
Zafarani-Moattar, M. T., Frouzesh, F., and Rafiee, H. R. (2014). The study of volumetric, acoustic and transport properties of ionic liquid, 1-butyl-3-methyl imidazolium chloride [bmim][cl] in aqueous lithium bromide solutions at T= 298.15–318.15 K. Fluid Phase Equilib. 376, 40–47. doi: 10.1016/j.fluid.2014.05.028
Zafarani-Moattar, M. T., and Sarmad, S. (2010). Effect of tri-potassium phosphate on volumetric, acoustic, and transport behaviour of aqueous solutions of 1-ethyl-3-methylimidazolium bromide at T=(298.15 to 318.15) K. J. Chem. Thermodyn. 42, 1213–1221. doi: 10.1016/j.jct.2010.04.025
Zafarani-Moattar, M. T., and Sarmad, S. (2011). Osmotic and activity coefficient of 1-ethyl-3-methylimidazolium chloride in aqueous solutions of tri-potassium phosphate, potassium carbonate, and potassium chloride at T= 298.15 K. Calphad 35, 331–341. doi: 10.1016/j.calphad.2011.04.008
Zafarani-Moattar, M. T., and Sarmad, S. (2012). Apparent molar volumes, apparent isentropic compressibilities, and viscosity B-coefficients of 1-ethyl-3-methylimidazolium bromide in aqueous di-potassium hydrogen phosphate and potassium di-hydrogen phosphate solutions at T=(298.15, 303.15, 308.15, 313.15, and 318.15) K. J. Chem. Thermodyn. 54, 192–203. doi: 10.1016/j.jct.2012.03.032
Zafarani-Moattar, M. T., Shekaari, H., and Jafari, P. (2020). Volumetric, acoustic and viscometric investigation of some choline amino acid ionic liquids in aqueous solutions of polypropylene glycol 400 and polyethylene glycol 400. J. Chem. Thermodyn. 142”106019. doi: 10.1016/j.jct.2019.106019
Keywords: ionic liquids, 1-ethyl-3-methyl-imidazolium chloride, apparent molar volume, apparent isentropic compressibility, viscosity B-coefficient
Citation: Sarmad S, Zafarani-Moattar MT, Nikjoo D and Mikkola J-P (2020) How Different Electrolytes Can Influence the Aqueous Solution Behavior of 1-Ethyl-3-Methylimidazolium Chloride: A Volumetric, Viscometric, and Infrared Spectroscopy Approach. Front. Chem. 8:593786. doi: 10.3389/fchem.2020.593786
Received: 11 August 2020; Accepted: 14 September 2020;
Published: 12 November 2020.
Edited by:
Ramesh L. Gardas, Indian Institute of Technology Madras, IndiaReviewed by:
Mohammad Tariq, Universidade NOVA de Lisboa, PortugalImran Khan, Sultan Qaboos University, Oman
Copyright © 2020 Sarmad, Zafarani-Moattar, Nikjoo and Mikkola. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shokat Sarmad, c2hva2F0LnNhcm1hZCYjeDAwMDQwO3VtdS5zZQ==; Dariush Nikjoo, RGFyaXVzaC5uaWtqb28mI3gwMDA0MDtsdHUuc2U=