We present new numerical methods, based on a Generalized Space Translation representation (GST), for solving the Time-Dependent Schroedinger equation (TDSE), to treat the nonlinear, nonperturbative interaction of molecules with intense laser pulses. Adopting a Lagrangian, moving coordinate system, we generalize the previous Space Translation method (ST), used in atom–laser interaction problem. In this representation (gauge), the reference system is moving with the laser pulse so the classical movement of free particles in the field, i.e. in the asymptotic region where electron-molecule potentials are negligible but the laser field is still present, is exactly described. As a consequence, the asymptotic quantum wave functions are exact in presence of the laser pulse. To solve numerically the GST–TDSE, all standard discrete propagators for the time discretization and all efficient space discretization methods can be applied. In order to illustrate different possibilities for the choice of discretization methods, we have tested several types of split-operator (SO) and alternating direction implicit (ADI) methods combined with adaptive finite difference methods for the 3D Born–Oppenheimer simulation of H 2 + in a short intense laser pulse field. Our comparison of convergence between the same discretization methods for different gauges have demonstrated the superiority of the GST method. As examples, we present for the first time the ionization of H 2 + exposed to a circularly or linearly polarized intense laser pulse perpendicular to the internuclear axis by an exact 3D simulation.