The imaginary part of the complex potential energy surface (CPES) provides the autoionization, Auger and ICD (Intermolecular Coulombic Decay) ionization decay rates as functions of the variations in the molecular electronic structure and geometry. We introduce here a novel approach which enables the calculations the molecular complex potential energy surfaces by the standard computational algorithms which were originally developed for calculating real PES of stable molecules. As an illustrative numerical example we have calculated the CPES for the molecular autoionization of hydrogen molecule, H-2[(1)Sigma(+)(g)] where the two electrons are doubly excited. (C) 2011 Elsevier B. V. All rights reserved.