Predicting the solubility of mixtures of sugars and their replacers using the Flory-Huggins theory.

In this paper we investigate whether the Flory-Huggins theory can describe the thermodynamics of solutions of simple carbohydrates, like sugars and polyols. In particular, we focus on the description of the solubility of the carbohydrates in water. This is investigated for both binary and ternary mixtures, having two types of these carbohydrates. This research question arises especially in the case of bakery products, where one seeks to replace sucrose with other simple carbohydrates - which are often polyols. Based on the model parameters obtained from fitting the theory to the experimental data of binary solutions, we show that the theory can predict (a) solubility data for ternary mixtures, over a broad range of concentrations and temperatures, and (b) the deliquescence point of binary mixtures of carbohydrate crystals as a function of temperature. The theory can even be applied to carbohydrates, which form hydrate crystals. Together with our earlier theories on the thermodynamics of complex food mixtures, we have now a complete thermodynamic framework to describe the phase and state transitions of food materials as confectionery and bakery products, where the question of sucrose replacement is urgent.


Introduction
In this paper we investigate whether the Flory-Huggins theory can predict the solubility of simple sugars and their replacers like polyols. Moreover, we are interested in the solubility of ternary mixtures of simple carbohydrates, being either polyols or sugars. This research is the initial step towards a thermodynamic description of bakery and confectionery products, where ( part of ) the sugars are replaced with other carbohydrates like polyols. The reformulation of confectionery and bakery products has proven to be difficult, because of the fact that sugars like sucrose have multiple functionality in these food products, apart from providing sweetness. [1][2][3] Sugars have a major impact on the phase transitions of biopolymers such as starch and gluten in bakery products, 4,5 the boiling point of water, and the sorption of water. 6,7 Often, the crystalline nature of sucrose contributes significantly to the texture of food products like biscuits, 2 and confectionery like fondants and creams. [8][9][10][11][12][13] During their dissolution in the biscuit dough mixing stage, the dissolving sugars compete for water with gluten. The degree of gluten hydration will determine the spread and rise of the biscuit. 14 During the cooling stage after baking the dissolved sugars can recrystallize, which is said to provide the fractured texture of the biscuit surface. 15 Research shows that sucrose cannot be simply replaced by a single compound. A mixture of compounds must be used to provide the various functions provided by sucrose in the original formulation. [16][17][18] As part of this mixture of sugar replacers one includes low molecular weight plasticizers like monosaccharides 14,19,20 or polyols. 21,22 Hence, for reformulation purposes it is important to understand the thermodynamics of complex mixtures including sugars and their replacers. The thermodynamics of complex food systems often helps also to understand food structuring processes. 23,24 This is exemplified by the use of the supplemented state diagram in the analysis of cookie baking, 25 bread baking, 26 and the expansion of starchy snacks. 27,28 In the supplemented state diagram one plots the glass transition, and other phase transitions as a function of moisture content. One can obtain a profound understanding of the food structuring process if one plots the processing path in the supplemented state diagram. Hereby, it is important to note that during processing physical systems will follow the phase transition lines like the boiling line. 25,27 At this moment little attention has been given to the thermodynamics of complex food matrices, comprising of at least three components: say a biopolymer, water and a carbohydrate. However, they are important food categories. In bakery products there are two biopolymers present: starch and glutenwhile in confectionery products a single biopolymer like gelatin is used. In bakery products also oils and fats are used, but they do not influence much the thermodynamics of hydrophilic compounds as biopolymers and carbohydrates. With respect to thermodynamics they can be regarded as inert fillers. In various empirical models predicting water activity from food composition this hypothesis has also been taken. [29][30][31][32] There it is clearly observed that water activity is determined only by compounds that are soluble in water.
In some of our recent papers we have shown that the freevolume-Flory-Huggins theory is able to describe the thermodynamics of complex food systems, such as (1) mixtures of proteins, salt and water, 33 and (2) mixtures of biopolymers, disaccharides and water. 7 In this paper we take a first step towards the description of complex mixtures as found in bakery and confectionery products containing sugar replacers. In particular we investigate ternary mixtures of two carbohydrates (either a sugar or a polyol) and water. Because of the important functionality of crystalline sugars, we give emphasis to the prediction of the solubility of these ternary mixtures.
To achieve the above goal, we extend first the Flory-Huggins theory towards polyols. Thereby, an important question we like to resolve is whether the Flory-Huggins interaction parameter for sugars and polyols is temperature dependent or not. Earlier, we have found that for maltodextrins and starch the interaction parameter is temperature independent. This we have assumed to hold also for low molecular weight carbohydrates. However, studies have stated that the water activity of sucrose is temperature dependent. 34 The knowledge on the temperature dependency of the Flory-Huggins interaction parameter is important for predicting the solubility of carbohydrates (mixtures) as a function of temperature. Furthermore, with this understanding we can construct the supplemented state diagram for sugars like sucrose, which is stated to be essential for understanding the structuring process of sweet cookie baking. 25 Subsequently, we will model the solubility for the binary system first, being aqueous solutions of carbohydrates, with the crystals being anhydrates or monohydrates. Furthermore, with the help of the theory we will construct the supplemented state diagram for sucrose and glucose. The latter forms monohydrate crystals, and thus the construction of the supplemented state diagram is more challenging than the case of sucrose. Only partial state diagrams of glucose are shown in earlier papers. 35,36 Finally, the Flory-Huggins theory will be applied to ternary mixtures of water, and two carbohydrates, where one of them may form monohydrate crystals. For the ternary system we can use most of the parameters already determined for binary systems. Only, the interaction parameter between the two carbohydrates has to be estimated.
Like sucrose, sugar replacers are often added to the food formulation in crystalline form (such as sweet biscuits 2 ). In the mixtures of crystals the phenomenon of deliquescence can occur. 37 At a defined water activity the crystals dissolve into liquid, which is absorbed from a water vapour phase at or above this defined water activity. In the mixtures of crystals the deliquescence occurs at the eutonic point. For the mixtures of sugars and polyols, these deliquescence or eutonic points have been determined. 38,39 These data will serve as an extra validation of our theory.

Flory-Huggins theory for ternary mixtures
Originally, the Flory-Huggins theory is developed to describe the thermodynamic interaction between the blends of synthetic polymers or synthetic polymer solutions. 40 Recently, we have shown in a multitude of papers that the theory also applies to food ingredients like polysaccharides (starch, dextran, and maltodextrins), proteins, and even low molecular weight sugars. 7,41,42 To describe these systems in their glassy state a free-volume extension of the Flory-Huggins theory is required, as developed by Vrentas and Vrentas. 43 For a detailed description of the application of these theories to food materials we refer the reader to our earlier papers. However, to make this paper more self-contained we have added a small introduction to the basics of the theory. As crystallization and dissolution occurs only if the temperature is above the glass temperature, we will not discuss the free-volume extension.
The Flory-Huggins theory is formulated in terms of a free energy function, from which chemical potentials and osmotic pressures can be derived. The free energy function is quite similar to that of regular solution theory. In both theories the free energy function has two contributions: an entropic, ideal mixing part and an enthalpic contribution, which is second order in composition (similar to the virial extension of the osmotic pressure). This second term is proportional with the so-called Flory-Huggins interaction parameter χ.
In the Flory-Huggins theory one calculates with volume fractions instead of molar fractions, because the long polymer molecule has multiple degrees of freedom. For synthetic polymers the interaction parameter accounts for the van der Waals interaction between polymers or the polymer and its solvent. Furthermore, it is assumed that the polymer is a flexible, linear chain.
Biopolymers and carbohydrates interact with their solvent (water) via hydrogen bonds, for which the original Flory-Huggins theory does not account for. For systems with hydrogen bonding fundamental theories do exist, but they are much more involving than the Flory-Huggins theory and the theory has only been applied to a few case studies, like PEG and PEO. 6,44,45 Later, it has been shown that these fundamental theories can be approximated by a composition dependent interaction parameter. Hence, in our application of the Flory-Huggins theory to food materials the interaction parameter has lost its fundamental meaning.
The Flory-Huggins theory is best introduced for a polymer solution. The polymer volume fraction is denoted by ϕ p = 1 − ϕ w , and N p is the ratio of the molar volumes of the polymer monomer and the solvent. The free energy density function is: where R is the gas constant, ν w is the molar volume of the solvent, and T is the absolute temperature. The term RT takes into account the primary effect of temperature. From the free energy density function one can derive chemical potentials, and thus equilibrium conditions, cf. Flory. 40 Below, we will state the chemical potentials for ternary mixtures of sugars, polyols and water using the Flory-Huggins theory. Note that crystallization or dissolution can only occur, if the system is not in the glassy state. Below the glass transition temperature the time scales for solute diffusion are so slow that crystallization or dissolution does not occur on practical time scales as that in laboratory experiments or industrial processing. Hence, it is commonly assumed in food science that crystallization practically does not occur. [46][47][48] The free energy density functional for ternary mixtures, according to the Flory-Huggins theory, is as follows: 49 where ν w is the molar volume of water, N i = ν i /ν w is the ratio of molar volumes, and χ ij is the Flory-Huggins interaction parameter. The chemical potential of component i, μ i , is derived using: 50μ where i ≠ j. ϕ i is the volume fraction of compound i, and μ i ¼ μ i =ν i , with ν i being the molar volume. The large extent of the free energy functional makes it easy to produce errors in the algebra for deriving the chemical potential. Henceforth, we have used Maple V (Maplesoft) to derive the equations for the chemical potentials. The fact that it is easy to make errors is evident by the inconsistencies between the expressions of various papers. 49,[51][52][53] We find that only the expressions in ref. 51 are in agreement with our derivation. Hence, the chemical potentials are: The index i = s indicates the primary carbohydrate (often sucrose), and index i = g indicates the secondary carbohydrate (often the sugar replacer).
The binary case follows for ϕ g = 0. The binary interaction parameters can be estimated from the sorption isotherm, rendering the water activity a w : where ϕ s = 1 − ϕ w .

Chemical potential of a crystalline phase
The chemical potential for the crystalline carbohydrate μ X follows from the general expression: where ΔH(T ) is the difference in enthalpy between the solid and the liquid phase, and ΔS(T ) is the difference in entropy between the solid and the liquid phase. For sugars and polyols the enthalpy and specific heat change significantly with temperature: where T x is the melting point of the carbohydrate crystal, and the specific heat difference between its solid and solubilized state is ΔC p,x , changing linearly with temperature as follows: Hence, the chemical potential of the crystal phase becomes: 36

Conditions for equilibrium
First, we use the thermodynamic theory for computing the solubility of simple carbohydrate forming anhydrate crystals. The solubility is governed by the thermodynamic equilibrium between the solution and crystalline phase. This is computed from the equality of the chemical potential of the carbohydrate in the solution and solid phase: μ s = μ X . For a carbohydrate dissolved in water the chemical potential for sugar follows the binary Flory-Huggins theory: χ sw is the interaction parameter between sucrose and water. For computing the solubility in ternary mixtures, containing water and two carbohydrateswith one of them in the crystalline phase, one still has the condition: μ s = μ X . However, now the chemical potential of the dissolved carbohydrate is computed using eqn (4).
The aqueous mixtures of two carbohydrates can also form an eutonic mixture, where two crystalline phases are in equili-brium with the solution phase. At a defined temperature the eutectic mixture is only at one specific composition. The thermodynamic condition for the eutonic point is: where s and g indicate the two carbohydrates. μ X,s and μ X,g are the chemical potentials of their crystalline phases. μ s and μ g are the chemical potentials of the two carbohydrates in the ternary solution, which are again computed using eqn (4). The eutonic point can be computed using the above conditions, and numerical techniques from the field of optimization. 54

Hydrate crystals
Various sugars form hydrate crystals, for which the thermodynamics is not yet described in the Flory-Huggins framework. In the framework of the UNIFAC theory the thermodynamics of hydrate crystals has been described, 55 taking a similar approach to that for cocrystals, salt hydrate and gas hydrate crystals. 56,57 However, these formulations are on a molar basis, instead of a volumetric basis as in the Flory-Huggins theory. We will derive a similar equation using the minimization procedure of the free energy density functional, cf. 58 The total free energy density functional is: For a system of an aqueous carbohydrate solution in equilibrium with a hydrate crystal, we can write it as: where X indicates the crystalline phase. We assume that the crystalline phase has n h water molecules per carbohydrate molecule. For monohydrates n h = 1. In experiments often only the total amounts of carbohydrates and water are given. Their total volume fractions are given by ϕ s,tot and ϕ w,tot (which add up to unity). The total carbohydrate and water will be partitioned over the crystalline and solution phase via: where ϕ s and ϕ w are the volume fractions in the solution phase. ϕ X is the total volume fraction of the crystalline phase. Hence, the free energy can be described in terms of ϕ X : Note that ϕ X is the only variable in the above expression of the free energy functional. The system will seek a minimum in the free energy density, which is given by: or, rather: the relative volume fraction the hydrate water occupies in the crystalline phase. We note that an identical expression for the equilibrium has been found for the solid-liquid transition of eye lens proteins, which also form hydrate crystals. 59 Hence, we are confident that the above derivation is valid.

Sorption isotherms of binary systems
In this section we investigate the water activity of various polyol solutions, as they have not yet been analysed by the Flory-Huggins theory. Furthermore, we revisit the water activity of sugar solutions, mainly because of the question whether the Flory-Huggins interaction parameter is temperature dependent or not.
For our analysis we have made use of data sources from the literature, which we have listed in Table 1. First, we have used only data in the range of 20-45 degrees Celsius in order to estimate the interaction parameters. We have fitted the FH-theory to the water activity data, as shown in Fig. 1. Here, the symbols represent the experimental data, and the lines with the corresponding color represent the fitting using the Flory-Huggins theory. The fitted values of the interaction parameter χ ws are listed in Table 4.
Furthermore, we show the interaction parameter as a function of the molar weight of the carbohydrate in Fig. 2. As one can observe χ ws is a smooth function of the molar weight. The sorption of polyols and sugars of similar molecular weights (like sorbitol and glucose) cannot be distinguished, as shown in Fig. 1. Thus, it appears that with respect to water sorption polyols and sugars can be treated as a single class of compounds, which are only characterized by their molar weight, similar to maltodextrins and dextrans. 6 However, the Flory-Huggins interaction parameter is composition dependent for maltodextrins with the degree of polymerization larger than 2: where χ 0 = 0.5 for all maltodextrins and polysaccharides, but χ 1 depending on the molar weight M s . For small sugars and polyols it holds that χ 0 = χ 1 ≤ 0.5. The data in Fig. 2 fit nicely with the curve of χ 1 as a function of molar weight. Here, we have also included the χ 1 values for maltodextrins.
The temperature dependence is investigated with extra datasets, spanning a wide range of temperatures. We have analysed the data on the basis of the excess chemical potential: The experimental data points are shown in Fig. 3, where we have used a similar color coding as in Fig. 1. We have compared the data with the expression for μ ex , cf. the Flory-Huggins theory, which is indicated with lines having similar color coding to the experimental datapoints. If water sorption of carbohydrates is independent of temperature, the data would coincide with the theoretical prediction. As one can observe, all experimental data for each compound fall onto a single curve, and follows the FH prediction. Hence, we conclude that the Flory-Huggins interaction parameter is indeed independent of temperature, for polyols, sugars and polysaccharides.

Binary systems
We have investigated the solubility of sucrose in water using the above theory. Because of its importance for understanding  Table 1. Data of simple sugars are duplicated from ref. 6. From top to bottom: ethyl glycol (magenta), glycerol (black), erythritol (yellow), xylitol (blue), sorbitol, mannitol, glucose, fructose (all red), sucrose, maltose (green), and maltotriose (cyan). Different symbols refer to different data sources.  Table 4. food structuring processes in complex food matrices containing sugar replacers, we will present our results in the supplemented state diagram of sucrose, Fig. 4. The phase transition lines are indicated using the available experimental data (symbols) from literature sources, 34,55,66,[74][75][76] and predictions made using the Flory-Huggins theory with the temperature independent interaction parameter of χ ws = 0.52. Earlier, we have modelled the glass transition data, 77 which are copied into the figure.
The solubility data are fitted using the following parameters: ΔH x = 57 kJ mol −1 , T x = 470 K, ΔC p,x = 260 J mol −1 K, and Δγ p,x = 0.008 J mol −1 K 2 . Their values are in a similar range to those in previous studies, 36,55,78 albeit that there is definitely some variation between these sources. Further investigation shows that there are multiple parameter sets that give reasonable fitting to the data. We note that the enthalpy, as defined in eqn (7), decreases very significantly, and approaches zero near the freezing line. Nevertheless, from the figure we observe that the Flory-Huggins theory can accurately predict all thermodynamic behaviour of a binary sucrose-water mixture.
Subsequently, we have focussed on the solubility of various sugar replacers, namely a multitude of polyols and xylose. We note that all investigated carbohydrates form anhydrate crystals. The experimental data from the literature are again compared with the fitted lines using the Flory-Huggins theory. The comparison between experimental data and theory can be observed in Fig. 5. The experimental data for polyols are obtained from ref. 79, while that for xylose are obtained from ref. 80. Most of the parameter values describing μ X are available in the literature. 79,81,82 Only Γ x has to be estimated, and in some cases ΔC p,x . All used values are listed in Table 4. As one can observe the values for all carbohydrates are in a similar range. Now, we turn the focus on binary systems, with the solute forming hydrate crystals. We will investigate glucose, lactose and maltose, all forming monohydrate crystals. As above, we will also list the estimated parameters in Table 4.
First, we investigate the case of glucose monohydrate crystals. The experimental data are obtained from ref. 35, 36, 55, and 66. This case was investigated earlier using group theory such as UNIFAC. 36,55,66 We have fitted the Flory-Huggins theory to solubility data for both glucose anhydrous and monohydrate crystals, with results shown in Fig. 6. Because of      Our values are quite comparable to other values in the literature. 36,55,66 However, it must be noted that there is a lot of variance in the values of the estimated parameters in these references. Similar to the above case of sucrose, there are multiple sets of parameter values that give similar fitting to the experimental data. This is not very surprising given the fact that there are four parameters involved.
Subsequently, we have tested the theory to the solubility of lactose, in the temperature range where it forms monohydrate crystals. The experimental data are obtained from ref. 78 and 83. The comparison between experimental data and theory is shown in Fig. 7. The fitted parameters are: T x = 487 K, ΔH x = 80 kJ mol −1 , ΔC p,x = 300 J mol −1 K −1 , and γ x = 0.015 J mol −1 K −2 .
Also, we have investigated the solubility of maltose in the temperature range, where it forms monohydrates. The experimental data are obtained from ref. 78. The experimental data and theory are compared in Fig. 8. The fitted parameters are: T x = 425 K, ΔH x = 73 kJ mol −1 , ΔC p,x = 400 J mol −1 K −1 , and γ x = 0.007 J mol −1 K −2 .

Ternary systems
In food materials, where one desires to lower the sucrose content, one often uses a mixture of sucrose and another small molecular weight carbohydrate, being a polyol or other sugar. The addition of a second carbohydrate substantially changes the solubility of sucrose, and the total solubility of all solutes. Here, we like to note that it is a direct consequence of thermodynamics that the total solubility of mixtures is always higher than that of the individual compounds. This follows directly from the Gibbs-Duhem relationship. 89,90 Here, we investigate whether the Flory-Huggins theory can predict the solubility of both solutes in the ternary mixtures of two carbohydrates and water. Experimental data on such ternary mixtures are quite scarce. The mixtures we investigate are indicated in Table 2, together with references to the experimental data, and the temperature range the measurements are performed. One can observe that the collection of experimental data also comprises of ternary mixtures containing monohydrate crystals. These datasets will provide a good validation of our theory for monohydrate crystals.
For the comparison of the theory and experiment, we will use most of the parameter values, as determined for binary   mixtures, as listed in Table 4. We have to estimate only a single extra parameter, namely the interaction parameter between the two carbohydrates, χ sg . We expect the value of χ sg to be small, 7 and that the prediction is not very sensitive to this value. We first investigate the solubility of sucrose in aqueous glycerol solutions. The results are shown in Fig. 9. Using parameters found for the binary systems, we have fitted the ternary Flory-Huggins theory to the data with χ 12 = 0.06. It is obvious that the mixture of polyol and water cannot be viewed as an effective solvent for sucrose, as is done in the case of biopolymers. 5 That is because sucrose hardly dissolves in pure glycerol, in the temperature range of 15-35. 91 At high volume fractions the saturated solution is very viscous, and it is hard to separate crystals from the solution. 91 We think that this explains the discrepancy between theory and experiment. Other experimental data points are well explained.
Subsequently, we have investigated the solubility of sucrose in invert sugar, which is an equimolar mixture of glucose and fructose (Fig. 10). Invert sugar does not crystallize. Given the fact that χ 1w = 0.35 for both glucose and fructose, one can treat invert sugar as a single non-crystallizing compound. We have fitted the ternary Flory-Huggins theory, resulting in χ 12 = 0.0. The experimental data are obtained from multiple literature sources, in the temperature range 20-70, with increments of 10 degree Celsius. The FH-theory predicts the experimental solubility data quite well for almost all temperatures, except for 20 degrees, which is probably due to experimental errors.
Another dataset comprises of solubility in the ternary mixture of sucrose, sorbitol and water. 87 Parameters for the crystalline state of sorbitol are obtained from ref. 92. The theory and the experimental data are compared in Fig. 11. They agree quite well using χ sg = 0.08. Observe that the figure shows two solubility lines, for sucrose (solid line), as well as for sorbitol (dashed line). The two solubility lines intersect at the eutonic point. Now, we turn our focus to ternary mixtures, where one of the carbohydrates can form a hydrate crystal. First, we investigate the ternary mixture of sucrose, glucose and water. Note that, as we have fitted χ sg already for the ternary mixture of sucrose and invert sugar, the comparison between theory and experiment can be done without any parameter fitting. All required parameter values are obtained above. The results are   shown in Fig. 12. One branch of the solubility line at T = 23°C relates to glucose monohydrate crystals, while this branch for T = 70°C relates to glucose anhydrate crystals, which are more stable at these temperatures, see the phase diagram of glucose - Fig. 6. One can observe that the theory gives quite good predictions of the solubility lines. Again, the two solubility lines intersect at the eutonic point. We observe that the eutonic point depends on the temperature. The final dataset we have analysed is that of the ternary mixture of sucrose, maltose and water. Solubility is measured at 30°C, where maltose monohydrate crystals are formed. We have obtained a good fit of the theory with experiment using χ sg = −0.22, as shown in Fig. 13. The non-linearity in the experimental data points in the sucrose branch is probably due to experimental errors. For other ternary mixtures the FHtheory predicts rather well the sucrose branches, which are all quite linear curves.

Eutonic points from deliquescence
As a final validation of our theory we compute the eutonic point of two carbohydrate crystals in equilibrium with their ternary aqueous solution. For the aqueous liquid that is in equilibrium with the two crystals, we can compute the water activity using the Flory-Huggins theory for ternary solutions, eqn (4). We compare the computed a w values with the experimental data on deliquescence of binary systems, a w,mix , taken from Mauer and coworkers. 39 Predictions and experimental values are presented in Table 3. Note that in the investigated temperature range glucose crystals are present in the monohydrate form, as follows from the glucose state diagram, see Fig. 6.
We have tuned the interaction parameter between the two carbohydrates in the range |χ sg | ≤ 0.1 to match the experi-mental data. For the predictions we have used the following values: χ sg = −0.05 for the sucrose/glucose mixture, χ sg = 0.1 for the sucrose/xylitol mixture, and χ sg = −0.05 for the xylitol/ glucose mixture. One can observe that these values are in a similar range to the values of χ sg obtained from solubility data, as presented in Table 2.
For the xylitol/glucose system we have obtained a very good prediction of a w at the eutonic point for all temperatures. For other systems there is a slight difference, but the trend with temperature is captured well. Concerning the experimental data Mauer and coworkers remark that for a dry mixture of binary crystals, poor contact between crystal types leads to deviations of the experimental observed values a w,mix , compared to the theoretical value. Another indication of the experimental error is obtained by comparing the value of a w,mix for sucrose and glucose at 25 degree Celsius from two different studies, 38,39 which are respectively 0.78 and 0.728. Hence, one may expect an experimental error of 0.05 in a w,mix .
Mauer and coworkers have predicted a w,mix using the deliquescence data of single components, and the Ross equation and Clausius-Clapeyron equation. The deviations of our predictions and experimental values are always less than or equal to 0.05. The majority of the predictions of Mauer and co- Fig. 12 Solubility data for the ternary mixture of sucrose, glucose and water, at T = 23°C (red symbols) and T = 70°C (yellow symbols). Lines are predictions based on the Flory-Huggins theory, using parameter sets obtained from fitting theory to the solubility data of the binary glucose/water mixture. Note that at room temperature glucose monohydrate crystals are formed, while at 70°C glucose anhydrate crystals are formed, as follows from the glucose state diagram, in Fig. 6. workers deviate more than the experimental error of 0.05 from the experimental observations. This is due to the fact that no interactions between solutes are assumed (i.e. χ sg = 0 is assumed), and the temperature dependence of their parameters is obtained from a limited temperature range 20 ≤ T ≤ 40°C, while the solubility data we have used are in the range 10 ≤ T ≤ 80°C.

Conclusions
In this paper we have investigated the solubility of various binary and ternary carbohydrate solutions with the Flory-Huggins theory. This investigation has been performed with the aim to be able to predict the solubility of sucrose in food matrices, where sucrose is partly replaced by another low molecular weight carbohydrate like polyols. In various sugar-rich food products the crystallinity of sucrose is an important functionality. Hence, it is advantageous if the solubility can be predicted from theory. In this paper we have shown that the Flory-Huggins theory is capable of doing that. If all parameters are obtained independently via analysis of binary mixtures, the Flory-Huggins theory can predict accurately the solubility of ternary mixtures containing sucrose. With the same theory one can also predict the deliquescence of binary mixtures of carbohydrate crystals, as we have shown above. We note that the enthalpy of carbohydrate crystals is strongly dependent on temperature, requiring three parameters to describe it (ΔH x , ΔC p,x and γ p,x ). These three parameters together with the melting temperature determine the chemical potential of crystals. Due to this number of parameters required to describe the crystalline system, and the variability in the experimental data, there are often multiple parameter sets that show comparable good fitting to the data. This is evident also in the variance of these parameters in the literature, where the theory is fitted to the same system. However, having obtained a parameter set that shows a reasonable fit with binary solutions, one can use this parameter set equally well for the application of the theory to ternary mixtures.
In this paper we have also resolved two additional research questions: (1) whether the Flory-Huggins interaction parameter is temperature dependent or not, and (2) whether the theory can also describe the solubility of sugars forming hydrate crystals. Our analysis shows that the interaction parameter between carbohydrates and water is indeed temperature independent. After minimization of the free energy density functional for the system with a hydrate crystal, we have obtained the condition for equilibrium described in terms of the chemical potentials on a volumetric basisμ i . With this condition we have been able to predict the solubility of sugars forming monohydrate crystals in both binary and ternary aqueous mixtures.
From the findings in this paper we can conclude that the Flory-Huggins theory is a highly valuable tool for the description of the thermodynamics of complex mixtures of carbohydrates. Such mixtures one encounters often in sweet food products, where sucrose is ( partially) replaced by alternative sweeteners like polyols or oligosaccharides. In these products other phase transitions like boiling and biopolymer melting are important, 4,93 which can also be described by the Flory-Huggins theory, as we have shown earlier. 5,6 Hence, the Flory-Huggins theory can describe much of the thermodynamics occurring during the manufacturing of sweet food products with reduced sucrose levels, such as confectionery and bakery products. Understanding of the thermodynamics of these complex food systems is often improved by the use of the supplemented state diagram. 25 As we have shown above for the case of sucrose and glucose, the Flory-Huggins theory can be used to predict all relevant phase transitions.