Prediction of an extremely large nonlinear refractive index for crystals at terahertz frequencies