In this work, the solubility of butyl paraben in 7 pure solvents and in 5 different ethanol-water mixtures has been determined from 1 ˚C to 50 ˚C. The solubility of ethyl paraben and propyl paraben in various solvents has been determined at 10 ˚C. The molar solubility of butyl paraben in pure solvents and its thermodynamic properties, measured by Differential Scanning Calorimetry, have been used to estimate the activity of the pure solid phase, and solution activity coefficients. More than 5000 nucleation experiments of ethyl paraben, propyl paraben and butyl paraben in ethyl acetate, acetone, methanol, ethanol, propanol and 70%, 90% ethanol aqueous solution have been performed. The induction time of each paraben has been determined at three different supersaturation levels in various solvents. The wide variation in induction time reveals the stochastic nature of nucleation. The solid-liquid interfacial energy, free energy of nucleation, nuclei critical radius and pre-exponential factor of parabens in these solvents have been determined according to the classical nucleation theory, and different methods of evaluation are compared. The interfacial energy of parabens in these solvents tends to increase with decreasing mole fraction solubility but the correlation is not very strong. The influence of solvent on nucleation of each paraben and nucleation behavior of parabens in each solvent is discussed. There is a trend in the data that the higher the boiling point of the solvent and the higher the melting point of the solute, the more difficult is the nucleation. This observation is paralleled by the fact that a metastable polymorph has a lower interfacial energy than the stable form, and that a solid compound with a higher melting point appears to have a higher solid-melt and solid-aqueous solution interfacial energy. It has been found that when a paraben is added to aqueous solutions with a certain proportion of ethanol, the solution separates into two immiscible liquid phases in equilibrium. The top layer is water-rich and the bottom layer is paraben-rich. The area in the ternary phase diagram of the liquid-liquid-phase separation region increases with increasing temperature. The area of the liquid-liquid-phase separation region decreases from butyl paraben, propyl paraben to ethyl paraben at the constant temperature. Cooling crystallization of solutions of different proportions of butyl paraben, water and ethanol have been carried out and recorded using the Focused Beam Reflectance Method, Particle Vision and Measurement, and in-situ Infrared Spectroscopy. The FBRM and IR curves and the PVM photos track the appearance of liquid-liquid phase separation and crystallization. The results suggest that the liquid-liquid phase separation has a negative influence on the crystal size distribution. The work illustrates how Process Analytical Technology (PAT) can be used to increase the understanding of complex crystallizations. By cooling crystallization of butyl paraben under conditions of liquid-liquid-phase separation, crystals consisting of a porous layer in between two solid layers have been produced. The outer layers are transparent and compact while the middle layer is full of pores. The thickness of the porous layer can reach more than half of the whole crystal. These sandwich crystals contain only one polymorph as determined by Confocal Raman Microscopy and single crystal X-Ray Diffraction. However, the middle layer material melts at lower temperature than outer layer material. / <p>QC 20130515</p> / investigate nucleation and crystallization of drug-like organic molecules