Thermodynamic Study of the Solubility of Acetaminophen in Propylene Glycol + Water Cosolvent Mixtures

Baseado nas equações de van’t Hoff e Gibbs, as funções termodinâmicas energia de Gibbs, entalpia e entropia de solução, mistura e solvatação de acetaminofen em misturas dos solventes propileno glicol + água (PG + W), foram avaliadas por medidas de solubilidade em diversas temperaturas. A solubilidade foi maior para 100% de PG em todas as temperaturas estudadas. A solvatação dessa droga nas misturas aumenta com o aumento da proporção de PG, atingindo um máximo em 70% de PG. De 0% até 20% de PG e de solução 70% até 100% de PG, foi encontrado um domínio da entropia sobre o processo, enquanto de água pura até 10% de PG e de 20% até 70% de PG, foi encontrado domínio da entropia. Estes fatos são explicados em termos de perda da estrutura da água, e uma diminuição na energia requerida para formação da cavidade no solvente, para misturas de 30% até 70% de PG.


Introduction
Acetaminophen is an analgesic and antipyretic drug widely used in modern therapeutics.This drug is specially indicated in the treatment of several minor diseases presented by pediatric patients. 1In the Colombian market, it is commercially available as tablets, syrups and concentrates, but it is not available as parenteral products.The later ones have been recently asked for by physicians and by other care practitioners.Injectable homogeneous liquid formulations supply relatively high doses of drug in small volumes.For this reason, some physicochemical properties such as the solubility and the occupied volumes by the drugs and other components in the solution are very important because they facilitate the design process of pharmaceutical dosage forms. 2 The solubility behavior of drugs in cosolvent mixtures takes great importance because cosolvent blends are frequently used in purification methods, preformulation studies, and pharmaceutical dosage forms design, among other applications. 3Nowadays several methods to calculate the solubility are available.However, these methods do not explain totally the mechanism of cosolvent action in mixtures.On the other hand, almost all of these methods in general do not consider the effect of temperature.For these reasons, it is important to determine, systematically, the solubility of drugs, in order to obtain complete information about physicochemical data of pharmaceutical systems.This information facilitates widely the labor of pharmacists associated to development and research of new products in pharmaceutical industry. 4Temperature-solubility dependence allows to carry out the respective thermodynamic analysis, which, at the same time, permits explain the molecular mechanisms, involved in the solution processes. 5. Braz.Chem.Soc.
The main objective of this study was to evaluate the effect of the cosolvent composition on solubility and solution thermodynamics of acetaminophen in propylene glycol + water cosolvent mixtures.The analysis was based on van't Hoff method, including the respective contributions by mixing and solvation of the drug on the solution processes.Ethanol and propylene glycol are probably the more widely used cosolvents in parenteral medications.This investigation expands the concepts developed for this drug in cosolvent systems by Pérez et al. 2 in ethanol + water, propylene glycol + water, and ethanol + propylene glycol mixtures at 25.0 °C, by Grant et al. 6 in water at several temperatures, by Etman and Naggar 7 in sugar aqueous solutions at 20.0 and 37.0 °C, by Bustamante and coworkers 8 in ethanol + water, ethanol + ethyl acetate, and dioxane + water mixtures at several temperatures, and Martínez 9 in propylene glycol + water mixtures at 25.0 °C, among others.

Solubility determinations
An excess of ACP was added to 20 cm 3 of each cosolvent mixture evaluated in glass flasks.The cosolvent mixtures were prepared by mass in quantities close to 100.0 g varying in 10%, m/m.Solid-liquid mixtures were stirred in a mechanical shaker for 1 hour.Samples were then allowed to stand in water baths kept at the appropriate temperature ± 0.05 ºC.All samples were maintained at least for 48 hours to reach the equilibrium.This equilibrium time was established in a previous investigation 11 about the dissolution rate and solubility of ACP in EtOH, PG and W at 20.0 °C.After this time the supernatant solutions were filtered (at isothermal conditions) to ensure that they were free of particulate matter before sampling.Concentrations were determined by measuring refractive indexes after appropriate dilution and interpolation from previously constructed calibration curves for ACP in each cosolvent mixture. 2 All the solubility experiments were repeated at least three times.In order to make the equivalence between molarity and mole fraction concentration scales, the density of the saturated solutions was determined with a digital density meter.

Results and Discussion
In Table 1, the molecular structure of ACP and some of their physicochemical properties are summarized. 8,12,13The melting point and enthalpy of fusion were reported by Bustamante and coworkers 8 while the enthalpy of sublimation was reported by Williams et al. 13 According to Romero et al. 14 this drug acts in solution mainly as a Lewis acid in order to establish hydrogen bonds with proton-acceptor groups in the solvents (oxygen in -OH groups).Dearden 15 demonstrated that both functional groups of this drug (-NH and -OH) were involved in complex formation with the carbonyl group of antipyrine.ACP could also act as a proton-acceptor compound by means of its carbonyl and -OH moieties.

Ideal and experimental solubility of ACP
The ideal solubility of a crystalline solute in a liquid solvent can be calculated by equation ( 1): (1) where X i 2 d is the ideal solubility of the solute as mole fraction, ∆H fus is the molar enthalpy of fusion of the pure solute (at the melting point), T fus is the absolute melting point, T is the absolute solution temperature, R is the gas constant (8.314J mol -1 K -1 ), and ∆C p is the difference between the molar heat capacity of the crystalline form and the molar heat capacity of the hypothetical supercooled liquid form, both at the solution temperature. 16ince ∆C p cannot be easily determined, one of the following assumptions has to be made: (a) ∆C p is negligible and can be considered zero or (b) ∆C p may be approximated to the entropy of fusion, ∆S fus .In this investigation the later consideration is assumed.
Table 2 summarizes the experimental solubilities of ACP, expressed as molarities and mole fractions, and the ideal solubilities calculated by means of equation ( 1) from ∆H fus , and T fus presented in Table 1.In all cases, the coefficients of variation for solubility were smaller than 2.0%.On the other hand, Figure 1 shows the solubility expressed in mole fraction at all the studied temperatures.In this cosolvent system a maximum in solubility is not obtained in contrast to that found in other cosolvent systems such as ethanol + water (EtOH + W). 8 The Hildebrand solubility parameter (δ) obtained for this drug in EtOH + W mixtures was 28.3 MPa 1/2 (13.8 cal 1/2 cm -3/2 ) at 25.0 °C. 8This value is outside of δ values obtained with PG + W mixtures, i.e., from 30.3 MPa 1/2 (14.8 cal 1/2 cm -3/2 ) up to 47.9 MPa 1/2 (23.0 cal 1/2 cm -3/2 ) .For this reason, the solubility obtained in EtOH + W mixtures is relatively larger in comparison with PG + W mixtures. 9 On the other hand, if molarity is considered, a maximum in solubility is obtained at 90% of PG at all temperatures.a In all cases the coefficients of variation (CV) were smaller than 2.0%. ) by means of van't Hoff equation (equation 2): (2) In more recent treatments, some corrections have been introduced to equation ( 2) in order to reduce the propagation of errors, and therefore, to separate the chemical effects from those due only to statistical treatments used in compensation plots.For this reason, the mean harmonic temperature (T hm ) is used in van' Hoff analysis.T hm is calculated as: 17 (3 where n is the number of tested temperatures.In our case the T hm value obtained was just 303 K.The corrected expression more widely used can be written as follows: 8 (4) As an example, Figure 2 shows the modified van't Hoff plot for ACP in mixtures having 80% and 90% of PG.Linear models with good correlation coefficients were obtained in all mixtures studied.For this reason, the equation ( 4) is useful to estimate the ∆H s 0 o a l p n p values for ACP in the PG + W cosolvent system.
For non-ideal solutions, the slope obtained in equation ( 4) does not give directly the heat of solution.Therefore, it is necessary to consider the variation of solute thermodynamic activity (a 2 ) with concentration at constant temperature and pressure. 8,18Then, the enthalpic change of solution is calculated as: (5)   in which, the second term of the right side is calculated by means of: 8,19 (6) The term "sat" indicates the saturation.In equation ( 6) the solute volumetric fraction (φ 2 ) is required.This property is calculated from the apparent specific volume of solute (ASV 2 ) at saturation, and the mixture composition.ASV 2 is calculated by means of: (7)   where, m 2 and m 1 are the masses of solute and solvent at saturation, respectively, SV 1 is the specific volume of solvent, and ρ is the solution density.Although in a more refined treatment, the partial specific volume of solute instead of ASV 2 should be used, the procedure proposed here is also adequate.
Since ACP is a solid, the thermodynamic activity at saturation equals the ideal solubility (X i 2 d ) 3 and therefore, it follows that: (8)   The term (X i 2 d / X s 2 at ) in equation ( 8) is equal to the solute activity coefficient in the solution (γ 2 ) and it is an indication of the deviation presented by this one in front to ideal behavior.Table 3 shows the experimental % (m/v) solubilities, saturated solution densities, cosolvent mixtures densities, solute volume fractions, solute activity coefficients, and correction factors at 30.0 °C.This temperature is the nearest to 303 K.In order to calculate the γ 2 and (∂ln a 2 /∂ln X 2 ) T,p values some methods for estimating propagation of errors were used. 20rom the γ 2 values presented in Table 3 a rough estimate of solute-solvent intermolecular interactions can be made by considering the following expression: (9) where w 11, w 22 and w 12 represent the solvent-solvent, solutesolute and solvent-solute interaction energies, respectively; V 2 is the molar volume of the supercooled liquid solute, and finally, φ 1 is the volume fraction of the solvent.In a first approach the term( V 2 φ 1 2 / RT) may be considered approximately constant at the same temperature, and then γ 2 depends almost exclusively on w 11, w 22 and w 12 . 21The w 11 and w 22 terms are unfavorable for solubility, while the w 12 term favors the solution process.
It can be seen in equation ( 9) that the contribution of w 22 represents the work necessary to move molecules from solid state to the vapor state, and therefore it is constant in all mixtures.On the other hand, Romero et al. 22 have demonstrate recently by using calorimetric, spectroscopic, and crystallographic techniques, that ACP solid phase in excess keeps its original crystalline properties in saturated solutions in several cosolvent mixtures varying in polarity and Lewis acid-base character.Although, an increase of 8 °C in the melting point has been reported for ACP solid phase at equilibrium with saturated solutions having cosolvent proportions greater than 50 % (v/v), 8 according to these authors, for practical purposes it may be considered that the contribution of solid phase toward the overall solution process is constant for this drug in the different saturated solutions studied.
The term w 11 is higher in water (δ = 47.9MPa 1/2 ) while it is comparatively smaller in PG (δ = 30.3MPa 1/2 ). 23The pure water and water-rich mixtures have larger γ 2 values, which means, high w 11 and low w 12 values.On the other hand, in PG-rich mixtures (with γ 2 values close to 1.0), the w 11 values are relatively low, whereas the w 12 values are higher.According to this fact, the solvation of ACP should be higher in PG-rich mixtures.
The apparent standard Gibbs energy change for the solution process (∆G s 0 o a l p n p ) has been traditionally calculated in literature as: 21 (10)   Nevertheless considering the approach proposed by Krug et al., 13 this property is more appropriately calculated by means of: (11)   in which, the intercept used is the one obtained from ln X 2 vs. 1/T -1/T hm plots (equation 4).This thermodynamic function is also corrected using the factor (∂ln a 2 /∂ln X 2 ) T,P in order to express it in terms of solute thermodynamic activity instead of solute concentration.
The standard entropic change for solution process (∆S s 0 oln ) is obtained from the respective ∆H s 0 oln and ∆G s 0 oln values by using: (12)   Table 4 summarizes the corrected standard thermodynamic functions for experimental solution process of ACP in all cosolvent mixtures including those functions for the ideal process.In order to calculate the thermodynamic magnitudes of experimental solution some methods for estimating propagation of errors were used. 20t was found that the standard Gibbs energy of solution was positive in all cases; i.e., the solution process apparently is not spontaneous, which may be explained in terms of the concentration scale used (mole fraction), where the reference state is the ideal solution having the unity as concentration of ACP, that is, the solid pure solute.
The enthalpy of solution is positive for all cases, therefore the process is always endothermic.The entropy of solution is also positive in all cases, indicating entropy driving on overall the solution processes.The With the aim to compare the relative contributions by enthalpy (%ζ H ) and by entropy (%ζ TS ) toward the solution process, equations ( 13) and ( 14) were employed respectively.
(13) (14)   From Table 4 it follows that in all cases the main contributor to standard free energy of solution process of ACP is the enthalpy (greater than 62% in all cases).

Thermodynamic functions of mixing
The solution process may be represented by the following hypothetic stages: 21 Solute (Solid) → Solute (Liquid) → Solute (Solution)   where, fusion and mixing are the respective partial processes toward the solution process at 303 K.This approximation permits to calculate the partial thermodynamic contributions to solution process by means of equations ( 15) and ( 16). .This value is coincident with the enthalpic change for ideal solution.In contrast, the entropy of fusion at 303 K (59.35 J mol -1 K -1 ) is not coincident with the entropy of ideal solution at this temperature (36.90 J mol -1 K -1 ).Nevertheless, for practical purposes, the ∆S s 0 o i l d n value was used instead of ∆S f 3 u 0 s 3 .In Table 5 the thermodynamic functions of mixing of ACP are summarized.
By analyzing the partial contributions by ideal solution (related to solute fusion process) and mixing processes to the enthalpy and entropy of solution, it is found that ∆H f  4).On the other hand, the contribution of the thermodynamic functions relative to mixing process toward the solution process is variable, that is, ∆H 0 m ix is positive in those mixtures with PG content equal or lower than 40% and negative for all other mixtures, while the entropy of mixing ( ∆S 0 m ix ) is positive in the mixture containing 20% of PG, but negative in all other mixtures.Therefore, the entropies of mixing are in general unfavorable (negative values: Table 5).However, considering the overall solution process (that is, data from Table 4), entropy change is the driven force (positive values: Table 4) because the solution process includes the favorable entropy of melting (positive value: Table 4).
The net variation in ∆H 0 m ix values results from the contribution of several kinds of interactions.The enthalpy of cavity formation is endothermic because energy must be supplied to overcome the cohesive forces of the solvent.This process decreases solubility.On the other hand, the enthalpy of solute-solvent interaction is exothermic and it is originated mainly from the van der Waals and Lewis acid-base a %ζ H and %ζ TS are the relative contributions by enthalpy and entropy toward Gibbs energy of solution.These values were calculated by means of equations ( 13) and ( 14 interactions.The structuring of water molecules around the nonpolar groups of solutes (hydrophobic hydration) contributes to decrease the net heat of mixing to small or even negative values in aqueous solutions.Nevertheless, this is not observed in the case of ACP in water (Table 5).As it was already said, the energy of cavity formation should be lower as the proportion of PG increases because the polarity of the medium decreases, which favors solute-solvent interactions.This fact is partially observed in Table 5, where ∆H 0 m ix is more negative over 50% of cosolvent with a maximum at 70% of PG.According to Romero et al. 14 in the initial portion of the solubility curve, the hydrogen bonding of ACP will increases with cosolvent concentration.At large cosolvent proportions, this interaction may be saturated, becoming a constant contribution.On the other hand, nonspecific and cavity effects are not saturated and vary with cosolvent concentration.
For comparative purposes, Figure 3 shows the thermodynamic functions of mixing, ∆G 0 m ix , ∆H 0 m ix , and Τ∆S 0 m ix .All functions vary nonlinearly with composition showing maxima for enthalpy and entropy at 20% of PG.
In order to verify the effect of cosolvent composition on the thermodynamic function driving the solution process Table 6 summarizes the thermodynamic functions of transfer of ACP from more polar solvents to those less polar solvents.These new functions were calculated as the differences in thermodynamic magnitudes of mixing between the less polar mixtures and the more polar mixtures.As a calculation example, in the case of transfer of ACP from pure water to 10% PG mixture (considering data of Table 5), the enthalpy of transfer (∆H 1 0 →2 ) corresponds to 5.6 kJ mol -1 (∆H 0 m ix in 10% PG mixture) minus 6.4 kJ mol -1 (∆H 0 m ix in pure water) obtaining the value: -0.8 kJ mol -1 (Table 6).All other thermodynamic magnitudes of transfer were calculated on the same way.
If the addition of PG to water is considered, it happens the following: At 10% of PG (∆G 1 0 →2 < 0 and ∆H 1 0 →2 < 0), the solution process is driven by enthalpy, which it is not easily explained.For ACP solubility in EtOH + W mixtures near to this composition the solution process was driven by entropy. 8t 20% of PG (∆G 1 0 →2 < 0 and ∆H 1 0 →2 > 0) the solubility process is driven by the increment of entropy (∆S 1 0 →2 >0).This suggests a decrease of the hydrophobic effect of water, produced by the addition of 20% of cosolvent.That is, the cosolvent disrupts the ordered structure of water (like "icebergs", waterwater hydrogen bonds) around de nonpolar groups of the drug, increasing the enthalpy and entropy of the system.From 30% up to 70% of PG (∆G 1 0 →2 < 0 and ∆H 1 0 →2 < 0) the solution processes are driven by enthalpy; this fact is probably due to solvation of ACP.In the mixture of 40% of PG in addition to enthalpy of mixing the entropy of mixing also contributes to solution process.Finally, from 80% up to 100% of PG (∆G 1 0 →2 < 0 and ∆H 1 0 →2 > 0) the entropy of mixing drives the solution process, because this property is positive and therefore it is favorable.

Thermodynamic functions of solvation
In addition to the hypothetic fusion-mixing stages previously exposed, the solution process may also be represented by the following hypothetic stages: 24  13) and ( 14), respectively.
where, the respective partial processes toward the solution process, are in this case, sublimation and solvation.This treatment permits calculate the partial thermodynamic contributions to solution process by means of equations ( 17) and ( 18), respectively, while the Gibbs energy of solvation is calculate by means of equation ( 19): where, ∆H s 0 ubl = 238.85kJ mol -1 was taken from Williams et al., 13  , and therefore ∆S s 0 ubl = 577.7 J mol -1 K -1 at the same temperature.In Table 7 the thermodynamic functions of solvation are presented, while on the other hand, with the aim to compare the relative contributions by enthalpy (%ζ H ) and entropy (%ζ TS ) toward the solvation process, two equations analogous to equations (13) and ( 14) were employed.
From the values of %ζ H and %ζ TS presented in Table 7 it follows that the main contributing force to standard Gibbs energy of the solvation process of ACP in all the cosolvent mixtures is the enthalpy (%ζ H are than 56% in all cases).
Because that not only the main driving force of solvation process of drug compounds is important, but also the balance between specific and non-specific solutesolvent interactions as well, therefore, parameters which describe the relative ratio of specific and non-specific solute-solvent interaction in terms of enthalpies (%ε H ) and in terms of entropies (%ε S ), were used according to the following definitions introduced by Perlovich and coworkers: 24  (25)   Cyclohexane was chosen as an "inert" solvent, which interacts with drug molecules solely by nonspecific interactions (dispersion forces), while the cosolvent mixtures interact with ACP by specific interactions such as hydrogen bonding.Benzene and hexane have also been used as inert solvents in the study of naproxen although important differences have been found between these two solvents, indicating some effect of π electrons and planar geometry of benzene on non-specific interactions of that drug. 24olubility data for ACP in cyclohexane taken from Baena et al. 25 were analyzed according to equations (4),  (CH) = 41.2J mol -1 K -1 .The apparent specific volume of ACP in cyclohexane obtained by using densities of solvent and saturated solutions was a negative value (due to very scarce solubility and uncertainty in density measurements).For this reason, in order to calculate the (∂ln a 2 /∂ln X 2 ) T,P for ACP in this solvent, the molar volume of drug was calculate by means of Fedors method 26 obtaining a value of 124.4 cm 3 mol -1 .From this value and the solubility at 303 K, the value obtained for (∂ln a 2 /∂ln X 2 ) T,P using equation (8) was 0.9994.Since this value is included into the uncertainty obtained in thermodynamic functions of solution, then, the apparent values were used instead of corrected values.
The %ε H and %ε S values for ACP solvation are also presented in Table 7.These values indicate that during dissolution of ACP in all mixtures studied, the specific solute-solvent interactions (hydrogen bonding, mainly) do not affect the entropic term of free energy with respect to non-specific interactions.With regard to the enthalpic term in all cases the non-specific solute-solvent interactions predominate.

Enthalpy-entropy compensation of solution
Bustamante et al. 8 have demonstrated some chemical compensation effects for the solubility of several drug compounds in aqueous cosolvent mixtures.This analysis was used in order to identify the mechanism of the cosolvent action.The making of weighted graphs of ∆H s o ln as a function of ∆G s 0 oln at mean harmonic temperature permits to observe similar mechanisms for the solution process according to tendencies obtained. 27r solubility of ACP in EtOH + W, Bustamante and coworkers 8 obtained a nonlinear trend using seven cosolvent compositions, including the pure solvents.Their data were adjusted to a parabolic regression model obtaining a maximum for 20% v/v of EtOH.From 0 up to 20% v/v of EtOH a negative slope was obtained while over this EtOH proportion a positive slope was obtained.According to these authors, this fact implies a change from entropy driving to enthalpy driving toward the solution process.
On the other hand, Figure 4 shows the corrected ∆H s 0 oln and ∆G s 0 oln values for solubility of ACP in PG + W mixtures obtaining a trend slightly more complex compared with that presented by Bustamante and coworkers 8 for EtOH + W mixtures.
Figure 4 shows fully that this solute-cosolvent system does not present linear ∆H s 0 oln -∆G s 0 oln compensation in all compositions studied.Nevertheless, if an interval from 20% up to 80% of PG is considered an apparent linear trend is observed with positive slope (without considering 30% of PG).According to this graph it follows that from pure water up to 10% of PG and from 20% up to 70% of PG the dominant mechanism for solubility is the enthalpy, while from 10% up to 20% of PG and from 70% up to 100% of PG the dominant mechanism is the entropy.

Conclusions
From all aspects discussed previously it can be concluded that the solution process of ACP in PG + W mixtures is very complex and highly dependent on cosolvent composition.The solvation of this drug is greater for PGrich mixtures especially at 70% of PG.In a similar way to that found for the solubility of this drug in EtOH + W mixtures, the solution process in PG + W mixtures does  13) and ( 14), respectively; b %ε H and%ε S are the relative ratio of specific and non specific solute-solvent interactions expressed in terms of enthalpy and entropy.These values were calculated by means of equations ( 20) to (25).J. Braz.Chem.Soc.
not follows linear enthalpy-entropy compensation using ∆H s 0 oln as a function of ∆G s 0 oln .This fact at its time is explained as entropy or enthalpy driving for solution process according to the cosolvent mixture composition.Finally, it can be said that all values presented amply the physicochemical information useful in the design of homogeneous liquid pharmaceutical dosage forms.

Figure 1 .
Figure 1.Solubility of ACP expressed in mole fraction in PG + W cosolvent mixtures at several temperatures
∆C p (T fus -T) using ∆S f M u P s instead of ∆C p obtaining a value of 17.98 kJ mol-1
and therefore, the function ∆H s 0 olv was calculated from ∆H s 0 oln values presented in Table 4.The respective entropy of sublimation was calculated as ∆S s T at 303 K, where ∆G s 0 ubl = -RT ln(p / p 0 ) with p = 1.05×10 -6 Pa at 303 K (calculated from some values presented by Williams et al. 9 ) and p 0 = 101325 Pa; then s 0 = 63.71kJ mol -1

Figure 4 .
Figure 4. Enthalpy-entropy compensation plot for solubility of ACP in PG + W cosolvent mixtures at 303 K.

Table 1 .
Some physicochemical properties of ACP a Taken from Ref. 12; b Taken from Ref. 8; c Taken from Ref. 13

Table 2 .
Experimental solubility of ACP in PG + W cosolvent mixtures expressed in mol L -1 and mole fraction including ideal solubility at several temperatures

Table 3 .
6olubility of ACP expressed in % (m/v), saturated solution and solvent densities, solute volumetric fraction, solute activity coefficient, and activity variation factor in PG + W cosolvent mixtures at 30.0 °C water is in good agreement with those presented by Grant et al.6and Bustamante and coworkers, 8 that is, 23.7 and 22.5 kJ mol -1 , respectively.The

Table 4 .
Thermodynamic functions relative to solution process of ACP in PG + W cosolvent mixtures including ideal process at 303 K

Table 5 .
Thermodynamic functions relative to mixing process of ACP in PG + W cosolvent mixtures at 303 K and %ζ TS are the relative contributions by enthalpy and entropy toward Gibbs energy of mixing.These values were calculated by means of equations (

Table 6 .
Thermodynamic functions of transfer of ACP from more polar solvents to less polar solvents in PG + W cosolvent mixtures at 303 K

Table 7 .
Thermodynamic functions relative to solvation process of ACP in PG + W cosolvent mixtures at 303 K PG / (%, m/m) ∆G 0 s olv / (kJ mol -1 ) ∆H 0 and%ζ TS are the relative contributions by enthalpy and entropy toward Gibbs energy of solvation.These values were calculated by means of equations (