In this work, we propose to correlate and predict the solubility in supercritical CO 2 of disubstituted aromatic isomers of hydroxybenzoic acid and methylbenzoic acid with a new methodology based on the expanded liquid theory, in which the solid-fluid equilibrium is modeled using the local composition model of UNIQUAC in which the interaction parameters are related to the solvent reduced density with an empiric exponential form equations. The experimental solubility of hydroxybenzoic acid isomers, methylbenzoic acid isomers and mixed isomers (mhydroxybenzoic acid+p-hydroxybenzoic acid) are used for evaluating the correlation and prediction capabilities of this new methodology. The results obtained using the proposed model show good agreement with the experimental data used.