In our previous study Int. J. Quant. Chem., to appear , we have developed a novel numerical calculation scheme for a dynamics of quantum network for linear molecular aggregates under intense time-dependent electric fields. In this approach, each molecule is assumed to be an electric dipole arranged linearly with an angle from the longitudinal axis, and the molecular interactions are taken into account by adding the radiations from these dipoles to the external electric fields. The effects of the radiations from all the dipoles involve the intermolecular distance, the speed of light, retarded polarization, and its first- and second-order time derivatives at the position of each dipole. The quantum dynamics is performed by solving coupled Liouville equations composed of the Liouville equation for each dipole. In the present study, we develop a calculation approach of nonperturbative second hyperpolarizability g in our novel approach and examine the g of dimer models composed of two-state molecules under the one-photon near resonant intense laser fields. Similar phase transition-like behavior in the field-intensity dependence of the g is observed. We also investigate the second hyperpolarizability spectra in the three-photon resonant region for dimers composed of three-state molecules, which mimic the electronic states of allyl cation. Contrary to the one-photon resonant case, phase transition-like behavior is not observed in the intensity dependence of g in the three-photon resonant region. Q 1999 John Wiley & Sons, Inc. Int J Quant Chem 71: 295)306, 1999