A simple correction to the infinite dilution activity coefficient computed via molecular simulation for a nonelectrolyte solid solute in solution is proposed. The methodology adopts the concept that the activity coefficient may be fundamentally interpreted as a product of a residual and combinatorial term. The residual contribution is assumed to be insensitive to concentration, and the combinatorial term is modeled using the athermal Flory-Huggins theory. The proposed method uses only properties for the solute computed at infinite dilution to estimate solution-phase properties at finite concentrations. Properties of the pure solid solute are estimated using the group contribution method of Gani and coworkers, allowing for efficient blind solubility predictions to be made. The method is applied to predict the solubility of solid phenanthrene in 17 different solvents. For all cases, the combinatorial correction lowers the predicted solubility relative to the infinite dilution approximation, and in general, improves agreement with experiment. (c) 2013 American Institute of Chemical Engineers AIChE J, 59: 2647-2661, 2013