Numerical simulations of self-focusing of ultrafast laser pulses

Simulation of nonlinear propagation of intense ultrafast laser pulses is a hard problem, because of the steep spatial gradients and the temporal shocks that form during the propagation. In this study we adapt the iterative grid distribution method of Ren and Wang [J. Comput. Phys. 159, 246 (2000)] to solve the two-dimensional nonlinear Schrödinger equation with normal time dispersion, space-time focusing, and self-steepening. Our simulations show that, after the asymmetric temporal pulse splitting, the rear peak self-focuses faster than the front one. As a result, the collapse of the rear peak is arrested before that of the front peak. Unlike what has sometimes been conjectured, however, collapse of the two peaks is not arrested through multiple splittings, but rather through temporal dispersion.