A new approach is introduced for computing probabilities of molecular transitions in time-dependent fields. The method is based on the stationary (t,t′) representation of the Schrödinger equation and is shown to be equivalent to infinite order time-dependent perturbation theory. Bound-to-bound (i.e., photoexcitation) and bound-to-continuum (i.e., photoreaction) transitions are regarded as reactive collisions with the ‘‘time coordinate’’ as the reaction coordinate in an extended Hilbert space. A numerical method based on imposing absorbing boundary conditions for the time coordinate in a discrete variable representation framework is introduced. A single operation of the Green’s operator provides all the state-specific transition probabilities as well as partial state-resolved (inclusive) reaction probabilities. Illustrative numerical applications are given for model systems.