We present a combined molecular dynamics/quantum chemical perturbation method for calculating the refractive index of liquid water at different temperatures. We compare results of this method with the refractive index obtained from other solvent models. The best agreement with the experimental refractive index of liquid water and its temperature dependence is obtained using correlated gas-phase polarizabilities in the classical Lorentz-Lorenz expression. Also, the iterative self-consistent reaction field approach in the semicontinuum implementation matches the experimental refractive index reasonably well.