The heat of fusion can also be used to predict
solubility for solids in liquids. Provided an
ideal solution is obtained the
mole fraction (x_2) of solute at saturation is a function of the heat of fusion, the
melting point of the solid (T_\text{fus}) and the
temperature (T) of the solution: \ln x_2 = - \frac {\Delta H^\circ_\text{fus}}{R} \left(\frac{1}{T}- \frac{1}{T_\text{fus}}\right) Here, R is the
gas constant. For example, the solubility of
paracetamol in water at 298
K is predicted to be: x_2 = \exp {\left[- \frac {28100 ~\text{J mol}^{-1}} {8.314 ~\text{J K}^{-1} ~\text{mol}^{-1}}\left(\frac{1}{298 ~\text{K}}- \frac{1}{442 ~\text{K}}\right)\right]} = 0.0248 Since the molar mass of water and paracetamol are and and the density of the solution is , an estimate of the solubility in grams per liter is: • \frac{0.0248 \times \frac{1000 ~\text{g L}^{-1}}{18.0153 ~\text{g mol}^{-1}}}{1-0.0248} \times 151.17 ~\text{g mol}^{-1} = 213.4 ~\text{g L}^{-1} • 1000 g/L * (mol/18.0153g) is an estimate of the number of moles of molecules in 1L solution, using water density as a reference; • 0.0248 * (1000 g/L * (mol/18.0153g)) is the molar fraction of substance in saturated solution with a unit of mol/L; • 0.0248 * (1000 g/L * (mol/18.0153g)) * 151.17g/mol is the solute's molar fraction equivalent mass conversion; • 1-0.0248 will be the fraction of the solution that is solvent. which is a deviation from the real solubility (240 g/L) of 11%. This error can be reduced when an additional
heat capacity parameter is taken into account.
Proof At
equilibrium the
chemical potentials for the solute in the solution and pure solid are identical: \mu^\circ_\text{solid} = \mu^\circ_\text{solute}\, or \mu^\circ_\text{solid} = \mu^\circ_\text{liquid} + RT\ln X_2\, with R\, the
gas constant and T\, the
temperature. Rearranging gives: RT\ln X_2 = -\left(\mu^\circ_\text{liquid} - \mu^\circ_\text{solid}\right)\, and since \Delta G^\circ_\text{fus} = \mu^\circ_\text{liquid} - \mu^\circ_\text{solid}\, the heat of fusion being the difference in chemical potential between the pure liquid and the pure solid, it follows that RT\ln X_2 = -\left(\Delta G^\circ_\text{fus}\right)\, Application of the
Gibbs–Helmholtz equation: \left( \frac{\partial \left( \frac{\Delta G^\circ_\text{fus} } {T} \right) } {\partial T} \right)_{p\,} = -\frac {\Delta H^\circ_\text{fus}} {T^2} ultimately gives: \left( \frac{\partial \left( \ln X_2 \right) } {\partial T} \right) = \frac {\Delta H^\circ_\text{fus}} {RT^2} or: \partial \ln X_2 = \frac {\Delta H^\circ_\text{fus}} {RT^2} \times \delta T and with
integration: \int^{X_2=x_2}_{X_2 = 1} \delta \ln X_2 = \ln x_2 = \int_{T_\text{fus}}^T \frac {\Delta H^\circ_\text{fus}} {RT^2} \times \Delta T the result is obtained: \ln x_2 = - \frac {\Delta H^\circ_\text{fus}} {R}\left(\frac{1}{T}- \frac{1}{T_\text{fus}}\right) ==See also==