The dipole polarizabilities, dipole moments, density of states, and second-order nonlinear optical (NLO) properties of organoimido derivatives of hexamolybdates have been investigated by using time-dependent density functional response theory. This class of organic-inorganic hybrid compounds possesses remarkably large and eye-catching molecular second-order NLO response, especially [Mo6O17(NC16H12NO2)(FeNC10H9)](2-) (7) and [Mo6O17(NC16H12NO2)(NC6H2(NH2)(3))](2-) (6) with static second-order polarizability (beta(vee)) computed to be 15766.27 x 10(-30) esu and 6299.59 x 10(-30) esu, respectively. Thus, these systems have the possibility to be excellent second-order nonlinear optical materials. Analysis of the major contributions to the beta(vec) value suggests that the charge transfer (CT) from polyanion to organic segment (D-A) along the z-axis plays the key role in NLO response; the polyanion acts as a donor (D) whereas organoimido acts as an acceptor (A) in all the studied systems. The computed beta(vec) values increase by incorporation of an electron acceptor (-NO2) at the end of the phenyl ring of the organoimido segment. Furthermore, substitution of amino (-NH2) or ferrocenyl (-FeC10H9) at the outer side of polyanion and an electron acceptor (-NO2) at the end of the phenyl ring in organoimido segment simultaneously is more important to enhance the optical nonlinearity. Orbital analysis shows that the degree of CT between the polyanion and organoimido segments was increased when ferrocenyl donor was introduced. The present investigation provides important insight into the remarkably large NLO properties of organoimido-substituted hexamolybdates.