Multiconfigurational second-order perturbation theory computations and reaction path mapping for the retinal protonated Schiff base models all-trans-nona-2,4,6,8-tetraeniminium and 2-cis-nona-2,4,6,8-tetraeniminium cation demonstrate that, in isolated conditions, retinal chromophores exhibit at least three competing excited-state double bond isomerization paths. These paths are associated with the photo-isomerization of the double bonds in positions 9, 11, and 13, respectively, and are controlled by barriers that favor the position 11. The computations provide a basis for the understanding of the observed excited-state lifetime in both naturally occurring and synthetic chromophores in solution and, tentatively, in the protein environment. In particular, we provide a rationalization of the excited-state lifetimes observed for a group of locked retinal chromophores which suggests that photoisomerization in bacteriorhodopsin is the result of simultaneous specific "catalysis" (all-trans --> 13-cis path) accompanied by specific "inhibition" (all-trans --> 11-cis path). The nature of the S-1-->S-0 decay channel associated with the three paths has also been investigated at the CASSCF level of theory. It is shown that the energy surfaces in the vicinity of the conical intersection for the photo isomerization about the central double bond of retinal (position 11) and the two corresponding lateral double bonds (positions 9 and 13) are structurally different.