CONTEXT: One of the prevalent methods for evaluating separation performance is to predict the interactions of solvent and solute molecules. The infinite dilution activity coefficient, Gibbs free energy, and sigma profiles provided insights into the solubilization of a solute and revealed the intensity of the solution's molecular interactions. The effective thermodynamic tools (infinite dilution activity coefficient, Gibbs free energy) were evaluated for predicting the efficiency of 18 polar and non-polar organic solvents in rubber seed oil (RSO) extraction. An infinite dilution activity coefficient was computed to evaluate the solubility of the rubber seed oil model compound (linoleic acid) in the organic solvents. Gibbs free energy was evaluated to show the energy change associated with the molecules mixing process and forecast the miscibility of linoleic acid molecules in the solvents. Moreover, the study examined the sigma profiles and sigma surfaces of organic solvents and linoleic acid to acquire a deeper insight into their similarities and how they interact molecularly. According to the computational prediction and experimental verification, the thermodynamic properties of Gibbs free energy and activity coefficient proved to be highly effective tools for screening polar and moderately polar solvents, predicting the molecular interactions with solute. Whereas the sigma profile and sigma surface were found to be the most efficient tools for evaluating the efficacy of non-polar solvents. Solvents with moderate polarity, such as tetrahydrofuran and diethyl ether, as well as non-polar solvents like pentane, heptane, and n-hexane, proved to be effective and favorable for oil extraction, resulting in the highest oil yields of approximately 27.0%. Overall, the COSMO-RS method demonstrates its utility in estimating the solubility of RSO in organic solvents, enabling early identification of the most effective solvent.METHODS: The initial geometry optimization of every component was conducted through density functional theory (DFT) using TmoleX software. A single-point density functional theory (DFT) computation using Becke Perdew 86 (BP86) and the Triple-Zeta Valence Potential (TZVPD) was performed to produce.cosmo files. COSMO-RS calculations were performed by applying the parameterization file BP_TZVPD_FINE_19.ctd using COSMOthermX software. The practical extraction of oil from plant seeds was performed using a sonicator bath to verify the accuracy of the COSMO-RS predictions.