Geometries, electronic states, and electron affinities of GamPm and GamPn- (m + n = 2-5) clusters have been examined using four hybrid and pure density functional theory (DFT) methods. Structural optimization and frequency analyses are performed with the basis of a 6-311+G(2df) one-particle basis set. The geometries are fully optimized with each DFT method independently. Three types of energy separations reported in this work are the adiabatic electron affinity (EA,,,), the vertical electron affinity (EA(vert)), and the vertical detachment energy (VDE). The calculation results show that the singlet structures have higher symmetry than that of doublet structures. The best method for predicting molecular structures was found to be BLYP, while other methods generally underestimated bond lengths. The most reliable adiabatic electron affinities and vertical detachment energy, obtained at the BP86 and B3LYP level of theory, are predicted to be 2.22 and 2.10 eV (GaP), 2.51 and 2.46 eV (Ga(2)p), 1.86 and 1.94 eV (GaP2), 1.96 and 2.27 eV (GaP3), 1.76 and 1.99 eV (Ga3P), 1.79 and 2.14 eV (Ga2P2), 2.85 and 3.67 eV (GaP4), 2.08 and 2.10 eV (Ga4P), 2.90 and 3.17 eV (Ga2P3), and 2.70 and 3.37 eV (Ga3P2), respectively. Those for Ga2P, Ga3P, Ga2P2, Ga4P, GaP4, Ga2P3, and Ga(3)p(2) are in good agreement with experiment, but the predicted EA(ad) values for GaP, Ga2p, GaP2, and GaP, are larger than the available experimental values. For the vibrational frequencies of the GamPn series, the B3LYP method produces good predictions with the average error only similar to 10 cm(-1) from available experimental and theoretical values. The other three methods overestimate or underestimate the vibrational frequencies, with the worst predictions given by the BLYP method. (c) 2005 Wiley Periodicals, Inc.