The static first hyperpolarizabilities and origin of nonlinear optical (NLO) properties of [(2-methylnaphthyl)imido]hexamolybdates derivatives have been investigated by density functional theory (DFT). The [(2-methylnaphthyl)imido]hexamolybdate has considerable large first hyperpolarizability, 6.780 x 10(-30) esu, and it is larger than that of [(2,6-dimethylphenyl)arylimido]hexamolybdate due to the double aromatic rings in the naphthylimido ligand. The naphthylimido ligand acts as an electron-donor and the polyanion acts as an electron-acceptor. The substituent position on the naphthylimido is a key factor to determine the first hyperpolarizability of (naphthylimido)hexamolybdate derivatives. The derivative, which the iodine atom locates on the para nitrogen on the naphthylimido ligand, has the largest beta(0) value among the iodine-substituted derivatives. It suggests that the iodine atom is quasi linear with nitrogen and Mo, which is bonded to the nitrogen atom, could generate a large static electronic field and give the large contribution to NLO response. The introducing of electron-donors significantly enhances the first hyperpolarizabilities of (naphthylimido)hexamolybdates comparing with the electron-acceptors as the electron-donating ability is significantly enhanced when the electron-donor is attached to the naphthylimido segment. The present investigation provides important insight into NLO properties of (arylimido)molybdate derivatives.