The (hyper)polarizabilities of GaN, GaP, and GaAs clusters have been calculated using an ab initio timedependent Hartree-Fock method and an even-tempered Gaussian basis set. The geometries of the clusters used in the study were optimized by ab initio Hartree-Fock calculations using the same basis set. The clusters used in this calculation are of the type Ga m X n (m ) 1, 3, 4 and n ) 1, 3, 4) where X ) N, P, or As. The Ga m X n clusters are in a charge neutral (q ) 0) state for m ) n and in an appropriately charged state for m * n. All of the clusters studied are stable in a nearly tetrahedral geometry. The magnitude of the calculated (hyper)polarizabilities appears to depend on the size of the cluster. The cluster-size dependence of calculated (hyper)polarizabilities is more pronounced for the first hyperpolarizability, β, than for the polarizability, R. The calculated frequency-dependent polarizability, R(ω), exhibits a small but finite dispersion with the frequency of the optical field. The magnitude of the calculated β(-ω σ ;ω 1 ,ω 2 ) corresponding to various second-order effects shows the expected trend of β(-2ω;ω,ω) > β(-ω;0,ω) > β(0;0,0).