Third order nonlinear response of liquid water is investigated by using a classical molecular dynamics simulation. The calculated response function yields a marked agreement with the recent experimental results. The physical origins of the third order spectrum is analyzed in terms of the translational and the librational motions. Its temperature dependence and isotope effect are also analyzed.