The theoretical diffraction efficiency upper limit of diffractive phase elements (DPE's) with finite apertures is investigated. A successful numerical method of evaluating the efficiency upper bound of DPE's is proposed. The method includes a hybrid optimization procedure that combines a genetic algorithm with the conjugate gradient method. This efficient global optimization technique can also be used to design DPE's. Simulation computations are detailed for rotationally symmetric beam shaping in which a Gaussian profile laser beam is converted into a uniform beam. Numerical results demonstrate that the estimated diffraction efficiency upper bound is consistent with the design results.