A general implementation of the lowest nonvanishing order perturbation theory for the calculation of molecular multiphoton ionization cross sections is proposed in the framework of density functional theory. Bound and scattering wave functions are expanded in a multicentric basis set and advantage is taken of the full molecular point group symmetry, thus enabling the application of the formalism to medium-size molecules. Multiphoton ionization cross sections and angular asymmetry parameters have been calculated for the two- and four-photon ionization of the H2+ molecule, for linear and circular light polarizations. Both fixed and random orientations of the target molecule have been considered. To demonstrate the efficiency of the proposed methodology, the two-photon cross section and angular asymmetry parameters for the HOMO and HOMO-1 orbital ionization of benzene are also presented.
