In this study, NZVI particles was prepared and studied for the removal of vat green 1 dye from aqueous solution. A four-factor central composite design (CCD) combined with response surface modeling (RSM) to evaluate the combined effects of variables as well as optimization was employed for maximizing the dye removal by prepared NZVI based on 30 different experimental data obtained in a batch study. Four independent variables, viz. NZVI dose (0.1-0.9 g/L), pH (1.5-9.5), contact time (20-100 s), and initial dye concentration (10-50 mg/L) were transform to coded values and quadratic model was built to predict the responses. The significant of independent variables and their interactions were tested by the analysis of variance (ANOVA). Adequacy of the model was tested by the correlation between experimental and predicted values of the response and enumeration of prediction errors. The ANOVA results indicated that the proposed model can be used to navigate the design space. Optimization of the variables for maximum adsorption of dye by NZVI particles was performed using quadratic model. The predicted maximum adsorption efficiency (96.97%) under the optimum conditions of the process variables (NZVI dose 0.5 g/L, pH 4, contact time 60 s, and initial dye concentration 30 mg/L) was very close to the experimental value (96.16%) determined in batch experiment. In the optimization, R2 and R2adj correlation coefficients for the model were evaluated as 0.95 and 0.90, respectively.