Solubility has been extensively investigated by the phase equilibria approach at the mesoscale level, but its origin on the molecular and electronic levels is poorly understood. This study explored the solubility behaviour of crystalline solid in selected pure solvents with various functional groups by using both phase equilibria and molecular modelling methods. The model compound tridecanedioic acid (TDDA) solubility in methanol, ethanol, acetic acid, acetone, and ethyl acetate was determined from T = (283.15 to 323.15) K by a static method. It was found that almost all solutions studied exhibit non-ideal behaviour and deviate positively from Raoult's law indicating the important role of homo-molecules interactions. Thermodynamic analyses of solution suggest that both enthalpy and entropy of solution govern the dissolution process. Computational studies on solubility behaviour were performed by using both density functional theory (DFT) calculations and molecular dynamic (MD) simulations. The results conclude that the (solute + solvent) interaction is not the only factor determining solubility, and (solvent + solvent) interaction also plays an important role. The simulated results are found to be qualitatively consistent with experimental values. Finally, solubility values were correlated by the empirically modified Apelblat equation and two local composition models of Wilson and NRTL. (C) 2015 Elsevier Ltd. All rights reserved.