In this paper we evaluate the sensitivity of using different solubility models on cooling crystallization. Specifically, the cooling crystallization of acetaminophen in ethanol is investigated. Empirical, correlative thermodynamic (namely van Laar, Wilson, and NRTL) and predictive thermodynamic (namely MOSCED, NRTL-SAC, and UNIFAC) models are considered. Equilibrium solubility model prediction determines the predicted supersaturation profile. The different solubility equations are used within a population balance model for prediction of crystal size properties. Incorrect prediction of the supersaturation profile results in incorrect prediction of crystal size distribution. The NRTL model was found to be more accurate at predicting equilibrium solubility and consequently crystal size. After the solubility sensitivity is evaluated, two methods are proposed to make the crystallization model more robust against solubility model errors. (C) 2010 Elsevier B.V. All rights reserved.