Next Article in Journal
Experimental Study of Entrainment and Mixing of Renewable Active Particles in Fluidized Beds
Previous Article in Journal
Reliability, Availability, and Maintainability (RAM) Study of an Ice Cream Industry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Accuracy of Molar Solubility Prediction from Hansen Parameters. An Exemplified Treatment of the Bioantioxidant l-Ascorbic Acid

by
Ralph J. Lehnert
1,2,* and
Richard Schilling
2,3
1
School of Chemistry, Reutlingen University of Applied Sciences, 72762 Reutlingen, Germany
2
Reutlingen Research Institute, Reutlingen University of Applied Sciences, 72762 Reutlingen, Germany
3
Lehr- und Forschungszentrum für Prozessanalyse und Technologie (LFZ-PA&T), Reutlingen University of Applied Sciences, 72762 Reutlingen, Germany
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(12), 4266; https://doi.org/10.3390/app10124266
Submission received: 11 April 2020 / Revised: 12 June 2020 / Accepted: 17 June 2020 / Published: 22 June 2020
(This article belongs to the Section Chemical and Molecular Sciences)

Abstract

:
Estimating molar solubility from the Hildebrand-Scott relation employing Hansen solubility parameters (HSP) is widely presumed a valid semi-quantitative approach. To test this presumption and to determine quantitatively the inherent accuracy of such a solubility prognosis, l-ascorbic acid (LAA) was treated as an example of a commercially important solute. Analytical calculus and Monte Carlo (MC) simulation were performed for 20 common solvents with total HSP ranging from 14.5 to 33.0 (MPa)0.5 utilizing validated material data. It was found that, due to the uncertainty of the material data used in the calculations, the solubility prediction had a large scattering and, thus, a low precision. Prediction power is most adversely affected by the uncertainty of the HSP estimates (solvent and solute), followed by the solute heat of fusion. The solute melting temperature and molar volume have minor effects. Computed and experimental solubilities show the same qualitative behavior, while quantitative discrepancies reach one to three orders of magnitude. Solubility estimates were found to provide, at best, rough guiding information but, with the quality of material data on LAA available, they cannot be rated semi-quantitative. It is assumed that these results generally apply at least to solute-solvent systems with a material data quality and solubility similar to LAA.

Graphical Abstract

1. Introduction

The production and processing of chemicals normally involves solute-solvent systems at some stage. Thus, solubility plays a key role in all sorts of industrial applications such as synthesis, extraction, or purification, as it also does, for instance, in packaging, painting, food technology, or pharmacology. In commercial environments, attention is placed on rational product design and process intensification. The development of solid amorphous dispersions in the pharmaceutical industry (e.g., drug containing polymers), the functionalization of food packaging films with antioxidants, or the formulation of multicomponent solvents for membrane production, separation purposes, or other applications, are examples of practical problems that benefit clearly from the possibility of reliable solubility prediction. Consequently, it is of great industrial value to predict or, at the least, reasonably estimate the solubility from basic physicochemical solute and solvent property data.
A frequently adopted approach using the solubility parameter concept has been proposed early on for nonpolar solvents by Hildebrand and Scott [1,2,3,4]. Despite its various assumptions, this approach is considered useful in many practical areas, not least because of its conceptual simplicity, even though its predictive accuracy has not been quantitatively assessed in the literature. This study evaluates the numerical uncertainty typically encountered when applying the Hildebrand-Scott approach to solubility prognosis of the commercially important natural antioxidant l-ascorbic acid (LAA). The equation giving the mole fraction solubility x1 of a solid solute reads in its basic form for a regular solution [5,6]:
ln ( x 1 ) = Δ H 1 R · (   1 T m , 1 1 T   ) V 1 · Φ 2 2 · ( δ 1 δ 2 ) 2 R T
where ΔH1 the heat of fusion in units of J/mol, Tm,1 the melting temperature, and V1 the molar volume; for the solute, Φ2 the solvent volume fraction, and the Hildebrand solubility parameters are δ1 (solute) and δ2 (solvent). R and T have their usual meaning. Throughout this paper, subscript “1” (“2”) designates the solute (solvent). The derivation of Equation (1) rests on assuming a vanishing term containing the molar heat capacity difference of the solid and supercooled liquid state of the solute.
Both regular solution theory as well as the Hildebrand solubility parameter concept originally treat the mixing of two liquids. To extend this approach to dissolving crystalline solids in liquids at T < Tm,1 the molar volume V1 of the (hypothetical) subcooled liquid solute needs to enter Equation (1). Generally, knowledge of the value of V1 necessitates approximations, for instance, by group contribution methods [7] or by the extrapolation of solute melt density data to the solution temperature. Applying Equation (1) to various solute-solvent systems, it was rated vaguely to provide correlations or predictions that were “quite good,” “fair,” or “unreliable” for nonpolar, polar aprotic, and, protic solvents, respectively [5] (p. 344). This rating stays unchanged if the Hildebrand solubility parameter δ is replaced by the totalHansen solubility parameter (HSP) as defined by δ t = ( δ d 2 + δ p 2 + δ h 2 ) 0.5 with, in extension to Hildebrand, now three contributions to the material cohesive energy accounting for nonpolar dispersive, polar, and hydrogen bonding interactions (subscripts ‘d’, ‘p’, and ‘h’). The HSP concept has been explained in detail in the literature [8,9]. In principle, the theoretical basis for Equation (1) limits its use to nonpolar systems, but any remedy for this drawback makes necessary a considerably more complex treatment leading to intricate equations. While genuine quantitative agreement with experimental solubility data is rather rare, the above relation is still widely believed to prove its practical worth.
Experimental values for the parameters V1 and δ1, commonly originate from solubility data in some mixtures and are available for comparatively few compounds only. Therefore, these values, instead, have to be estimated with inevitable but usually unassessed error. Outlines of estimation methods as well as widely accepted δ2 values for numerous solvents can be found in the literature [5,7,8,9]. While Equation (1) is still relatively simple, there are significant data requirements for this approach relating to thermal (ΔH1, Tm,1), volumetric (V1), and cohesive properties (δ1, δ2). Besides its widespread use and economical bearing, LAA was chosen in this study as a paradigmatic compound for its low solubility (x1 ≤ 10−3) in numerous solvents. Such low solubilities ensure that the solvent volume fraction Φ2 can safely be approximated by unity. Thus, iterative numerical approximations (e.g., the method of successive approximations [10]) do not need to be performed so that x1 can straightforwardly be computed from Equation (1). Elimination of iterative optimization procedures in the calculation reduces numerical error and makes it possible to easily apply Monte Carlo (MC) methods, as described later, without the need for specific programming or dedicated software. The literature known to the authors does not provide quantitative information for LAA or other compounds on the error margins that occur when deriving solubility from HSP.
The present paper aims to deliver solubility predictions for LAA, including a comprehensive, exemplified treatment of the aforementioned accuracy issue. Another objective is to ascertain the fundamental uncertainty of estimating molar solubility following the concepts underlying the ISO Guide to the Expression of Uncertainty in Measurement (GUM) [11]. Available material data on LAA were scrutinized and thoroughly consolidated to provide values and their respective uncertainties for all relevant material quantities. From these, confidence intervals (CI) for the resulting solubility estimates of LAA are derived, the influences of the engaged quantities on the prediction error are investigated and, finally, the results are compared with available experimental data.

2. Methods

2.1. Determination of Molar Volume

Pure melt, solution, particle, and crystallographic literature data of LAA were assessed to derive an estimate for the molar volume V1 of the (hypothetical) supercooled liquid at 25 °C, the usual reference temperature for HSP data. The LAA molecular weight was taken as 176.124 g.
In addition, because the solubility parameter is a function of molar volume ( δ V 0.5 ) , the molar volumes, V2, of all solvents investigated here were calculated from experimental densities listed in reference [12] and compared with the values compiled in HSP standard works [5,8,9]. The purpose was to check to which temperatures the data in [5,8,9] actually refer to. In all cases where the solvent molar volume did not relate to measurements at 25 °C or 20 °C, the literature was searched to obtain data that are more appropriate.

2.2. Determination of Solubility Parameter

In line with common practice, the extension of the Hildebrand parameter, the total Hansen solubility parameter δt is taken to enter Equation (1) throughout. For the solute, δ1 is calculated from additive group contributions according to the scheme by Hoftyzer and van Krevelen (HVK) using the molar volume as determined in this study and referencing the SMILE (simplified molecular-input line-entry) structure OC[C@H](O)[C@H]1OC(=O)C(O)=C1O [7,9]. For the solvents, δ2 is computed from the Hansen parameters compiled in [9]. In case significant modification to the molar volume is applicable, δ2 is additionally calculated by the HVK method to verify plausibility. Solvent data used for the calculations are compiled in Table S1.

2.3. Determination of Solute Melting Parameters

Published data on the heat of fusion, ΔH1, and the melting temperature, Tm,1, of LAA were evaluated to deduce the most probable values and minimum-maximum ranges for both quantities.

2.4. Solubility Calculations

Mole fraction solubilities x1 were computed from Equation (1) using Φ2 = 1. In addition, GUM-compliant simulations were carried out within the respective uncertainty ranges of the material variables relevant here [11]. While maintaining the usual pattern of the Monte Carlo method, the following four consecutive steps were performed: (1) an input domain of five input variables (ΔH1, Tm,1, V1, δ1, δ2) was defined, (2) values for each variable were randomly generated from a uniform probability distribution (boxcar) over the respective variable range (mid-range value ± ½∙(maximum-minimum)), (3) molar solubilities x1 were deterministically calculated according to Equation (1) and using the input values from step (2), and, finally, (4) the calculated solubility values were aggregated and analyzed. For each solvent, the quartiles and 95% confidence limits (LCL, UCL) of x1 were numerically derived from 1000 random variable quintuples. The solution temperature was set to T = 25 °C for all calculations. The total solubility parameter and molar volume of LAA were both used as determined in this work. All computations and simulations were performed using a commercial spreadsheet (Excel, v. 2016, using the solver and analytical functions add-ins, Microsoft Corp., Redmond, WA, USA). Specific coding was not required. It is emphasized, that no molecular dynamics calculations were performed and that the term “Monte Carlo simulation” refers to the determination of the uncertainty of the molar solubilities calculated within the GUM framework as above.

3. Results and Discussion

3.1. Molar Volumes

Solute. Literature in the HSP context persistently quotes a molar volume of 124.85 cm3 for LAA. Tracing the source, it becomes evident that this value originates from melt density data that were interpolated to 192 °C [13,14]. Apparently, this volume refers to the amorphous-liquid state at a temperature far beyond 25 °C, suggesting substantial overestimation. To arrive at a less questionable figure, more reliable data were searched. From X-ray crystallographic unit cell dimensions by Hvoslef, a value of 103.67 cm3 at room temperature was derived [15]. Particle density at 25 °C gave 106.74 cm3 with the particle purity and crystallinity both unknown [12,14]. As HSP values generally relate to the liquid state, it is instructive to look also at the LAA partial molar volume (apparent molar volume at infinite dilution) in the solution. By linearly extrapolating aqueous LAA solution density to zero molality at 25 °C, one obtains partial molar volumes ranging from 102.38 to 106.49 cm3 (n = 9) with a median (mean) of 105.40 (105.04) cm3 [16,17,18,19]. All these data point at a molar volume of about 104 cm3 for highly ordered solid LAA at 25 °C.
In [13,14], the temperature dependence of melt density above 192 °C is approximated by a modified Rackett equation. Extending this relation way down to 25 °C yields an estimate for the molar volume of supercooled liquid LAA of 109.85 cm3, which is about 6% larger compared with the crystalline state. Since the coefficients of the approximate formula lack any indication as to their accuracy, this value should be considered tentative albeit, for instance, analogous polymer data suggest it to be reasonable. Plausibility is corroborated by the fact that, at 25 °C, amorphous (rubbery or glassy) polymers possess molar volumes that exceed their crystalline counterparts typically by up to 10% [7] (pp. 77–85). The foregoing results are summarized in Table 1, and it is proposed that 109.9 cm3/mol might be taken as the most probable value for V1 (25 °C), while presuming an error of about ±4% seems reasonable.
Solvents. All experimental solvent data published [5,8,9] relate to liquid phases at temperatures of either 20 °C or 25 °C, except methyl hydroperoxide. Values at 20 °C were also considered adequate since any discrepancies that might arise were rated insignificant within the scope of this study. The only experimental information on methyl hydroperoxide relates to its liquid density at 15 °C. Reference [12] lists 1.9967 g/cm3, whereas the traced original source [20,21] specifies 0.9967 g/cm3, suggesting a typographic error on data transfer in the past. This drastically impacts V2 (15 °C), which is 48.2 cm3 instead of the commonly cited 24.1 cm3. For 20 °C, only the calculated theoretical density (0.982 g/cm3) is available, leading to V2 (20 °C) = 48.9 cm3 [21]. In the absence of further data, this value was adopted for subsequent evaluation.

3.2. Solubility Parameters

Solute. The LAA solubility parameter of the subcooled liquid state was computed for two different molar volumes (temperatures) according to the HVK scheme using the commercial software HSPiP [9]. At 20 °C and 25 °C, the resulting HSP components, as well as δ1, were the same within less than 0.3% error, so that δ1 (20 °C) ≈ δ1 (25 °C) = 35.3 (MPa)0.5 will be used further. To derive an estimate for the uncertainty, the ±4% error of V1 (25 °C) was accounted for so that one arrives at an imprecision of ±0.5 (MPa)0.5 for the individual HSP components and about ±1 (MPa)0.5 for δ1. These values are congruent with precision claims by Abbott [9] and Stefanis [22] on determining HSP from experimental data. Therefore, a margin of ±1 (MPa)0.5 is viewed a realistic lower error boundary for δ1. Table 2 lists calculation results and literature data.
While δ1 by Hansen agrees well with the calculations in this study, it remains unclear which thermodynamic conditions Abbott et al.’s value is associated with. Park et al. obtained results by HVK for the solid particle state at 25 °C, which appear doubtful since we arrive at (δd, δp, δh) = (18.6, 12.4, 28.2) when using their molar volume (V1 = 106.7 cm3/mol), suggesting some error with respect to δp in [24]. This seems likely in view of Park et al. citing work by Ravindra et al. that had been relevant to their calculus but was later shown to be substantially erroneous [25].
Solvents. The solubility parameters δ2 of all solvents relate to 20 °C or 25 °C. For methyl hydro peroxide, the Hansen components were not taken from [9] but, instead, calculated by the HVK method yielding at 20 °C (δd, δp, δh) = (14.9, 13.1, 21.7) and δ2 = 29.4, each in units of (MPa)0.5. With δ2 = 36.7 (MPa)0.5, according to [9], the difference stays roughly within the δ V 0.5 dependency, so that there appear to be no major discrepancies other than those arising from double the molar volume. Like δ1, an error margin of Δδ2 = ±1 (MPa)0.5 is assumed a reasonable figure for the solvents. All solvent values are presented in the overall results table further below.

3.3. Solute Fusion Parameters

LAA is known to possess low thermal stability and to commence decomposing near the melting point [26,27,28,29,30]. Consequently, the assessment of fusion parameters is not only dependent on the measurement method, but also on experimental conditions, such as heating rates, causing discrepancies in the literature data, which should be used with due care. Corvis et al. [30] reference eight Differential Scanning Calorimetry (DSC) investigations, reporting melting points from 190.4 to 197.0 °C with a median of 192.9 °C. Kofler quoted fusion temperatures of 191 °C (Kofler bench), 185–190 °C (hot-stage microscope), 189–194 °C (capillary tube), and 188 °C (Differential Thermal Analysis, DTA) [31]. Error estimates for the melting temperature were attempted by the Thermodynamic Research Centre [32] using the NIST Thermodata Engine [12] and by the Design Institute for Physical Property Data [14], all aiming to provide comprehensive measures of the overall data reliability. These evaluated data for Tm,1 are compiled in Table 3. Weighting these individual estimates with the reciprocal uncertainty yields a mean of 190.6 °C with a presumed most probable error of ±5 °C. The latter value is justified by allowing for a safety margin in view of an average compound error of ±4.25 °C, and an uncertainty of ±4 °C (95% CI) for the most comprehensive error estimation published [12] (pp. 3–1, 3–30).
Because LAA degradation and melting phenomena appear to superimpose the exact effects experimental conditions have on fusion enthalpy, measurements are not yet fully understood. In consequence, there still is considerable uncertainty in the literature as regards ΔH1. Corvis [30] and Ziderman [33] quote seven values from six studies ranging from 37.0 to 47.1 kJ/mol. While the 95% CI (42.7 ± 4.0) kJ/mol agrees well with the mid-range value (42.1 kJ/mol), the median of 45.3 kJ/mol hints at an asymmetric, and thus, non-Gaussian distribution, rendering the use of mean and standard deviation doubtful. For the lack of more consistent data, the rounded mid-range value and half the (max-min) span are taken as “true” heat of fusion and its experimental uncertainty, respectively.

3.4. Solubilities

Calculated values of x1 are compiled in Table 4 together with the quartiles Qi and the span ratio (Q3/Q1) as determined from the MC simulation. Because symmetric boxcar probability distributions were used and x1 = f(ΔH1, Tm,1, V1, δ1, δ2) is a monotonous function with regard to each variable, the median solubility from simulation is expected to meet Q2x1 within statistical scatter, which was confirmed by the results. Also, as is expected from Equation (1), increasing solubility is matched by a decreasing absolute solubility parameter difference Δδ = |δ1δ2|. The numerical effect the material parameter uncertainties have on solubility will be analyzed in the following.
Because all calculated solubilities meet x1 < 0.002, the solvent volume fraction is well approximated by Φ2 = 1 so that Equation (1) is simplified and the partially factorized expression for x1 reads:
x 1 = e Δ H 1 R · ( 1 T m , 1 1 T ) · e V 1 R T · ( δ 1 δ 2 ) 2 .
Keeping T at 25 °C, both the sensitivity of x1 to the five material parameters and their individual contributions to the uncertainty of x1 can approximately be assessed by computing the solubility for the maximal and minimal value of any one variable while the other variables stay fixed at their respective mid-range values. Since x1 is a continuous and monotonous function in all variables, the span ratio R = (maximum x1)/(minimum x1) rates the influence a material parameter has on the target quantity solubility. Values for this uncertainty measure or, in the wording of the GUM, sensitivity measure are listed in Table 5. The exact expressions that were used for the calculation of the span ratio R are compiled in Table S2.
Apparently, the solubility parameter errors Δδi cause the principal largest uncertainty equally in x1. Their influence on solubility prognosis is the more dominant it is, then the larger the solvent-solute difference Δδ. Since Δδi = ±1 (MPa)0.5 already is a rather optimistic error estimate, it is unlikely that solubility uncertainty can be improved to this end. The second strongest influence is exerted by the melting enthalpy. Its effect is independent of the solvent, and it is the only prominent factor for small Δδ. Reducing the relative error of the fusion enthalpy to ±4% would substantially lower the corresponding (max/min) ratio from 17.8 to 1.6, offering a promising route towards a more precise solubility prediction. In the following, the MC results and experimental literature data presented in Figure 1 and Figure 2 are treated. The solid lines in Figure 1 correspond, from top to bottom, to the 100%, 97.5%, 75%, 25%, 2.5%, and 0% percentile from simulation. The dotted line represents the analytical computation results with nodes (filled squares) for each solvent.
Besides the solubility as calculated from Equation (1), Figure 1 also depicts the 100%, 95%, and 50% CI as derived from the MC simulation. The diagram reveals expected parabolic behavior for allCIs, as well as decreasing solubility uncertainty (narrowing CIs) as δ1 and δ2 become more similar. The latter point is taken up again in more detail in Figure 2.
The finite number of simulations makes extreme values improbable to occur and, in line with the GUM approach, the 95% CI or interquartile range should be considered to provide a meaningful measure of uncertainty. The quotients of conjugate confidence limits range from 2.1 to 11 (50% CI) and 2.8 to 483 (95% CI), respectively, indicating substantial imprecision.
Additionally, in Figure 1, experimental solubilities at 20 °C and 25 °C from [34,35,36] are plotted (See Table S3 for numerical data). A difference of 5 °C in solution temperature results in an 8% to 30% change in solubility dependent on solvent and study. The inter-study variability of x1 was found to be marked, reaching even 435% (!) for some solvents for unresolved reasons. Here, despite their magnitude, these discrepancies can be disregarded when compared with the ‘intrinsic’ uncertainty as evaluated above or with the systematic offset of the measured data. Compared with the 100% (50%) percentile the experimental values are larger by a factor 5 to 103 (20 to 2∙104). Generally, higher solubility can be caused by a decrease of ΔH1, Tm,1, V1, or Δδ = |δ1δ2|, or any combination thereof. The first three causes are ruled out on the grounds of uncertainties of these material parameters being already bracketed by the 100% CI as displayed in Figure 1. Considering the solubility parameters, it is evident that shifting all δ2 values equally by some −5 (MPa)0.5 could make the estimated and experimental data match. This would correspond to a systematic and equally large error of the same sign for all solvents, which is considered extremely improbable. On the other hand, Δδ also decreases if δ1 (solute) becomes smaller. A non-linear least-squares fitting of Equation (1) to the experimental values optimizing δ1, and leaving all other parameters fixed at their respective mid-range values, consistently arrives at δ1 ≈ 27 (MPa)0.5, some 24% below the value of 35.3 (MPa)0.5 and derived according to HVK. Apparently, a too high solute solubility parameter could, in principle, explain the observed differences, but other reasons, such as a general insufficiency of the underlying theory or oversimplification due to dropping the heat capacity term in deriving Equation (1), could also play a role. Normally, this term is positive and effectively increases solubility to some extent making, in principle, predicted and experimental data match better. Because, (a) there is no indication as to the actual value of this term, and (b) it is considered very unlikely that a difference of 3 to 4 log(x1)-units can be consistently explained along this line, this point needs to be left unaddressed. Another explanation for systematically predicting too low solubilities could be a solute activity coefficient that is consistently overestimated by the model for all solvents. With due caution, this could possibly hint at some specific solute-solvent interactions playing a role in this context [37]. Exploring this issue further would be beyond the scope of this study.
Overall, the results show that the scatter of the deterministically calculated solubility values is strong due to the inherent inaccuracies of the material data. The quality of the material data of LAA is well within the usual range and can be classified as typical. Furthermore, the accuracy of the determination of the HSP values by group contribution methods is approximately the same for all compounds. If one adds that the concrete way of calculating the solubilities (purely deterministic versus numerical-iterative) could have an influence, one can assume that the statements about LAA are also valid for compounds with similar material data uncertainty and solubility.

4. Conclusions

Solubility estimates for l-ascorbic acid utilizing the Hildebrand-Scott relation are subject to severe uncertainty. Quantifying uncertainty for solvents with solubility parameters between 14.5 and 33.0 (MPa)0.5 yields upper and lower confidence limits that are up to a factor of 10 (50% CI) and 500 (95% CI) apart, respectively. With the physicochemical data on LAA available to date, solubility uncertainty originates mainly from the HSP of both solvent and solute, and the solute heat of fusion. The most promising route towards achieving better accuracy is to obtain values for the heat of fusion that are more precise. Computed and experimental solubilities show the same qualitative behavior, although quantitative discrepancies are distinct, reaching one to three orders of magnitude for exact reasons still to be resolved. Further, it is assumed that the prediction uncertainty deduced for LAA in this work is generally to the same extent for any HSP-based solubility estimates of other compounds that possess comparable material data quality and solubility. In consequence, it appears that the notion of a semi-quantitative solubility prediction using the HSP methodology should be abandoned.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-3417/10/12/4266/s1, Table S1: Data used for computation of solvent solubility parameters (δ2). Table S2: Expressions used for the calculation of the solubility maximum-minimum ratio. Table S3: Experimental and calculated mole fraction solubilities x1 of l-ascorbic acid.

Author Contributions

Conceptualization, Methodology, Formal analysis, Data curation, Writing—review & editing, R.J.L. and R.S.; Software (simulation), R.J.L.; Writing—original draft, R.J.L. All authors have read and agreed to the submitted manuscript version.

Funding

This research received no external funding.

Acknowledgments

The article processing charge was funded in part by the Baden-Wuerttemberg Ministry of Science, Research and Culture within the Open Access Publishing programme.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hildebrand, J.H.; Scott, R.L. The Solubility of Nonelectrolytes, 3rd ed.; American Chemical Society Monograph No. 17; Reinhold Publishing Corp.: New York, NY, USA, 1950; pp. 27, 122, 129, 271. [Google Scholar]
  2. Prausnitz, J.M. Molecular Thermodynamics of Fluid-Phase Equilibria; Prentice-Hall Inc.: Englewood Cliffs, NJ, USA, 1969. [Google Scholar]
  3. Hildebrand, J.H.; Prausnitz, J.M.; Scott, R.L. Regular and Related Solutions; Van Nostrand Reinhold Co.: New York, NY, USA, 1970. [Google Scholar]
  4. Shinoda, K. Principles of Solution and Solubility; Marcel Dekker: New York, NY, USA, 1978. [Google Scholar]
  5. Barton, A.F.M. CRC Handbook of Solubility Parameters and Other Cohesion Parameters, 2nd ed.; CRC Press: Boca Raton, FL, USA, 1991; pp. 95–107, 177–193, 341–367. [Google Scholar]
  6. Jouyban, A. Handbook of Solubility Data for Pharmaceuticals; CRC Press: Boca Raton, FL, USA, 2010; pp. 30–31. [Google Scholar]
  7. Van Krevelen, D.W.; te Nijenhuis, K. Properties of Polymers. Their Correlation with Chemical Structure, Their Numerical Estimation and Prediction from Additive Group Contributions, 4th ed.; Elsevier: Amsterdam, The Netherlands, 2009; pp. 71–108, 189–227. [Google Scholar]
  8. Hansen, C.M. Hansen Solubility Parameters. A User’s Handbook, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2007; pp. 1–43, 345–483. [Google Scholar]
  9. Abbott, S.; Hansen, C.M.; Yamamoto, H. Hansen Solubility Parameters in Practice, 5th ed.; Hansen-Solubility: Boca Raton, FL, USA, 2015; pp. 1–38, 227–239. ISBN 978-0-9551220-2-6. [Google Scholar]
  10. Chertkoff, M.J.; Martin, A.N. The solubility of Benzoic Acid in Mixed Solvents. J. Am. Pharmaceut. Assoc. 1960, 49, 444–447. [Google Scholar] [CrossRef]
  11. ISO/IEC Guide 98-3:2008, Uncertainty of Measurement. Part 3: Guide to the Expression of Uncertainty in Measurement, Revision 2010, and Supplement 1: Propagation of Distributions Using a Monte Carlo Method; International Organization for Standardization: Geneva, Switzerland, 2008.
  12. CRC Handbook of Chemistry and Physics, 95th ed.; CRC Press: Boca Raton, FL, USA, 2014; pp. 3-1–3-553.
  13. Yaws, C.L. Chemical Properties Handbook; Mc-Graw Hill: New York, NY, USA, 1998; pp. 13, 195. [Google Scholar]
  14. Daubert, T.E.; Danner, R.P. Part 3, extent 1993. In Physical and Thermodynamic Properties of Pure Chemicals: Data Compilation; Taylor and Francis: New York, NY, USA, 1989; pp. II-1–V-1, Table C6H8O6. [Google Scholar]
  15. Hvoslef, J. The Crystal Structure of L-Ascorbic Acid, ‘Vitamin C’. I. The X-ray Analysis. J. Acta Cryst. 1968, B24, 23–35. [Google Scholar] [CrossRef] [PubMed]
  16. Banipal, T.S.; Singh, H.; Banipal, P.K.; Singh, V. Volumetric and viscometric studies on L-ascorbic acid, nicotinic acid, thiamine hydrochloride and pyridoxine mcmcloride in water at temperatures (288.15–318.15) K and at atmospheric pressure. Thermochim. Acta 2013, 553, 31–39. [Google Scholar] [CrossRef]
  17. Dhondge, S.S.; Deshmukh, D.W.; Paliwal, L.J.; Dahasahasra, P.N. The study of molecular interactions of ascorbic acid and sodium ascorbate with water at temperatures (278.15, 288.15 and 298.15) K. J. Chem. Thermodyn. 2013, 67, 217–226. [Google Scholar] [CrossRef]
  18. Singh, D.; Sharma, G.; Gardas, R.L. Exploration of the solvation behaviour of ascorbic acid in aqueous solutions of 1,2,4-triazolium based ionic liquid. J. Mol. Liquids 2017, 244, 55–64. [Google Scholar] [CrossRef]
  19. Sharma, R.; Thakur, R.C.; Sani, B.; Kumar, H. Properties of L-Ascorbic Acid in Water and Binary Aqueous Mixtures of D-Glucose and D-Fructose at Different Temperatures. Russ. J. Phys. Chem. A 2017, 91, 2389–2396. [Google Scholar] [CrossRef]
  20. Rieche, A.; Hitz, F. Über Monomethyl-hydroperoxyd. (III. Mitteilung über Alkylperoxyde). Ber. Dtsch. Chem. Ges. 1929, 62, 2458–2474. [Google Scholar] [CrossRef]
  21. SciFinder (Copyright American Chemical Society). Available online: https://scifinder.cas.org (accessed on 21 March 2020).
  22. Stefanis, E.; Panayiotou, C. A new expanded solubility parameter approach. Int. J. Pharm. 2012, 426, 29–43. [Google Scholar] [CrossRef] [PubMed]
  23. Hansen, C.M. Polymer additives and solubility parameters. Prog. Org. Coat. 2004, 51, 109–112. [Google Scholar] [CrossRef]
  24. Park, K.A.; Lee, H.J.; Hong, I.K. Solubility prediction of bioantioxidants for functional solvent by groupd contribution method. J. Ind. Eng. Chem. 2010, 16, 490–495. [Google Scholar] [CrossRef]
  25. Lehnert, R.J.; Kandelbauer, A. Comments on “Solubility parameter of chitin and chitosan” Carbohydrate Polymers 36 (1998) 121–127. Carbohydr. Polym. 2017, 175, 601–602. [Google Scholar] [CrossRef] [PubMed]
  26. Lerdkanchanaporn, S.; Dollimore, D.; Alexander, K.S. A thermogravimetric study of ascorbic acid and its exipients in pharmaceutical formulations. Thermochim. Acta 1996, 284, 115–126. [Google Scholar] [CrossRef]
  27. Vemin, G.; Chakib, S.; Rogacheva, S.M.; Obretenov, T.D.; Parkanyi, C. Thermal decomposition of ascorbic acid. Carbohydr. Res. 1998, 305, 1–15. [Google Scholar]
  28. Tan, Y. Study on Thermal Stability and Kinetics of VC by Using DSC. J. Wuxi Univ. Light Ind. 2003, 22, 102–105. [Google Scholar]
  29. Jingyan, S.; Yuwen, L.; Zhiyong, W.; Cunxin, W. Investigation of thermal decomposition of ascorbic acid by TG-FTIR and thermal kinetics analysis. J. Pharm. Biomed. Anal. 2013, 77, 116–119. [Google Scholar] [CrossRef]
  30. Corvis, Y.; Menet, M.-C.; Ne´grier, P.; Lazergesa, M.; Espeau, P. The role of stearic acid in ascorbic acid protection from degradation: A heterogeneous system for homogeneous thermodynamic data. New J. Chem. 2013, 37, 761–768. [Google Scholar] [CrossRef] [Green Version]
  31. Kofler, L.; Sitte, H. Zur Schmelzpunktbestimmung von Substanzen, die unter Zersetzung schmelzen. Mon. Chem. 1950, 81, 619–626. [Google Scholar] [CrossRef]
  32. NIST Chemistry WebBook, Standard Reference Database Number 69. Available online: https://webbook.nist.gov/chemistry/ (accessed on 26 March 2020). [CrossRef]
  33. Ziderman, I.I.; Gregorski, K.S.; Lopez, S.V.; Friedman, M. Thermal Interaction of Ascorbic Acid and Sodium Ascorbate with Proteins in relation to Nonenzymatic Browning and Maillard Reactions of Foods. J. Agric. Food Chem. 1989, 37, 1480–1486. [Google Scholar] [CrossRef]
  34. Shalmashi, A.; Eliassi, A. Solubility of L-(+)-Ascorbic Acid in Water, Ethanol, Methanol, Propan-2-ol, Acetone, Acetonitrile, Ethyl Acetate, and Tetrahydrofuran from (293 to 323) K. J. Chem. Eng. Data 2008, 53, 1332–1334. [Google Scholar] [CrossRef]
  35. Neto, A.C.R.; Pires, R.F.; Malagoni, R.A.; Franco, M.R. Solubility of Vitamin C in Water, Ethanol, Propanol-1-ol, Water+Ethanol, and Water+Propan-1-ol at (298.15 and 308.15) K. J. Chem. Eng. Data 2010, 55, 1718–1721. [Google Scholar] [CrossRef]
  36. Galvão, A.C.; Robazza, W.S.; Bianchi, A.D.; Matiello, J.A.; Paludo, A.R.; Thomas, R. Solubility and thermodynamics of vitamin C in binary liquid mixtures involving water, methanol, ethanol and isopropanol at different temperatures. J. Chem. Thermodyn. 2018, 121, 8–16. [Google Scholar] [CrossRef]
  37. Acree, W.E. Thermodynamic Properties of Nonelectrolyte Solutions; Academic Press Inc.: Orlando, FL, USA, 1984; pp. 197–239. [Google Scholar]
Figure 1. Molar solubility x1 of l-ascorbic acid at 25 °C from simulation (solid lines) and analytical computations (■). Experimental data at 20 °C (○) and 25 °C (●) from references [34,35,36]. Dashed vertical line indicates the total solubility parameter of l-ascorbic acid δ1 = 35.3 (MPa)0.5.
Figure 1. Molar solubility x1 of l-ascorbic acid at 25 °C from simulation (solid lines) and analytical computations (■). Experimental data at 20 °C (○) and 25 °C (●) from references [34,35,36]. Dashed vertical line indicates the total solubility parameter of l-ascorbic acid δ1 = 35.3 (MPa)0.5.
Applsci 10 04266 g001
Figure 2. Relative width of the confidence intervals (CI) of l-ascorbic acid solubilities at 25 °C from the Monte Carlo simulation. The width is expressed as ratio = UCL/LCL of the respective upper and lower confidence limit of the numerical simulation results.
Figure 2. Relative width of the confidence intervals (CI) of l-ascorbic acid solubilities at 25 °C from the Monte Carlo simulation. The width is expressed as ratio = UCL/LCL of the respective upper and lower confidence limit of the numerical simulation results.
Applsci 10 04266 g002
Table 1. Estimates for the molar volume V1 of l-ascorbic acid.
Table 1. Estimates for the molar volume V1 of l-ascorbic acid.
Method V1 (cm3/mol)RemarkT (°C)Ref.
X-ray103.7crystallographic unit cell volumer.t.[15]
aq. solution density105.4partial molar volume (median), values from 102.38 to 106.4925[16,17,18,19]
particle density 106.7solid state: purity and crystallinity unknown25[12,13,14]
melt density109.5 ± 4subcooled liquid state: values20[13,14]
109.9 ± 4extrapolated from melt data, presumed probable error ± 4 %25
Table 2. Hansen solubility parameters (δd, δp, δh) and total solubility parameter δ1 of l-ascorbic acid computed according to Hoftyzer and Krevelen.
Table 2. Hansen solubility parameters (δd, δp, δh) and total solubility parameter δ1 of l-ascorbic acid computed according to Hoftyzer and Krevelen.
State of Soluteδdδpδhδ1ReferenceV1T
(cm3/mol)(°C)
subcooled liq.18.112.027.835.3this work109.925
subcooled liq.18.212.027.935.4this work109.520
subcooled liq.18.015.727.436.4Hansen [23]109.925
solid particle18.620.428.239.5Park et al. [24]106.725
unspecified18.011.725.533.3Abbott et al. [9]124.8unspecified
Units for all solubility parameters are in (MPa)0.5. Calculations based on ascorbic acid molar volume V1 (T). Presumed uncertainty Δδ1 = ±1 (MPa)0.5.
Table 3. Thermal fusion properties of l-ascorbic acid.
Table 3. Thermal fusion properties of l-ascorbic acid.
ParameterRangeRemarkRef.
Fusion temperature Tm,1 (°C)188 ± 5DTA measurements[31,32]
191 ± 2Kofler bench measurements[31,32]
191 ± 4compound value, various sources[12]
192 ± 6none[14]
190.6 ± 5weighted mean ± probable errorthis study
Heat of fusion37.0–47.1DSC measurements from 7 studies[30,33]
ΔH1 (kJ/mol)42.1 ± 5mid-range value ±½∙(max-min)this study
Table 4. Calculated mole fraction solubility x1 of l-ascorbic acid at 25 °C for various solvent solubility parameters δ2.
Table 4. Calculated mole fraction solubility x1 of l-ascorbic acid at 25 °C for various solvent solubility parameters δ2.
Solventδ2 (MPa)0.5x1aInterquartile Range of x1 b95% CI b
Q1Q3 Q 3 Q 1 U C L L C L
Pentane14.501.13 × 10−113.30 × 10−123.63 × 10−1111.0483.2
n-Hexane14.902.34 × 10−118.46 × 10−127.22 × 10−118.5424.9
Heptane15.304.79 × 10−111.67 × 10−111.37 × 10−108.2313.9
Diethyl ether15.496.65 × 10−112.27 × 10−111.96 × 10−108.6363.4
Cyclohexane16.806.19 × 10−102.23 × 10−101.74 × 10−97.8282.9
n-Butyl acetate17.411.65 × 10−95.52 × 10−104.19 × 10−97.6209.0
Xylene18.104.79 × 10−91.74 × 10−91.25 × 10−87.2171.9
Ethyl acetate18.155.24 × 10−91.88 × 10−91.37 × 10−87.3174.1
Toluene18.165.33 × 10−91.93 × 10−91.24 × 10−86.4165.0
Benzene18.518.92 × 10−93.63 × 10−92.44 × 10−86.7161.9
Tetrahydrofuran19.463.53 × 10−81.36 × 10−88.00 × 10−85.9128.3
Acetone19.946.80 × 10−82.82 × 10−81.61 × 10−75.7106.7
Acetic acid21.374.35 × 10−72.08 × 10−79.37 × 10−74.562.2
Propan-2-ol23.585.35 × 10−62.75 × 10−61.03 × 10−53.738.9
Acetonitrile24.401.22 × 10−56.66 × 10−62.32 × 10−53.534.8
1-Propanol24.601.48 × 10−58.17 × 10−62.96 × 10−53.630.4
Ethanol26.527.76 × 10−54.55 × 10−51.27 × 10−42.819.8
Methyl hydroperoxide29.405.03 × 10−43.46 × 10−47.76 × 10−42.29.9
Methanol29.415.03 × 10−43.38 × 10−47.26 × 10−42.110.4
Ethylene glycol32.951.83 × 10−31.20 × 10−32.49 × 10−32.12.8
a: Solubility was calculated utilizing solvent solubility parameters δ2 from Table S1 and material data from Table 2 and Table 3. For all calculations, LAA solubility parameter δ1 = 35.3 (MPa)0.5 was used. b: Upper and lower confidence limits (UCL, LCL) of the 95% confidence interval (CI) and quartiles (Q1, Q3) of the distribution of x1 values were obtained by Monte Carlo calculations. Results are presented in ascending order of solubility.
Table 5. Influence of variable uncertainty on molar solubility x1 of l-ascorbic acid at 25 °C.
Table 5. Influence of variable uncertainty on molar solubility x1 of l-ascorbic acid at 25 °C.
VariableUncertaintySpan Ratio
AbsoluteRelativeR = max(x1)/min(x1)
ΔH1±5 kJ/mol(±11.9%)17.8
Tm,1±5 K(±1.1%)1.3
V1±4 cm3/mol(±3.6%)1.6 to 4.0 b
δ1±1 (MPa)0.5(±2.8%)1.5 to 40.0 b
δ2±1 (MPa)0.5(≤±6.8%) a1.5 to 40.0 b
a: error ranges from ±2.7% for δ2 = 32.95 (ethylene glycol) to ±6.8% for δ2 = 14.50 (pentane). b: ratio ranges from lowest for ethylene glycol to largest for pentane.

Share and Cite

MDPI and ACS Style

Lehnert, R.J.; Schilling, R. Accuracy of Molar Solubility Prediction from Hansen Parameters. An Exemplified Treatment of the Bioantioxidant l-Ascorbic Acid. Appl. Sci. 2020, 10, 4266. https://doi.org/10.3390/app10124266

AMA Style

Lehnert RJ, Schilling R. Accuracy of Molar Solubility Prediction from Hansen Parameters. An Exemplified Treatment of the Bioantioxidant l-Ascorbic Acid. Applied Sciences. 2020; 10(12):4266. https://doi.org/10.3390/app10124266

Chicago/Turabian Style

Lehnert, Ralph J., and Richard Schilling. 2020. "Accuracy of Molar Solubility Prediction from Hansen Parameters. An Exemplified Treatment of the Bioantioxidant l-Ascorbic Acid" Applied Sciences 10, no. 12: 4266. https://doi.org/10.3390/app10124266

APA Style

Lehnert, R. J., & Schilling, R. (2020). Accuracy of Molar Solubility Prediction from Hansen Parameters. An Exemplified Treatment of the Bioantioxidant l-Ascorbic Acid. Applied Sciences, 10(12), 4266. https://doi.org/10.3390/app10124266

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop