Solvent usage is a major source of environmental waste in pharmaceutical industry. The current paradigm shift toward continuous manufacturing in pharmaceutical industry has renewed the interest in continuous crystallization, which offers the prospect of easy solvent recycling. However, the selection of solvents for an integrated crystallization processes is nontrivial due to the likely trade-off between optimal solvent properties for crystallization and solvent separation and recycling. A systematic approach for the simultaneous optimization of process conditions and solvent selection for continuous crystallization including solvent recycling is presented. A unified perturbed-chain statistical associating fluid theory model framework is applied to predict thermodynamic properties related to solubility and vapor-liquid equilibrium, which is integrated with a process model. A continuous mapping procedure is adopted to solve the optimization problem effectively. A case study based on continuous antisolvent crystallization of paracetamol with solvent separation via flash demonstrates the approach. (c) 2017 American Institute of Chemical Engineers AIChE J, 64: 1205-1216, 2018