Time-dependent perturbation: a recursive generation of the transition amplitudes between bound states

We present a new expansion of the solution to the time-dependent Schrödinger equation it ∂U/∂t = {H₀ + V(t)}U. A complete set of eigenvectors of H₀ spanning the Hilbert space in which H₀ and V operate is partitioned in two subsets. Transition amplitudes from the first subspace to the second one are calculated by building an adequate series of intermediate representations with respect to the couplings which produce the transition from the model space into its orthogonal complement. These expansions yield a coupling matrix series V⁽ⁿ⁾ for which an iterative solution V⁽ⁿ⁾ → V⁽ⁿ⁺¹⁾ is derived. This solution leads to a recursive numerical treatment of the calculation of transition amplitudes. A simple example, concerning a harmonic oscillator under an intense laser field, is considered using a Fourier analysis of the perturbation.