Solubility data in different solvents are an important issue for separation processes involving complex
molecules such as natural products and pharmaceutical drugs. Nonetheless, solubility data are in general
scarce and difficult to obtain, and so models are important tools to generate the necessary estimates.
Different correlative, statistical and thermodynamic models have been proposed to evaluate solubilities.
From these, the more theoretically sound thermodynamic models allow to generate estimates at
broader temperature, pressure and composition conditions while using a smaller amount of experimental
information. Among these, the cubic-plus-association equation of state that combines the simplicity
and robustness of a cubic equation of state with the Wertheim’s association contribution has been under
attention. In this work, this EoS is for the first time proposed to model organic phase solubilities of
drug-like molecules in a wide range of temperatures.
Solubilities of acetanilide, acetylsalicylic acid, adipic acid, ascorbic acid, hydroquinone, ibuprofen,
paracetamol and stearic acid were estimated in alcohols, ketones, alkanes, esters, acids, aromatics,
chlorinated solvents, as well as in other common solvents. The hydrogen bonding behaviour was explicitly
accounted for with each associating group being treated individually, as well as multiple group substitutions.
Accurate correlations were obtained using a single binary interaction parameter (global AAD of 24.2%), while considering the complexity of the studied systems predictions were generally also satisfactory.