Two-temperature model for pulsed-laser-induced subsurface modifications in Si