Determination of Energy Interaction Parameters for the UNIFAC Model Based on Solvent Activity Coefficients in Benzene–D2EHPA and Toluene–D2EHPA Systems

This study examines the activity coefficients of benzene, toluene, and di-(2-ethylhexyl)phosphoric acid (D2EHPA) in binary benzene–D2EHPA and toluene–D2EHPA systems, as well as the ternary n-hexane–toluene–D2EHPA system, using gas chromatography at 293.0 K. The primary objective was to determine UNIFAC model interaction parameters and validate their accuracy for predicting thermodynamic behavior in these systems. Experimental measurements revealed activity coefficient maxima for benzene and toluene at mole fractions of 0.8–0.9, decreasing to 0.46–0.67 in dilute solutions. The UNIFAC interaction parameters were calculated as follows: ACH–HPO4 (−334, 4605), ACCH3–HPO4 (680, 467), and refined CH2–HPO4 (54, 1199). The UNIFAC model achieved deviations of less than 2% from experimental data in both binary and ternary systems. A novel methodology incorporating intermediate standards for gas chromatography was developed to overcome challenges in measuring volatile solvent concentrations, enhancing measurement precision. These findings enable accurate prediction of activity coefficients in mixtures of alkanes, cycloalkanes, and monoaromatic hydrocarbons with D2EHPA, offering significant implications for optimizing metal liquid–liquid extraction processes.