Thermodynamic investigation on the aqueous mixtures of choline chloride/propylene glycol deep eutectic solvent at T = (293.15 to 313.15) K

This work focuses on physicochemical features of the choline chloride/propylene glycol deep eutectic solvent (DES) + water mixtures by determining their density values in mono- and mixed- states over various temperatures in the range of 293.15–318.15 K. The density data obtained from measurements were utilized for the computation of various quantities such as excess molar volumes, molar volume, apparent molar volume, limiting apparent molar expansibility, and isobaric thermal expansion coefficient. Furthermore, the experimental densities were fitted to some mathematical equations such as Jouyban-Acree, Jouyban-Acreevan’t Hoff, modified Jouyban-Acree-van’t Hoff, Redlich–Kister and Emmerling. Studies of this nature can provide useful insights into solute–solvent interactions in aqueous solutions of DES, especially about to their novel application in drug solubilization. Supplementary Information The online version contains supplementary material available at 10.1186/s13065-024-01153-y.


Introduction
Deep eutectic solvents (DESs) as green solvents are an alternatives class for ionic liquids and are synthesized using a mixture of two or more organic substances.The combination of components in DESs results in a homogeneous mixture with a melting point that is lower than that of each component [1].DESs are structured based on hydrogen bonding, where a hydrogen bond acceptor and hydrogen bond donor components interact with each other [2].The main benefits of DESs are their biodegradability, cost-effectively, non-flammability and no further purification is required for their preparations [3,4].Moreover, DESs have a significant effect on the solute solubilization and enhance its aqueous solubility [5].DESs are mostly based on choline chloride (ChCl), carboxylic acids, and other hydrogen bond donors e.g.succinic acid, urea, glycerol, and citric acid [6].Before a solvent selection for any purpose, it needs to obtain some essential information on the targeted solvent including its volumetric features.Knowing these properties for DESs and their mixtures with water provides some helpful information on the solute-solvent interactions.Excess molar volume as a volumetric property demonstrates the non-ideality of a system and explains the intermolecular interactions [7][8][9].
Conducting such studies can yield useful insights regarding solute-solvent interactions in DESs' aqueous solutions, which have recently been employed for drug solubilization purposes.To date, there was some data on the thermodynamic behavior of an aqueous mixture of DESs including ChCl/ethylene glycol [10], ChCl /urea [11], ChCl/malonic acid [12], ChCl/phenol [13], ChCl/ glycerol [14], ChCl/glucose [15], ChCl/oxalic acid [16], ChCl/triethylene glycol [17], halide salts/ethylene glycol [18], and amino acids/ lactic acid [19].However, there is no data for the physicochemical properties of ChCl/ propylene glycol (PG) at various temperatures in the literature.The goal of this work is to expand the database for physicochemical and thermodynamic properties of DESs + water mixtures at different temperatures for use in the solubilization procedures.It should be noted that the exact choice of solvent system will depend on various factors such as drug solubility, compatibility, and stability.The ChCl/PG + water system may not be suitable for all drugs and formulations, and alternative solvent systems may be required in some cases.However, regarding the comparison with other glycols, such as ethylene glycol and polyethylene glycols (PEGs), choosing the PG in the structure of DES has some advantages including: (i) ethylene glycol is toxic and not suitable for pharmaceutical applications.It is primarily used as an industrial solvent and antifreeze.Therefore, it is not a preferred choice for pharmaceutical formulations due to its potential health risks [20].(ii) PEGs are polymers of ethylene glycol and are commonly used in pharmaceutical formulations as a solubilizing agent, stabilizer, and viscosity modifier [21].However, compared to PG, PEGs may have higher viscosity and can be less efficient in solubilizing certain drugs.The specific ratio of 1 part ChCl to 3 parts PG was chosen for this study.The reason is that this composition has been reported as an effective DES for drug solubilization purposes in previous studies [22].The DES system is formed by mixing ChCl and PG in a specific molar ratio, resulting in a lower melting point than either component alone.This lower melting point allows for the solubilization of drugs that are poorly soluble in traditional organic solvents.The 1:3 ChCl:PG ratio was selected as it provides a balance between solubilization efficiency and solvent viscosity, which is important for drug formulation and processing.

Materials
This work utilized PG sourced from Scharlau Chemie (Spain), with a mass fraction purity greater than 0.995, ChCl obtained from Daejung (Korea) with a mass fraction purity greater than 0.999, and double-distilled water produced in the laboratory (Table 1).

Preparation of ChCl/PG DES
Drying of ChCl was carried out in an oven at a temperature of 50 °C for 8 h.Subsequently, ChCl and PG were weighed in the desired amounts for preparing the solvent mixture, at a molar ratio of 1:3.The two components were then mixed and heated on a hot plate at a temperature of 80 ͦ C while stirring with a magnetic stirrer until a clear liquid was formed [23].

Solvent mixtures preparation and determination of their densities
After DESs preparation, the eleven mixtures i.e. two mono-solvents and nine solvent mixtures in the DES ratio of 0.1-0.9 with a mass fraction interval of 0.1 were prepared and incubated at 293.2-313.2K in a Nabziran Industrial Group incubator (Iran) until to reach temperature equilibrium.A 2 mL pycnometer was utilized for density measuring for each solution.Calibration of the pycnometer involves using the densities of distilled water at the respective temperatures used [24].The pycnometer was filled with each prepared solution at desired temperatures and measured using a balance (model AB204-S, Mettler Toledo, Switzerland).With considering the weight of the empty pycnometer, the weight of the solvent mixture was obtained and used for density computing.The reported data were the mean of three replications at each temperature and composition.

Density data and mathematical modeling
The measured densities (ρ m ) at a temperature range of 293.15 to 313.15 K for the ChCl/PG DES and water mixtures in all compositions were given in Table 2. Figure 1 illustrated the density plot of ChCl/PG DES + water as a function of temperature, enabling clarity in visualization.Reliability of measured data using a pycnometer was investigated by repeating measurements using a density meter (Anton Paar DSA 5000 M).Density values for ChCl/PG DES + water mixtures at 298.2 K were determined using the density meter and results were given in electronic supplementary material (ESM) in Additional file 1:Table S1.As can be seen, there was no significant difference between the results of both measurements which confirmed the data provided by the pycnometer.Table 2 and Fig. 1 demonstrated that the densities rise with an increase in the mass fraction of DESs, while they decreased with an increase in temperature.The decrease in density can be attributed to the thermal expansion of the liquid volume with a temperature increase, leading to a less dense mixture at higher temperatures.
Equation (1) presents the temperature dependence of density values in the studied solutions, using a polynomial expression.
Table 3 showed the model constants (a, b, and c) obtained from the least squares analysis for temperature T/K.Furthermore, the corresponding molar volume was calculated using the density of neat DES with the following equation: where M DES and ρ DES are the molar mass of DES and its density.The molar mass is computed by Eq. ( 3) [25].
where x HBA , x HBD , M HBA and M HBD are molar ratio and molar masses of ChCl as HBA and PG as HBD, respectively. (1) x HBA M HBA + x HBD M HBD x HBA + x HBD Table 2 The density data (g•cm -3 ) for a binary mixtures of ChCl/PG DES + water at different temperatures a   Table 4 showed that the variation of V m with tem- perature was comparable across the DES samples, with V m values exhibiting a slight increase as temperature increases.
Jouyban-Acree [26,27] (Eq.( 4)), and Jouyban-Acreevan't Hoff [28] (Eq.5) models were used to represent the densities of pseudo-binary mixtures of ChCl/PG DES at different temperatures and each model constants together with the ARDs% and SDs were collected in Table 5. in these equations, ρ m , ρ DES , ρ w are the density of mix- ture, neat DES and neat water, respectively.x DES , x w correspond to the mole fraction of DES and water, respectively.J i is the model constants.
To assess the precision of each model in depicting density values in the mixtures, the average relative deviation (ARD%) (Eq.6) and standard deviation (SD) (Eq.10) were calculated.
in Eqs. ( 6) and (7) m and N are the experimental densities, the mean experimentally densities, the back-computed densities from different equations and the data point numbers, respectively.

Excess molar volume and data correlation
To evaluate the non-ideality of (ChCl/PG DES + water) pseudo-binary mixtures, the V E values were calculated using Eq. ( 8).
Tables 6 reported the V E values found for the mixtures investigated.These values were fitted with Eq. ( 9) and the corresponding results were presented in Table 7.Additionally, Table 7 displayed the correlation outcomes of V m values with Eq. ( 4).
Figure 2 displayed the experimentally and calculated V E values as a function of composition and temperature, allowing for a visual interpretation of their behavior.The negative trend observed throughout the range of composition and temperature indicated that the volume of mixtures was lower than that of an ideal mixture due to the strong interactions between DES and water molecules, likely attributed to hydrogen bonding between ChCl and PG with water molecules.This observation was consistent with previous findings for other pseudo-binary mixtures of (DES + water) [11-13, 19, 33, 34].Additionally, the V E values became less negative at higher tempera- tures in the measured system, indicating that the intermolecular interactions between DES and water were weakened at elevated temperatures.This may result from an increase in hydrogen bond strength due to decreased ( 8)

Apparent properties
The density values acquired for the aforementioned mixtures were utilized to calculate the apparent molar volume of DES ( V ϕ,DES ) in a water medium using the fol- lowing equation: The calculated values of ( V ϕ,DES ) were collected in Table 8.It was evident that the values of V ϕ,DES were affected by the DES mole fraction and temperature.As both DES mole fraction and temperature increase, the values of V ϕ,DES also increase.The dependency of V ϕ,DES values on molality under isothermal conditions was fitted to the Redlich-Mayer equation [35,36], which can be expressed as: Equation ( 11) expresses V 0 ϕ,DES (cm 3 mol −1 ) and A υ (cm 3 kg 1/2 mol −3/2 ) as the apparent molar volume in infinite dilution and the Debye-Huckel limiting slopes for the apparent molar volume, respectively.In addition,.B υ (cm 3 kg mol −2 ) and C υ (cm 3 kg 3/2 mol −5/2 ) are considered as empirical constants that depend on the solute, solvent, and temperature.m DES represents the molality of DES present in the mix- tures.The parameters of Eq. ( 11) were calculated for DES in its aqueous pseudo-binary mixtures using the least square fit (10) Table 6 The dependence of V E values on DES mole fractions at operational temperatures for pseudo-binary mixtures of ChCl/PG DES (molar ratio of 1:3) in water a Standard uncertainties (u) are u(T) = 0.1 K, u(p) = 0.5 kPa, u(x DES ) = 0.005 and the average combined expanded uncertainties   )  )) method, and the outcomes were depicted in Table 9.According to the presented results in Table 9, it was evident that the values for DES in water were positive, indicating the existence of robust interactions between DES and water.Alternatively, the rise in V 0 ϕ,DES values at higher temperatures can be explained by examining the primary and secondary solvation layers.An augmentation in temperature results in the release of solvent molecules from the secondary hydration layer of solutes into the bulk of the solvent, leading to the expansion of the primary layer.This was reported in recent studies for numerous systems [19,37,38].In comparison to the B υ values, the higher V 0 ϕ,DES values suggest that the interactions between DES and water were stronger than those of the self-interactions of molecules (such as DES-DES or water-water) across all temperatures.
Apparent molar volume was utilized to determine the partial molar volume ( V DES ) of ChCl/PG DES in water medium by the following equation [39]: Table 8 provided the computed V DES values, which exhibit the same pattern as the V ϕ,DES values.
Two crucial parameters were calculated using the first and second derivatives of V 0 ϕ,DES values with regard to temperature, i.e. the limiting apparent molar expansibility ( φ 0 E,DES ) of neat DES and Hepler's constant ( ) with the Eqs.( 14) and (15).The results were also reported in Table 9.
Table 9 displayed the favorable φ 0 E,DES values for ChCl/ PG DES aqueous pseudo-binary solutions at the operating temperature, which indicate solvation and electrostriction of solutes in an aqueous environment [40].This observation could be explained by the speedy transfer of water molecules from the DES hydration (12) layers to the neat water.Additionally, the positive values of Hepler's constant revealed that ChCl/PG DES has a structure-making effect on the water medium at the operating temperatures, as mentioned in reference [41].
Understanding the impact of temperature on liquid expansion is crucial for equipment design.As the temperature of a liquid increases, its volume also increases due to the kinetic energy of its molecules.This phenomenon is known as thermal expansion.The amount of expansion is proportional to the temperature change and the initial volume of the liquid.This relationship is expressed by the coefficient of thermal expansion, which is a measure of how much a material expands per unit length for a given temperature change.This information can be obtained through the isobaric thermal expansion coefficient ( α p,m ) calculated using Eq. ( 16).Table 9 indicates that the α p,m values declined as the temperature increased.

Conclusion
The densities of ChCl and PG (at a ratio of 1:3) deep eutectic solvents were experimentally determined over a range of temperatures, and based on these measurements, certain physicochemical properties such as V E , φ 0 E,DES , φ 0 E,DES , α p,m , and were calculated.The fact that V E values were negative for the entire range of mole fractions and temperatures in the aqueous pseudobinary DES mixtures suggests that DES and water molecules interact favorably with each other.Additionally, an increase in the temperature resulted in an increase in the limiting apparent molar volume ( V 0 ϕ,DES ), which further supports the idea that favorable DES-water interactions occur in these mixtures.In addition, the fact that Helper's constant had positive values indicates that DES functions as a structure-maker in a water-based medium.Furthermore, several models, including Jouyban-Acree, Jouyban-Acree-van't Hoff, modified Jouyban-Acree-van't Hoff, Redlich-Kister, and Emmerling were utilized to represent the experimental densities.Based on the results, the order of the models' effectiveness in density correlation was determined as follows:

a
Standard uncertainty (u) for pressure, temperature, mole fraction of DESs and density is u (P) = 0.5 kPa, u (T) = 0.1 K, u (x DES ) = 0.005, u (ρ m ) = 0.05 kg m −3 , respectively b x DES is the mole fraction of ChCl/PG DESs dissolved in water c Standard uncertainty (u) for DES composition was estimated to be less than 0.05 mol ratio x DES b,c

Fig. 1 Table 3
Fig. 1 Density-temperature plot for a binary mixture of ChCl/PG DES and water

Table 1
Some details of the purity and chemical structure of the employed materials a a The purity

of the employed chemicals was provided by the suppliers Material Mass fraction purity a Source Chemical formula Molar mass/g mol −1
H O H

Table 4
The molar masses and volumes of DES measured at various temperatures

Table 5
Model constants and the ARDs% along with SD for the densities of pseudo-binary (ChCl/PG DES + water) mixtures

Table 7
The model constants of Redlich-Kister obtained by fitting the V E values for the investigated solutions at different temperatures

Table 9
Limiting apparent molar volume of DES ( V 0