Thermophysical and Transport Properties of Deep Eutectic Solvent Mixtures Containing Pure Choline Chloride/Ethylene Glycol with a Mole Ratio of 1:4 and Mixed with LiPF6 for Application in Redox-flow Batteries

1. Introduction

The rise in the world’s population, coupled with apprehensions about the overconsumption of fossil fuels and the threat of global warming, has necessitated a collective endeavor toward creating eco-friendly battery systems that can efficiently store electrical energy[1]. An effective power grid requires a certain degree of flexibility to maintain a balance between consumers’' electricity demand and its supply. Due to its intermittent nature, maintaining this balance is difficult, and as a result, traditional methods of electricity generation are still needed to fill the necessary gap[2]. Using energy storage to prevent such interruption and provide electricity even during times of low production and high demand is a promising solution. Thus far, various energy storage technologies have existed, including mechanical, electrical, chemical, and electrochemical storage.

Batteries are a family of energy storage technologies that cover a wide range of applications and have created great hopes for themselves[3]. Batteries take advantage of their ability to convert reversible electrical energy into chemical energy through redox reactions that occur inside the cell between positive and negative electrodes. Different materials and configurations of batteries enhance the efficiency of these energy storage technologies[3]. Current research trends focus on developing newer, more stable, and efficient materials, and methods for their production, which predict lower costs, higher capacities, and better performance systems[4]. This promising outlook is an important factor for continuous development and cost reduction, which must be considered when choosing a suitable storage technology for network scale use.

Redox flow batteries (RFBs) are an energy storage technology that stores large amounts of energy for long periods using two electrolyte solutions that undergo redox reactions[5]. During charging, electrical energy is used to oxidize the chemical species in one electrolyte while reducing the chemical species in the other. This process accumulates energy in the electrolyte, which is released as electrical energy during discharge[6]. RFBs consist of large external tanks containing positive and negative oxidizing and reducing materials in liquid form, which circulate through a central cell where redox reactions occur on carbon-based electrodes. A membrane separates the positive and negative materials to prevent cross-contamination while allowing only the passage of charge-balancing ions.

Deep eutectic solvents (DESs) are a new and less well-known class of liquid mixtures that have high conductivity, low vapor pressure, strong solvency power, and are made from inexpensive and environmentally friendly materials. These solvents have properties that make them suitable candidates for use in redox flow batteries[7].

The majority of studies on DESs as electrolytes for RFBs have focused on inorganic metal salts[8,9,10]. DESs consisting of choline chloride and urea have high solubility of vanadium and iron, resulting in an increased capacity compared with aqueous electrolytes [11]. However, the high viscosity of the DESs presented sluggish redox kinetics and large internal resistance, thus resulting in low voltage efficiency. Elevating the operating temperature has been shown to be lower viscosity and improve energy storage efficiency[11]. Recent studies have sought to build hybrid RFBs with DES electrolytes using metal chlorides as hydrogen bond acceptors (HBAs) coupled with hydrogen bond donors (HBDs) such as urea and acetamide[12,13,14]. These hybrid RFBs provide high concentrations of the redox species with high metal ion activities. These studies demonstrate that the H-bonding network in DESs increases viscosity and reduces conductivity, but it can also be modified through additives to weaken Coulombic interactions.

Metal ion coordination was investigated by Miller et al. to understand its impact on the electrochemical performance of an all-iron battery using ChlCl:Eg electrolytes[15]. The study found that the Fe⁺² /Fe⁺³ redox reaction was more robust compared with Fe⁺² /Fe⁰. The coordination structure of Fe, which depends on the Cl:Fe ratio in the electrolyte, plays a crucial role. At low concentrations of Cl⁻ (Cl:Fe < 4:1), the formation of Fe- Eg complexes dominated over that of tetrahedral iron complexes. These complexes improved electrochemical stability, reduced side reactions, and enhanced plating efficiency. However, the conductivity was lower in the 3:1 mixture than in the 4:1 mixture (Cl:Fe). Similarly, Shen et al. studied the effect of complexation on Cu⁺² /Cu⁺¹ redox kinetics in ChlCl:Eg electrolytes[16]. They found that the exchange current density (i₀) for Cu²⁺ to Cu¹⁺ reduction reaction was significantly smaller in ChlCl:Eg electrolytes than in aqueous electrolytes. The complexation of Cu⁺² and Cu⁺¹ with Cl⁻ was believed to be responsible for this difference. Further research is needed to understand the complexation and speciation of metal ions in DESs to gain insights into the reaction mechanisms and facilitate the design of redox flow batteries (RFBs).

To date, relatively little is known about the molecular properties of DESs that give rise to their unique macroscopic properties, and most of their applications have been specific to the aforementioned fields. Therefore, the design of DESs is mainly a trial-and-error process, and very little information is available on the molecular underpinnings that justify the performance of deep eutectic solvents or the methods for designing them to improve their macroscopic behavior. Using classical molecular dynamics (CMD) simulations and a multidisciplinary team of experimental and computational experts, models for DESs can be constructed and their properties validated. Using these models, the fundamental structural and dynamic properties of these solvents are finally revealed, and the way to exploit them to improve their performance as electrolytes for redox flow batteries is determined.

The current study focuses on experimental measurement and MD simulation of the density, viscosity, and ionic conductivity of ChlCl/Eg with a mole ratio of 1:4 DES doped with LiPF₆ salt. The aim of this study was to examine how the addition of LiPF₆ affects the density, viscosity, ion self-diffusion coefficients, and ionic conductivity of the solution from the perspective of an RFB electrolyte, using an experimental and MD simulation approach. A schematic representation of the imaginary RFB and the electrolyte used in this study are shown in Figure 1a and b, respectively.

2. Methods

2.1. Experimental Section

Materials, Apparatus, and Procedure 

To determine the transport properties of the ChlCl/Eg-LiPF₆ salt mixture, experimental techniques were used. Initially, the ChlCl/Eg mixture was produced by subjecting the components to a 1:4 molar ratio and, heating and stirring them at 351 K for 1 h until a transparent and homogeneous solution was obtained. Subsequently, LiPF₆ was dissolved in dry methanol under dry nitrogen gas and added to the ChlCl/Eg mixture. The resulting solution was stirred before the methanol was removed under vacuum. The chemicals employed in the experiment are listed in Table 1, accompanied by their source, CAS number, purity, and structure.

An analytical balance (AW 220, GR220, Shimadzu, Japan) was employed to prepare the DES with a precision of 10^-4 g. The experimental density of the mixture was measured using an Anton Para DSA 5000 at a frequency of approximately 3 MHz after calibration with distilled deionized and degassed water and dry air to maintain a constant temperature during the measurement. The Peltier device integrated within the densimeter provided a precision of 0.01 K.

The viscosity of the ChlCl/4Eg-LiPF₆ salt mixture was measured in a temperature-controlled water bath using a Ubbelohde-type viscometer. The experiment involved measuring the flow time and density of the mixture to obtain the viscosity data, which was calculated using the following equation:

$$\frac{\eta }{d}=kt-\frac{l}{t}$$

(1)

Where η is dynamic viscosity, d is density, k is the viscometer constant, and t is flow time. To determine the viscometer constants, the calibration of the viscometers was conducted at different temperatures using the density and viscosity of ethylene glycol. The average of three flow times was obtained using a stopwatch with a resolution of ±0.01 s. Precise temperature control was ensured throughout the measurement process.

The measurement of specific electrical conductance, represented by $\kappa$, was conducted using a Metrohm model 712 conductivity meter from Switzerland with an accuracy of ±0.005 μ S cm⁻¹. To calibrate the apparatus, the conductivity meter values were adjusted to the KCl solution with a concentration of 0.01 mol·kg⁻¹. The mixtures were created by adding pre-determined quantities of API-IL to the cell of the conductivity meter and stirring it with a magnetic stirrer. The temperature of the mixture was maintained using a thermostatic water bath, with an uncertainty of ±0.02 K. Using the formula $\sigma =\frac{1000\kappa }{C}$, the molar conductivities ($\sigma$) of the mixtures were computed given the molar concentration of API-IL and the specific conductance of the (glycine + API-IL + water) mixture.

2.2. Computational Section

The GROMACS 21.4 software package was utilized for performing calculations on four systems with a ChlCl/Eg mole ratio of 1:4 and LiPF₆ salt concentrations of 0. 0111, 0.0553, 0.1013, and 0.1502 M [17]. Table 2 provides information on the number of molecules, box size, and LiPF₆ concentration in each mixture. The module used for inserting molecules was employed to pack the molecules within the box. To observe the impact of solvents on thermophysical and transport properties, simulations were performed at different temperatures, including 298.15, 308.15, 318.15, 328.15, 348.15, and 398.15 K, and simulation times were selected to obtain sufficient sampling of thermophysical properties. A modified Berendsen thermostat and Parrinello−Rahman barostat were employed to control temperature and pressure, with coupling constants of 0.4 and 2.0 ps, respectively[18]. The LINC algorithm constrained all bond lengths[19], whereas the PME method calculated long-range electrostatic interactions [20]. In addition, a short-range cutoff of 1.1 nm was used for VDW and electrostatic interactions. The Verlet leapfrog algorithm [21] was used to integrate equations of motion with a time step of 2 fs, and periodic boundary conditions were employed in three directions.

Calculation of the Physicochemical Properties 

Density. Molecular dynamic simulations were performed using GROMACS software. The DES mixture containing choline chloride and ethylene glycol in a 1:4 molar ratio in a pure form and LiPF₆ was modeled in a cubic box. The system was energy-minimized to reduce any bad steric interactions. Temperature equilibration simulations were then performed in the NPT ensemble at given temperatures to obtain the required temperature. The equilibration protocol involved simulated annealing with an initial 2 ns at 500 K followed by cooling to the target temperature for 2 ns. The system was then equilibrated for an additional 6 ns at constant temperature and pressure. A time step of 2 fs was used for the integration of the equations of motion. Production runs were then performed for 40 ns in the NPT ensemble at constant temperature and pressure to calculate the density. Configurations were saved every 2 ps for data analysis.

Self-Diffusion Coefficients. MD simulations were performed using the GROMACS software, following a protocol similar to that for density calculations. The system was energy-minimized and then equilibrated at the desired temperature through simulated annealing. However, to calculate transport properties such as self-diffusion coefficients, production runs were performed in the NVT ensemble at constant number of particles, volume, and temperature, rather than in the NPT ensemble. This is a more suitable ensemble for transport property calculations. Production runs of 40 ns duration were performed in the NVT ensemble at the equilibrated temperature, with configurations saved every 2 ps. The mean-square displacements of the molecules were then determined from the trajectories, and the self-diffusion coefficients were calculated using the Einstein relation: [22]:

$$D=\frac{1}{6}\underset{t\to \infty }{lim}\frac{d}{dt}⟨|r_i(t)-r_0(t){|}^2⟩$$

(2)

where r_i is the position vector of particle i at time t. D was obtained from the slope of the linear region of the mean-square displacement (MSD) curves.

Shear Viscosity. To compute the viscosities of the DES mixtures, the same simulation methodology for diffusion coefficient calculations, as described in the previous section, was executed. In brief, the process entailed energy minimization, followed by simulated annealing for a duration of 4 ns in the NPT ensemble, and subsequently, a 40 ns production run in the NVT ensemble to calculate the transverse-current correlation function (TCCF) [23] via the transverse momentum fields [24]. Using this approach, the viscosity η_∞ was determined by fitting the 16 acquired values of η(k) corresponding to the 16 distinct k vectors, in line with the following equation:

$$\eta \left(k\right)={\eta }_{\infty }(1-a{k}^2)$$

(3)

In this equation a is the fitting parameter. The system size does not affect the viscosity in this method [25].

Ionic conductivity. The production runs performed in the NVT ensemble for calculating the self-diffusion coefficients and viscosity were also used to determine the ionic conductivity of the mixtures. The Nernst-Einstein equation was employed, which relates conductivity to self-diffusion coefficients and other system properties. The equation takes following form: [26]:

$${\sigma }_{NE}={\sigma }_{+}+{\sigma }_{-}$$

(4)

$${\sigma }_{NE}=\frac{q_{+}^2\rho D_{+}}{kT}+\frac{q_{-}^2\rho D_{-}}{kT}$$

(5)

where σ_NE is the Nernst-Einstein conductivity, D₊ and D₋ are the self-diffusion coefficients of the cation and anion, respectively, q₊ and q_- are their respective charges, ρ is the number density of the salt, k is the Boltzmann constant, and T is the temperature.

The Nernst-Einstein equation, while computationally convenient, overestimates ionic conductivity because of its inability to account for ion-ion correlations. However, it still provides a valuable estimate of conductivity based on the available simulation data.

3. Results and Discussion

To go deep inside the solvent structure at the atomic level, one needs to choose an appropriate forcefield to portrait the interaction of particles in the simulation box. Previous studies [27,28] suggested the FFM3 model from the work of Ferreira et al. [29] and found it to be the most suitable forcefield for pure ChlCl:Eg in 1:2 and 1:3 molar ratios and for their mixture with LiPF₆ and discussed the thermophysical, transport, and structural properties of the mixtures mentioned above. These works encourage us to investigate the same DES with a 1:4 molar ratio through experiment and MD simulation studies. Variable properties such as density, viscosity, diffusion coefficient, and ionic conductivity were calculated from the MD simulation and compared with experimental data for ChlCl: 4Eg:LiPF₆ mixtures to illustrate the applicability of the applied atomic model.

After careful consideration of the findings reported by Moradi et al. [27,28], we performed a charge scaling in this case study. However, our results did not yield satisfactory results for certain thermophysical properties, particularly the diffusion coefficient. Therefore, we have decided not to consider charge scaling in our study.

3.1. Thermophysical Properties

Primarily, our focus lies on the experimental density data presented in Table 3 for the ChlCl/4Eg + LiPF₆ mixtures at various temperatures and LiPF₆ concentrations. These experimental values can be compared with the density data from MD simulations to evaluate the accuracy of the simulation (Table 4). Upon scrutiny of these tables, it can be observed that the density of the mixtures increases with increasing LiPF₆ concentration for both experimental and MD simulation data. This trend is expected because LiPF₆ is a salt that imparts more ions to the mixture, consequently leading to a rise in the density of the solution.

Regarding temperature, the density of the mixtures decreases with increasing temperature for both experimental and MD simulation data. Figure 2(a) illustrates this trend for the simulation data while incorporating the uncertainties for each data point. This trend is also expected because augmenting the temperature results in an expansion of the solution's volume, subsequently decreasing density. However, a comparison of the simulation and experimental densities, as plotted in Figure 2(b), in terms of the percentage of relative error, indicates that the MD simulation data tend to underestimate the density compared with the experimental data. This discrepancy could be attributed to several factors, such as the force field employed and the simulation conditions.

Table 5 displays the results from MD simulations of the diffusion coefficients (D) for choline, chloride, ethylene glycol, lithium, and PF₆ in the deep eutectic solvent mixture. It is noteworthy that the force field parameters in the simulation did not use a reduced charge scaling of 0.8. To accurately estimate the mean square displacement (MSD) from molecular dynamics simulation trajectories, Del Popolo and Voth [30] proposed determining the onset of the diffusive regime by calculating the β(t) parameter. The β(t) parameter is defined as

$$\beta \left(t\right)=\frac{d{log}_{10}〈∆r\left(\left.t^2\right)\right.〉}{d{log}_{10}}$$

(6)

where <Δr(t)² >is the average MSD of the species and t is the time. The optimized value for β (t) indicates the best region for calculating the diffusion coefficient.

In a study by Moradi et al., it was observed that the diffusion coefficient of ethylene glycol (Eg) molecules in ChlCl:2Eg and ChlCl:3Eg mixtures peaked and rose as the number of Eg molecules in the solutions increased [27]. From their study, it is understandable that choline [Chl+] has the lowest diffusion coefficient because of its larger size than ethylene glycol (Eg) and chloride [Cl-] ions. However, the fact that chloride [Cl-] has a lower diffusion coefficient than ethylene glycol molecules despite its smaller size can be explained by considering that chloride [Cl-] did not diffuse independently as a single ion. Rather, the chloride [Cl-] ions likely formed complexes or were associated with the choline [Chl+] ions in the mixture, reducing their mobility and diffusion rate. This indicates that the chloride [Cl-] ions did not diffuse as freely as the ethylene glycol molecules but instead were coupled to the diffusion of the choline [Chl+] ions to some extent.

In this study, we can analyze the trends in the diffusion coefficient data and compare them with expected behavior based on previous knowledge of the system. First, we examine the overall trend in the data for each component. For example, we can observe that the diffusion coefficients for choline, ethylene glycol, and chloride tend to increase with increasing temperature (Figure 3), whereas the diffusion coefficients for [Li+] and [PF₆-] exhibit a drop when the temperature shifts from 298.15 K to 303.15 K and then increase with increasing temperature (Figure 3b-e). This behavior aligns with the expected behavior of ions in solution, where higher temperatures lead to increased thermal motion of the solvent molecules, enabling the ions to move more freely and increasing their diffusion coefficients.

Next, we examined the impact of LiPF₆ concentration on the diffusion coefficients. Figure 4 illustrates this trend in the diffusion coefficients. According to this figure, for choline, ethylene glycol, and chloride, the diffusion coefficients tend to marginally decrease with increasing LiPF₆ concentration. This is likely due to the increased ionic strength of the solution at higher LiPF₆ concentrations, which can lead to increased electrostatic interactions between ions and reduced mobility. However, based on Figure 4, we are not observing regular and systematic behavior for the diffusion coefficients of [Li+] and [PF₆-] ions in the mixtures at different temperatures by increasing the LiPF₆ concentration.

In general, viscosity is a measure of a liquid’s resistance to flow, and it depends on various factors such as the size, shape, and interaction of molecules, concentration of solutes, and temperature. Electrostatic interactions between ions and the formation of ion clusters can also influence viscosity. The experimental viscosity data of ChlCl/4Eg mixtures with LiPF₆ at different concentrations and temperatures are presented in Table 6. The data show that viscosity generally increases with concentration and temperature, as expected for a solute added to a solvent. Increasing the concentration of solutes generally increases the viscosity of the solution.

The effect of temperature on viscosity can be explained by the fact that increasing temperature usually reduces the strength of intermolecular forces, making molecules more mobile and decreasing viscosity. In addition, we calculated the viscosity of the aforementioned mixtures using computational approaches, and the results are presented in Table 7. The viscosities obtained from the transverse-current correlation function (TCCF) are compared with the experimental data. The simulation data show a similar trend of increasing viscosity with increasing concentration and temperature. However, MD simulations tend to overestimate viscosity, especially at lower concentrations and temperatures. This indicates that MD simulations may overestimate the strength of the interactions between molecules in the mixture. Figure 5 shows the percentage of relative error between the simulation and experimental data.

Both experimental and MD simulation data show a decrease in viscosity with increasing temperature, as higher temperatures lead to increased molecular motion and reduced intermolecular interactions. The rate of decrease in viscosity with temperature seems to be faster for the experimental data than for the MD simulation data, especially at lower concentrations.

In terms of concentration dependance, both experimental and MD simulation data show an increase in viscosity with increasing concentration of LiPF₆. The addition of LiPF₆ leads to increased ionic strength and stronger ion-ion interactions, which in turn increase the viscosity. The rate of increase in viscosity with concentration seems to be faster for the experimental data than for the MD simulation data, especially at lower temperatures.

Ionic conductivity is a measure of a solution’s ability to conduct electricity, which depends on the mobility of ions in the solution and is influenced by factors such as ion concentration, temperature, and solvent properties. As ion concentration and/or temperature increase, ionic conductivity typically increases. The experimental ionic conductivity data for mixtures of ChlCl/Eg + LiPF₆ at different temperatures and compositions, are presented in Table 8.

The Nernst-Einstein equation relates the ionic conductivity of a solution to the diffusion coefficient and charge of the ions in the solution (equation 6). The equation was used to calculate the ionic conductivity of these mixtures, and the results are presented in Table 9. To use this equation to calculate the ionic conductivity, we used the diffusion coefficient data described above. Then, the number density of the ions is obtained from the concentration of the ions and the volume of the simulation box. The MD simulation data shows similar trends to the experimental data, but the magnitudes of the ionic conductivity values are different. The MD simulation results tend to underestimate the ionic conductivity compared with the experimental data, especially at lower LiPF₆ concentrations. This could be due to various factors such as force field parameters, simulation time, or treatment of long-range electrostatic interactions.

Figure 6 illustrates the trend of these data by variation in temperature (Figure 6a) and concentration of LiPF₆ (Figure 6b). The data show that the ionic conductivity decreases with increasing LiPF₆ concentration at all temperatures, indicating that the presence of LiPF₆ reduces the ionic conductivity of the mixture. This result can be attributed to the increase in ion-pairing and decrease in ion mobility due to the presence of LiPF₆. Higher temperatures at the same LiPF₆ concentration result in higher ionic conductivity due to increased thermal energy, leading to greater ion mobility (Figure 6a). A plot of the percentage of relative error between the simulation and experimental ionic conductivity data for ChlCl:4Eg versus the concentration of LiPF₆ in different temperature also represented in Figure 6c.

3.1. Structural Analysis

Determining the structural configuration of each constituent in a mixture can provide valuable information regarding the main internal forces between the components of the liquid. To achieve this, we used atom-atom radial distribution functions (RDF) analysis at a temperature of 298 K. We first tried to find the arrangement of choline toward the Li⁺ ion. Because both sides of choline, the nitrogen (N) and oxygen (O) atom, are of interest, we calculate the RDF of these atoms with Li⁺. Figure 6a shows the RDF for the N atom in choline with Li⁺ at different concentrations. According to this figure, at any LiPF₆ concentrations (0.01-0.15 M), there is essentially no coordination between nitrogen and lithium, as indicated by RDF values 0 up to approximately 4 angstroms. Then, an RDF peak begins to emerge around 4-4.5 angstroms, indicating the formation of the first solvation shell with lithium ions coordinated to nitrogen atoms. The height and width of the first peak increase with LiPF₆ concentration, showing that more lithium ionsinteract with nitrogen and the distances of these interactions increase. As we can see in Figure 6b, the same RDF peak begins around 4-4.5 angstroms, but the intensity for 0.01M LiPF₆ is higher than that for the other concentrations indicating that the Li ion tends to arrange toward both sides of the choline molecules at lower concentrations. The interaction between the nitrogen atom in choline and chloride ion is represented as the RDF curve in Figure 6c. This indicates that the first solvation shell with Cl ions coordinated to nitrogen atoms begins to appear around 4-4.5 angstroms. According to Figure 6d, the first coordination shell, which is the distance at which the first peak occurs, is approximately 2.9 angstroms, indicating that the chloride ion and oxygen atom in choline are on average separated by this distance in the DES mixture. According to this figure, increasing the Li ion does not change the intensity of correlation between these two atoms, indicating that the hydrogen bonds between chloride ion and oxygen atom has been not affect by introducing more LiPF₆ into the mixture. The RDF values are generally higher at higher LiPF₆ concentrations, indicating greater correlation and interaction between the nitrogen and the hydrogen atoms. This shows that as the salt concentration increases, the local structure around these atoms is modified and they tend to come closer together. Figure 6e shows the RDF data of the nitrogen atom in choline and the hydrogen atoms in ethylene glycol in 0, 0.01, 0.05, 0.1 and 0.15 molar of LiPF₆. There is a slight difference in RDF values between the 0.01M, 0.1M, and 0.15M concentrations, which slightly increases as the distance decreases. the peaks in the RDF data around 4 Angstroms and 6 Angstroms indicate that the nitrogen and hydrogen atoms exhibit preferential correlations or interactions at those distances. The 4 Angstrom correlations are much stronger and become more prominent at higher LiPF₆ concentrations, indicating that they may be facilitated by the salt ions. The 6 Angstrom correlations were weaker and showed a smaller increase with concentration. We investigated the correlation between oxygen atoms in choline and hydrogen atoms in ethylene glycol in 0.01, 0.05, 0.1, and 0.15 molar of LiPF₆ and the results are depicted in Figure 6f. It shows that the oxygen and hydrogen atoms have a stronger propensity to interact, and correlate at all distances, as seen from the generally higher RDF values. The local environment around the atoms is perturbed more substantially, allowing for stronger interactions, especially at closer lengths (<4 Angstroms). Preferential correlations emerge at approximately 3–4 Angstroms and are enhanced at higher concentrations, indicating that the LiPF₆ ions facilitate these interactions. Figure 6g illustrates the interaction between oxygen in ethylene glycol and Li ions. We can see a strong interaction between these two atoms, indicating that there is a large coordination of Li ions around the ethylene glycol molecules. Increasing the LiPF₆ concentration makes this interaction stronger, as we expected, because more Li ions are accessible to the ethylene glycol in the mixture. The same correlation observed in the RDF plot of chloride ions with hydrogen atoms in ethylene glycol, which is depicted in Figure 6h. the distance and intensity of the peak in this figure indicates that strong interaction takes place between these two atoms in the mixture. It is also interesting to evaluate the arrangement of Li ion toward chloride ion in the mixture. We have such a plot in Figure 6i, which shows that strong coordination between Li and chloride ions occurs at a distance of 2 angstroms and increases with increasing the concentration of LiPF₆.

4. Conclusions

In this study, we synthesized DES mixtures with varying concentrations of LiPF6 salt while maintaining a constant constituent ratio of ChlCl/Eg as 1:4 mole ratio. The mixtures remained stable at temperatures between 298 and 398 K. The experimental density, viscosity, and ionic cunductivity of the mixtures increased with salt concentration. The study also used MD simulation to validate the forcefield used for other transport and structural properties. The ionic conductivities of the mixtures decreased with increasing salt concentration, with the highest conductivity observed in the mixture containing 0.011 M LiPF₆. The diffusion coefficients for all cations and anions increased with temperature, but decreased with increasing salt concentration due to electrostatic interactions between counter ions. The RDFs showed that the Li+ ion was in close contact with Cl and strongly connected to oxygen atoms in the Eg molecules.