Julio I. Urzúa*a,
María Luisa Valenzuelab,
Jenifer Cavieresc and
María José Inestrosa-Izurietaa
aCentro de Materiales para la Transición y Sostenibilidad Energética, Comisión Chilena de Energía Nuclear, Ruta 68, km 20, Pudahuel, Santiago, Chile. E-mail: julio.urzua@cchen.cl
bUniversidad Autónoma de Chile, Facultad de Ingeniería, Instituto de Ciencias Químicas Aplicadas, Grupo de Investigación en Energía y Procesos Sustentables, Av. El Llano Subercaseaux 2801, Santiago, Chile
cDepartamento de Ingeniería Química y Ambiental, Universidad Técnica Federico Santa María, Av. Vicuña Mackenna 3939, San Joaquín, Santiago, Chile
First published on 27th June 2023
Fifteen equimolar binary mixtures are synthesized and thermophysically evaluated in this study. These mixtures are derived from six ionic liquids (ILs) based on methylimidazolium and 2,3-dimethylimidazolium cations with butyl chains. The objective is to compare and elucidate the impact of small structural changes on the thermal properties. The preliminary results are compared to previously obtained results with mixtures containing longer eight-carbon chains. The study demonstrates that certain mixtures exhibit an increase in their heat capacity. Additionally, due to their higher densities, these mixtures achieve a thermal storage density equivalent to that of mixtures with longer chains. Moreover, their thermal storage density surpasses that of some conventional materials commonly used for energy storage.
These storage technologies encompass various categories, each relying on different mechanisms.1–3 Chemical storage involves energy conversion through chemical reactions, exemplified by fuel cells, lithium-ion and post-lithium batteries.4–6 Mechanical storage harnesses energy using mechanisms like flywheels and air compression.7,8 Electrochemical storage utilizes supercapacitors or electrochemical capacitors to store energy as electrical charge.9,10 And finally, thermal storage uses heat for later use, by heating materials such as molten salts or water, or by storing latent heat in phase-change materials.
This latest technology has aroused great scientific and technological interest due to its simple and practical application. When applied on an industrial scale, known as concentrating solar power (CSP), it demonstrates an impressive ability to generate electricity from solar energy. CSP power plants consist of numerous solar collector mirrors that redirect the sun's rays to a specific area. Within this zone, a heat transfer fluid (HTF) is stored or housed in a thermal energy storage (TES) system.11–14
Nitrate salts (40% KNO3 + 60% NaNO3), eutectic mixture of diphenyl (DPO) and biphenyl oxide (Therminol VP-1), and mineral or silicone oil are the most commonly used materials as HTF and TES in CSP solar plants. However, each has aspects that hinder the efficiency of the process. For example, nitrated salts have a high melting point, which requires high temperatures, and run the risk of solidification in cold weather; additionally, they generate a high degree of corrosion. Therminol VP-1 has a high vapor pressure at high temperature, which could be harmful to pipelines and increase the risk of ignition; mineral oils and silicones, on the other hand, have a limited working range since they decompose at low temperatures.15 In view of the above, the development of new storage materials is key. In this context, ionic liquids (ILs) are materials that have been proposed in a wide range of applications,16 and their use as heat transfer and thermal storage fluids has been proposed based on their unique properties, such as low vapor pressure, non-flammability, chemical and thermal stability, and a wide temperature range in the liquid state, properties that offer significant advantages over molten salts and thermal oils.
Considerable studies indicate that methylimidazolium-based ionic liquids exhibit interesting thermal properties; however, many of them do not outperform the thermal properties of conventional materials, such as thermal storage density, a key property in HTF and TES that relates density (ρ), heat capacity (Cp), and operating temperature. This has prompted to explore the properties of ionic liquid mixtures, where systems formed by coulombic, van der Waals, and hydrogen bonding interactions of more than two types of ions are considered as an alternative to improve thermal properties.17,18
In addition to the studies performed by other researchers evaluating the thermophysical properties of ionic liquid mixtures,19–23 our group has previously demonstrated that certain binary mixtures of ionic liquids based on 1-octyl-3-methylimidazolium and 1-octyl-2,3-dimethyl imidazolium increased their heat capacity and, therefore, their thermal storage density, which allowed preliminary comparison with materials conventionally used in solar energy storage.24
In order to explore the behaviour of mixtures of similar ILs, with small structural differences, and to compare their properties with previously reported mixtures and conventional materials, in this work we carried out the synthesis and thermophysical characterization of six ionic liquids with butyl chains (Table 1), and their equimolar mixtures, in attempt to determine the variation of their properties from those previously defined by the pure components and, in turn, to evaluate these new compounds as potential heat transfer fluids.
No. | Compound | Abbreviation | State at r.t. |
---|---|---|---|
IL1 | 1-Butyl-3-methylimidazolium tetrafluoroborate | [Bmim][BF4] | Liquid |
IL2 | 1-Butyl-3-methylimidazolium hexafluorophosphate | [Bmim][PF6] | Liquid |
IL3 | 1-Butyl-3-methylimidazolium bis(trifluoromethyl sulfony) imide | [Bmim][Tf2N] | Liquid |
IL4 | 1-Butyl-2,3-dimethylimidazolium tetrafluoroborate | [Bdmim][BF4] | Liquid |
IL5 | 1-Butyl-2,3-dimethylimidazolium hexafluorophosphate | [Bdmim][PF6] | Solid |
IL6 | 1-Butyl-2,3-dimethylimidazolium bis(trifluoromethyl sulfony)imide | [Bdmim][Tf2N] | Liquid |
ρ = A0 + A1T + A2T2 | (1) |
As reported in the literature on ILs,25,26 the density of mixtures and ILs decreases with increasing temperature. In order to compare the results obtained, the density values at 373 K are shown in Table 2. The density values obtained for the mixtures, in most cases, correspond to average values of the densities of their pure components. The measured density range for all species is from 1137.5 kg m−3 for [Bdmim][BF4] to 1354.9 kg m−3 corresponding to mixture 9 [Bmim][Tf2N] + [Bdmim][Tf2N], in this case slightly increasing the density respect to the pure systems.
No. | IL | Density |
---|---|---|
IL1 | [Bmim][BF4] | 1143.7 |
IL2 | [Bmim][PF6] | 1301.2 |
IL3 | [Bmim][Tf2N] | 1350.6 |
IL4 | [Bdmim][BF4] | 1137.5 |
IL5 | [Bdmim][PF6] | 1282.2 |
IL6 | [Bdmim][Tf2N] | 1340.5 |
No. | Mixture | Density |
---|---|---|
M1 | [Bmim][BF4] + [Bmim][Tf2N] | 1287.4 |
M2 | [Bmim][BF4] + [Bmim][PF6] | 1235.8 |
M3 | [Bmim][PF6] + [Bmim][Tf2N] | 1332.9 |
M4 | [Bdmim][BF4] + [Bdmim][Tf2N] | 1268.5 |
M5 | [Bdmim][BF4] + [Bdmim][PF6] | 1211.8 |
M6 | [Bdmim][PF6] + [Bdmim][Tf2N] | 1316.6 |
M7 | [Bmim][BF4] + [Bdmim][BF4] | 1151.9 |
M8 | [Bmim][PF6] + [Bdmim][PF6] | 1283.6 |
M9 | [Bmim][Tf2N] + [Bdmim][Tf2N] | 1354.9 |
M10 | [Bmim][BF4] + [Bdmim][PF6] | 1226.1 |
M11 | [Bmim][BF4] + [Bdmim][Tf2N] | 1268.3 |
M12 | [Bmim][PF6] + [Bdmim][BF4] | 1215.9 |
M13 | [Bmim][PF6] + [Bdmim][Tf2N] | 1323.9 |
M14 | [Bmim][Tf2N] + [Bdmim][BF4] | 1271.6 |
M15 | [Bmim][Tf2N] + [Bdmim][PF6] | 1326.2 |
On the other hand, the mixtures previously obtained with longer alkyl chains (Octyl) have a range of densities from 1059.8 kg m−3 of [Omim][BF4] to 1255.1 kg m−3 of the [Omim][Tf2N] + [Odmim][Tf2N] mixtures. The variation in density values is in agreement with that reported in the literature27 regarding the increase in density with decreasing chain length of imidazolium cation-based ionic liquids.
No. | IL | Tonset | Cp | No. | Mixture | Tonset | Cp | No. | Mixture | Tonset | Cp |
---|---|---|---|---|---|---|---|---|---|---|---|
IL1 | [Bmim][BF4] | 675 | 1.72 | M1 | [Bmim][BF4] + [Bmim][Tf2N] | 663 | 1.55 | M10 | [Bmim][BF4] + [Bdmim][PF6] | 658 | 1.39 |
IL2 | [Bmim][PF6] | 649 | 1.61 | M2 | [Bmim][BF4] + [Bmim][PF6] | 640 | 1.11 | M11 | [Bmim][BF4] + [Bdmim][Tf2N] | 678 | 2.62 |
IL3 | [Bmim][Tf2N] | 704 | 1.45 | M3 | [Bmim][PF6] + [Bmim][Tf2N] | 650 | 1.37 | M12 | [Bmim][PF6] + [Bdmim][BF4] | 659 | 1.80 |
IL4 | [Bdmim][BF4] | 767 | 1.81 | M4 | [Bdmim][BF4] + [Bdmim][Tf2N] | 685 | 1.02 | M13 | [Bmim][PF6] + [Bdmim][Tf2N] | 669 | 1.68 |
IL5 | [Bdmim][PF6] | 656 | 1.65 | M5 | [Bdmim][BF4] + [Bdmim][PF6] | 656 | 1.38 | M14 | [Bmim][Tf2N] + [Bdmim][BF4] | 683 | 2.23 |
IL6 | [Bdmim][Tf2N] | 707 | 1.51 | M6 | [Bdmim][PF6] + [Bdmim][Tf2N] | 671 | 1.36 | M15 | [Bmim][Tf2N] + [Bdmim][PF6] | 639 | 2.13 |
M7 | [Bmim][BF4] + [Bdmim][BF4] | 652 | 1.75 | ||||||||
M8 | [Bmim][PF6] + [Bdmim][PF6] | 647 | 1.50 | ||||||||
M9 | [Bmim][Tf2N] + [Bdmim][Tf2N] | 690 | 1.22 |
In some cases, the properties of the mixtures can be simply deduced from the neat systems, as is the case for density, but in others the interactions do not show the ideally expected behaviour. Most of the time, the thermal properties must be understood by studying the intermolecular interactions and charge distribution of the ionic liquids.
The six initial ILs have an average decomposition temperature (Tonset) ranging between 664 K and 703 K, but the average decomposition temperatures of the mixtures are at lower values with respect to their pure systems, i.e., when two ionic liquids are mixed, their thermal stability decreases.
On the other hand, compared to the values obtained previously with octyl chains, mixtures with shorter chains show slightly lower decomposition values, contrary to those reported in the literature.28
In the case of binary mixtures where two anions or two cations are equal, or in other words, three ions, it is observed in the thermogram DTA as a trend, a single peak of maximum decay is similar to one of the pure ILs. This occurs in mixtures 2, 5, 7, 8, and 9. This could indicate that there is an exchange of ions forming an equilibrium.
In M2 (IL1 + IL2), the thermogram DTA acquires the same thermal behaviour as IL2; in the case of M5 (IL4 + IL6), it behaves like IL4; M7 (IL1 + IL4) like IL4; M8 (IL2 + IL5) like IL5, and M9 (IL3 + IL6) like IL6.
In these cases, there is clearly a dominance of the thermal behaviour of ILs based on 1-butyl-2,3-dimethylimidazolium. This could be due to the influence exerted by the methyl group located at the C(2) position of the imidazolium cation (Fig. 1), since this site represents a predominant site for interionic hydrogen bonding. The C(2)–H bond possesses a slightly acidic nature, being substituted by a methyl, it favour charge distribution by forcing the anion to interact to a greater extent via hydrogen bonds with C(4) and C(5).29
![]() | ||
Fig. 1 Structures of imidazolium cations. The dotted arrow indicates the positions of possible hydrogen bonds. |
This domain can also be observed in M4 (IL4 + IL6) and M6 (IL5 + IL6) where the behaviour of each of the base ILs prevails and two decay peaks are observed. On the other hand, in M1 (IL1 + IL3) and M3 (IL2 + IL3) the behaviour is exclusively due to the influence of the anions.30
In binary reciprocal mixtures, composed of four ions, they also reflect the influence of the methyl group located at the C(2) position of the imidazolium cation, since two peaks corresponding to each of the base ILs are clearly defined. This could be one of the reasons for the higher thermal stability of the 1-alkyl-2,3-dimethyl imidazolium cations,31 and unlike other properties measured in reciprocal binary mixtures,32–34 in this study it was observed that within a mixture with four different ions, some prior organization could persist that would give rise to different thermal behaviours, even when there are no stability differences between the ions involved.
In the case of M10 (IL1 + IL5), in the thermogram DTA, two peaks of thermal degradation are clearly and well defined, the first peak being IL1 and the second peak IL5; in the case of M11 (IL1 + IL6), M12 (IL2 + IL4), M13 (IL2 + IL6), M14 (IL3 + IL4), and M15 (IL3 + IL5), all two peaks of thermal degradation in the order described.
Finally, in this study, the TGA analysis shows a high thermal stability of the mixtures, and the average Tonset value for all of them is around 662 K.
Cp was calculated from its definition using the following eqn (2).35
![]() | (2) |
![]() | (3) |
T0 fit the heat capacity data, the quadratic polynomial eqn (4) was used. The fitting parameters, correlation factor and other data are available in the ESI.†
Cp = A0 + A1T + A2T2 | (4) |
For comparison, heat capacities data at 373 K are presented in Table 3.
Part of the aim of this study is to corroborate the increase in Cp after mixing two ionic liquids. In the group of mixtures under study, only four of the fifteen mixtures (mixtures 11, 13, 14, and 15) increased their Cp values from neat systems, three of them (1, 7, and 12) showed an intermediate value, and the remaining eight decreased their heat capacity compared to the initial ILs.
Compared to previously reported13 values for mixtures with longer alkyl chains, the Cp values are lower, but with an increase reasonably proportional to the basal ionic liquids. A comparison of the Cp values obtained in this study and those obtained previously is shown in Table 4. In addition, the values of the starting ionic liquids are incorporated. In the case of the four mixtures, reported in this work, which increased their Cp values from the neat systems, they are all binary reciprocal mixtures.
Mixture | Cp IL | Cp | Mixture | Cp IL | Cp |
---|---|---|---|---|---|
[Bmim][BF4] + [Bdmim][Tf2N] | 1.72 | 2.62 | [Odmim][BF4] + [Odmim][Tf2N] | 1.83 | 2.01 |
1.51 | 1.21 | ||||
[Bmim][PF6] + [Bdmim][Tf2N] | 1.61 | 1.68 | [Odmim][BF4] + [Odmim][PF6] | 1.83 | 2.83 |
1.51 | 2.61 | ||||
[Bmim][Tf2N] + [Bdmim][BF4] | 1.45 | 2.23 | [Omim][PF6] + [Odmim][PF6] | 1.59 | 2.97 |
1.81 | 2.61 | ||||
[Bmim][Tf2N] + [Bdmim][PF6] | 1.45 | 2.13 | [Omim][Tf2N] + [Odmim][Tf2N] | 1.23 | 1.99 |
1.65 | 1.21 | ||||
[Omim][BF4] + [Odmim][Tf2N] | 1.63 | 1.70 | |||
1.21 | |||||
[Omim][PF6] + [Odmim][BF4] | 1.59 | 2.28 | |||
1.83 | |||||
[Omim][PF6] + [Odmim][Tf2N] | 1.59 | 2.31 | |||
1.21 |
There is no structural pattern that can be analysed for each of the mixtures, so the variation in heat capacity is difficult to predict. Ion pairs in an IL are organized into polar and non-polar domains, resulting in structural, transition and relative size segregation of high and low charge regions. These aspects can be influenced by alkyl chain length, and have a direct impact on fluid properties. Thus, the study of intermolecular interactions between anions and cations could give some insight into the molecular behaviour that influences each of the thermophysical properties, with greater interest when studying mixtures, which involve a larger number of ions.
Looking at the cation charges by molecular dynamics, the hydrogen atoms in C(2), C(4), and C(5) apparently carry most of the positive charge, while the negative charge is distributed in N(1) and N(2), as well as between the C(4) and C(5) atoms. Thus, the C(4) and C(5) atoms appear to be almost neutral, whereas the C(2) atom possesses exclusively a positive charge due to electron deficit, which results in higher acidity and a tendency to form strong hydrogen bonds, while, on the contrary, C(4) and C(5) possess electron-donating properties.36,37
The impact of chain length is another aspect to consider. Kasemi et al. conducted molecular dynamics simulations to investigate the influence of chain length in [Cnmim][Tf2N] (n = 2, 8) ionic liquids. Their findings revealed that with a longer alkyl chain (n = 8), there is a greater likelihood of locating the C(2)–H groups of the cation and the anionic oxygen atoms within the polar domain. This means that there is a higher probability that hydrogen bonds will form, leading to higher vibrational energy transfer.38 Garaga et al.39 showed that, as the alkyl chain length increases, the hydrogen bonds and van der Waals interactions become stronger, while the coulombic strength decreases, and the study of the imidazolium cation and the [Tf2N]− anion in [Cnmim[[Tf2N]] revealed that the overall rotational mobility decreases with n, while different dynamic properties are observed, with the polar domains being the most dynamic and the imidazolium ring and its nearest chain segment being the most rigid. Concluding that the formation of segregated nanodomains in [Cnmim][Tf2N] is accompanied by the formation of nanodynamic heterogeneities, therefore, the cation slows down its rotational mobility more than the anion, contrary to what is observed for shorter chains.40 This reinforces the idea put forward by Rocha et al.,41 who show the increase of the heat capacity in ionic liquids [Cnmim][Tf2N] as the chain length increases. This analysis could respond to the difference in Cp observed between mixtures with octyl and butyl substituents.
These different analyses could help to understand the heat capacity increase of the mixtures under study from their structures and components. However, for a deeper understanding, it is necessary to further analyse the intermolecular effects as a set because these are complex systems where the binary reciprocal equimolar mixtures differ in their structures and the properties do not follow a deducible pattern.
According to Hunt et al.,42 the effects caused by the loss of hydrogen bonds could be compensated by the loss of entropy. The disorder in the system is reduced by the elimination of ion pair conformers that are stable for 1-alkyl-3-methyl imidazolium cations, but not for 1-alkyl-2,3-dimethyl imidazolium cations and an increase in the rotational barrier of the butyl chain, which limits free rotation and facilitates the association of the alkyl chain. This analysis could explain that all the mixtures where Cp increases involve a 1-alkyl-2,3-dimethyl imidazolium cation.
For the calculation of E in mixtures and ILs, eqn (5) is used:43,44
E = ρCp(T0 − Ti) | (5) |
Fluid | Reference | Cp (kJ kg−1 K−1) | ρ (kg m−3) | E (MJ m−3) |
---|---|---|---|---|
[Bmim][BF4] | This work | 1.7 | 1143 | 194 |
[Bmim][PF6] | This work | 1.6 | 1301 | 208 |
[Bmim][Tf2N] | This work | 1.4 | 1351 | 189 |
[Bdmim][BF4] | This work | 1.8 | 1138 | 205 |
[Bdmim][PF6] | This work | 1.6 | 1282 | 205 |
[Omim][PF6] | 24 | 1.6 | 1181 | 189 |
[Odmim][BF4] | 24 | 1.8 | 1190 | 214 |
[Odmim][PF6] | 24 | 2.6 | 1088 | 283 |
[Omim][PF6] + [Odmim][PF6] | 24 | 2.9 | 1116 | 324 |
[Omim][PF6] + [Odmim][BF4] | 24 | 2.3 | 1199 | 276 |
[Bmim][BF4] + [Bdmim][Tf2N] | This work | 2.6 | 1268 | 330 |
[Bmim][PF6] + [Bdmim][Tf2N] | This work | 1.7 | 1324 | 225 |
[Bmim][Tf2N] + [Bdmim][BF4] | This work | 2.2 | 1272 | 280 |
[Bmim][Tf2N] + [Bdmim][PF6] | This work | 2.1 | 1326 | 278 |
In Table 6, it is possible to compare the E values of the four mixtures that increased their Cp in this study vs. the mixtures with higher Cp with eight-carbon alkyl chain obtained in a previous report.
Fluid | Reference | ρ (kg m−3) | CP (kJ kg−1 K−1) | T (K) | E (MJ m−3) |
---|---|---|---|---|---|
Tinlet − Toulet | |||||
Mineral oil | 43 | 770 | 2.60 | 473.15 − 573.15 | 200 |
Silicone oil | 43 | 900 | 2.10 | 573.15 − 673.15 | 189 |
Synthetic oil | 43 | 900 | 2.30 | 523.15 − 533.15 | 207 |
Therminol VP-1 | 45 | 991 | 1.77 | 523.15 − 663.15 | 246 |
[Omim][PF6] + [Odmim][PF6] | 24 | 1116 | 2.97 | 523.15 − 663.15 | 464 |
[Omim][PF6] + [Odmim][BF4] | 24 | 1199 | 2.28 | 523.15 − 663.15 | 383 |
[Bmim][BF4] + [Bdmim][Tf2N] | This work | 1268 | 2.62 | 523.15 − 663.15 | 465 |
[Bmim][Tf2N] + [Bdmim][BF4] | This work | 1272 | 2.23 | 523.15 − 663.15 | 397 |
Nitrate salts | 46 | 1870 | 1.54 | 538.15 − 838.15 | 864 |
The results show the dependence of E on density and Cp, since the [Omim][PF6] + [Odmim][PF6] blend has the highest Cp = 2.9 kJ kg−1 K−1 and a density of 1116 kg m−3, resulting in E = 324 MJ m−3, while the mixture [Bmim][BF4] + [Bdmim][Tf2N] has a lower value of Cp = 2.6, but a higher density, 1268 kg m−3, which results in a slightly higher value of E, 330 MJ m−3.
It is necessary to consider the working temperatures to make an objective comparison of the thermal storage density values between the mixtures studied and the conventional materials, these values can be observed in Table 6.
Table 6 also shows that the thermal storage densities of the thermal oils are substantially lower than the mixtures, mainly due to their lower density and heat capacity values.
In comparison to the nitrate salts, there is an interesting relationship because, while the density is higher in the nitrate salts, the mixtures have higher Cp values. These differences could compensate the final E value; however, the working temperature (defined mainly by thermal stability) finally determines higher values for nitrated salts.
It is important to highlight that the results obtained in this article demonstrate high reproducibility. The experiments were conducted independently and repeated in triplicate, ensuring the robustness of the data. The results showed significant consistency, with an average variation of less than 5%. Furthermore, comparisons were made with previous research in the field, and it was found that the results are consistent with the trends reported in the literature.
1H and 13C NMR spectra were recorded on Bruker AMX-400 (11.74 T, 400 MHz to 1H and 126 MHz to 13C) NMR spectrometers using the residual proton or the carbon signal of the deuterated solvent as an internal standard.
Thermogravimetric analysis (TGA-DTA) was conducted using a PerkinElmer TGA400S model, operating at a voltage range of 100–240 V and frequency of 50–60 Hz. The experiments were carried out in the temperature range of 303 K to 873 K, with a flow rate of 50 ml min−1 and a heating rate of 10 °C min−1. The samples, weighing between 2.5 to 5 mg, were placed in platinum crucibles, and the analysis was performed under a nitrogen atmosphere as an inert gas. The equipment has a temperature uncertainty of ±1 °C and a mass uncertainty of 0.005%, ensuring accurate and precise measurements during the analysis.
Differential scanning calorimetry (DSC) was carried out using a Mettler Toledo calorimeter (model 822e, USA) at a heating rate of 10 K min−1 and 5 to 10 mg samples in aluminum pots under nitrogen as an inert gas. Calibration for temperature and enthalpy measurements was performed using certified indium and zinc standards as a reference, with an uncertainty of ±0.2 K for temperature measurement and 2% for enthalpy change. As a comparison test, the procedure was applied to the base fluids, and the results were consistent with Cp values previously reported in the literature.
Densities were measured with a densimeter Anton Paar DMA 4500 M with an accuracy of 0.00005 g cm−3 and repeatability of 0.00001 g cm−3 in the temperature range of 293.15 to 363.15 K.
The water content of ionic liquids was determined in triplicate by Karl-Fischer titration (model 870 KF Titrino Plus, Metrohm, Switzerland), after drying it at 70 °C under vacuum in the presence of phosphorus pentoxide conditions for 36 h, giving a water content less than 100 ppm.
The mixtures were prepared by combining equimolar amounts of each IL and were stirred for 12 h at room temperature to ensure complete homogeneity. Nine binary mixtures ([A][X][Y] or [A][B][X]) and six reciprocal binary mixtures ([A][B][X][X][Y]) were obtained, according to the classification of Welton et al.47
While density may be a predictable property, thermal stability and heat capacity depend directly on the intermolecular interactions between cations and anions of the mixtures, and being complex systems, further studies and modeling must be done to understand the variation of their properties.
Compared to a previous report with mixtures consisting of eight-carbon alkyl chains, the mixtures with butyl substituents presented, on average, higher density values, but lower decomposition temperatures and lower heat capacity. Also, unlike the previous mixtures, only four mixtures increased their heat capacity (mixtures 11, 13, 14 and 15), all of them binary reciprocals.
Despite the decrease in Cp and thermal stability, the thermal storage density of mixture 11 ([Bmim][BF4] + [Bdmim][Tf2N]), was slightly higher than the previous mixture with higher Cp ([Omim][PF6] + [Odmim][PF6]), mainly due to the relationship with density.
Finally, as with previous mixtures, two mixtures (11 and 14) demonstrated higher thermal storage density than thermal oils used for energy storage, while compared to nitrated salts, the need for increased density or operating temperatures despite higher heat capacity values is evident. Although this study has successfully demonstrated that equimolar mixing of some ionic liquids significantly increases the thermal storage density, further studies are needed to clearly define a possible use as a thermal fluid.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ra02901h |
This journal is © The Royal Society of Chemistry 2023 |