The accurate determination of solubility of an active pharmaceutical ingredient in a solvent is essential for the design and development of the purification process. The solubility of chlorzoxazone (CHZ) in the pure solvents such as isopropanol, ethyl acetate, 2-ethoxyethyl acetate, and toluene are determined in the temperature range of 278.15 to 333.15 K by shake-flask method at a pressure of 101.2 kPa. The solubility of CHZ in selected solvents is found to increase with temperature. The solubility of CHZ is maximum in 2-ethoxyethyl acetate and minimum in toluene. The mole fraction (x1) solubility of CHZ is correlated with four semi-empirical temperature dependent equations such as the modified Apelblat equation, λh equation, van’t Hoff equation, and non-random two liquid (NRTL) equation. The maximum root-mean-square deviation (RMSD) and relative average deviation (RAD) predicted via four thermodynamic models are 528.80 ×10−4 and 1.146 ×10−2 respectively. The modified Apelblat equation shows a better fit with experimental values for solubility of CHZ. The R2 value obtained from the modified Apelblat equation is >0.99 and the values of RMSD and RAD are smaller for all the selected solvents. The thermodynamic energy properties of the solution are calculated using the van't Hoff equation. The obtained values ΔHm, ΔSm and ΔGm are found to be positive which indicates the dissolution process is an endothermic and entropy-driven process.