One of many steps in palm oil refining process is deacidification which aims to separate free fatty acids and other compounds from the oil. The deacidification process was using a green solvent, known as NADES, that consisted of betaine monohydrate and propionic acid at molar ratio of 1:8. In this study, the process conditions were optimized using the response surface method (RSM) through central composite design in order to predict the maximum distribution coefficient of palmitic acid. The obtained regression equation of the basic model for optimization was: y = 0.717 + 0.003x1 + 0.043 x2 + 0.148x3 - 0.005 x1x1 - 0.030 x2x2 + 0.047 x3x3 - 0.008 x1x2 + 0.008 x1x3 + 0.033 x2x3. The independent variables are x1 ≡ temperature (40, 60, 80°C), x2 ≡ amount of palmitic acid in the palm oil (2, 5, 8 %) and x3 ≡ mass ratios of oil to NADES (1:2, 1:1, 2:1). The optimum process condition found was temperature of 62.3°C, palmitic acid content of 8%, and NADES to palm oil mass ratio of 1:2, resulting in the maximum distribution coefficient of 0.96.